Identifying Optimal Loci for the Molecular Diagnosis of Microsatellite Instability

Background Microsatellite instability (MSI) predicts oncological response to checkpoint blockade immunotherapies. Although microsatellite mutation is pathognomonic for the condition, loci have unequal diagnostic value for predicting MSI within and across cancer types. Methods To better inform molecular diagnosis of MSI, we examined 9438 tumor-normal exome pairs and 901 whole genome sequence pairs from 32 different cancer types and cataloged genome-wide microsatellite instability events. Using a statistical framework, we identified microsatellite mutations that were predictive of MSI within and across cancer types. The diagnostic accuracy of different subsets of maximally informative markers was estimated computationally using a dedicated validation set. Results Twenty-five cancer types exhibited hypermutated states consistent with MSI. Recurrently mutated microsatellites associated with MSI were identifiable in 15 cancer types, but were largely specific to individual cancer types. Cancer-specific microsatellite panels of 1 to 7 loci were needed to attain ≥95% diagnostic sensitivity and specificity for 11 cancer types, and in 8 of the cancer types, 100% sensitivity and specificity were achieved. Breast cancer required 800 loci to achieve comparable performance. We were unable to identify recurrent microsatellite mutations supporting reliable MSI diagnosis in ovarian tumors. Features associated with informative microsatellites were cataloged. Conclusions Most microsatellites informative for MSI are specific to particular cancer types, requiring the use of tissue-specific loci for optimal diagnosis. Limited numbers of markers are needed to provide accurate MSI diagnosis in most tumor types, but it is challenging to diagnose breast and ovarian cancers using predefined microsatellite locus panels.

Engines of change: Transposable element mutation rates are high and vary widely among genotypes and populations of Daphnia magna

Transposable elements (TEs) represent a large and dynamic portion of most eukaryotic genomes, yet little is known about their mutation rates or the correspondence between rates and long-term patterns of accrual. We compare TE activity over long and short time periods by quantifying TE profiles and mutation rates (with and without minimizing selection) among 9 genotypes from three populations of Daphnia magna sampled along a latitudinal gradient. The patterns of genome-wide variation observed in nature mirror direct estimates of rates and spectra observed in a multi-year laboratory mutation accumulation experiment, where net rates range from -11.98 to 12.79 x 10-5 per copy per generation across genotypes. Overall, gains outnumber losses and both types of events are highly deleterious based on comparing lines with and without selection minimized. The rate and spectrum of TE mutations vary widely among genotypes and across TE families/types, even within the same population. We compare TE mutation rates to previously published rates of base substitution, microsatellite mutation, and gene conversion for the same genotypes, and show a correlation only with the latter. Our study provides strong evidence for the notion that TEs represent a highly mutagenic force in the genome. Furthermore, the variation we observe underscores the need to expand the repertoire of mutations studied to include a wider array of mutation types with different underlying mechanisms in order to better understand the evolution of the mutation rate and the ways in which genetic variation is generated genome wide.

A family study to examine clonal diversity in unisexual salamanders (genus Ambystoma)

Unisexual Ambystoma are the oldest known unisexual vertebrates and comprise a lineage of eastern North American all female salamanders that reproduce by stealing sperm from as many as five normally bisexual congeneric species. The sperm may be used to only stimulate egg development by gynogenesis but can be incorporated in the zygote to elevate the ploidy level or to replace one of the female’s haploid genomes. This flexible and unique reproductive system, termed kleptogenesis, is investigated using a microsatellite examination of 988 offspring from 14 unisexual mothers. All mothers produced clonal and ploidy-elevated offspring. Genome replacement and multiple paternity are confirmed for the first time in unisexual Ambystoma. Microsatellite mutations were found in all five microsatellite loci and the estimated microsatellite mutation rate varied by locus and by genome. Clonal variation is attributed to the inclusion of sperm donors’ haploid genomes for ploidy elevation, genome replacement, mutations, and natural selection.

