Estimation of genotyping error rate from repeat genotyping,unintentional recaptures and known parent–offspring comparisons in 16 microsatellite loci for brown rockfish (Sebastes auriculatus)

Genotyping errors are present in almost all genetic data and can affect biological conclusions of a study, particularly for studies based on individual identification and parentage. Many statistical approaches can incorporate genotyping errors, but usually need accurate estimates of error rates. Here, we used a new microsatellite data set developed for brown rockfish (Sebastes auriculatus) to estimate genotyping error using three approaches: (i) repeat genotyping 5% of samples, (ii) comparing unintentionally recaptured individuals and (iii) Mendelian inheritance error checking for known parent–offspring pairs. In each data set, we quantified genotyping error rate per allele due to allele drop‐out and false alleles. Genotyping error rate per locus revealed an average overall genotyping error rate by direct count of 0.3%, 1.5% and 1.7% (0.002, 0.007 and 0.008 per allele error rate) from replicate genotypes, known parent–offspring pairs and unintentionally recaptured individuals, respectively. By direct‐count error estimates, the recapture and known parent–offspring data sets revealed an error rate four times greater than estimated using repeat genotypes. There was no evidence of correlation between error rates and locus variability for all three data sets, and errors appeared to occur randomly over loci in the repeat genotypes, but not in recaptures and parent–offspring comparisons. Furthermore, there was no correlation in locus‐specific error rates between any two of the three data sets. Our data suggest that repeat genotyping may underestimate true error rates and may not estimate locus‐specific error rates accurately. We therefore suggest using methods for error estimation that correspond to the overall aim of the study (e.g. known parent–offspring comparisons in parentage studies).     

Plucked hair samples as a source of DNA: reliability of dinucleotide microsatellite genotyping

To test whether plucked hairs are a reliable source of DNA for genotyping microsatellite loci, we carried out experiments using one, three, or 10 hairs per extract for 50 alpine marmots. For each extract, seven independent genotypings were performed for the same locus (multiple-tubes approach). Two types of genotyping errors were recorded: a false homozygote defined as the detection of only one allele of a true heterozygote, and a false allele defined as a PCR-generated allele that was not one of the alleles of the true genotype. Using DNA extracted from one, three, or 10 hairs, the overall error rate was 14.00%, 4.86%, and 0.29%, respectively. Based on our results, we conclude that 10 hairs should be used to obtain consistently reliable genotypings using the single-tube approach, and that a single plucked hair could represent a reliable source of DNA if the multiple-tubes approach is used. For future studies of dinucleotide repeat diversity using DNA extracted from one to three shed or plucked hairs, we strongly recommend initiating an appropriate pilot study to quantify the error rate and to determine the reliability of the single-tube approach.     

A novel trinucleotide repeat expansion at chromosome 3q26.2 identified by a CAG/CTG repeat expansion detection array

CAG/CTG repeat expansions cause at least 12 different neurological disorders, and additional disorders of this type probably exist. Using the repeat expansion detection (RED) assay, we identified an expanded CAG/CTG repeat in a 50-year-old woman with an autosomal dominant syndrome with prominent progressive sensory neuropathy. The expansion could not be accounted for by any of the CAG/CTG repeats known to undergo expansion. To identify the locus of the expansion, we created a PCR array to assess the repeat length of all repeats of eight or more CAG or CTG triplets in the human genome. The expansion was localized to a repeat contained in an intron of a Genscan-predicted gene, 185 nt downstream of a predicted exon that is conserved through mouse. The closest experimentally verified gene in the region (TNIK, encoding a serine/threonine kinase) occurs approximately 63 Kb downstream from the repeat. The length of the expansion in the proband is 98 triplets. This repeat is not expanded in the proband’s cousin (the only other affected family member for whom DNA is currently available) and no expansions were detected in a set of 230 patients with movement disorders of unknown cause. An expanded allele containing 58 triplets was detected in a single control individual, and no other expansions were detected in a set of 255 controls. The normal repeat length ranges from 5 to 30 triplets, with 8 triplets the most common allele. Our results suggest that this new repeat expansion is probably not the direct cause of the phenotype in the proband. Whether the repeat contributes to the patient’s phenotype, or is associated with another phenotype, remains to be determined.Electronic Supplementary Material Supplementary material is available for this article at .     

