Microsatellites, from molecules to populations and back

Development and polymorphism of simple sequence repeat DNA markers for Shorea curtisii and other Dipterocarpaceae species

Nine simple sequence repeat (SSR) markers were developed from Shorea curtisii using two different methods. One SSR locus was isolated by the commonly used method of screening by colony hybridization, and the other eight loci were isolated by a vectorette PCR method. Primer pairs were designed based on the sequences of all these SSR loci. Analysis of 40 individuals of S. curtisii from natural forest in Malaysia revealed that all SSR loci were polymorphic. Four SSR markers, Shc01, Shc04, Shc07 and Shc09, were highly polymorphic. We have also tested the applicability of these SSR printers to other species of Dipterocarpaceae using PCR amplification. Because the flanking region sequences of the S. curtisii SSRs were well conserved within this family, the SSR primers for S. curtisii can be applied to almost all species of Dipterocarpaceae.

Genealogical inference from microsatellite data

Ease and accuracy of typing, together with high levels of polymorphism and widespread distribution in the genome, make microsatellite (or short tandem repeat) loci an attractive potential source of information about both population histories and evolutionary processes. However, microsatellite data are difficult to interpret, in particular because of the frequency of back-mutations. Stochastic models for the underlying genetic processes can be specified, but in the past they have been too complicated for direct analysis. Recent developments in stochastic simulation methodology now allow direct inference about both historical events, such as genealogical coalescence times, and evolutionary parameters, such as mutation rates. A feature of the Markov chain Monte Carlo (MCMC) algorithm that we propose here is that the likelihood computations are simplified by treating the (unknown) ancestral allelic states as auxiliary parameters. We illustrate the algorithm by analyzing microsatellite samples simulated under the model. Our results suggest that a single microsatellite usually does not provide enough information for useful inferences, but that several completely linked microsatellites can be informative about some aspects of genealogical history and evolutionary processes. We also reanalyze data from a previously published human Y chromosome microsatellite study, finding evidence for an effective population size for human Y chromosomes in the low thousands and a recent time since their most recent common ancestor: the 95% interval runs from approximately 15, 000 to 130,000 years, with most likely values around 30,000 years.

Density and variability of dinucleotide microsatellites in the parthenogenetic polyploid snail Melanoides tuberculata

Characterization of microsatellites in the parthenogenetic polyploid snail Melanoides tuberculata revealed an unusual high density of dinucleotide repeats. Multiple banding patterns were obtained at these loci, and interpreted as a consequence of polyploidy. Microsatellite variability was low within, but high between, shell morphotypes. Genotypes were wholly transmitted from mothers to offspring. These results suggest that reproduction is strictly apomictic, and that shell morphotypes are genetic clones.

Cloning and characterization of a microsatellite in the mitochondrial control region of the African side-necked turtle, Pelomedusa subrufa

The nucleotide sequence of the African side-necked turtle mitochondrial control region and its flanking tRNA genes was determined. This 73% A+T-rich region is 1194 bp long. Several conserved motifs involved in the regulation of the mitochondrial genome replication process, including one conserved sequence block (CSB1), and three termination-associated sequences were identified. The most remarkable feature found in this control region was the presence of six microsatellite-containing tandem repeats between the CSB1 motif and the tRNAPhe gene. The potential usefulness of this microsatellite sequence for population-level studies is enhanced by its unique localization in the maternally inherited mitochondrial molecule.

Diploid Musa acuminata genetic diversity assayed with sequence-tagged microsatellite sites

The sequence-tagged microsatellite site (STMS) discrimination potential was explored using nine microsatellite primer pairs. STMS polymorphism was assayed by nonradioactive urea-polyacrylamide gel electrophoresis. Genetic relationships were examined among 59 genotypes of wild or cultivated accessions of diploid Musa acuminata. The organization of the subspecies was confirmed and some clone relationships were clarified.

Population structure of nuclear and mitochondrial DNA variation among humpback whales in the North Pacific

