Microsatellites: consensus and controversy

Detecting microsatellites in genome data: variance in definitions and bioinformatic approaches cause systematic bias

Microsatellites are currently one of the most commonly used genetic markers. The application of bioinformatic tools has become common practice in the study of these short tandem repeats (STR). However, in silico studies can suffer from study bias. Using a meta-analysis on microsatellite distribution in yeast we show that estimates of numbers of repeats reported by different studies can differ in the order of several magnitudes, even within a single genome. These differences arise because varying definitions of microsatellites, spanning repeat size, array length and array composition, are used in different search paradigms, with minimum array length being the main influencing factor. Structural differences in the implemented search algorithm additionally contribute to variation in the number of repeats detected. We suggest that for future studies a consistent approach to STR searches is adopted in order to improve the power of intra- and interspecific comparisons.

DNA polymerase kappa produces interrupted mutations and displays polar pausing within mononucleotide microsatellite sequences

Microsatellites are ubiquitously present in eukaryotic genomes and are implicated as positive factors in evolution. At the nucleotide level, microsatellites undergo slippage events that alter allele length and base changes that interrupt the repetitive tract. We examined DNA polymerase errors within a [T]₁₁ microsatellite using an in vitro assay that preferentially detects mutations other than unit changes. We observed that human DNA polymerase kappa (Pol κ) inserts dGMP and dCMP within the [T]₁₁ mononucleotide repeat, producing an interrupted 12-bp allele. Polymerase β produced such interruptions at a lower frequency. These data demonstrate that DNA polymerases are capable of directly producing base interruptions within microsatellites. At the molecular level, expanded microsatellites have been implicated in DNA replication fork stalling. Using an in vitro primer extension assay, we observed sequence-specific synthesis termination by DNA polymerases within mononucleotides. Quantitatively, intense, polar pausing was observed for both pol κ and polymerase α-primase within a [T]₁₁ allele. A mechanism is proposed in which pausing results from DNA bending within the duplex stem of the nascent DNA. Our data support the concept of a microsatellite life-cycle, and are consistent with the models in which DNA sequence or secondary structures contributes to non-uniform rates of replication fork progression.

Using principle component analysis to compare genetic diversity across polyploidy levels within plant complexes: an example from British Restharrows (Ononis spinosa and Ononis repens)

The investigation of genetic diversity between related plant populations which differ in ploidy levels is problematic, with common statistical methods developed for diploids being inappropriate for polyploid species. Studies into gene flow in such complexes are critical and can shed light on the mechanisms that generate and maintain populations of different polyploidy levels. We have investigated the use of principle component (PCO) analysis as one approach to elucidate population structure within British Restharrows (Leguminsoae, Ononis spp). Restharrows were common agricultural weed species until the advent of mechanical ploughing and both diploid (2n=2x=30; O. spinosa and O. intermedia) and tetraploid (2n=4x=60; O. repens and O. maritima) taxa exist. Patterns of genetic diversity were investigated among British Restharrows using 10 microsatellite loci with 21 Restharrow populations analysed (411 individual plants) from Central and Eastern Britain. PCO analysis revealed clear genetic differentiation of the sampled plants into two groups, corresponding to O. spinosa/O. intermedia (diploid) and O. repens/O. maritima (tetraploid) plants. Evidence of genetic differentiation by distance was also revealed for O. repens/O. maritima, but not for O. spinosa/O. intermedia. The data suggest the presence of strong reproductive barriers between diploid and tetraploid Restharrows in Britain, but not within ploidy levels. This genetic isolation between ploidy levels is confirmed by a detailed analysis of a sympatric site (Harton Down Hill). These results demonstrate that PCA analysis is a suitable general tool for comparing related species of different ploidy levels.

Spatial analysis of nuclear and cytoplasmic DNA diversity in wild sea beet (Beta vulgaris ssp. maritima) populations: do marine currents shape the genetic structure?

Patterns of seed dispersal in the wild sea beet (Beta vulgaris ssp. maritima) are predicted to be influenced by marine currents because populations are widely distributed along the European Atlantic coast. We investigated the potential influence of marine currents on the pattern of spatial genetic structuring in natural populations of sea beet. Populations were located along the French coasts of the Anglo-Norman gulf that features peculiar marine currents in the Channel. Thirty-three populations were sampled, among which 23 were continental and 10 were insular populations located in Jersey, Guernsey and Chausey, for a total of 1224 plants genotyped. To validate the coastal topography influence and the possibility of marine current orientated gene flow on the genetic features of sea beet populations, we assessed patterns of genetic structuring of cytoplasmic and nuclear diversity by: (i) searching for an isolation-by-distance (IBD) pattern using spatial autocorrelation tools; (ii) using the Monmonier algorithm to identify genetic boundaries in the area studied; and (iii) performing assignment tests that are based on multilocus genotype information to ascertain population membership of individuals. Our results showed a highly contrasted cytoplasmic and nuclear genetic differentiation and highlighted the peculiar situation of island populations. Beyond a classical isolation-by-distance due to short-range dispersal, genetic barriers fitting the orientation of marine currents were clearly identified. This suggests the occurrence of long-distance seed dispersal events and an asymmetrical gene flow separating the eastern and western part of the Anglo-Norman gulf.

Frequency, type, and distribution of EST-SSRs from three genotypes of Lolium perenne, and their conservation across orthologous sequences of Festuca arundinacea, Brachypodium distachyon, and Oryza sativa

BACKGROUND