The impact of library preparation protocols on the consistency of allele frequency estimates in Pool‐Seq data

Sequencing pools of individuals (Pool‐Seq) is a cost‐effective method to determine genome‐wide allele frequency estimates. Given the importance of meta‐analyses combining data sets, we determined the influence of different genomic library preparation protocols on the consistency of allele frequency estimates. We found that typically no more than 1% of the variation in allele frequency estimates could be attributed to differences in library preparation. Also read length had only a minor effect on the consistency of allele frequency estimates. By far, the most pronounced influence could be attributed to sequence coverage. Increasing the coverage from 30‐ to 50‐fold improved the consistency of allele frequency estimates by at least 27%. We conclude that Pool‐Seq data can be easily combined across different library preparation methods, but sufficient sequence coverage is key to reliable results.     

Ligase-based detection of mononucleotide repeat sequences

Up to 15% of all colorectal cancers are considered to be replication error positive (RER(+)) and contain mutations at hundreds of thousands of microsatellite repeat sequences. Recently, a number of intragenic mononucleotide repeat sequences have been demonstrated to be targets for inactivating genes in RER(+)colorectal tumors. In this study, thermostable DNA ligases were tested for the ability to detect alterations in microsatellite sequences in colon tumor samples. Ligation profiles on mononucleotide repeat sequences were determined for four related thermostable DNA ligases, Thermus thermophilus ( Tth ) ligase, Thermus sp. AK16D ligase, Aquifex aeolicus ligase and the K294R mutant of the Tth ligase. While the limit of detection for point mutations was one mutation in 1000 wild-type sequences, the ability to detect a single base deletion in a 10 base mononucleotide repeat was one mutation in 100 wild-type sequences. Furthermore, the misligation error increased exponentially as the length of the mono-nucleotide repeat increased, and was 10% of the correct signal for a 19 base mononucleotide repeat. A fluorescent ligase-based assay [polymerase chain reaction/ligase detection reaction (PCR/LDR)] correlated with results obtained using a radioactive assay to detect instability within the TGF-beta Type II receptor gene. PCR/LDR was also used to detect the APCI1307K mononucleotide repeat allele which has a carrier frequency of 6.1% in Ashkenazi Jewish individuals. In a blind study, 30 samples that had been typed for the presence of the APCI1307K allele were tested. The PCR/LDR results correlated with those obtained using sequencing and allele-specific oligonucleotide hybridization for 16 samples carrying the mutation and 13 wild-type samples. Ligation assays that characterize mononucleotide repeats can be used to rapidly detect somatic mutations in tumors, and to screen for individuals who have a hereditary predisposition to develop colon cancer.     

A prevalent POLG CAG microsatellite length allele in humans and African great apes

The human nuclear gene for the catalytic subunit of mitochondrial DNA polymerase (POLG) contains within its coding region a CAG microsatellite encoding a polyglutamine repeat. Previous studies demonstrated an association between length variation at this repeat and male infertility, suggesting a mechanism whereby the prevalent (CAG)₁₀ allele, which occurs at a frequency of >80% in different populations, could be maintained by selection. Sequence analysis of the POLG CAG microsatellite region of more than 1000 human chromosomes reveals that virtually all allelic variation at the locus is accounted for by length variation of the CAG repeat. Analysis of POLG from African great apes shows that a prevalent length allele is present in each species, although its exact length is species-specific. In common chimpanzee (Pan troglodytes) a number of different sequence variants contribute to the prevalent length allele, strongly supporting the idea that the length of the POLG microsatellite region, rather than its exact nucleotide or amino acid sequence, is what is maintained. Analysis of POLG in other primates indicates that the repeat has expanded from a shorter, glutamine-rich sequence, present in the common ancestor of Old and New World monkeys.     