Mode and Tempo of Microsatellite Length Change in a Malaria Parasite Mutation Accumulation Experiment

Malaria parasites have small extremely AT-rich genomes: microsatellite repeats (1–9 bp) comprise 11% of the genome and genetic variation in natural populations is dominated by repeat changes in microsatellites rather than point mutations. This experiment was designed to quantify microsatellite mutation patterns in Plasmodium falciparum. We established 31 parasite cultures derived from a single parasite cell and maintained these for 114–267 days with frequent reductions to a single cell, so parasites accumulated mutations during ∼13,207 cell divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Microsatellite calls were 99.89% (GATK), 99.99% (freeBayes), and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries, whereas introduction of microsatellite mutations into the reference genome revealed a low false negative calling rate (0.68%). We observed 98 microsatellite mutations. We highlight several conclusions: microsatellite mutation rates (3.12 × 10−7 to 2.16 × 10−8/cell division) are associated with both repeat number and repeat motif like other organisms studied. However, 41% of changes resulted from loss or gain of more than one repeat: this was particularly true for long repeat arrays. Unlike other eukaryotes, we found no insertions or deletions that were not associated with repeats or homology regions. Overall, microsatellite mutation rates are among the lowest recorded and comparable to those in another AT-rich protozoan (Dictyostelium). However, a single infection (>1011 parasites) will still contain over 2.16 × 103 to 3.12 × 104 independent mutations at any single microsatellite locus.

Quantifying Microsatellite Mutation Rates from Intestinal Stem Cell Dynamics in Msh2-Deficient Murine Epithelium

Microsatellite sequences have an enhanced susceptibility to mutation, and can act as sentinels indicating elevated mutation rates and increased risk of cancer. The probability of mutant fixation within the intestinal epithelium is dictated by a combination of stem cell dynamics and mutation rate. Here, we exploit this relationship to infer microsatellite mutation rates. First a sensitive, multiplexed, and quantitative method for detecting somatic changes in microsatellite length was developed that allowed the parallel detection of mutant [CA]n sequences from hundreds of low-input tissue samples at up to 14 loci. The method was applied to colonic crypts in Mus musculus, and enabled detection of mutant subclones down to 20% of the cellularity of the crypt (∼50 of 250 cells). By quantifying age-related increases in clone frequencies for multiple loci, microsatellite mutation rates in wild-type and Msh2-deficient epithelium were established. An average 388-fold increase in mutation per mitosis rate was observed in Msh2-deficient epithelium (2.4 × 10−2) compared to wild-type epithelium (6.2 × 10−5).

Intraspecific Variation in Microsatellite Mutation Profiles in Daphnia magna

Microsatellite loci (tandem repeats of short nucleotide motifs) are highly abundant in eukaryotic genomes and often used as genetic markers because they can exhibit variation both within and between populations. Although widely recognized for their mutability and utility, the mutation rates of microsatellites have only been empirically estimated in a few species, and have rarely been compared across genotypes and populations within a species. Here, we investigate the dynamics of microsatellite mutation over long- and short-time periods by quantifying the starting abundance and mutation rates for microsatellites for six different genotypes of Daphnia magna, an aquatic microcrustacean, collected from three populations (Finland, Germany, and Israel). Using whole-genome sequences of these six starting genotypes, descendent mutation accumulation (MA) lines, and large population controls (non-MA lines), we find each genotype exhibits a distinctive initial microsatellite profile which clusters according to the population-of-origin. During the period of MA, we observe motif-specific, highly variable, and rapid microsatellite mutation rates across genotypes of D. magna, the average of which is order of magnitude greater than the recently reported rate observed in a single genotype of the congener, Daphnia pulex. In our experiment, genotypes with more microsatellites starting out exhibit greater losses and those with fewer microsatellites starting out exhibit greater gains—a context-dependent mutation bias that has not been reported previously. We discuss how genotype-specific mutation rates and spectra, in conjunction with evolutionary forces, can shape both the differential accumulation of repeat content in the genome and the evolution of mutation rates.