Simple sequence repeat (SSR) markers are highly informative and widely used for genetic and breeding studies in several plant species. They are used for cultivar identification, variety protection, as anchor markers in genetic mapping, and in marker-assisted breeding. Currently, a limited number of SSR markers are publicly available for perennial ryegrass (Lolium perenne). We report on the exploitation of a comprehensive EST collection in L. perenne for SSR identification. The objectives of this study were 1) to analyse the frequency, type, and distribution of SSR motifs in ESTs derived from three genotypes of L. perenne, 2) to perform a comparative analysis of SSR motif polymorphisms between allelic sequences, 3) to conduct a comparative analysis of SSR motif polymorphisms between orthologous sequences of L. perenne, Festuca arundinacea, Brachypodium distachyon, and O. sativa, 4) to identify functionally associated EST-SSR markers for application in comparative genomics and breeding.

RESULTS

From 25,744 ESTs, representing 8.53 megabases of nucleotide information from three genotypes of L. perenne, 1,458 ESTs (5.7%) contained one or more SSRs. Of these SSRs, 955 (3.7%) were non-redundant. Tri-nucleotide repeats were the most abundant type of repeats followed by di- and tetra-nucleotide repeats. The EST-SSRs from the three genotypes were analysed for allelic- and/or genotypic SSR motif polymorphisms. Most of the SSR motifs (97.7%) showed no polymorphisms, whereas 22 EST-SSRs showed allelic- and/or genotypic polymorphisms. All polymorphisms identified were changes in the number of repeat units. Comparative analysis of the L. perenne EST-SSRs with sequences of Festuca arundinacea, Brachypodium distachyon, and Oryza sativa identified 19 clusters of orthologous sequences between these four species. Analysis of the clusters showed that the SSR motif generally is conserved in the closely related species F. arundinacea, but often differs in length of the SSR motif. In contrast, SSR motifs are often lost in the more distant related species B. distachyon and O. sativa.

CONCLUSION

The results indicate that the L. perenne EST-SSR markers are a valuable resource for genetic mapping, as well as evaluation of co-location between QTLs and functionally associated markers.

A set of multiplex panels of microsatellite markers for rapid molecular characterization of rice accessions

BACKGROUND

This study aimed to analyze the efficiency of three new microsatellite multiplex panels, which were designed to evaluate a total of 16 loci of the rice genome, based on single PCR reactions of each panel. A sample of 548 accessions of traditional upland rice landraces collected in Brazil in the last 25 years was genotyped, a database of allelic frequencies was established, estimates of genetic parameters were performed and analysis of genetic structure of the collection was developed.

RESULTS

The three panels yielded a combined matching probability of 6.4 x 10-21, polymorphism information content (PIC) of 0.637, and a combined power of exclusion greater than 99.99%. A few samples presented a genetic background of indica rice. The 16 SSR loci produced a total of 229 alleles. Gene diversity values averaged 0.667, and PIC values averaged 0.637. Genetic structure analysis of the collection using a Bayesian approach detected three possible major clusters, with an overall FST value of 0.177. Important inputs on the knowledge about upland rice germplasm differentiations which happened in Brazil in the last few centuries were also achieved and are discussed.

CONCLUSION

The three multiplex panels described here represent a powerful tool for rice genetic analysis, offering a rapid and efficient option for rice germplasm characterization. The data gathered demonstrates the feasibility of genotyping extensive germplasm collections using panels of multiplexed microsatellite markers. It contributes to the advancement of research on large scale characterization and management of germplasm banks, as well as identification, protection and assessments of genetic relationship of rice germplasm.

RAPD with microsatellite as a tool for differentiation of Candida genus yeasts isolated in brewing

Fifteen wild yeast strains were isolated in two factories of a lager brewing company in Poland. Their identification with API 32C system showed mainly the presence of Candida sake species (7/15). To differentiate the isolates, randomly amplified polymorphic DNA (RAPD) with (GTG)(5), (GAC)(5), (GACA)(4) microsatellite primers and M13 core sequence (5'-GAG GGT GGC GGT TCT-3') were chosen. The results of patterns similarity are presented as dendrograms for each RAPD analysis and for overall patterns. On the overall patterns, all isolates identified as C. sake, except Strain No. 1, were regrouped in one cluster. Collection strain C. sake CBS 617 was similar in 46% to the cluster with six isolates (Strain Nos. 3, 6, 8, 11, 13, 14). The second reference strain C. sake CBS 159 and the Strain No. 1 were regrouped with other Candida species (collection strains) showing, respectively, only 20% and 42% of similarity to other C. sake strains. The similarity based on the overall dendrogram between isolate Nos. 3, 6, 8, 11, 13, 14 and C. sake CBS 617 was 49%. Between those strains and other Candida, the similarity was only 37%.

Detecting microsatellites within genomes: significant variation among algorithms

Background

Microsatellites are short, tandemly-repeated DNA sequences which are widely distributed among genomes. Their structure, role and evolution can be analyzed based on exhaustive extraction from sequenced genomes. Several dedicated algorithms have been developed for this purpose. Here, we compared the detection efficiency of five of them (TRF, Mreps, Sputnik, STAR, and RepeatMasker).

Results

Our analysis was first conducted on the human X chromosome, and microsatellite distributions were characterized by microsatellite number, length, and divergence from a pure motif. The algorithms work with user-defined parameters, and we demonstrate that the parameter values chosen can strongly influence microsatellite distributions. The five algorithms were then compared by fixing parameters settings, and the analysis was extended to three other genomes (Saccharomyces cerevisiae, Neurospora crassa and Drosophila melanogaster) spanning a wide range of size and structure. Significant differences for all characteristics of microsatellites were observed among algorithms, but not among genomes, for both perfect and imperfect microsatellites. Striking differences were detected for short microsatellites (below 20 bp), regardless of motif.

Conclusion