The population structure of variation in a nuclear actin intron and the control region of mitochondrial DNA is described for humpback whales from eight regions in the North Pacific Ocean: central California, Baja Peninsula, nearshore Mexico (Bahia Banderas), offshore Mexico (Socorro Island), southeastern Alaska, central Alaska (Prince Williams Sound), Hawaii and Japan (Ogasawara Islands). Primary mtDNA haplotypes and intron alleles were identified using selected restriction fragment length polymorphisms of target sequences amplified by the polymerase chain reaction (PCR-RFLP). There was little evidence of heterogeneity in the frequencies of mtDNA haplotypes or actin intron alleles due to the year or sex composition of the sample. However, frequencies of four mtDNA haplotypes showed marked regional differences in their distributions (phi ST = 0.277; P < 0.001; n = 205 individuals) while the two alleles showed significant, but less marked, regional differences (phi ST = 0.033; P < 0.013; n = 400 chromosomes). An hierarchical analysis of variance in frequencies of haplotypes and alleles supported the grouping of six regions into a central and eastern stock with further partitioning of variance among regions within stocks for haplotypes but not for alleles. Based on available genetic and demographic evidence, the southeastern Alaska and central California feeding grounds were selected for additional analyses of nuclear differentiation using allelic variation at four microsatellite loci. All four loci showed significant differences in allele frequencies (overall FST = 0.043; P < 0.001; average n = 139 chromosomes per locus), indicating at least partial reproductive isolation between the two regions as well as the segregation of mtDNA lineages. Although the two feeding grounds were not panmictic for nuclear or mitochondrial loci, estimates of long-term migration rates suggested that male-mediated gene flow was several-fold greater than female gene flow. These results include and extend the range and sample size of previously published work, providing additional evidence for the significance of genetic management units within oceanic populations of humpback whales.

The identification and characterization of microsatellites in the compact genome of the Japanese pufferfish, Fugu rubripes: perspectives in functional and comparative genomic analyses

Fugu rubripes (Fugu) has one of the smallest recorded vertebrate genomes and is an economic tool for comparative DNA sequence analysis. Initial characterization of 128 kb of Fugu DNA attributed the compactness of this genome, in part, to a sparseness of repetitive DNA sequence compared with mammalian genomic sequences. This paper describes a new and comprehensive analysis in which 501 theoretically possible microsatellites with a repeat unit of one to six bases were used to query two orders of magnitude more Fugu DNA (i.e. 11.338 Mb). A total of 6042 microsatellites were identified and categorized. In decreasing order, the 20 most frequently occurring microsatellites are AC, A, C, AGG, AG, AGC, AAT, AAAT, ACAG, ACGC, ATCC, AAC, ATC, AGGG, AAAG, AAG, AAAC, AT, CCG and TTAGGG. The 20 most frequently occurring microsatellites represent 81.79% of all microsatellites identified. Our results indicate that one microsatellite occurs every 1.876 kb of DNA in Fugu, 11.55% of the microsatellites are detected in open reading frames that are predicted protein coding regions. With respect to the proportion of microsatellites present in open reading frames and the total abundance (bp) of all microsatellites, the genome of Fugu is similar to the genome of many other vertebrate species. Previous estimates performed indicate that approximately 1% of many vertebrate genomes are comprized of microsatellite sequences. However, many differences prevail in the abundance and frequency of the individual microsatellite classes. Many of the frequently occurring microsatellites in Fugu are known to code in other species for regions in proteins such as transcription factors, whilst others are associated with known functions, such as transcription factor binding sites and form part of promoter regions in DNA sequences of genes. Therefore, it is likely that such repeats in genomes have a role in the evolution of genes, regulation of gene expression and consequently the evolution of species.

Microsatellites reveal heterosis in red deer

The fitness consequences of inbreeding and outbreeding are poorly understood in natural populations. We explore two microsatellite-based variables, individual heterozygosity (likely to correlate with recent inbreeding) and a new individual-specific internal distance measure, mean d2 (focusing on events deeper in the pedigree), in relation to two measures of fitness expressed early in life, birth weight and neonatal survival, in 670 red deer calves (Cervus elaphus) born on the Isle of Rum between 1982 and 1996. For comparison, we also analyse inbreeding coefficients derived from pedigrees in which paternity was inferred by molecular methods. Only 14 out of 231 calves (6.1%) had non-zero inbreeding coefficients, and neither inbreeding coefficient nor individual heterozygosity was consistently related to birth weight or neonatal survival. However, mean d2 was consistently related to both fitness measures. Low mean d2 was associated with low birth weight, especially following cold Aprils, in which foetal growth is reduced. Low mean d2 was also associated with low neonatal survival, but this effect was probably mediated by birth weight because fitting birth weight to the neonatal survival model displaced mean d2 as an explanatory variable. We conclude that in the deer population fitness measures expressed early in life do not show evidence of inbreeding depression, but they do show evidence of heterosis, possibly as a result of population mixing. We also demonstrate the practical problems of estimating inbreeding via pedigrees compared with a direct marker-based estimate of individual heterozygosity. We suggest that, together, individual heterozygosity and mean d2, estimated using microsatellites, are useful tools for exploring inbreeding and outbreeding in natural population.

The distribution and frequency of microsatellite loci in Drosophila melanogaster