Mode and tempo of microsatellite length change in a malaria parasite mutation accumulation experiment

Microsatellite sequences are widely assumed to evolve neutrally, but also play an important role in bacterial pathogenesis, human disease and transcript abundance. The malaria parasite Plasmodium falciparum genome is extraordinarily AT-rich, containing 132,449 microsatellites-stretches of perfect 1-9 bp repeats between 10-1000bp, which comprise 10.74% of the 23 Mb genome. This project was designed to determine the mode and tempo of microsatellite mutations in malaria parasites. We maintained 31 parasite lines derived from a single 3D7 parasite cell for 114-267 days, with frequent bottlenecking to a single cell to minimize effective population size, allowing us to measure mutations accumulated over ~13,207 mitotic divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Calls were 99.89% (GATK), 99.99% (freeBayes) and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries. We observed 98 microsatellite mutations, giving rates of 2.11 × 10-7 - 1.46 × 10-8 /cell division that were strongly influenced by repeat motif and array length. Mutation rate was low relative to other organisms. However, despite this, in a single infection (1011 parasites) there will be 1.46 × 103 - 2.11 × 104 independent mutations at any single microsatellite locus. Given that many microsatellites are found in promotors, introns, within or close to coding sequences, we suggest that they may be important regulators of transcriptional and phenotypic variation in this pathogen.Author summaryMutation is central to evolution: in pathogens, the rate of mutation may determine how rapidly drug resistance evolves or how effectively pathogens can escape immune attack. Malaria parasites have small extremely AT-rich genomes, and genetic variation in natural populations is dominated by repeat number changes in short tandem repeats (microsatellites) rather than point mutations. We therefore focused on quantifying microsatellite mutation. We established 31 parasite cultures in the laboratory all derived from a single parasite cell. These were maintained for 114-267 days with frequent reductions to a single cell, so parasites accumulated mutations during ~13,207 cell divisions. We sequenced the parasite genomes at the end of the experiment to count the mutations. We highlight several conclusions: like other organisms studied, microsatellite mutation rates are associated with both repeat number and repeat motif. However, 41% of changes resulted from loss or gain of more than one repeat: this was particularly true for long repeat arrays. Unlike other eukaryotes, we found no insertions or deletions that were not associated with repeats or homology regions. Overall, we found that microsatellite mutation rates in malaria were amongst the lowest recorded and comparable to those in another AT-rich protozoan (the slime mold Dictyostelium).

Intraspecific Variation in Microsatellite Mutation Profiles in Daphnia magna

Microsatellite loci (tandem repeats of short nucleotide motifs) are highly abundant in eukaryotic genomes and are often used as genetic markers because they can exhibit variation both within and between populations. Although widely recognized for their mutability and utility, the mutation rates of microsatellites have only been empirically estimated in a few species and have rarely been compared across genotypes and populations and intraspecific differences in overall microsatellite content have rarely been explored. To investigate the accumulation of microsatellite DNA over long-and short-time periods, we quantified the abundance and genome-wide mutation rates in whole-genome sequences of 47 mutation accumulation (MA) lines and 12 non-MA lines derived from six different genotypes of the crustacean Daphnia magna collected from three populations (Finland, Germany, and Israel). Each genotype possessed a distinctive microsatellite profile and clustered according to their population of origin. During the period of mutation accumulation, we observed very high microsatellite mutation rates (a net change of −0.19 to 0.33 per copy per generation), which surpass rates reported from a closely-related congener, D. pulex, by an order of magnitude. Rates vary between microsatellite motifs and among genotypes, with those starting with high microsatellite content exhibiting greater losses and those with low microsatellite content exhibiting greater gains. Our results show that microsatellite mutation rates depend both on characteristics of the microsatellites and the genomic background. These context-dependent mutation dynamics may, in conjunction with other evolutionary forces that may differ among populations, explain the differential accumulation of repeat content in the genome over long time periods.