Since the algorithm used strongly influences empirical distributions, studies analyzing microsatellite evolution based on a comparison between empirical and theoretical size distributions should therefore be considered with caution. We also discuss why a typological definition of microsatellites limits our capacity to capture their genomic distributions.

IMEx: Imperfect Microsatellite Extractor

MOTIVATION

Microsatellites, also known as simple sequence repeats, are the tandem repeats of nucleotide motifs of the size 1-6 bp found in every genome known so far. Their importance in genomes is well known. Microsatellites are associated with various disease genes, have been used as molecular markers in linkage analysis and DNA fingerprinting studies, and also seem to play an important role in the genome evolution. Therefore, it is of importance to study distribution, enrichment and polymorphism of microsatellites in the genomes of interest. For this, the prerequisite is the availability of a computational tool for extraction of microsatellites (perfect as well as imperfect) and their related information from whole genome sequences. Examination of available tools revealed certain lacunae in them and prompted us to develop a new tool.

RESULTS

In order to efficiently screen genome sequences for microsatellites (perfect as well as imperfect), we developed a new tool called IMEx (Imperfect Microsatellite Extractor). IMEx uses simple string-matching algorithm with sliding window approach to screen DNA sequences for microsatellites and reports the motif, copy number, genomic location, nearby genes, mutational events and many other features useful for in-depth studies. IMEx is more sensitive, efficient and useful than the available widely used tools. IMEx is available in the form of a stand-alone program as well as in the form of a web-server.

AVAILABILITY

A World Wide Web server and the stand-alone program are available for free access at http://203.197.254.154/IMEX/ or http://www.cdfd.org.in/imex.

Profiling and authentication of human cell lines using short tandem repeat (STR) loci: Report from the National Cell Bank of Iran

Assurance of cell line homogeneity and capability of cell contamination detection are among the most essential steps of cell based research. Due to high discriminatory efficiency, low cost and reliability, analysis of short tandem repeats (STR) has been introduced as a method of choice for human cell line authentication. In the present study 13 Combined DNA Index System (CODIS) based STRs along with the gender determination (Amelogenin) gene were utilized to establish a reproducible approach for the authentication of 100 human cell lines deposited in the National Cell Bank of Iran (NCBI), using the polymerase chain reaction (PCR) method. PCR products were subsequently analyzed by polyacrylamide gel electrophoresis (PAGE) and visualized by silver staining followed by gel documentation and software analysis. STR profiles obtained were compared with those of the American Type Culture Collection (ATCC) and the Japanese Collection of Research Bioresource (JCRB) as STR references. We detected 18.8% cross contamination among the NCBI human cell lines. To our knowledge, this is the first report of authentication of human cell lines using the 13 CODIS core STRs combined with Amelogenin.

The rise, fall and renaissance of microsatellites in eukaryotic genomes

Microsatellites are among the most versatile of genetic markers, being used in an impressive number of biological applications. However, the evolutionary dynamics of these markers remain a source of contention. Almost 20 years after the discovery of these ubiquitous simple sequences, new genomic data are clarifying our understanding of the structure, distribution and variability of microsatellites in genomes, especially for the eukaryotes. While these new data provide a great deal of descriptive information about the nature and abundance of microsatellite sequences within eukaryotic genomes, there have been few attempts to synthesise this information to develop a global concept of evolution. This review provides an up-to-date account of the mutational processes, biases and constraints believed to be involved in the evolution of microsatellites, particularly with respect to the creation and degeneration of microsatellites, which we assert may be broadly viewed as a life cycle. In addition, we identify areas of contention that require further research and propose some possible directions for future investigation.

Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers

Recent improvements in genetic analysis and genotyping methods have resulted in a rapid expansion of the power of molecular markers to address ecological questions. Microsatellites have emerged as the most popular and versatile marker type for ecological applications. The rise of commercial services that can isolate microsatellites for new study species and genotype samples at reasonable prices presents ecologists with the unprecedented ability to employ genetic approaches without heavy investment in specialized equipment. Nevertheless, the lack of accessible, synthesized information on the practicalities and pitfalls of using genetic tools impedes ecologists' ability to make informed decisions on using molecular approaches and creates the risk that some will use microsatellites without understanding the steps needed to evaluate the quality of a genetic data set. The first goal of this synthesis is to provide an overview of the strengths and limitations of microsatellite markers and the risks, cost and time requirements of isolating and using microsatellites with the aid of commercial services. The second goal is to encourage the use and consistent reporting of thorough marker screening to ensure high quality data. To that end, we present a multistep screening process to evaluate candidate loci for inclusion in a genetic study that is broadly targeted to both novice and experienced geneticists alike.

Microsatellite analysis reveals genetic structure of Leishmania tropica

The current rapid spread of leishmaniases caused by Leishmania tropica and the complexity of its clinical spectrum call for this parasite's epidemiological and evolutionary investigation. Evaluation of its population structure by isoenzyme electrophoresis and previous molecular biological analysis has proved difficult. In this study, we used 21 microsatellite loci to type 117 strains from different African and Asian locations. Eighty-one different genotypes were found. A genetic bottleneck supported by a gradient in the number of alleles and consistent with the geographical structure of the Middle East suggests an African origin of this species. A Bayesian approach identified 10 genetic clusters that correlated predominantly with geographical origin. The strains in the 'Asia' cluster form a very heterogeneous sub-population, with a varied but inter-related genotype that is geographically very widely dispersed and consistent with anthroponotic transmission of the parasite. The other nine clusters were more homogenous. The propagation of L. tropica appears to be predominantly clonal. In Africa and the Middle East, anthroponotic and zoonotic systems of distribution may contribute to the development of overlapping, genetically distinct populations of L. tropica.

