High-throughput microsatellite isolation through 454 GS-FLX Titanium pyrosequencing of enriched DNA libraries

Microsatellites (or SSRs: simple sequence repeats) are among the most frequently used DNA markers in many areas of research. The use of microsatellite markers is limited by the difficulties involved in their de novo isolation from species for which no genomic resources are available. We describe here a high-throughput method for isolating microsatellite markers based on coupling multiplex microsatellite enrichment and next-generation sequencing on 454 GS-FLX Titanium platforms. The procedure was calibrated on a model species (Apis mellifera) and validated on 13 other species from various taxonomic groups (animals, plants and fungi), including taxa for which severe difficulties were previously encountered using traditional methods. We obtained from 11,497 to 34,483 sequences depending on the species and the number of detected microsatellite loci ranged from 199 to 5791. We thus demonstrated that this procedure can be readily and successfully applied to a large variety of taxonomic groups, at much lower cost than would have been possible with traditional protocols. This method is expected to speed up the acquisition of high-quality genetic markers for nonmodel organisms.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

Development and characterization of microsatellite markers in a clonal plant,Dioscorea japonica Thunb.

We developed 10 microsatellite loci from genomic DNA of a dioecious clonal plant, Dioscorea japonica. Out of 384 clones, 148 contained microsatellite repeats. Polymerase chain reaction primer pairs were designed for 95 of these clones from their sequence data, of which, 10 pairs produced successful amplification. Thirty‐eight individuals were genotyped for allelic diversity. We detected three to nine alleles per locus, and the expected heterozygosity ranged from 0.461 to 0.851.     

豚鼠基因组26个多态性微卫星标记的筛选

目的筛选豚鼠基因组的多态性微卫星标记,为豚鼠遗传质量控制及基因定位等工作奠定基础。方法采用磁珠富集法和豚鼠基因组数据库筛选法获取微卫星位点序列,通过分析和初步筛选,挑选部分候选位点,根据其序列设计引物,对5种不同来源的豚鼠基因组DNA标本进行PCR扩增,以期获得多态性分子标记。结果本实验采用磁珠富集法共获得微卫星序列304个,设计引物125对,最终获得多态性位点1个,暂未发现多态性的特异性位点17个;用数据库筛选法共获得微卫星序列292个,设计并合成相应引物178对,最终发现多态性位点25个,暂未发现多态性的特异性位点28个。结论本实验获得26个多态性微卫星标记,45个潜在的候选标记,为微卫星标记在豚鼠遗传质量监测及突变基因定位等工作的应用奠定了基础。     

Polymorphic microsatellite loci for the common marmoset (Callithrix jacchus) designed using a cost- and time-efficient method

We describe a cost- and time-efficient method for designing new microsatellite markers in any species with substantial genomic DNA sequence data available. Using this technique, we report 14 new polymorphic dinucleotide microsatellite loci isolated from the common marmoset. The relative yield of new polymorphisms was higher with less labor than described in previous marmoset studies. Of 20 loci initially evaluated, 14 were polymorphic and amplified reliably (70% success rate). The number of alleles ranged from 3 to 9 with heterozygosity varying from 0.48 to 0.83.     

An integrated genetic linkage map of avocado

 An avocado genomic library was screened with various microsatellite repeats. (A/T)_n and (TC/AG)_n sequences were found to be the most frequent repeats. One hundred and seventy-two positive clones were sequenced successfully of which 113 were found to contain simple sequence repeats (SSR). Polymerase chain reaction primers were designed to the regions flanking the SSR in 62 clones. A GenBank search of avocado DNA sequences revealed 1 sequence containing a (CT)₁₀ repeat. A total of 92 avocado-specific SSR markers were screened for polymorphism using 50 offspring of a cross between the avocado cultivars ‘Pinkerton’ and ‘Ettinger’. Both are standard avocado cultivars which are normally outcrossed and highly heterozygous. Fifty polymorphic SSR loci, 17 random amplified polymorphic DNA (RAPD) and 23 minisatellite DNA Fingerprint (DFP) bands were used to construct the avocado genetic map. The resulting data were analyzed with various mapping programs in order to assess which program best accommodated data from progeny of heterozygous parents. The analyses resulted in 12 linkage groups with 34 markers (25 SSRs, 3 RAPDs and 6 DFP bands) covering 352.6 cM. This initial map can serve as a basis for developing a detailed genomic map and for detection of linkage between markers and quantitative trait loci. Received: 2 April 1996 / Accepted: 28 February 1997     

Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L.)

