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.     

Characterization of chloroplast DNA microsatellites from Saccharum spp and related species

Microsatellites, or simple sequence repeats (SSRs), and their flanking regions in chloroplast genomes (plastomes) of some species of the family Poaceae were analyzed in silico to look for DNA sequence variations. Comparison of the complete chloroplast DNA sequences (cpDNAs) of sugarcane (Saccharum hybrid cv. SP-80-3280 and S. officinarum cv. NCo310) and related species, Agrostis stolonifera, Brachypodium distachyon, Hordeum vulgare subsp vulgare, Lolium perenne, Oryza nivara, O. sativa subsp indica, O. sativa subsp japonica, Sorghum bicolor, Triticum aestivum, Zea mays, and Z. mays cv. B73, allowed us to examine the organization of chloroplast SSRs (cpSSRs) in genic and intergenic regions. We identified 204 cpSSRs in the sugarcane cpDNA; 22.5% were in genic regions. The ndh, rps, trn, and rpl gene clusters of the chloroplasts had the most repeats. Mononucleotide repeats were the most abundant cpSSRs in these species; however, di-, tri-, tetra-, penta-, and hexanucleotide repeats were also identified. Many base substitutions and deletions/insertions were identified in the cpSSR loci and their flanking regions. Multiple alignments of all cpSSR sequences of Poaceae species made identification of nucleotide variability possible; repeat motifs are not uniformly distributed across the Poaceae plastomes, but are mostly confined to intergenic regions. Phylogeny was determined by maximum parsimony and neighbor-joining inference methods. The cpSSRs of these species were found to be polymorphic. It appears that individual cpSSRs in the Poaceae are stable, at least over short periods of evolutionary time. We conclude that the plastome database can be exploited for phylogenetic analysis and biotechnological development.     

Evolution of the primate beta-globin gene region. High rate of variation in CpG dinucleotides and in short repeated sequences between man and chimpanzee

A 5500 base-pair fragment including the beta-globin gene downstream from codon 122 and about 4000 base-pairs of its 5' flanking sequence was cloned from chimpanzee DNA and thoroughly sequenced before being compared with the corresponding human sequence: 88 point differences (83 substitutions and 5 deletions or insertions of 1 base-pair) were detected as well as seven more important deletion/insertion events. These changes occur preferentially in two kinds of structure. First, 40% of the CpG dinucleotides present in either human or chimpanzee sequences are affected by nucleotide variations. This corresponds to a divergence level considerably higher than that expected. Second, most short repeated sequences found in the 5' extragenic sequence are involved in mutational events (amplification or contraction of the number of basic motifs as well as point substitutions or deletions/insertions of 1 base-pair). Considering the very low level of nucleotide sequence divergence between these two closely related species, our data provide direct evidence for CpG and tandem array instability.     

Isolation, cloning and characterisation of motifs containing (GA/TC)n repeats isolated from vetch, Vicia bithynica

Microsatellites are widely distributed in plant genomes and comprise unstable regions that undergo mutational changes at rates much greater than that observed for non-repetitive sequences. They demonstrate intrinsic genetic instability, manifested as frequent length changes due to insertions or deletions of repeat units. Detailed analysis of 1600 clones containing genomic sequences of Vicia bithynica revealed the presence of microsatellite repeats in its genome. Based on the screening of a partial DNA library of plasmids, 13 clones harbouring (GA/TC)n tracts of various lengths of repeated motif were identified for further analysis of their internal sequence organization. Sequence analyses revealed the precise length, number of repeats, interruptions within tracts, as well as sequence composition flanking the repeat motifs. Representative plasmids containing different lengths of (GA/TC)n embedded in their original flanking sequence were used to investigate the genetic stability of the repeats. In the study presented herein, we employed a well characterised and tractable bacterial genetic system. Recultivations of Escherichia coli harbouring plasmids containing (GA/TC)n inserts demonstrated that the genetic instability of (GA/TC)n microsatellites depends highly on their length (number of repeats). These observations are in agreement with similar studies performed on repetitive sequences from humans and other organisms.     

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.     

