SSR allelic variation in almond (Prunus dulcis Mill.)

Sixteen SSR markers including eight EST-SSR and eight genomic SSRs were used for genetic diversity analysis of 23 Chinese and 15 international almond cultivars. EST- and genomic SSR markers previously reported in species of Prunus, mainly peach, proved to be useful for almond genetic analysis. DNA sequences of 117 alleles of six of the 16 SSR loci were analysed to reveal sequence variation among the 38 almond accessions. For the four SSR loci with AG/CT repeats, no insertions or deletions were observed in the flanking regions of the 98 alleles sequenced. Allelic size variation of these loci resulted exclusively from differences in the structures of repeat motifs, which involved interruptions or occurrences of new motif repeats in addition to varying number of AG/CT repeats. Some alleles had a high number of uninterrupted repeat motifs, indicating that SSR mutational patterns differ among alleles at a given SSR locus within the almond species. Allelic homoplasy was observed in the SSR loci because of base substitutions, interruptions or compound repeat motifs. Substitutions in the repeat regions were found at two SSR loci, suggesting that point mutations operate on SSRs and hinder the further SSR expansion by introducing repeat interruptions to stabilize SSR loci. Furthermore, it was shown that some potential point mutations in the flanking regions are linked with new SSR repeat motif variation in almond and peach. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Conservation of Nuclear SSR Loci Reveals High Affinity of <Emphasis Type="Italic">Quercus infectoria</Emphasis> ssp. <Emphasis Type="Italic">veneris</Emphasis> A. Kern (Fagaceae) to Section <Emphasis Type="Italic">Robur</Emphasis>

Conservation of 16 nuclear microsatellite loci, originally developed for Quercus macrocarpa (section Albae), Q. petraea, Q. robur (section Robur), and Q. myrsinifolia, (subgenus Cyclobalanopsis) was tested in a Q. infectoria ssp. veneris population from Cyprus. All loci could be amplified successfully and displayed allele size and diversity patterns that match those of oak species belonging to the section Robur. At least in one case, limited amplification and high levels of homozygosity support the occurrence of “null alleles” caused by a possible mutation in the highly conserved primer areas, thus hindering PCR. The sampled population exhibited high levels of diversity despite the very limited distribution of this species in Cyprus and extended population fragmentation. Allele sizes of Q. infectoria at locus QpZAG9 partially match those of Q. alnifolia and Q. coccifera from neighboring populations. However, sequencing showed homoplasy, excluding a case of interspecific introgression with the latter, phylogenetically remote species. Q. infectoria ssp. veneris sequences at this locus were concordant to those of other species of section Robur, while sequences of Quercus alnifolia and Quercus coccifera were almost identical to Q. cerris.     

The Complex Genetic Structure of Sugarcane Limits Identification of Additional SNP-Defined Simplex Alleles in Microsatellite Loci

Cultivated sugarcane possesses a large number of chromosomes (typically >80) which can be organised into eight homology groups, each containing approximately 12 homo(eo)logous chromosomes. Currently, microsatellite (SSR) markers are the most easily used markers for marker-assisted selection and other genetic applications. However, only SSR alleles that segregate as simplex and duplex markers can be incorporated into sugarcane maps using populations of ～300 progeny. Consequently only a subset of possible alleles have been mapped for a given SSR locus and are available for subsequent QTL analyses. Three sugarcane SSR loci, mSSCIR8, mSSCIR17 and mSSCIR18, were amplified, cloned and sequenced from the parents of an Australian sugarcane mapping population, IJ76-514 and . The sequences were examined to identify nucleotide sequence polymorphisms in the flanking regions that could provide additional simplex SSR allele markers. Alignment of the sequences revealed SNP, indel and repeat length variation and the pattern of sequence variation suggested multiple alleles for each SSR, including all of the alleles previously scored by fragment length. While the flanking regions of the SSR loci contained numerous SNPs and indels, none defined new simplex alleles within multiplex fragments and no new SSR simplex alleles were mapped. Furthermore, many of the sequence-defined alleles appear to be spurious and may have arisen from PCR-mediated recombination. These results confirm the difficulties associated with characterising allelic diversity in a polyploid species and the complexity of mapping in sugarcane.     

Evolution of Microsatellite Alleles in Four Species of Mice (Genus Apodemus)

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or ``true' alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species. Received: 3 November 1999 / Accepted: 2 May 2000     

Insertions–Deletions in a Microsatellite Flanking Region May Be Resolved by Variation in Stuttering Patterns