Genetic integrity of somaclonal variants in tea (Camellia sinensis (L.) O Kuntze) as revealed by inter simple sequence repeats

Adoption of inter simple sequence repeats (ISSR) technique to analyze the genetic variability of somatic embryo derived tea plants was evaluated. Morphological characterisation of the field grown plants revealed no identical character aligning with the parent, UPASI-10. Out of 40 primers, 15 exhibited concurrent polymorphism were selected for the study. Genetic variability of somaclones derived from single line cotyledonary culture ranged from 33.0 to 55.0%. A unique fragment of 1.2Kb was visible in majority of the accessions whereas the fragments below the length of 0.6Kb were noticed only in 50% of the variants. Out of 120 interactions attempted using Pearson's coefficient correlation, only 9.2% of somaclones exhibited significant similarity at genetic level. Dendrogram constructed based on simple matching coefficient revealed a distance of 2.257-3.317 between the final clusters. This strengthens the existence of wide genetic variation among the somaclones.

TRAP: automated classification, quantification and annotation of tandemly repeated sequences

TRAP, the Tandem Repeats Analysis Program, is a Perl program that provides a unified set of analyses for the selection, classification, quantification and automated annotation of tandemly repeated sequences. TRAP uses the results of the Tandem Repeats Finder program to perform a global analysis of the satellite content of DNA sequences, permitting researchers to easily assess the tandem repeat content for both individual sequences and whole genomes. The results can be generated in convenient formats such as HTML and comma-separated values. TRAP can also be used to automatically generate annotation data in the format of feature table and GFF files.

A molecular identification system for grasses: a novel technology for forensic botany

Our present inability to rapidly, accurately and cost-effectively identify trace botanical evidence remains the major impediment to the routine application of forensic botany. Grasses are amongst the most likely plant species encountered as forensic trace evidence and have the potential to provide links between crime scenes and individuals or other vital crime scene information. We are designing a molecular DNA-based identification system for grasses consisting of several PCR assays that, like a traditional morphological taxonomic key, provide criteria that progressively identify an unknown grass sample to a given taxonomic rank. In a prior study of DNA sequences across 20 phylogenetically representative grass species, we identified a series of potentially informative indels in the grass mitochondrial genome. In this study we designed and tested five PCR assays spanning these indels and assessed the feasibility of these assays to aid identification of unknown grass samples. We confirmed that for our control set of 20 samples, on which the design of the PCR assays was based, the five primer combinations produced the expected results. Using these PCR assays in a 'blind test', we were able to identify 25 unknown grass samples with some restrictions. Species belonging to genera represented in our control set were all correctly identified to genus with one exception. Similarly, genera belonging to tribes in the control set were correctly identified to the tribal level. Finally, for those samples for which neither the tribal or genus specific PCR assays were designed, we could confidently exclude these samples from belonging to certain tribes and genera. The results confirmed the utility of the PCR assays and the feasibility of developing a robust full-scale usable grass identification system for forensic purposes.

Microsatellite mutation directed by an external stimulus

Microsatellites are regions of DNA containing tandem repeats of a core 2-6 bp nucleotide sequence. To test the hypothesis that microsatellite mutation can be directed by exposure to specific external cues, control and treatment groups of resistant and susceptible wheat varieties were grown under controlled conditions and genotyped at a number of microsatellite loci that map to chromosomes known to contain Fusarium head blight (FHB) resistance/susceptibility loci. Genotyping was undertaken both prior to and following exposure to Fusarium graminearum, the FHB pathogen. Within a month of inoculation of inflorescences, 58% of experimental plants, and no control plants, had acquired a novel allele at the locus Xgwm112.1. This allele was detected only in head blight affected tissue. Uninoculated control plants, and leaf samples from inoculated plants, showed no mutation. Cloning and sequencing of PCR products indicates that the new allele was generated by contraction of the (CT)(n) repeat motif. Observation of the same deletion-based mutation in all varieties, its absence in control plants not exposed to the head blight pathogen, and the detection of no similar mutational events in a control panel of loci not expected to show mutation, indicates that this example of microsatellite mutation is induced and/or caused by FHB infection.

Genetic structure of a population of the fungus Leptosphaeria maculans in a disease nursery of Brassica napus in Australia

Microsatellite, minisatellite and mating type markers were used to determine the genetic structure of the fungus Leptosphaeria maculans within a disease nursery, where Brassica napus lines were screened for resistance to blackleg disease under high inoculum pressure. Fungal isolates were collected from pseudothecia in infected stubble and pycnidia within cotyledon lesions on seedlings within the nursery. Genetic diversity was high with gene diversity at H=0.700 across four polymorphic loci, and genotypic diversity at D=0.993. Among the 159 isolates analysed, 102 multilocus genotypes were identified. The even distribution of mating type idiomorphs MAT1-1 and MAT1-2 and gametic equilibrium within the population provided further evidence of random mating. Genetic diversity was distributed on a very fine scale in the disease nursery. The majority of genetic diversity (67%) was distributed among conidia within a lesion or among ascospores from a piece of stubble, while the remainder (33%) was distributed within lesions on seedlings or different stubble pieces. There were no among-group differences between samples from stubble and seedlings. This is consistent with the low level of genetic differentiation between the ascospore and conidia samples (F (ST)=0.017) indicating that all isolates of L. maculans from the disease nursery most likely belong to one population, and that ascospores form the primary inoculum in the disease nursery.

Population structure and host-plant specialization in two Scaptodrosophila flower-breeding species