High incidence of nonslippage mechanisms generating variability and complexity in eurasian badger microsatellites

The use of microsatellites in population genetics is hindered by a lack of understanding of the pattern and origin of mutations, the need to develop more specific and better computational models, and a paucity of information about specific taxa and loci. We analyzed between 4 and 10 allele sequences from 10 different microsatellites in Eurasian badgers in order to determine the compliance of the sequences with stepwise mutation models and the origin of that variability which cannot be detected through standard genotyping procedures. All microsatellite loci exhibited imperfections and/or substitutions and indels in the flanking region, as well as additions or deletions of repeat units. Our data set of sequences showed a higher number of imperfect repeats than other published badger and carnivore sequences. This could be attributed to the process of loci isolation because when genetic variability is low, researchers may be more likely to use imperfect loci if these are variable in the population being studied. Locus Mel15 had 2 repetitive arrays: one was part of a polypyrimidine region of a carnivoran short interspersed nuclear element (CAN-SINE) and the other was located in an A-rich region typical of these insertions. In spite of this complexity, heterozygosity was correlated with the maximum number of repeats. Thus, although new theoretical models are being evolved to cover complex patterns of microsatellite mutation, sequencing electromorphs is needed to identify microsatellites or portions of them whose evolution can be modeled under simple models.     

Identification and molecular characterization of a large insertion within the repetitive domain of a high-molecular-weight glutenin subunit gene from hexaploid wheat

High-molecular-weight (HMW) glutenin subunits are a particular class of wheat endosperm proteins containing a large repetitive domain flanked by two short N- and C-terminal non-repetitive regions. Deletions and insertions within the central repetitive domain has been suggested to be mainly responsible for the length variations observed for this class of proteins. Nucleotide sequence comparison of a number of HMW glutenin genes allowed the identification of small insertions or deletions within the repetitive domain. However, only indirect evidence has been produced which suggests the occurrence of substantial insertions or deletions within this region when a large variation in molecular size is present between different HMW glutenin subunits. This paper represents the first report on the molecular characterization of an unusually large insertion within the repetitive domain of a functional HMW glutenin gene. This gene is located at the Glu-D1 locus of a hexaploid wheat genotype and contains an insertion of 561 base pairs that codes for 187 amino acids corresponding to the repetitive domain of a HMW glutenin subunit encoded at the same locus. The precise location of the insertion has been identified and the molecular processes underlying such mutational events are discussed.     

Evolution of a B2 tagged sequence from a long-range repeat family in the genus Mus

A long-range repeat family of more than 50 kb repeat size is clustered in Chromosomes (Chr) 1 of Mus musculus and M. spretus. In M. musculus this long-range repeat family shows considerable variation of copy-number frequency and contains coding regions for at least two genes. In an intron of a gene, which is part of the repeat, a B2 small interspersed repetitive element (SINE) is inserted at identical positions. The B2 element is present in all copies of the long-range repeat family; it was presumably a component of the ancestral single-copy precursor sequence that gave rise by amplification to the repeat family. Copies of the long-range repeat family vary with respect to the number of TAAA tandem repeats in the A-rich 3 end region of the B2 element. As inferred from polymerase chain reaction (PCR) data, presence and frequency of repeat number variants in the (TAAA)_n block are strain and species specific. The B2 element and its flanking regions were sequenced from two copies of the long-range repeat family. Sequence divergence between the two copies (only non-CG base substitutions and deletions/insertions) was determined to be 2.6%. Based on the drift rate in human Alu elements and a correction for the higher drift rates in rodents, and estimate for the divergence time of 1.7 million years was calculated. Since the long-range repeat family is present in M. musculus and M. spretus, it must have evolved by amplification before the separation of the two species about 1–4 million years ago.     

Comparative sequencing of a microsatellite locus reveals size homoplasy within and between european oak species (<Emphasis Type="Italic">Quercus</Emphasis> spp.)