The growing number of rice microsatellite markers warrants a comprehensive comparison of allelic variability between the markers developed using different methods, with various sequence repeat motifs, and from coding and non-coding portions of the genome. We have performed such a comparison over a set of 323 microsatellite markers; 194 were derived from genomic library screening and 129 were derived from the analysis of rice-expressed sequence tags (ESTs) available in public DNA databases. We have evaluated the frequency of polymorphism between parental pairs of six inter- subspecific crosses and one inter-specific cross widely used for mapping in rice. Microsatellites derived from genomic libraries detected a higher level of polymorphism than those derived from ESTs contained in the GenBank database (83.8% versus 54.0%). Similarly, the other measures of genetic variability [the number of alleles per locus, polymorphism information content (PIC), and allele size ranges] were all higher in genomic library-derived microsatellites than in their EST-database counterparts. The highest overall degree of genetic diversity was seen in GA-containing microsatellites of genomic library origin, while the most conserved markers contained CCG- or CAG-trinucleotide motifs and were developed from GenBank sequences. Preferential location of specific motifs in coding versus non-coding regions of known genes was related to observed levels of microsatellite diversity. A strong positive correlation was observed between the maximum length of a microsatellite motif and the standard deviation of the molecular-weight of amplified fragments. The reliability of molecular weight standard deviation (SDmw) as an indicator of genetic variability of microsatellite loci is discussed. Received: 5 May 1999 / Accepted: 16 August 1999     

Sequence-tagged microsatellite profiling (STMP): a rapid technique for developing SSR markers

We describe a technique, sequence-tagged microsatellite profiling (STMP), to rapidly generate large numbers of simple sequence repeat (SSR) markers from genomic or cDNA. This technique eliminates the need for library screening to identify SSR-containing clones and provides an ~25-fold increase in sequencing throughput compared to traditional methods. STMP generates short but characteristic nucleotide sequence tags for fragments that are present within a pool of SSR amplicons. These tags are then ligated together to form concatemers for cloning and sequencing. The analysis of thousands of tags gives rise to a representational profile of the abundance and frequency of SSRs within the DNA pool, from which low copy sequences can be identified. As each tag contains sufficient nucleotide sequence for primer design, their conversion into PCR primers allows the amplification of corresponding full-length fragments from the pool of SSR amplicons. These fragments permit the full characterisation of a SSR locus and provide flanking sequence for the development of a microsatellite marker. Alternatively, sequence tag primers can be used to directly amplify corresponding SSR loci from genomic DNA, thereby reducing the cost of developing a microsatellite marker to the synthesis of just one sequence-specific primer. We demonstrate the utility of STMP by the development of SSR markers in bread wheat.     

Sequence-tagged microsatellite profiling (STMP): a rapid technique for developing SSR markers