Microsatellites or simple sequence repeats (SSRs) may display polymerase-chain-reaction-amplified fragment lengths mismatching the patterns expected from repeat copy number variation. We sequenced alleles of a nuclear dinucleotide SSR locus in two oak species which showed 2- and 1-bp length differences between alleles and three types of stuttering patterns in fragment length analysis. In accordance with the variation in stuttering, we identified three allele classes characterized by insertions–deletions in the flanking regions and overlapping repeat copy number ranges. Different alleles could thus only be safely separated when considering these stuttering patterns. Our results raise the question of how to adequately delimit alleles when such size homoplasy is present. We advise to thoroughly characterize SSR sequence variation during marker development and to carefully place primer sites along flanking regions to facilitate automated allele scoring and to minimize labor-intensive visual inspection.     

Exon-primed intron-crossing (EPIC) PCR markers of Helicoverpa armigera (Lepidoptera: Noctuidae)

Applying microsatellite DNA markers in population genetic studies of the pest moth Helicoverpa armigera is subject to numerous technical problems, such as the high frequency of null alleles, occurrence of size homoplasy, presence of multiple copies of flanking sequence in the genome and the lack of PCR amplification robustness between populations. To overcome these difficulties, we developed exon-primed intron-crossing (EPIC) nuclear DNA markers for H. armigera based on ribosomal protein (Rp) and the Dopa Decarboxylase (DDC) genes and sequenced alleles showing length polymorphisms. Allele length polymorphisms were usually from random indels (insertions or deletions) within introns, although variation of short dinucleotide DNA repeat units was also detected. Mapping crosses demonstrated Mendelian inheritance patterns for these EPIC markers and the absence of both null alleles and allele 'dropouts'. Three examples of allele size homoplasies due to indels were detected in EPIC markers RpL3, RpS6 and DDC, while sequencing of multiple individuals across 11 randomly selected alleles did not detect indel size homoplasies. The robustness of the EPIC-PCR markers was demonstrated by PCR amplification in the related species, H. zea, H. assulta and H. punctigera.     

What phylogeny and gene genealogy analyses reveal about homoplasy in citrus microsatellite alleles

Sixty-five microsatellite alleles amplified from ancestral citrus accessions classified in three separate genera were evaluated for sequence polymorphism to establish the basis of inter- and intra-allelic genetic variation, evaluate the extent of size homoplasy, and determine an appropriate model (stepwise or infinite allele) for analysis of citrus microsatellite alleles. Sequences for each locus were aligned and subsequently used to determine relationships between alleles of different taxa via parsimony. Interallelic size variation at each SSR locus examined was due to changes in repeat copy number with one exception. Sequencing these alleles uncovered new distinct point mutations in the microsatellite region and the region flanking the microsatellite. Several of the point mutations were found to be genus, species, or allele specific, and some mutations were informative about the inferred evolutionary relationships among alleles. Overall, homoplasy was observed in alleles from all three loci, where the core microsatellite repeat was changed causing alleles of the same size class to be identical in state but not identical by descent. Because nearly all changes in allele size (with one exception) were due to expansion or contraction of the repeat motif, this suggests that a stepwise mutation model, which assumes homoplasy may occur, would be the most appropriate for analyzing Citrus SSR data. The collected data indicate that microsatellites can be a useful tool for evaluating Citrus species and two related genera since repeat motifs were reasonably well retained. However, this work also demonstrated that the number of microsatellite alleles is clearly an underestimate of the number of sequence variants present.     

Cross-species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: implications for the transferability of SSRs in plants

We investigated the transferability of 31 soybean (Glycine max) simple sequence repeat (SSR) loci to wild congeners and to other legume genera. Up to 65% of the soybean primer pairs amplified SSRs within Glycine, but frequently, the SSRs were short and interrupted compared with those of soybeans. Nevertheless, 85% of the loci were polymorphic within G. clandestina. Cross-species amplification outside of the genus was much lower (3%-13%), with polymorphism restricted to one primer pair, AG81. AG81 amplified loci in Glycine, Kennedia, and Vigna (Phaseoleae), Vicia (Vicieae), Trifolium (Trifolieae), and Lupinus (Genisteae) within the Papilionoideae, and in Albizia within the Mimosoideae. The primer conservation at AG81 may be explained by its apparent proximity to the seryl-tRNA synthetase gene. Interspecific differences in allele size at AG81 loci reflected repeat length variation within the SSR region and indels in the flanking region. Alleles of identical size with different underlying sequences (size homoplasy) were observed. Our findings and the emerging patterns in other plant studies suggest that in contrast to animals, successful cross-species amplification of SSRs in plants is largely restricted to congeners or closely related genera. Because mutations in both the SSR region and the flanking region contribute to variation in allele size among species, knowledge of DNA sequence is essential before SSR loci can be meaningfully used to address applied and evolutionary questions.      