C9ORF72 Repeat Expansion in Australian and Spanish Frontotemporal Dementia Patients

A hexanucleotide repeat expansion in C9ORF72 has been established as a common cause of frontotemporal dementia (FTD). However, the minimum repeat number necessary for disease pathogenesis is not known. The aims of our study were to determine the frequency of the C9ORF72 repeat expansion in two FTD patient collections (one Australian and one Spanish, combined n = 190), to examine C9ORF72 expansion allele length in a subset of FTD patients, and to examine C9ORF72 allele length in ‘non-expansion’ patients (those with <30 repeats). The C9ORF72 repeat expansion was detected in 5–17% of patients (21–41% of familial FTD patients). For one family, the expansion was present in the proband but absent in the mother, who was diagnosed with dementia at age 68. No association was found between C9ORF72 non-expanded allele length and age of onset and in the Spanish sample mean allele length was shorter in cases than in controls. Southern blotting analysis revealed that one of the nine ‘expansion-positive’ patients examined, who had neuropathologically confirmed frontotemporal lobar degeneration with TDP-43 pathology, harboured an ‘intermediate’ allele with a mean size of only ∼65 repeats. Our study indicates that the C9ORF72 repeat expansion accounts for a significant proportion of Australian and Spanish FTD cases. However, C9ORF72 allele length does not influence the age at onset of ‘non-expansion’ FTD patients in the series examined. Expansion of the C9ORF72 allele to as little as ∼65 repeats may be sufficient to cause disease.     

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

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

Multiple levels of single-strand slippage at cetacean tri- and tetranucleotide repeat microsatellite loci.

Between three and six tri- and tetranucleotide repeat microsatellite loci were analyzed in 3720 samples collected from four different species of baleen whales. Ten of the 18 species/locus combinations had imperfect allele arrays, i.e., some alleles differed in length by other than simple integer multiples of the basic repeat length. The estimate of the average number of alleles and heterozygosity was higher at loci with imperfect allele arrays relative to those with perfect allele arrays. Nucleotide sequences of 23 different alleles at one tetranucleotide repeat microsatellite locus in fin whales, Balaenoptera physalus, and humpback whales, Megaptera novaeangliae, revealed sequence changes including perfect repeats only, multiple repeats, and partial repeats. The relative rate of the latter two categories of mutation was estimated at 0.024 of the mutation rate involving perfect repeats only. It is hypothesized that single-strand slippage of partial repeats may provide a mechanism for counteracting the continuous expansion of microsatellite loci, which is the logical consequence of recent reports demonstrating directional mutations. Partial-repeat mutations introduce imperfections in the repeat array, which subsequently could reduce the rate of single-strand slippage. Limited computer simulations confirmed this predicted effect of partial-repeat mutations.     

Analysis of somatic mutations at human minisatellite loci in tumors and cell lines

Hypervariable human minisatellite loci show a substantial level of germline instability, and spontaneous mutation rates to new length alleles have been measured directly by pedigree analysis. We now show that mutation events altering the number of minisatellite repeat units are not restricted to the germline, but also arise in other tissues. Mutant alleles can be detected at a very low frequency in lymphoblastoid cell lines and at much higher frequencies in clonal tumor cell populations, most particularly in gastrointestinal adenocarcinomas. Mutant alleles in these tumors are usually present at a dosage equal to or greater than that of the progenitor allele, indicating that most or all of the tumor cells carry the same clonally derived mutant allele. As with germline mutation, the incidence of somatic mutations in tumors varies from locus to locus, with the same locus showing the highest level of germline and somatic instability. Most length changes, as those in the germline, are of only a few repeat units; however, very large changes are also observed, implying that such mutations can occur in the absence of meiosis.     

The influence of sequence divergence between alleles of the human MS205 minisatellite incorporated into the yeast genome on length-mutation rates and lethal recombination events during meiosis