We report the results of a comprehensive search of Drosophila melanogaster DNA sequences in GenBank for di-, tri-, and tetranucleotide repeats of more than four repeat units, and a DNA library screen for dinucleotide repeats. Dinucleotide repeats are more abundant (66%) than tri- (30%) or tetranucleotide (4%) repeats. We estimate that 1917 dinucleotide repeats with 10 or more repeat units are present in the euchromatic D. melanogaster genome and, on average, they occur once every 60 kb. Relative to many other animals, dinucleotide repeats in D. melanogaster are short. Tri- and tetranucleotide repeats have even fewer repeat units on average than dinucleotide repeats. Our WorldWide Web site (http://www.bio.cornell.edu/genetics/aquadro/+ ++aquadro.html) posts the complete list of 1298 microsatellites (> or = five repeat units) identified from the GenBank search. We also summarize assay conditions for 70 D. melanogaster microsatellites characterized in previous studies and an additional 56 newly characterized markers.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR-single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Genetic analysis of the freshwater crayfish Cherax tenuimanus

The marron (Cherax tenuimanus) is one of the few species of freshwater crayfish native to Australia that is suitable for aquaculture and occurs only in the southwest of Western Australia. This study describes polymorphic microsatellite markers which differentiate marron populations from several geographically distinct regions (including rivers and streams, dams, and commercial marron farms) throughout Western Australia. Twenty microsatellite loci, primarily of the (CA)n. (GT)n type, were isolated and sequenced from a marron cosmid library. Three of these loci were characterised further. Two loci exhibited extensive polymorphism and one was monomorphic. The polymorphic loci exhibited Mendelian codominant inheritance in the family group comprising two individual parents and approximately 100 offspring bred for this study. These loci permitted differentiation between the five geographically distinct populations studied and thus provide a basis for genetic characterisation of marron stock in Western Australia.

Assessment of genetic diversity by DNA profiling and its significance in silkworm, Bombyx mori

Silkworm genetic resources that are being maintained in different countries are yet to be adequately tapped to develop elite varieties that are suited to different agro-eco-climatic conditions of countries like India. This is mostly due to unavailability of efficient protocols that could uncover usable genetic variability in silkworms. Molecular markers are known to provide unambiguous estimates of genetic variability of populations since they are independent of confounding effects of environment. The DNA fingerprinting assays, based on random amplified polymorphic DNA (RAPD) and banded krait minor satellite DNA (Bkm) 2(8) multilocus probes, which successfully characterise the diverse silkworm genotypes at their DNA level, are described. The use of these two DNA fingerprinting assays in estimation of within- and between-population genetic diversity is discussed.

Low frequency of microsatellites in the avian genome

A better insight into the occurrence of microsatellites in a range of taxa may help to understand the evolution of simple repeats. Previous studies have found the relative abundance of several repeat motifs to differ among mammals, invertebrates, and plants. Absolute numbers of microsatellites also tend to correlate positively with genome size. We analyzed the occurrence, frequency, and distribution of microsatellites in birds, a taxon with one of the smallest known genome sizes among vertebrates. Dot-blot hybridization revealed that about half of 22 different di-, tri-, and tetranucleotide repeat motifs were clearly more common in human than in three species of birds: chicken, woodpecker, and swallow. For the remaining motifs no clear difference was found. From searching avian database sequences we estimated there to be 30,000-70,000 microsatellites longer than 20 bp in the avian genome. The number of (CA) > or = 10 would be around 7000-9000 and the number of (CA) > or = 14 about 3000. The calculated density of avian microsatellites (total, one every 20-39 kb; (CA) > or = 10, one every 136-150 kb) is much lower than that estimated for the human genome (one every 6 and 30 kb, respectively). This may be explained by the fact that the avian genome contains relatively less noncoding DNA than most mammals and that avian SINE/LINE elements do not terminate in poly(A) tails, which are known to provide a resource for the evolution of simple repeats in mammals. We found no association between microsatellites and SINEs in birds. Primed in situ labeling suggested fairly even distribution of (CA)n repeats over chicken macrochromosomes and intermediate chromosomes, whereas the microchromosomes, a large part of the Z and W chromosomes, and most telomeres and centromeres had very low concentrations of (CA)n microsatellites. The scarcity of microsatellites on the microchromosomes is compatible to these regions likely being unusually rich in coding sequences. The low microsatellite density in the genome in general and on the microchromosomes in particular imposes an obstacle for the development of marker-rich genetic maps of chicken and other birds, and for the localization of quantitative trait genes.

Simple sequence repeats as a source of quantitative genetic variation

Most traits in biological populations appear to be under stabilizing selection, which acts to eliminate quantitative genetic variation. Yet, virtually all measured traits in biological populations continue to show significant quantitative genetic variation. The paradox can be resolved by postulating the existence of an abundant, though unspecified, source of mutations that has quantitative effects on phenotype, but does not reduce fitness. Does such a source actually exist? We propose that it does, in the form of repeat-number variation in SSRs (simple sequence repeats, of which the triplet repeats of human neurodegenerative diseases are a special case). Viewing SSRs as a major source of quantitative mutation has broad implications for understanding molecular processes of evolutionary adaptation, including the evolutionary control of the mutation process itself.