Mutational bias in penguin microsatellite DNA

Analysis of nucleotide sequence variation at a microsatellite DNA locus revealed extensive size homoplasy of alleles in Adélie penguins (Pygoscelis adeliae). Variation in the flanking regions at this locus allowed discrimination between mechanisms proposed for length changes in microsatellite DNA alleles. We further examined the structure of alleles for the same microsatellite DNA locus across 11 additional species of penguin (Spheniscidae) by mapping allele sequences onto an independent penguin phylogeny. Our analysis indicated that the repeat motifs appear to have evolved independently on several occasions. We observed sequence instability in the region bordering the repeat tract with a transversional bias predominating. We propose that this bias results from inaccurate DNA replication owing to the sequence context of this repeat tract. Because we show that regions flanking repeat sequences exhibit this mutational bias, this cautions against the use of such regions for phylogeny reconstruction.     

Patterns and relative rates of nucleotide and insertion/deletion evolution at six chloroplast intergenic regions in new world species of the Lecythidaceae

Insertions and deletions (indels) in chloroplast noncoding regions are common genetic markers to estimate population structure and gene flow, although relatively little is known about indel evolution among recently diverged lineages such as within plant families. Because indel events tend to occur nonrandomly along DNA sequences, recurrent mutations may generate homoplasy for indel haplotypes. This is a potential problem for population studies, because indel haplotypes may be shared among populations after recurrent mutation as well as gene flow. Furthermore, indel haplotypes may differ in fitness and therefore be subject to natural selection detectable as rate heterogeneity among lineages. Such selection could contribute to the spatial patterning of cpDNA haplotypes, greatly complicating the interpretation of cpDNA population structure. This study examined both nucleotide and indel cpDNA variation and divergence at six noncoding regions (psbB-psbH, atpB-rbcL, trnL-trnH, rpl20-5'rps12, trnS-trnG, and trnH-psbA) in 16 individuals from eight species in the Lecythidaceae and a Sapotaceae outgroup. We described patterns of cpDNA changes, assessed the level of indel homoplasy, and tested for rate heterogeneity among lineages and regions. Although regression analysis of branch lengths suggested some degree of indel homoplasy among the most divergent lineages, there was little evidence for indel homoplasy within the Lecythidaceae. Likelihood ratio tests applied to the entire phylogenetic tree revealed a consistent pattern rejecting a molecular clock. Tajima's 1D and 2D tests revealed two taxa with consistent rate heterogeneity, one showing relatively more and one relatively fewer changes than other taxa. In general, nucleotide changes showed more evidence of rate heterogeneity than did indel changes. The rate of evolution was highly variable among the six cpDNA regions examined, with the trnS-trnG and trnH-psbA regions showing as much as 10% and 15% divergence within the Lecythidaceae. Deviations from rate homogeneity in the two taxa were constant across cpDNA regions, consistent with lineage-specific rates of evolution rather than cpDNA region-specific natural selection. There is no evidence that indels are more likely than nucleotide changes to experience homoplasy within the Lecythidaceae. These results support a neutral interpretation of cpDNA indel and nucleotide variation in population studies within species such as Corythophora alta.     

Repetitive Indel Markers within the ALMT1 Gene Conditioning Aluminium Tolerance in Wheat (Triticum aestivum L.)

ALMT1 gene encoding a membrane protein that facilitates an aluminium stimulated malate efflux has been characterised and mapped in wheat (Triticum aestivum L.). Here, we have identified molecular markers targeting insertion/deletion (indel) and SSR repeats within intron 3 region of the ALMT1 gene. Both the markers: ALMT1-SSR3a and ALMT1-SSR3b based on repetitive indels, exhibited complete cosegregation with Al tolerance, malate efflux, and a CAPS marker discriminating ALMT1-1 and ALMT1-2 alleles, in a doubled haploid population derived from Diamondbird (Al-tolerant)/Janz (Al-sensitive). A parental screen of 20 diverse wheat genotypes with repetitive indel markers indicated that six allele variants exist at the ALMT1SSR3 locus. Sequence analysis confirmed that these variations were due to indels, copy number of SSR repeats, and base substitution within SSR repeats. The higher level of variation in intron 3 suggests that this genomic region has been constrained by indels, SSR and single nucleotide polymorphisms. Results have proven that repetitive indel markers cosegregating with the Al tolerance locus will be useful for marker assisted selection and population and evolution studies.     

Development of genome-wide insertion and deletion markers for maize,based on next-generation sequencing data

Background

Insertions and deletions (indels) are the most abundant form of structural variation in all genomes. Indels have been increasingly recognized as an important source of molecular markers due to high-density occurrence, cost-effectiveness, and ease of genotyping. Coupled with developments in bioinformatics, next-generation sequencing (NGS) platforms enable the discovery of millions of indel polymorphisms by comparing the whole genome sequences of individuals within a species.

Results

A total of 1,973,746 unique indels were identified in 345 maize genomes, with an overall density of 958.79 indels/Mbp, and an average allele number of 2.76, ranging from 2 to 107. There were 264,214 indels with polymorphism information content (PIC) values greater than or equal to 0.5, accounting for 13.39 % of overall indels. Of these highly polymorphic indels, we designed primer pairs for 83,481 and 29,403 indels with major allele differences (i.e. the size difference between the most and second most frequent alleles) greater than or equal to 3 and 8 bp, respectively, based on the differing resolution capabilities of gel electrophoresis. The accuracy of our indel markers was experimentally validated, and among 100 indel markers, average accuracy was approximately 90 %. In addition, we also validated the polymorphism of the indel markers. Of 100 highly polymorphic indel markers, all had polymorphisms with average PIC values of 0.54.

Conclusions

The maize genome is rich in indel polymorphisms. Intriguingly, the level of polymorphism in genic regions of the maize genome was higher than that in intergenic regions. The polymorphic indel markers developed from this study may enhance the efficiency of genetic research and marker-assisted breeding in maize.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1797-5) contains supplementary material, which is available to authorized users.     