Certain minisatellites exhibit hypervariability with respect to the number of repeat units and, thus, allele length. Such polymorphism is generated by germline-specific recombinational events that occur at high frequencies and lead to the gain or loss of repeat units. In order to elucidate the molecular details of mutagenesis in minisatellites, we have integrated human minisatellites into the yeast genome in the vicinity of a hotspot for meiotic double-strand breaks (DSBs). Here, we describe the results of tetrad analyses of mutations in the human MS205 minisatellite in yeast strains heterozygous for alleles composed of 51 and 31 repeat units, as well as in a strain homozygous for the same 51 repeat unit allele. The length-mutation rate was twice as high in the heterozygous strain as in the homozygous strain, suggesting that sequence divergence between alleles enhances the generation of length mutations. In the case of heterozygotes, the frequency of length mutants resulting from inter-allelic exchange was significantly higher in tetrads with three viable spores than in tetrads with four viable spores, indicating that there is a higher probability for spore mortality in tetrads originating from meioses during which inter-allelic exchange of repeat units occurs. In an attempt to explain these findings, we propose a model for minisatellite mutation involving recombination, in which sequence divergence between alleles results in a heteroduplex containing numerous mismatches. We suggest that convergent mismatch-repair tracts in this heteroduplex give rise to a DSB that may be repaired by an additional round of recombination resulting in mutation of a third allele, or be lethal if such recombination fails. It appears probable that the formation of such additional mutants is the major explanation for the difference in meiotic length-mutation rates between the heterozygous and homozygous yeast strains, and that this phenomenon contributes to high germline length-mutation frequencies at minisatellites in humans.     

Reduction of stability of arabidopsis genomic and transgenic DNA-repeat sequences (microsatellites) by inactivation of AtMSH2 mismatch-repair function

Highly conserved mismatch repair (MMR) systems promote genomic stability by correcting DNA replication errors, antagonizing homeologous recombination, and responding to various DNA lesions. Arabidopsis and other plants encode a suite of MMR protein orthologs, including MSH2, the constant component of various specialized eukaryotic mismatch recognition heterodimers. To study MMR roles in plant genomic stability, we used Arabidopsis AtMSH2::TDNA mutant SALK_002708 and AtMSH2 RNA-interference (RNAi) lines. AtMSH2::TDNA and RNAi lines show normal growth, development, and fertility. To analyze AtMSH2 effects on germ line DNA fidelity, we measured insertion-deletion mutation of dinucleotide-repeat sequences (microsatellite instability) at nine loci in 16 or more progeny of two to four different wild-type or AtMSH2-deficient plants. Scoring 992 total alleles revealed 23 (2.3%) unique and 51 (5.1%) total repeat length shifts ([+2], [-2], [+4], or [-4] bp). For the six longest repeat loci, the corresponding frequencies were 22/608 and 50/608. Two of four AtMSH2-RNAi plants showed similar microsatellite instability. In wild-type progeny, only one unique repeat length allele was found in 576 alleles tested. This endogenous microsatellite instability, shown for the first time in MMR-defective plants, is similar to that seen in MMR-defective yeast and mice, indicating that plants also use MMR to promote germ line fidelity. We used a frameshifted reporter transgene, (G)(7)GUS, to measure insertion-deletion reversion as blue-staining beta-glucuronidase-positive leaf spots. Reversion rates increased only 5-fold in AtMSH2::TDNA plants, considerably less than increases in MSH2-deficient yeast or mammalian cells for similar mononucleotide repeats. Thus, MMR-dependent error correction may be less stringent in differentiated leaf cells than in plant equivalents of germ line tissue.     

Noninvasive genotyping and Mendelian analysis of microsatellites in African savannah elephants