Microsatellites (simple sequence repeats [SSRs]) are highly variable molecular markers that are a rich and readily assayed source of variation for population genetic studies. Cross-amplification between closely related species is possible when there are no (or few) sequence differences in the primer binding sites. The occurrence of nonhomologous fragments of the same size (size homoplasy) is a contraint of microsatellites. Size homoplasy can be caused by insertions/deletions (indels) in SSR flanking regions. We found that size variation in locus ssrQZAG9 is due to different repeat numbers of the SSR motifs but also to indels in SSR flanking regions. Indels were found within species belonging to sectionsRobur andCerris of genusQuercus and also between species of the 2 sections. In sectionRobur (Quercis robur L.,Quercus petraea [Matt.] Liebl.,Quercus pubescens Willd.), we detected rare alleles with an indel of 57 bp or 62 bp followed by a smaller indel of 12 bp in the SSR flanking regions. These alleles show a size range overlapping with that of alleles amplified inQuercus cerris L. (sectionCerris). Multiple alignments with sequences of sectionRobur revealed the same SSR repeat motif but multiple indels in SSR flanking regions inQ. cerris. We discuss the effects of size homoplasy of SSR loci for the study of interspecific gene flow and on estimates of population differentiation.     

Mutation and evolution of microsatellite loci in Neurospora

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.     

Comparative analysis of the intergenic spacer regions and population structure of the species complex of the pathogenic yeast Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic basidiomycete responsible for the high incidence of cryptococcosis in patients with AIDS and in other immune-compromised individuals. This study, which focused on the molecular structure and genetic variability of the two varieties in the C. neoformans and Cryptococcus gattii species complex, employed sequence analysis of the intergenic spacer regions, IGSI and IGSII. The IGS region is the most rapidly evolving region of the rDNA families. The IGSI displayed the most genetic variability represented by nucleotide base substitutions and the presence of long insertions/deletions (indels). In contrast, the IGSII region exhibited less heterogeneity and the indels were not as extensive as those displayed in the IGSI region. Both intergenic spacers contained short, interspersed repeat motifs, which can be related to length polymorphisms observed between sequences. Phylogenetic analysis undertaken in the IGSI, IGSII and IGSI +5S rRNA + IGSII regions revealed the presence of six major phylogenetic lineages, some of which segregated into subgroups. The major lineages are represented by genotypes 1 (C. neoformans var. grubii), genotype 2 (C. neoformans var. neoformans), and genotypes 3, 4, 5 and 6 represented by C. gattii. Genotype 6 is a newly described IGS genotypic group within the C. neoformans species complex. With the inclusion of IGS subgenotypic groups, our sequence analysis distinguished 12 different lineages. Sequencing of clones, which was performed to determine the presence of multiple alleles at the IGS locus in several hybrid strains, yielded a single IGS sequence type per isolate, thus suggesting that the selected group of cloned strains was mono-allelic at this locus. IGS sequence analyses proved to be a powerful technique for the delineation of the varieties of C. neoformans and C. gattii at genotypic and subgenotypic levels.     

Complex mutation at a microsatellite locus in sturgeons: Acipenser sinensis, A. schrenckii, A. gueldenstaedtii and A. baerii

Cross‐species amplifications of microsatellite locus Spl‐106, which was originally screened from the genome of shovelnose sturgeon (Scaphirhynchus platorynchus) with a perfect TAGA repeat motif, were carried out in four other species of the genera Acipenser. A total of 34 polymerase chain reaction (PCR) products representing 16 different alleles of this locus was sequenced. Sequence analysis results showed that besides the number changes of repeat units, many mutational events, such as single‐base substitutions and various insertion/deletion (indels) occurred not only at species level but also at individual level, even among the different alleles within the same individual. The repeat motifs varied from perfect (TAGA)n array to perfect compound (TAAA)m (GAAA)n and perfect or imperfect compound (TAAA)m (TAGA)n (TAAA)x arrays in different species and different individuals. The evolution dynamics of this locus in sturgeons was inferred in that it may evolve from a single perfect to different perfect or imperfect compounds.     

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.     