In contrast to phytophagous insect species, little attention has been paid to the possibility of host races in the Drosophilidae, although flower-breeding species, where courtship and mating take place on the flowers, are likely candidates. Two species of Scaptodrosophila, S. hibisci and S. aclinata, are restricted to flowers of Hibiscus species (section Furcaria), and the Furcaria specialization likely predated the separation of S. hibisci and S. aclinata. In all, 20 microsatellite loci were analysed in nine populations of S. hibisci and five of S. aclinata. For two pairs of S. hibisci populations in close proximity, but breeding on different Hibiscus species, differentiation between the populations of each of these pairs was similar to that between the populations that were from the same Hibiscus species, but geographically distant, suggesting the early stages of host-race formation. Genetic variability was significantly less in S. aclinata than in S. hibisci, suggesting greater drift effects in the former. However, of 253 alleles detected, 82 were present in both species, 160 in S. hibisci only and 11 in S. aclinata only, indicating that S. aclinata was derived from S. hibisci, following a strong bottleneck at the time of separation--possibly 40,000 years BP. Analyses and interpretation of Hardy-Weinberg equilibrium and F statistics needed to account for null alleles known to be present at eight loci in S. hibisci, and possibly present at other loci. The results emphasize the need for caution in studies where the presence of null alleles is inferred only from population data.

Linkage maps of microsatellite DNA markers for the Pacific oyster Crassostrea gigas

We constructed male and female consensus linkage maps for the Pacific oyster Crassostrea gigas, using a total of 102 microsatellite DNA markers typed in 11-day-old larvae from three families. We identified 11 and 12 linkage groups in the male and female consensus maps, respectively. Alignment of these separate maps, however, suggests 10 linkage groups, which agrees with the haploid chromosome number. The male linkage map comprises 88 loci and spans 616.1 cM, while the female map comprises 86 loci and spans 770.5 cM. The male and the female maps share 74 loci; 2 markers remain unlinked. The estimated coverages for the consensus linkage maps are 79% for the male and 70-75% for the female, on the basis of two estimates of genome length. Ninety-five percent of the genome is expected to lie within 16 and 21 cM of markers on the male and female maps, respectively, while 95% of simulated minimum distances to the male and female maps are within 10.1 and 13.6 cM, respectively. Females have significantly more recombination than males, across 118 pairs of linked markers in common to the parents of the three families. Significant differences in recombination and orders of markers are also evident among same-sex parents of different families as well as sibling parents of opposite sex. These observations suggest that polymorphism for chromosomal rearrangements may exist in natural populations, which could have profound implications for interpreting the evolutionary genetics of the oyster. These are the first linkage maps for a bivalve mollusc that use microsatellite DNA markers, which should enable them to be transferred to other families and to be useful for further genetic analyses such as QTL mapping.

Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects

Malaria control projects based on the introduction and spread of transgenes into mosquito populations depend on the extent of isolation between those populations. On the basis of the distribution of paracentric inversions, Anopheles gambiae has been subdivided into five subspecific chromosomal forms. Estimating gene flow between and within these forms of An. gambiae presents a number of challenges. We compared patterns of genetic divergence (F(ST)) between sympatric populations of the Bamako and Mopti forms at five sites. We used microsatellite loci within the j inversion on chromosome 2, which is fixed in the Bamako form but absent in the Mopti form, and microsatellites on chromosome 3, a region void of inversions. Estimates of genetic diversity and F(ST)'s suggest genetic exchanges between forms for the third chromosome but little for the j inversion. These results suggest a role for the inversion in speciation. Extensive gene flow within forms among sites resulted in populations clustering according to form despite substantial gene flow between forms. These patterns underscore the low levels of current gene flow between chromosomal forms in this area of sympatry. Introducing refractoriness genes in areas of the genome void of inversions may facilitate their spread within forms but their passage between forms may prove more difficult than previously thought.

The impact of microsatellite electromorph size homoplasy on multilocus population structure estimates in a tropical tree (Corythophora alta) and an anadromous fish (Morone saxatilis)

Microsatellite allelic states are determined by electrophoretic sizing of polymerase chain reaction fragments to define electromorphs. Numerous studies have documented that identical microsatellite electromorphs are potentially heterogeneous at the DNA sequence level, a phenomenon called electromorph size homoplasy. Few studies have examined the impact of electromorph size homoplasy on estimates of population genetic parameters. We investigated the frequency of microsatellite electromorph size homoplasy for 12 loci in the tropical tree Corythophora alta and 11 loci in the anadromous fish Morone saxatilis by sequencing 14-23 homozygotes per locus sampled from multiple populations for a total of 453 sequences. Sequencing revealed no homoplasy for M. saxatilis loci. Seven C. alta loci exhibited homoplasy, including single and compound repeat motifs both with and without interruptions. Between 12.5 and 42.9% of electromorphs sampled per locus showed size homoplasy. Two methods of correction for homoplasy in C. alta generally produced little or no change in single-locus estimates of RST, except for two loci in which some additional differentiation among populations was revealed. Twelve-locus estimates of RST (including the seven loci corrected for homoplasy) were slightly greater than estimates from uncorrected data, although the 95% confidence intervals overlapped. The frequency of methodological errors such as clerical mistakes or sample mislabelling per genotype scored was estimated at 5.4 and 7.3% for C. alta and M. saxatilis, respectively. Simulations showed that the increase in RST produced by homoplasy correction was only slightly larger than variation in RST estimates expected to be caused by methodological errors.

TRbase: a database relating tandem repeats to disease genes for the human genome

MOTIVATION

Tandem repeats are associated with disease genes, play an important role in evolution and are important in genomic organization and function. Although much research has been done on short perfect patterns of repeats, there has been less focus on imperfect repeats. Thus, there is an acute need for a tandem repeats database that provides reliable and up to date information on both perfect and imperfect tandem repeats in the human genome and relates these to disease genes.