We obtained fresh dung samples from 202 (133 mother-offspring pairs) savannah elephants (Loxodonta africana) in Samburu, Kenya, and genotyped them at 20 microsatellite loci to assess genotyping success and errors. A total of 98.6% consensus genotypes was successfully obtained, with allelic dropout and false allele rates at 1.6% (n = 46) and 0.9% (n = 37) of heterozygous and total consensus genotypes, respectively, and an overall genotyping error rate of 2.5% based on repeat typing. Mendelian analysis revealed consistent inheritance in all but 38 allelic pairs from mother-offspring, giving an average mismatch error rate of 2.06%, a possible result of null alleles, mutations, genotyping errors, or inaccuracy in maternity assignment. We detected no evidence for large allele dropout, stuttering, or scoring error in the dataset and significant Hardy-Weinberg deviations at only two loci due to heterozygosity deficiency. Across loci, null allele frequencies were low (range: 0.000-0.042) and below the 0.20 threshold that would significantly bias individual-based studies. The high genotyping success and low errors observed in this study demonstrate reliability of the method employed and underscore the application of simple pedigrees in noninvasive studies. Since none of the sires were included in this study, the error rates presented are just estimates.     

Estimation of allele frequencies for VNTR loci

VNTR loci provide valuable information for a number of fields of study involving human genetics, ranging from forensics (DNA fingerprinting and paternity testing) to linkage analysis and population genetics. Alleles of a VNTR locus are simply fragments obtained from a particular portion of the DNA molecule and are defined in terms of their length. The essential element of a VNTR fragment is the repeat, which is a short sequence of basepairs. The core of the fragment is composed of a variable number of identical repeats that are linked in tandem. A sample of fragments from a population of individuals exhibits substantial variation in length because of variation in the number of repeats. Each distinct fragment length defines an allele, but any given fragment is measured with error. Therefore the observed distribution of fragment lengths is not discrete but is continuous, and determination of distinct allele classes is not straightforward. A mixture model is the natural statistical method for estimating the allele frequencies of VNTR loci. In this article we develop nonparametric methods for obtaining the distribution of allele sizes and estimates of their frequencies. Methods for obtaining maximum-likelihood estimates are developed. In addition, we suggest an empirical Bayes method to improve the maximum-likelihood estimates of the gene frequencies; the empirical Bayes procedure effects a local smoothing. The latter method works particularly well when measurement error is large relative to the repeat size, because the estimated distribution of allele frequencies when maximum likelihood is used is unreliable because of an alternating pattern of over- and underestimation. We define alleles and estimate the allele frequencies for two VNTR loci from the human genome (D17S79 and D2S44), from data obtained from Lifecodes, Inc.     

MultiplexSSR: A pipeline for developing multiplex SSR‐PCR assays from resequencing data

Next‐generation sequencing has greatly promoted the investigation of single nucleotide polymorphisms, while studies of simple sequence repeats are sharply decreasing. However, simple sequence repeats still present some advantages in conservation genetics. In this study, an end‐to‐end pipeline referred to as MultiplexSSR was established to develop multiplex PCR assays in batches with highly polymorphic simple sequence repeats for capillary platforms from resequencing data. The distribution of single sequence repeats in the genome, the error profiles of genotypes and allelotypes, and the increase in the allele length range depending on the number of individuals were investigated. A total of 98% of single sequence repeats presented lengths of less than 100 bp. The error rate of the genotyping and allelotyping of dimeric patterns was ten times higher than those for other patterns. The error rate of allelotyping was less than that of genotyping. The allele length range reached approximate saturation with 10 individuals. This pipeline uses allele numbers to select highly polymorphic loci, masks loci with variation, and applies in silico PCR to improve primer specificity. The application of the developed multiplex SSR‐PCR assays validated the pipeline's robustness, showing higher polymorphism and stability for the developed simple sequence repeats and a lower cost for genotyping and providing low‐depth resequencing data from less than a dozen individuals for the development of markers. This pipeline fills the gap between next‐generation sequencing and multiplex SSR‐PCR.     

A high frequency of length polymorphisms in repeated sequences adjacent to Alu sequences