Characterisation and analysis of new HMW-glutenin alleles encoded by the Glu-R1 locus of Secale cereale

This work reports the molecular characterisation of new alleles of the previously reported Glu-R1 locus. Wheat lines carrying the chromosome substitution 1R(1D), rye cultivars and related wild species were analysed. Five new x-type and four y-type Glu-R1 glutenin subunits were isolated and characterised. The coding region of the sequences shows the typical structure of the HMW glutenin genes previously described in wheat, with the N and C-terminal domains flanking the central repetitive region. Tri-, hexa- and nona-peptides found in the central repetitive region of wheat glutenin genes were also present in the rye genes. Duplications and deletions of these motifs are responsible for allelic variation at the Glu-R1 locus. Orthologous genes (from different genomes) were more closely related than paralogous genes (x- and y-type), supporting the hypothesis of gene duplication before Triticeae speciation. Differences in the number and position of cysteine residues identified alleles which in wheat are associated with good dough quality. SDS proteins encoded by some characterised alleles were presumptively identified.     

The unstable CCTG repeat responsible for myotonic dystrophy type 2 originates from an AluSx element insertion into an early primate genome

Myotonic dystrophy type 2 (DM2) is a subtype of the myotonic dystrophies, caused by expansion of a tetranucleotide CCTG repeat in intron 1 of the zinc finger protein 9 (ZNF9) gene. The expansions are extremely unstable and variable, ranging from 75-11,000 CCTG repeats. This unprecedented repeat size and somatic heterogeneity make molecular diagnosis of DM2 difficult, and yield variable clinical phenotypes. To better understand the mutational origin and instability of the ZNF9 CCTG repeat, we analyzed the repeat configuration and flanking regions in 26 primate species. The 3'-end of an AluSx element, flanked by target site duplications (5'-ACTRCCAR-3'or 5'-ACTRCCARTTA-3'), followed the CCTG repeat, suggesting that the repeat was originally derived from the Alu element insertion. In addition, our results revealed lineage-specific repetitive motifs: pyrimidine (CT)-rich repeat motifs in New World monkeys, dinucleotide (TG) repeat motifs in Old World monkeys and gibbons, and dinucleotide (TG) and tetranucleotide (TCTG and/or CCTG) repeat motifs in great apes and humans. Moreover, these di- and tetra-nucleotide repeat motifs arose from the poly (A) tail of the AluSx element, and evolved into unstable CCTG repeats during primate evolution. Alu elements are known to be the source of microsatellite repeats responsible for two other repeat expansion disorders: Friedreich ataxia and spinocerebellar ataxia type 10. Taken together, these findings raise questions as to the mechanism(s) by which Alu-mediated repeats developed into the large, extremely unstable expansions common to these three disorders.     

Interdependent MHC-DRB exon-plus-intron evolution in artiodactyls

Exon 2 sequences of an expressed MHC-DRB locus from sheep were examined for polymorphisms in both the antigen-binding regions and the adjacent intronic mixed simple tandem repeat. Twenty-one novel exon 2 Ovar-DRB alleles were identified. Short nucleotide motifs are extensively shared between certain exon 2 regions of Ovar-DRB alleles. The simple repeat variations, the number of different amino acids at usually polymorphic sites, and the number of silent substitutions were reduced in the intraspecies analyses of sheep DRB sequences, compared with those of cattle and goats. It was paradoxical that the abundance of different sheep alleles was similar to that of cattle and goats. This paradox may be explained by postulating a relatively small number of "ancient" alleles, with the present-day Ovar-DRB alleles being generated by reciprocal exchange of nucleotide motifs. At the antigen-binding sites, new combinations of amino acids were maintained in Ovar-DRB alleles by strong positive selection. In sheep--and less pronounced in goats and cattle--the DRB alleles can be divided into two groups. In one group, silent substitutions are increased when compared with the other. This suggests separate evolutionary pathways for certain groups of DRB alleles within a species. The simple repetitive sequences are also discussed with respect to the evolution of DRB alleles.      