The pattern of insertion/deletion polymorphism in Arabidopsis thaliana

Little is known about variation of nucleotide insertion/deletions (indels) within species. In Arabidopsis thaliana, we investigated indel polymorphism patterns between two genome sequences and among 96 accessions at 1215 loci. Our study identified patterns in the variation of indel density, size, GC content and distribution, and a correlation between indels and substitutions. We found that the GC content in indel sequences was lower than that in non-indel sequences and that indels typically occur in regions with lower GC content. Patterns of indel frequency distribution among populations were more consistent with neutral expectation than substitution patterns. We also found that the local level of substitutions is positively correlated with indel density and negatively correlated with their distance to the closed indel, suggesting that indels play an important role in nucleotide variation.     

Microsatellite Allelic Homoplasy Due to Variable Flanking Sequences

Microsatellite DNA sequences have become the dominant source of nuclear genetic markers for most applications. It is important to investigate the basis of variation between alleles and to know if current assumptions about the mechanisms of microsatellite mutation (that is to say, variations involving simple changes in the number of repeat) are correct. We have characterized, by DNA sequencing, the human alleles of a new highly informative (CA)n repeat localized approximately 20 kb centromeric to the HLA-B gene. Although 12 alleles were identified based on conventional length criteria, sequencing of the alleles demonstrated that differences between alleles were found to be more complex than previously assumed: A high degree of microsatellite variability is due to variation in the region immediately flanking the repeat. These data indicate that the mutational process which generates polymorphism in this region has involved not only simple changes in the number of dinucleotide CA repeats but also perturbations in the nonrepeated 5′ and 3′ flanking sequences. Three families of alleles (not visible from the overall length of the alleles), with presumably separate evolutionary histories, exist and can yield to homoplasy of size. Effectively, we can observe alleles of the same size with different internal structures which are separated by a significant amount of variation. Although allelic homoplasy for noninterrupted microsatellite loci has been suggested between different species, it has not been unequivocally demonstrated within species. A strong association is noted between alleles defined at the sequence level and HLA-B alleles. The observation of several families of alleles at the population level provides information about the evolutionary history and mutation processes of microsatellites and may have implications for the use of these markers in phylogenetic, linkage disequilibrium studies, and gene mapping. Received: 14 May 1996 / Accepted: 9 September 1996     

Efficiency of microsatellite enrichment in<Emphasis Type="Italic">Prosopis chilensis</Emphasis> using magnetic capture