RESULTS

This paper presents a web-accessible relational tandem repeats database that relates tandem repeats to gene locations and disease genes of the human genome. In contrast to other available databases, this database identifies both perfect and imperfect repeats of 1-2000 bp unit lengths. The utility of this database has been illustrated by analysing these repeats for their distribution and frequencies across chromosomes and genomic locations and between protein-coding and non-coding regions. The applicability of this database to identify diseases associated with previously uncharacterized tandem repeats is demonstrated.

Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species

We studied microsatellite frequency and distribution in 21.76-Mb random genomic sequences, 0.67-Mb BAC sequences from the Z chromosome, and 6.3-Mb EST sequences of Bombyx mori. We mined microsatellites of >/=15 bases of mononucleotide repeats and >/=5 repeat units of other classes of repeats. We estimated that microsatellites account for 0.31% of the genome of B. mori. Microsatellite tracts of A, AT, and ATT were the most abundant whereas their number drastically decreased as the length of the repeat motif increased. In general, tri- and hexanucleotide repeats were overrepresented in the transcribed sequences except TAA, GTA, and TGA, which were in excess in genomic sequences. The Z chromosome sequences contained shorter repeat types than the rest of the chromosomes in addition to a higher abundance of AT-rich repeats. Our results showed that base composition of the flanking sequence has an influence on the origin and evolution of microsatellites. Transitions/transversions were high in microsatellites of ESTs, whereas the genomic sequence had an equal number of substitutions and indels. The average heterozygosity value for 23 polymorphic microsatellite loci surveyed in 13 diverse silkmoth strains having 2-14 alleles was 0.54. Only 36 (18.2%) of 198 microsatellite loci were polymorphic between the two divergent silkworm populations and 10 (5%) loci revealed null alleles. The microsatellite map generated using these polymorphic markers resulted in 8 linkage groups. B. mori microsatellite loci were the most conserved in its immediate ancestor, B. mandarina, followed by the wild saturniid silkmoth, Antheraea assama.

Haplotypic relationship between SNP and microsatellite markers at the NOS2A locus in two populations

The density of genetic markers required for successful association mapping of complex diseases depends on linkage disequilibrium (LD) between non-functional markers and functional variants. The haplotypic relationship between stable markers and potentially unstable but highly informative markers (e.g. microsatellites) indicates that LD might be maintained over considerable genetic distance in non-African populations, supporting the use of such 'mixed marker haplotypes' in LD-based mapping, and allowing inferences to be drawn about human origins. We investigated sequence variation in the proximal 2.6 kb of the inducible nitric oxide synthase (NOS2A) promoter and the relationship between SNP haplotypes and a pentanucleotide microsatellite (the 'NOS2A(-2.6) microsatellite') in Gambians and UK Caucasians. UK Caucasians exhibited a subset of sequence diversity observed in Gambians, sharing four of 11 SNPs and a similar haplotypic structure. Five SNPs were found in the sequence of interspersed repetitive DNA elements. In both populations, there was dramatic loss of LD between SNP haplotypes and microsatellite alleles across a very short physical distance, suggesting a high intrinsic mutation rate of the NOS2A(-2.6) microsatellite, the SNP haplotypes are relatively ancient, or that this was a region of frequent recombination. Understanding locus- and population-specific LD is essential when designing and interpreting genetic association studies.

Microsatellite analysis reveals a progressive widening of the genetic basis in the elite durum wheat germplasm

It has been argued that the level of genetic diversity in the modern durum wheat ( Triticum turgidum L. var. durum) elite germplasm may have declined due to the high selection pressure applied in breeding programs. In this study, 58 accessions covering a wide spectrum of genetic diversity of the cultivated durum wheat gene pool were characterized with 70 microsatellite loci (or simple sequence repeats, SSRs). On average, SSRs detected 5.6 different allelic variants per locus, with a mean diversity index (DI) equal to 0.56, thus revealing a diversity content comparable to those previously observed with SSRs in other small-grain cereal gene pools. The mean genetic similarity value was equal to 0.44. A highly diagnostic SSR set has been identified. A high variation in allele size was detected among SSR loci, suggesting a different suitability of these loci for estimating genetic diversity. The B genome was characterized by an overall polymorphism significantly higher than that of the A genome. Genetic diversity is organised in well-distinct sub-groups identified by the corresponding foundation-genotypes. A large portion (92.7%) of the molecular variation detected within the group of 45 modern cvs was accounted for by SSR alleles tracing back to ten foundation-genotypes; among those, the most recent CIMMYT-derived founders were genetically distant from the old Mediterranean ones. On the other hand, rare alleles were abundant, suggesting that a large number of genetic introgressions contributed to the foundation of the well-diversified germplasm herein considered. The profiles of recently released varieties indicate that the level of genetic diversity present in the modern durum wheat germplasm has actually increased over time.

Comparative genetic study confirms exceptionally low genetic variation in the ancient and endangered relictual conifer, Wollemia nobilis (Araucariaceae)

The Wollemi pine, Wollemia nobilis (Araucariaceae), was discovered in 1994 as the only extant member of the genus, previously known only from the fossil record. With fewer than 100 trees known from an inaccessible canyon in southeastern Australia, it is one of the most endangered tree species in the world. We conducted a comparative population genetic survey at allozyme, amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) loci in W. nobilis, Araucaria cunninghamii and Agathis robusta - representatives of the two sister genera. No polymorphism was detected at 13 allozyme loci, more than 800 AFLP loci or the 20 SSR loci screened in W. nobilis. In Ag. robusta only one of 12 allozyme loci, five of 800 AFLP loci and none of the 15 SSR loci were variable. For A. cunninghamii, 10 of > 800 AFLP loci and five of 20 SSR loci were variable. Thus low genetic diversity characterizes all three species. While not ruling out the existence of genetic variation, we conclude that genetic diversity is exceptionally low in the Wollemi pine. To our knowledge this is the most extreme case known in plants. We conclude that the combination of small population effects, clonality and below-average genetic variation in the family are probable contributing factors to the low diversity. The exceptionally low genetic diversity of the Wollemi pine, combined with its known susceptibility to exotic fungal pathogens, reinforces current management policies of strict control of access to the pines and secrecy of the pine locations.