We describe a technique, sequence-tagged microsatellite profiling (STMP), to rapidly generate large numbers of simple sequence repeat (SSR) markers from genomic or cDNA. This technique eliminates the need for library screening to identify SSR-containing clones and provides an approximately 25-fold increase in sequencing throughput compared to traditional methods. STMP generates short but characteristic nucleotide sequence tags for fragments that are present within a pool of SSR amplicons. These tags are then ligated together to form concatemers for cloning and sequencing. The analysis of thousands of tags gives rise to a representational profile of the abundance and frequency of SSRs within the DNA pool, from which low copy sequences can be identified. As each tag contains sufficient nucleotide sequence for primer design, their conversion into PCR primers allows the amplification of corresponding full-length fragments from the pool of SSR amplicons. These fragments permit the full characterisation of a SSR locus and provide flanking sequence for the development of a microsatellite marker. Alternatively, sequence tag primers can be used to directly amplify corresponding SSR loci from genomic DNA, thereby reducing the cost of developing a microsatellite marker to the synthesis of just one sequence-specific primer. We demonstrate the utility of STMP by the development of SSR markers in bread wheat.     

Development of a maize molecular evolutionary genomic database

PANZEA is the first public database for studying maize genomic diversity. It was initiated as a repository of genomic diversity for an NSF Plant Genome project on 'Maize Evolutionary Genomics'. PANZEA is hosted at the Bioinformatics Research Center, North Carolina State University, and is open to the public (http://statgen.ncsu.edu/panzea). PANZEA is designed to capture the interrelationships between germplasm, molecular diversity, phenotypic diversity and genome structure. It has the ability to store, integrate and visualize DNA sequence, enzymatic, SSR (simple sequence repeat) marker, germplasm and phenotypic data. The relational data model is selected and implemented in Oracle. An automated DNA sequence data submission tool has been created that allows project researchers to remotely submit their DNA sequence data directly to PANZEA. On-line database search forms and reports have been created to allow users to search or download germplasm, DNA sequence, gene/locus data and much more, directly from the web.     

Lessons learned from microsatellite development for nonmodel organisms using 454 pyrosequencing

Microsatellites, also known as simple sequence repeats (SSRs), are among the most commonly used marker types in evolutionary and ecological studies. Next Generation Sequencing techniques such as 454 pyrosequencing allow the rapid development of microsatellite markers in nonmodel organisms. 454 pyrosequencing is a straightforward approach to develop a high number of microsatellite markers. Therefore, developing microsatellites using 454 pyrosequencing has become the method of choice for marker development. Here, we describe a user friendly way of microsatellite development from 454 pyrosequencing data and analyse data sets of 17 nonmodel species (plants, fungi, invertebrates, birds and a mammal) for microsatellite repeats and flanking regions suitable for primer development. We then compare the numbers of successfully lab‐tested microsatellite markers for the various species and furthermore describe diverse challenges that might arise in different study species, for example, large genome size or nonpure extraction of genomic DNA. Successful primer identification was feasible for all species. We found that in species for which large repeat numbers are uncommon, such as fungi, polymorphic markers can nevertheless be developed from 454 pyrosequencing reads containing small repeat numbers (five to six repeats). Furthermore, the development of microsatellite markers for species with large genomes was also with Next Generation Sequencing techniques more cost and time‐consuming than for species with smaller genomes. In this study, we showed that depending on the species, a different amount of 454 pyrosequencing data might be required for successful identification of a sufficient number of microsatellite markers for ecological genetic studies.     

Distribution of dinucleotide microsatellites in the Drosophila melanogaster genome.

Microsatellites, a special class of repetitive DNA, have become one of the most popular genetic markers. The progress of various genome projects has made it possible to study the genomic distribution of microsatellites and to evaluate the potential influence of several parameters on their genesis. We report the distribution of dinucleotide microsatellites in the genome of Drosophila melanogaster. When considering only microsatellites with five or more repeat units, the average length of dinucleotide repeats in D. melanogaster is 6.7 repeats. We tested a wide range of parameters which could potentially influence microsatellite density, and we did not detect a significant influence of recombination rate, number of exons, or total length of coding sequence. In concordance with the neutral expectation for the origin of microsatellites, a significant positive correlation between AT content and (AT/TA)n microsatellite density was detected. While this pattern may indicate that microsatellite genesis is a random process, we also found evidence for a nonrandom distribution of microsatellites. Average microsatellite density was higher on the X chromosome, but extreme heterogeneity was observed between different genomic regions. Such a clumping of microsatellites was also evident on a more local scale, as 38.9% of the contiguous sequences analyzed showed a deviation from a random distribution of microsatellites.     

Characterization of new Schistosoma mansoni microsatellite loci in sequences obtained from public DNA databases and microsatellite enriched genomic libraries

In the last decade microsatellites have become one of the most useful genetic markers used in a large number of organisms due to their abundance and high level of polymorphism. Microsatellites have been used for individual identification, paternity tests, forensic studies and population genetics. Data on microsatellite abundance comes preferentially from microsatellite enriched libraries and DNA sequence databases. We have conducted a search in GenBank of more than 16,000 Schistosoma mansoni ESTs and 42,000 BAC sequences. In addition, we obtained 300 sequences from CA and AT microsatellite enriched genomic libraries. The sequences were searched for simple repeats using the RepeatMasker software. Of 16,022 ESTs, we detected 481 (3%) sequences that contained 622 microsatellites (434 perfect, 164 imperfect and 24 compounds). Of the 481 ESTs, 194 were grouped in 63 clusters containing 2 to 15 ESTs per cluster. Polymorphisms were observed in 16 clusters. The 287 remaining ESTs were orphan sequences. Of the 42,017 BAC end sequences, 1,598 (3.8%) contained microsatellites (2,335 perfect, 287 imperfect and 79 compounds). The 1,598 BAC end sequences 80 were grouped into 17 clusters containing 3 to 17 BAC end sequences per cluster. Microsatellites were present in 67 out of 300 sequences from microsatellite enriched libraries (55 perfect, 38 imperfect and 15 compounds). From all of the observed loci 55 were selected for having the longest perfect repeats and flanking regions that allowed the design of primers for PCR amplification. Additionally we describe two new polymorphic microsatellite loci.     

A rapid method for detection of plant genomic instability using unanchored-microsatellite primers

In order to assess the feasibility of microsatellite primers as markers for genomic instability, we conducted a study of DNA stability in cauliflower callus. A protocol is described for the rapid screening of a large number of putative variant calli and plants. Genomic DNA is isolated and screened by microsatellite primers. We believe that inter-simple sequence repeat PCRs can conveniently detect and measure common genetic events underlying plant genomic instability. These include deletions, amplifications, translocations, insertions, recombination or chemical alterations. Our results indicate that instability occurred in an early step in the process of callogenesis. The technique is fast, reproducible, and is a new application for ISSR markers.     

Mining new crystal protein genes from Bacillus thuringiensis on the basis of mixed plasmid-enriched genome sequencing and a computational pipeline

We have designed a high-throughput system for the identification of novel crystal protein genes (cry) from Bacillus thuringiensis strains. The system was developed with two goals: (i) to acquire the mixed plasmid-enriched genomic sequence of B. thuringiensis using next-generation sequencing biotechnology, and (ii) to identify cry genes with a computational pipeline (using BtToxin_scanner). In our pipeline method, we employed three different kinds of well-developed prediction methods, BLAST, hidden Markov model (HMM), and support vector machine (SVM), to predict the presence of Cry toxin genes. The pipeline proved to be fast (average speed, 1.02 Mb/min for proteins and open reading frames [ORFs] and 1.80 Mb/min for nucleotide sequences), sensitive (it detected 40% more protein toxin genes than a keyword extraction method using genomic sequences downloaded from GenBank), and highly specific. Twenty-one strains from our laboratory's collection were selected based on their plasmid pattern and/or crystal morphology. The plasmid-enriched genomic DNA was extracted from these strains and mixed for Illumina sequencing. The sequencing data were de novo assembled, and a total of 113 candidate cry sequences were identified using the computational pipeline. Twenty-seven candidate sequences were selected on the basis of their low level of sequence identity to known cry genes, and eight full-length genes were obtained with PCR. Finally, three new cry-type genes (primary ranks) and five cry holotypes, which were designated cry8Ac1, cry7Ha1, cry21Ca1, cry32Fa1, and cry21Da1 by the B. thuringiensis Toxin Nomenclature Committee, were identified. The system described here is both efficient and cost-effective and can greatly accelerate the discovery of novel cry genes.     

Development and polymorphism of simple sequence repeat DNA markers for the evergreen shrub Ilex leucoclada M.

We have developed 13 microsatellite loci from an enrichment library of genomic DNA in the evergreen shrub Ilex leucoclada. One hundred and eighty‐nine out of 432 clones were found to contain microsatellite repeats. Primer pairs were designed for 92 of these clones according to their sequence data. Thirteen of these primer pairs revealed polymorphism among 36 individuals sampled from 12 populations. Three to 27 alleles per locus were detected, and the expected heterozygosity ranged from 0.133 to 0.971. Because these 13 microsatellite markers showed high degrees of genetic variation, they should be useful tools for studying population and ecological genetics of I. leucoclada.     

Silene tatarica microsatellites are frequently located in repetitive DNA

The genomic distribution of microsatellites can be explained by DNA slippage, slippage like processes and base substitutions. Nevertheless, microsatellites are also frequently associated with repetitive DNA, raising the question of the relative contributions of these processes to microsatellite genesis. We show that in Silene tatarica about 50% of the microsatellites isolated by an enrichment cloning protocol are associated with repetitive DNA. Based on the flanking sequences, we distinguished seven different classes of repetitive DNA. PCR primers designed for the flanking sequences of an individual clone amplified a heterogeneous family of repetitive DNA. Despite considerable variation in the flanking sequence (pi = 0.108), the microsatellite repeats did not show any evidence for decay. Rather, we observed the emergence of a new repeat type that probably arose by mutation and was spread by replication slippage. In fact, a complete repeat type switch could be observed among the analysed clones. We propose that the analysis of microsatellite sequences embedded in repetitive DNA provides a hitherto largely unexplored tool to study microsatellite evolution.     

Determination of sex and chimaerism in the domestic sheep by DNA amplification using HMG-box and microsatellite sequences

We describe a rapid, reliable method for the sexing of the domestic sheep (Ovis aries) by amplification of Y-chromosome-specific sequences in male genomic DNA using the polymerase-chain reaction (PCR). Oligonucleotide primers were selected from a conserved sequence, the HMG box, in the sequence of ovine Sry, permitting amplification of a defined 161 bp fragment only from male-specific genomic DNA. As a control, microsatellite primers also were used in PCR reactions, recognising a sequence that is amplifiable in genomic DNA from both males and females. In addition, we demonstrate the feasibility of using this technique for the detection of Y-specific sequences in foetal biopsies (specifically small numbers of foetal germ cells), and in reconstruction mixtures of male and female genomic DNA to simulate the analysis of intersex chimaeras which would be produced when pluripotent cells have been established for this species.     

Mobile elements and the genesis of microsatellites in dipterans

Factors that influence the genesis and genomic distribution of microsatellite DNA are poorly understood. We have identified a novel class of Dipteran mobile elements, mini-me elements, which help elucidate both of these issues. These retroposons contain two internal proto-microsatellite regions that commonly expand into lengthy microsatellite repeats. These elements are highly abundant, accounting for approximately 1.2% of the Drosophila melanogaster genome, giving them the potential to be a prolific source of microsatellite DNA variation. They also give us the opportunity to observe the outcomes of multiple microsatellite genesis events (initiating from the same proto-microsatellite) at separate mini-me loci. Based on these observations, we determined that the genesis of microsatellites within mini-me elements occurs through two separate mutational processes: the expansion of preexisting tandem repeats and the conversion of sequence with high cryptic simplicity into tandemly repetitive DNA. These modes of microsatellite genesis can be generalized beyond the case of mini-me elements and help to explain the genesis of microsatellites in any sequence region that is not constrained by selection.