We describe a new class of DNA length polymorphism that is due to a variation in the number of tandem repeats associated with Alu sequences (Alu sequence-related polymorphisms). The polymerase chain reaction was used to selectively amplify a (TTA)n repeat identified in the 3-hydroxy-3-methylglutaryl coenzyme A (HMG CoA) reductase gene from genomic DNA of 41 human subjects, and the size of the amplified products was determined by gel electrophoresis. Seven alleles were found that differed in size by integrals of three nucleotides. The allele frequencies ranged from 1.5% to 52%, and the overall heterozygosity index was 62%. The polymorphic TTA repeat was located adjacent to a repetitive sequence of the Alu family. A homology search of human genomic DNA sequences for the trinucleotide TTA (at least five members in length) revealed tandem repeats in six other genes. Three of the six (TTA)n repeats were located adjacent to Alu sequences, and two of the three (in the genes for beta-tubulin and interleukin-1 alpha) were found to be polymorphic in length. Tandemly repetitive sequences found in association with Alu sequences may be frequent sites of length polymorphism that can be used as genetic markers for gene mapping or linkage analysis.     

Data on the CGG repeat at the fragile X site in the non-retarded Japanese population and family suggest the presence of a subgroup of normal alleles predisposing to mutate

The fragile X mutation is the result of amplification in the repeat number of p(CGG) _n in FMR-1; alleles with more than 52 repeats have been shown to be so unstable as to mutate in the repeat number in almost every transmission. To improve our understanding of mutations in normal alleles of FMR-1, the following studies were carried out in the Japanese population: a study on length variation in the repeat to determine the allele distribution of the repeat length in a non-retarded population, family studies to observe new mutations in normal allele, and haplotype analyses with microsatellite markers flanking the repeat to confirm estimated mutation rates and founder chromosomes in the fragile X syndrome. Analysis of the p(CGG) _n in 370 unrelated males detected 24 distinct alleles with repeats of 18–44. A comparison with previously reported data suggests the presence of racial/ethnic differences in the allele distribution. No premutation allele was found in 824 unrelated X chromosomes examined by the polymerase chain reaction and Southern blot analysis. Family studies detected one new mutation in a total of 303 meioses. However, the mutation rate was not in accordance with the expected or observed heterozygosities in the population or with linkage disequilibrium observed between the repeat numbers and the haplotypes of the markers flanking the CGG. The haplotype in the chromosome in which the new mutation was found was the same as that frequently found in the Japanese fragile X chromosomes, and the variance in the CGG repeat number was wider in chromosomes with the haplotypes frequently found in the fragile X chromosome than in those with the other haplotypes. These observations suggest that a subgroup is present in normal alleles and that this subgroup is more liable to mutate than others.     

The influence of evolutionary distance between cross-species microsatellites and primer base-pair composition on allelic dropout rates

Allelic dropouts (ADO) are an important source of genotyping error and because of their negative impact on non-invasive sampling techniques, have become the focus of considerable attention. Previous studies have noted that ADO rates are greater with increasing allele size and in tetranucleotides. It has also been suggested, but not tested, that ADO rates may be higher in studies using cross-species microsatellites and that mutations may play a role in ADO rates. Here we examine the relationship between ADO rates and the relationship between evolutionary distance since divergence time between species for which the microsatellite was designed for and species on which it was used (divergence times), and how this may interact with median allele size. In addition, as the adenosine (A) and thymine (T) content of the primer may increase mutation rates, we also included total % AT content of the primer in the analyses. Finally, we examined whether other commonly associated causes of ADO (i.e. repeat motif length, median allele size and allele number) co-varied. We found that ADO rates were positively associated to divergence time and median allele size. Repeat motif length, median allele size and allele number positively covaried suggesting a link between mutability and these parameters. Results from previous studies that did not correct for co-variation among these parameters may have been confounded. AT content of the primer was positively associated with ADO rates. The best linear regression model contained divergence time, median allele size and total % AT content, explaining 21% of the variation in ADO rates. The available evidence suggests that mutations partly cause ADO and that studies using cross-species microsatellites may be at higher risk of ADO. Based on our results we highlight some important considerations in the selection of microsatellites for all conservation genetic studies.