Molecular analysis of natural populations of Populus nigra L intermingled with cultivated hybrids

In this study six simple sequence repeats (SSR or microsatellites) were selected for their ability to fingerprint a total of 60 commercial clones of Populus deltoides Marsh. and Populus x canadensis Moench (typically derived from crosses between Populus nigra L and P. deltoides) and to characterize a natural population of P. nigra growing along the Ticino river in the North of Italy. Out of six SSRs used, four microsatellite loci were found to have alleles which were species-specific to P. deltoides and could therefore be used as markers for introgression of P. deltoides into P. nigra. In the studied region hybrid poplars and P. deltoides commercial clones are cultivated as monoclonal stands close to the area where black poplar has its natural habitat. SSR analysis was performed to investigate whether there was evidence of introgression between the natural population and the monoclonal plantations of hybrids and P. deltoides clones cultivated in the surrounding area. Three stages of the natural population were analysed: a group of old trees about a hundred years old, a younger population (aged 2-30 years) and the seedlings of three females of this population. Alleles specific to P. deltoides were detected only in the old cohort of the natural population, while no introgression was observed in the younger individuals and their progenies. These results were also confirmed by isozyme analysis of loci PGI-B, PGM and LAP-A, which were previously identified as diagnostic for P. nigra, P. deltoides and P.xcanadensis.

Cross-species amplification from crop soybean Glycine max provides informative microsatellite markers for the study of inbreeding wild relatives

The development of microsatellite markers through transfer of primers from related species (cross-species amplification) remains a little-explored alternative to the de novo method in plants. In this study of 100 microsatellite loci from Glycine max, we examined two aspects of primer transfer. First, we tested if source locus properties can predict primer transfer and polymorphism in Glycine cyrtoloba and Glycine clandestina. We transferred 23 primers to G. cyrtoloba and 42 to G. clandestina, with 19 loci polymorphic within G. clandestina. However, we could not predict transfer or polymorphism from the source locus properties. Second, we evaluated the subset of 11 polymorphic loci for study in G. clandestina populations representing two local morphotypes. All loci were informative within populations (population mean He +/- SE = 0.58 +/- 0.04). We directly sequenced 28 alleles at 4 representative loci. The allelic patterns and sequencing results established that 8 of 11 loci were typical microsatellites, confirming the utility of primer transfer as an alternative to de novo development. Additionally, we found that morphotypic differentiation between populations was paralleled by changes in polymorphism level at six loci and size homoplasy at one locus. We interpret these patterns as being a product of selfing in G. clandestina. Our results demonstrate the value of allele sequence knowledge for the most effective use of microsatellites.

Molecular genetic evidence that dinoflagellates belonging to the genus Symbiodinium freudenthal are haploid

Microscopic and cytological evidence suggest that many dinoflagellates possess a haploid nuclear phase. However, the ploidy of a number of dinoflagellates remains unknown, and molecular genetic support for haploidy in this group has been lacking. To elucidate the ploidy of symbiotic dinoflagellates belonging to the genus Symbiodinium, we used five polymorphic microsatellites to examine populations harbored by the Caribbean gorgonians Plexaura kuna and Pseudopterogorgia elisabethae; we also studied a series of Symbiodinium cultures. In 690 out of 728 Symbiodinium samples in hospite (95% of the cases) and in all 45 Symbiodinium cultures, only a single allele was recovered per locus. Statistical testing of the Symbiodinium populations harbored by P. elisabethae revealed that the observed genotype frequencies deviate significantly from those expected under Hardy-Weinberg equilibrium. Taken together, our results confirm that, in the vegetative life stage, members of Symbiodinium, both cultured and in hospite, are haploid. Furthermore, based on the phylogenetics of the dinoflagellates, haploidy in vegetative cells appears to be an ancestral trait that extends to all 2,000 extant species of these important unicellular protists.

Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review

Microsatellites, or tandem simple sequence repeats (SSR), are abundant across genomes and show high levels of polymorphism. SSR genetic and evolutionary mechanisms remain controversial. Here we attempt to summarize the available data related to SSR distribution in coding and noncoding regions of genomes and SSR functional importance. Numerous lines of evidence demonstrate that SSR genomic distribution is nonrandom. Random expansions or contractions appear to be selected against for at least part of SSR loci, presumably because of their effect on chromatin organization, regulation of gene activity, recombination, DNA replication, cell cycle, mismatch repair system, etc. This review also discusses the role of two putative mutational mechanisms, replication slippage and recombination, and their interaction in SSR variation.

Microsatellite-centromere mapping in the loach, Misgurnus anguillicaudatus

Primer sets for 15 polymorphic microsatellite loci were developed in the loach, Misgurnus anguillicaudatus (Cobitidae) by molecular cloning and sequencing techniques. Mendelian inheritance was confirmed for the 15 loci by examining the genotypic segregation produced with the primer sets in two full-sib families. The loci were mapped in relation to their centromere in four gynogenetic diploid lines, which were induced by inhibition of the second meiotic division after fertilization with genetically inert sperm. Microsatellite-centromere recombination rates ranged between 0.06 and 0.95 under the assumption of complete interference. Thus, these loci are distributed from the centromeres to the telomeres of their respective chromosomes. The success of mitotic gynogenesis, produced by suppression of the first cleavage, was verified by homozygosity at three diagnostic microsatellite loci that exhibited high gene-centromere meiotic recombination rates in the same family. The differences in heterozygosity levels observed with these markers were attributed to differences in the temporal application of heat shock following inert sperm activation.