Microsatellites (i.e., simple sequence repeats [SSRs]) are highly variable genetic markers that are widely used at an intraspecific level in population genetic studies. Here we employed an enrichment strategy for microsatellite isolation by using microsatellite oligoprobes and magnetic capture of the fragments (Fischer and Bachmann, 1998) inProsopis chilensis (Mol.) Stuntz (Fabaceae). We analyzed the obtained level of enrichment by sequencing 120 enriched genomic fragments. A total of 521 SSR motives were detected. According to specific search criteria (SSR motifs ≥3 repeat units and ≥6 bp length), 95.8% of the clones contained SSR motifs. Of these, 7.8% showed homology to chloroplast sequences and 92.2% to nuclear sequences. When regarding only nuclear SSRs with 5 or more repeat units and a minimum length of 10 bp, the level of enrichment was 30.8%. A FASTA search against the European Molecular Biology Laboratory (EMBL) database univocally revealed 4 clones in transcribed regions, 102 clones in genomic regions with unknown function, and 9 clones in chloroplast regions. Among the loci with longer repeat units (≥10 bp, ≥5 repeat units), 3 were in transcribed regions and 65 were in other genomic regions. We discuss the applicability of these markers for population genetic studies.     

Causes of Size Homoplasy Among Chloroplast Microsatellites in Closely Related <Emphasis Type="Italic">Clusia</Emphasis> Species

Chloroplast DNA sequences and microsatellites are useful tools for phylogenetic as well as population genetic analyses of plants. Chloroplast microsatellites tend to be less variable than nuclear microsatellites and therefore they may not be as powerful as nuclear microsatellites for within-species population analysis. However, chloroplast microsatellites may be useful for phylogenetic analysis between closely related taxa when more conventional loci, such as ITS or chloroplast sequence data, are not variable enough to resolve phylogenetic relationships in all clades. To determine the limits of chloroplast microsatellites as tools in phylogenetic analyses, we need to understand their evolution. Thus, we examined and compared phylogenetic relationships of species within the genus Clusia, using both chloroplast sequence data and variation at seven chloroplast microsatellite loci. Neither ITS nor chloroplast sequences were variable enough to resolve relationships within some sections of the genus, yet chloroplast microsatellite loci were too variable to provide any useful phylogenetic information. Size homoplasy was apparent, caused by base substitutions within the microsatellite, base substitutions in the flanking regions, indels in the flanking regions, multiple microsatellites within a fragment, and forward/reverse mutations of repeat length resulting in microsatellites of identical base composition that were not identical by descent.     

Extensive size homoplasy at a microsatellite locus in the Japanese bumblebee, Bombus diversus

An extensive size homoplasy was found at microsatellite locus B11 of the bumblebee, Bombus diversus, in northern to central Honshu, Japan. A total of 16 alleles of different nucleotide sequences in five length morphs was obtained at B11 for this species. Of these alleles, five were 141 base pairs (bp) in length, five were 137 bp and four were 133 bp. Allele diversity in each length morph was high compared with previous studies. It is noteworthy that this extensive size homoplasy was found in a relatively small geographic area, in contrast to results from previous studies. Reconstruction of a median‐joining network revealed the complicated evolutionary process of the locus, involving insertion/deletion and point mutations. Preliminary estimation of the mutation rate of the B11 locus in B. diversus gives a value comparable to those estimated from experimental Drosophila populations. Effects of the extensive size homoplasy in population genetic studies is discussed.     

Microsatellites in Zea - variability,patterns of mutations,and use for evolutionary studies

To evaluate the performance of microsatellites or simple sequence repeats (SSRs) for evolutionary studies in Zea, 46 microsatellite loci originally derived from maize were applied to diverse arrays of populations that represent all the diploid species of Zea and 101 maize inbreds. Although null phenotypes and amplification of more than two alleles per plant were observed at modest rates, no practical obstacle was encountered for applying maize microsatellites to other Zea species. Sequencing of microsatellite alleles revealed complex patterns of mutation including frequent indels in the regions flanking microsatellite repeats. In one case, all variation at a microsatellite locus came from indels in the flanking region rather than in the repeat motif. Maize microsatellites show great variability within populations and provide a reliable means to measure intraspecific variation. Phylogeographic relationships of Zea populations were successfully reconstructed with good resolution using a genetic distance based on the infinite allele model, indicating that microsatellite loci are useful in evolutionary studies in Zea. Microsatellite loci show a principal division between tropical and temperate inbred lines, and group inbreds within these two broad germplasm groups in a manner that is largely consistent with their known pedigrees. Received: 10 February 2001 / Accepted: 21 May 2001