DNA indels in coding regions reveal selective constraints on protein evolution in the human lineage

Background  

Insertions and deletions of DNA segments (indels) are together with substitutions the major mutational processes that generate genetic variation. Here we focus on recent DNA insertions and deletions in protein coding regions of the human genome to investigate selective constraints on indels in protein evolution.     

Characterization of the Sperm Molecule Bindin in the Sea Urchin Genus <Emphasis Type="Italic">Paracentrotus</Emphasis>

Bindin is a sea urchin gamete-recognition protein that plays an essential role in the specificity of egg–sperm interactions and thus may be evolving under sexual selection and be related to speciation. Bindin has been found to evolve under strong selection in some sea urchin genera and neutrally in others. In this study, we characterized bindin in the two extant species of the genus Paracentrotus: P. lividus from the Atlanto-Mediterranean region and P. gaimardi from Brazil. The structure of the bindin molecule in Paracentrotus is similar to that of other genera studied thus far, consisting of a conserved core flanked by two variable regions and an intron of variable length located at the same conserved position as in other genera. Polymorphism in P. lividus is caused mainly by point substitutions and insertions/deletions, and length variations are caused mainly by the number of repeated motifs in the flanking regions. There is no evidence of recombination. Positive selection is acting on amino acid sites located in two regions flanking the conserved core.     

CpG/CpNpG motifs in the coding region are preferred sites for mutagenesis in the breast cancer susceptibility genes

The range of BRCA1/BRCA2 gene mutations is diverse and the mechanism accounting for this heterogeneity is obscure. To gain insight into the endogenous mutational mechanisms involved, we evaluated the association of specific sequences (i.e. CpG/CpNpG motifs, homonucleotides, short repeats) and mutations within the genes. We classified 1337 published mutations in BRCA1 (1765 BRCA2 mutations) for each specific sequence, and employed computer simulation combined with mathematical calculations to estimate the true underlying tendency of mutation occurrence. Interestingly, we found no mutational bias to homonucleotides and repeats in deletions/insertions and substitutions but striking bias to CpG/CpNpG in substitutions in both genes. This suggests that methylation-dependent DNA alterations would be a major mechanism for mutagenesis.     

Correlations among oligonucleotide repeats,nucleotide substitutions,and insertion–deletion mutations in chloroplast genomes of plant family Malvaceae

The co‐occurrence of mutational events including substitutions and insertions–deletions (InDels) with oligonucleotide repeats has previously been reported for a limited number of prokaryotic, eukaryotic, and organelle genomes. In this study, the correlations among these mutational events in chloroplast genomes of species in the eudicot family Malvaceae were investigated. This study also reported chloroplast genome sequences of Hibiscus mutabilis, Malva parviflora, and Malvastrum coromandelianum. These three genomes and 16 other publicly available chloroplast genomes from 12 genera of Malvaceae were used to calculate the correlation coefficients among the mutational events at family, subfamily, and genus levels. In these comparisons, chloroplast genomes were pairwise aligned to record the substitutions and the InDels in mutually exclusive, 250nucleotide long bins. Taking one among the two genomes as a reference, the coordinate positions of oligonucleotide repeats in the reference genome were recorded. The extent of correlations among repeats, substitutions, and InDels was calculated and categorized as follows: very weak (0.1–0.19), weak (0.20–0.29), moderate (0.30–0.39), and strong (0.4–0.69). The extent of correlations ranged 0.201–0.6 between “InDels and single‐nucleotide polymorphism (SNP)”, 0.182–0.513 between “InDels and repeat” and 0.055–0.403 between “SNPs and repeats”. At family‐ and subfamily‐level comparisons, 88%–96% of the repeats showed co‐occurrence with SNPs, whereas at the genus level, 23%–86% of the repeats co‐occurred with SNPs in same bins. Our findings support the previous hypothesis suggesting the use of oligonucleotide repeats as a proxy for finding the mutational hotspots.