Mutation rate and specificity analysis of tetranucleotide microsatellite DNA alleles in somatic human cells

We have systematically varied microsatellite sequence composition to determine the effects of repeat unit size, G+C content, and DNA secondary structure on microsatellite stability in human cells. The microsatellites were inserted in frame within the 5' region of the herpes simplex virus thymidine kinase (HSV-tk) gene. The polypyrimidine/polypurine microsatellites displayed enhanced S1 nuclease sensitivity in vitro, consistent with the formation of non-B-form DNA structures. Microsatellite mutagenesis studies were performed with a shuttle vector system in which inactivating HSV-tk mutations are measured after replication in a nontumorigenic cell line. A significant increase in the HSV-tk mutation frequency per cell generation was observed after insertion of [TTCC/AAGG]9, [TTTC/AAAG]9, or [TCTA/AGAT]9 sequences (P <or= 0.0002), relative to the HSV-tk gene control. We observed that the G + C content of the microsatellite may affect mutagenesis, as the mean microsatellite mutation rates of the [TTTC/AAAG]9 and [TCTA/AGAT]9 alleles were sevenfold and 11-fold higher, respectively, than the [TTCC/AAGG]9 allele. A bias toward expansion mutations was noted for the majority of clones bearing the [TTCC/AAGG]9 allele as well as a [TC/AG]17 microsatellite of similar allele length. The mean microsatellite mutation rate of the [TTCC/AAGG]9 allele did not differ significantly from that for a [TC/AG]11 allele, demonstrating that these tetranucleotide and dinucleotide alleles are of equivalent stability. It is known that microsatellite mutagenesis is affected by the number of repeat units within an allele. Our data suggest that additional biochemical factors may regulate both the rate and specificity of somatic cell microsatellite mutagenesis.

Towards selective breeding of Atlantic salmon for sea louse resistance: approaches to identify trait markers

Several species of sea lice, in particular Lepeophtheirus salmonis (Krøyer), affect the welfare and condition of farmed salmon, with an estimated annual cost to the Scottish industry of 15-30 Pounds million. In Atlantic salmon, some stocks show resistance to L salmonis. Such natural resistance could be utilized for stock improvement using molecular genetic technologies. The development of molecular markers linked to resistance genes, allowing the identification of resistant fish, could increase the efficacy of selective breeding programmes. Various approaches to achieve this goal are described. One way to identify genes conferring resistance is to develop screens for salmon genes that are activated upon louse infection. One such screen--differential display--requires no previous knowledge of gene sequences, involves no preconceptions about which gene families are involved and can therefore identify novel genes. Preliminary results of comparative gene expression in sea-louse-challenged and control fish illustrate the application of differential display.

Towards a standard battery of microsatellite markers for the analysis of the Leishmania donovani complex

The investigation of microsatellite markers has recently superseded that of isoenzymes for many population-biology applications. Microsatellites have the advantages of being dominant, neutral, highly polymorphic and easily scored by high-throughput methods. However, it is necessary to develop a new panel of markers for each group of organisms of interest. Previously, only about 5% of the markers that amplify Leishmania major microsatellite loci were also found to amplify L. donovani loci. A panel of 20 microsatellite markers that are polymorphic in L. donovani and L. infantum has now been developed, using a rapid-enrichment method that will be suitable for developing libraries of markers for other trypanosomatid species. This is the first panel of polymorphic microsatellite markers, to be isolated de novo from any species of Leishmania, that is large enough for population-biology applications.

The gene of chlamysin, a marine invertebrate-type lysozyme, is organized similar to vertebrate but different from invertebrate chicken-type lysozyme genes

In a recent publication we reported the protein purification, characterization, and the gene isolation of a cDNA encoding the antibacterial cold-active lysozyme-like protein chlamysin from the marine bivalve Chlamys islandica. A 4.2 kb genomic chlamysin gene has now been amplified and sequence-analyzed. By comparison to the cDNA sequence and its translation product, the coding region was found separated in four exons of 38-252 bp. The introns range in size from 0.8 to 1.5 kb, and have traditional spliceosomal intron 5'-GT donor and 3'-AG acceptor sites for splicing. Two of the introns contain multiple copies of three sequence motifs not found repeated in other published genes. The over-all gene organization of chlamysin resembles chicken-type (c-type) lysozyme genes in vertebrates, but is different from the three-exon structure in invertebrate c-type lysozyme genes. A phylogenetic analysis of invertebrate-type (i-type) and c-type lysozyme proteins demonstrated a large evolutionary distance between the i-type and the c-type enzyme classes. Exons of the i-type genes are not equally organized according to their homolog protein domains.

Review: proteins with repeated sequence--structural prediction and modeling

The relationship between the amino acid sequence and the three-dimensional structure of proteins with internal repeats is discussed. In particular, correlations between the amino acid composition and the ability to fold in a unique structure, as well as classification of the structures based on their repeat length, are described. This analysis suggests rules that can be used for the structural prediction of repeat-containing proteins. The paper is focused on prediction and modeling of solenoid-like proteins with the repeat length ranging between 5 and 40 residues. The models of leucine-rich repeat proteins and bacterial proteins with pentapeptide repeats are examined in light of the recently solved structures of the related molecules.