A brief review of short tandem repeat mutation

Short tandem repeats （STRs） are short tandemly repeated DNA sequences that involve a repetitive unit of 1-6 bp. Because of their polymorphisms and high mutation rates, STRs are widely used in biological research. Strand-slippage replication is the predominant mutation mechanism of STRs, and the stepwise mutation model is regarded as the main mutation model. STR mutation rates can be influenced by many factors. Moreover, some trinucleotide repeats are associated with human neurodegenerative diseases. In order to deepen our knowledge of these diseases and broaden STR application, it is essential to understand the STR mutation process in detail. In this review, we focus on the current known information about STR mutation.     

The Conserved PFT1 Tandem Repeat Is Crucial for Proper Flowering in Arabidopsis thaliana

It is widely appreciated that short tandem repeat (STR) variation underlies substantial phenotypic variation in organisms. Some propose that the high mutation rates of STRs in functional genomic regions facilitate evolutionary adaptation. Despite their high mutation rate, some STRs show little to no variation in populations. One such STR occurs in the Arabidopsis thaliana gene PFT1 (MED25), where it encodes an interrupted polyglutamine tract. Although the PFT1 STR is large (∼270 bp), and thus expected to be extremely variable, it shows only minuscule variation across A. thaliana strains. We hypothesized that the PFT1 STR is under selective constraint, due to previously undescribed roles in PFT1 function. We investigated this hypothesis using plants expressing transgenic PFT1 constructs with either an endogenous STR or synthetic STRs of varying length. Transgenic plants carrying the endogenous PFT1 STR generally performed best in complementing a pft1 null mutant across adult PFT1-dependent traits. In stark contrast, transgenic plants carrying a PFT1 transgene lacking the STR phenocopied a pft1 loss-of-function mutant for flowering time phenotypes and were generally hypomorphic for other traits, establishing the functional importance of this domain. Transgenic plants carrying various synthetic constructs occupied the phenotypic space between wild-type and pft1 loss-of-function mutants. By varying PFT1 STR length, we discovered that PFT1 can act as either an activator or repressor of flowering in a photoperiod-dependent manner. We conclude that the PFT1 STR is constrained to its approximate wild-type length by its various functional requirements. Our study implies that there is strong selection on STRs not only to generate allelic diversity, but also to maintain certain lengths pursuant to optimal molecular function.     

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.     

STR polymorphism of mtDNA D-loop in rhesus macaques of Bangladesh

Molecular variation of mitochondrial DNA (mtDNA) was investigated for rhesus macaques (Macaca mulatta) of Bangladesh. A partial sequence (583–599 bp) of mtDNA containing the second variable region of the D-loop was compared for 39 individuals from five localities in the country. A total of seven haplotypes were detected with substitutional or insertion/deletion mutations. They contained a unique polymorphism of pentanucleotide STRs (short tandem repeats). There were at least four different length types, from two to five repeats of the unit nucleotide. One site of substitution and one site of single nucleotide insertion/deletion were also involved in the polymorphism. The mutation hot spots of the STR polymorphism were located between the first and second conserved sequence blocks (CSB1 and CSB2), as observed previously in some other mammals. The geographical distribution of the STR polymorphism revealed local differences; the northeastern population was polymorphic with three STR haplotypes, but other local populations were simply monomorphic with a single STR haplotype. Molecular phylogenetic analysis with reported sequences from outside Bangladesh indicated a low substitution diversity of mtDNA in Bangladesh. Clustering results suggested a close relationship to India and divergence from Laos and China.     

Mitochondrial DNA variation and reticulate evolution of the genus Abies

The sequences of three regions of mitochondrial DNA (mtDNA) of a total length of 5226 bp were used to study the phylogeography of the genus Abies. The mtDNA haplotype network, comprising 36 studied Abies taxa, consisted of two branches; the first represented all American species plus two Asian, and the second included the remaining Eurasian species. Within these clusters, the haplotypes formed nine major groups, generally corresponding to the clades of the previously obtained phylogeny based on chloroplast DNA (cpDNA), but the relationships of these groups were significantly different; species assignment to the particular mtDNA haplotype group was more in line with its geographical distribution. In addition, the mtDNA haplotype network contains cycles indicating the recombination. It is assumed that the incongruence of cpDNA and mtDNA phylogenies is caused by the introgression capture of alien mtDNA during species hybridization and thus contains information about past migrations. The cases of incongruence of mitochondrial and chloroplast DNA suggesting a migration of Abies between Asia and North America are discussed.     

Mutation analysis of 28 autosomal short tandem repeats in the Chinese Han population

Short tandem repeats (STRs) have been extensively used in forensic genetics. However, according to previous studies, the mutation rates of STRs are relatively high and are affected by many factors. Therefore, it is important to analyze STR mutations and determine the influence of underlying factors on STR mutation rates. Mutation rates of 28 autosomal STRs were determined from 8708 paternity testing cases in the Chinese Han population, and the relationships between STR mutation rates and population, sex, age, allele length and heterozygosity were investigated. A total of 279 mutations were observed at 27 loci in a total of 233,530 meiosis cases, including 273 (97.8%) one-step, 5 (1.8%) two-step and 1 (0.4%) three-step mutations. The overall average mutation rate was 1.19?×?10^–3 (95% CI 1.06?×?10^–3???1.34?×?10^–3) ranging from 0 (TPOX) to 2.79?×?10^–3 (D13S325). Mutation rate comparisons revealed statistically significant differences at several STRs among populations. Paternal mutations occurred more frequently than maternal mutations, at a ratio of 6.04:1, and the mutation rate tended to increase with paternal age. Moreover, our study revealed a bias towards contraction mutations for long alleles and expansion mutations for short alleles. No obvious bias was observed in the overall mutation direction. In addition, STR loci with higher expected heterozygosity (H_exp) tended to have higher mutation rates. This work revealed the relationships between STR mutation rates and several influencing factors, providing useful data and information for further research on STR mutations in forensic genetics.

     

Lineage-specific group II intron gains and losses of the mitochondrial rps3 gene in gymnosperms

According to PCR assays and sequencing, we now report the shared presence of two rps3 introns, namely the rps3i74 and the rps3i249, in the mitochondria of all the classes representing the surviving lineages of gymnosperms, and unveil several lineages experiencing intron loss.Interestingly, the rps3 intron gains and losses within the four groups of gymnosperms let us sort out the Pinaceae and the non-Pinaceae into intron (+)- and intron (?)-lineages, respectively. Worthy of mention is also the finding that only Gnetum within the Gnetales harbours both the rps3 introns.This intron distribution pattern is consistent with the hypothesis that the two rps3 introns were likely present in the common ancestor of the seed plants and, then, independently lost in the non-Pinaceae during gymnosperm evolution.The derived secondary structural model of the novel group IIA intron improves our understanding of the significance and origin of the extraordinary length polymorphisms observed among rps3i249 orthologs.Despite the remarkable structural plasticity to adopt and reject introns, the rps3 mRNAs undergo accurate processing by splicing and extensive editing in gymnosperm mitochondria.This study provides additional insights into the evolutionarily high dynamics of mitochondrial introns which may come and go in closely related plant species.The turnover of the mitochondrial rps3 group II introns seen among lineages of seed plants further suggests that these introns might be an additional signature to discriminate between particularly cryptical taxonomic groups for which there is a need of a further evaluation of their evolutionary affiliation.     

Mitochondrial DNA exhibits resistance to induced point and deletion mutations

The accumulation of somatic mitochondrial DNA (mtDNA) mutations contributes to the pathogenesis of human disease. Currently, mitochondrial mutations are largely considered results of inaccurate processing of its heavily damaged genome. However, mainly from a lack of methods to monitor mtDNA mutations with sufficient sensitivity and accuracy, a link between mtDNA damage and mutation has not been established. To test the hypothesis that mtDNA-damaging agents induce mtDNA mutations, we exposed Muta^TMMouse mice to benzo[a]pyrene (B[a]P) or N-ethyl-N-nitrosourea (ENU), daily for 28 consecutive days, and quantified mtDNA point and deletion mutations in bone marrow and liver using our newly developed Digital Random Mutation Capture (dRMC) and Digital Deletion Detection (3D) assays. Surprisingly, our results demonstrate mutagen treatment did not increase mitochondrial point or deletion mutation frequencies, despite evidence both compounds increase nuclear DNA mutations and demonstrated B[a]P adduct formation in mtDNA. These findings contradict models of mtDNA mutagenesis that assert the elevated rate of mtDNA mutation stems from damage sensitivity and abridged repair capacity. Rather, our results demonstrate induced mtDNA damage does not readily convert into mutation. These findings suggest robust mitochondrial damage responses repress induced mutations after mutagen exposure.     

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1–6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.     

Discordant mtDNA and cpDNA phylogenies indicate geographic speciation and reticulation as driving factors for the diversification of the genus Picea

The phylogeny of the genus Picea was investigated by sequencing three loci from the paternally inherited chloroplast genome (trnK, rbcL and trnTLF) and the intron 2 of the maternally transmitted mitochondrial gene nad1 for 35 species. Significant topological differences were found between the trnK tree and the rbcL and trnTLF phylogenetic trees, and between cpDNA and mtDNA phylogenies. None of the phylogenies matched morphological classifications. The mtDNA phylogeny was geographically more structured than cpDNA phylogenies, reflecting the different inheritance of the two cytoplasmic genomes in the Pinaceae and their differential dispersion by seed only and seed and pollen, respectively. Most North American taxa formed a monophyletic group on the mtDNA tree, with topological patterns suggesting geographic speciation by range fragmentation or by dispersal and isolation. Similar patterns were also found among Asian taxa. Such a trend towards geographic speciation is anticipated in other Pinaceae genera with similar life history, autecology and reproductive system. Incongruences between organelle phylogenies suggested the occurrence of mtDNA capture by invading cpDNA. Incongruences between cpDNA partitions further suggested heterologous recombination presumably also linked to ancient reticulate evolution. Whilst cpDNA appears potentially valuable for molecular taxonomy and systematics purposes, these results emphasize the reduced value of cpDNA to infer vertical descent and the speciation history for plants with paternal transmission and high dispersal of their chloroplast genome.     

Oxidative Stress Is Not a Major Contributor to Somatic Mitochondrial DNA Mutations

The accumulation of somatic mitochondrial DNA (mtDNA) mutations is implicated in aging and common diseases of the elderly, including cancer and neurodegenerative disease. However, the mechanisms that influence the frequency of somatic mtDNA mutations are poorly understood. To develop a simple invertebrate model system to address this matter, we used the Random Mutation Capture (RMC) assay to characterize the age-dependent frequency and distribution of mtDNA mutations in the fruit fly Drosophila melanogaster. Because oxidative stress is a major suspect in the age-dependent accumulation of somatic mtDNA mutations, we also used the RMC assay to explore the influence of oxidative stress on the somatic mtDNA mutation frequency. We found that many of the features associated with mtDNA mutations in vertebrates are conserved in Drosophila, including a comparable somatic mtDNA mutation frequency (∼10⁻⁵), an increased frequency of mtDNA mutations with age, and a prevalence of transition mutations. Only a small fraction of the mtDNA mutations detected in young or old animals were G∶C to T∶A transversions, a signature of oxidative damage, and loss-of-function mutations in the mitochondrial superoxide dismutase, Sod2, had no detectable influence on the somatic mtDNA mutation frequency. Moreover, a loss-of-function mutation in Ogg1, which encodes a DNA repair enzyme that removes oxidatively damaged deoxyguanosine residues (8-hydroxy-2′-deoxyguanosine), did not significantly influence the somatic mtDNA mutation frequency of Sod2 mutants. Together, these findings indicate that oxidative stress is not a major cause of somatic mtDNA mutations. Our data instead suggests that somatic mtDNA mutations arise primarily from errors that occur during mtDNA replication. Further studies using Drosophila should aid in the identification of factors that influence the frequency of somatic mtDNA mutations.     

Ancient DNA analysis of human neolithic remains found in northeastern Siberia

We successfully extracted DNA from a bone sample of a Neolithic skeleton (dated 3,600 +/- 60 years BP) excavated in northeastern Yakutia (east Siberia). Ancient DNA was analyzed by autosomal STRs (short tandem repeats) and by sequencing of the hypervariable region I (HV1) of the mitochondrial DNA (mtDNA) control region. The STR profile, the mitochondrial haplotype, and the haplogroup determined were compared with those of modern Eurasian and Native American populations. The results showed the affinity of this ancient skeleton with both east Siberian/Asian and Native American populations.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

Fast evolution of the retroprocessed mitochondrial rps3 gene in Conifer II and further evidence for the phylogeny of gymnosperms

The popular view that plant mitochondrial genome evolves slowly in sequence has been recently challenged by the extraordinarily high substitution rates of mtDNA documented mainly from several angiosperm genera, but high substitution rate acceleration accompanied with great length variation has been very rarely reported in plant mitochondrial genes. Here, we studied evolution of the mitochondrial rps3 gene that encodes the ribosomal small subunit protein 3 and found a dramatically high variation in both length and sequence of an exon region of it in Conifer II. A sequence comparison between cDNA and genomic DNA showed that there are no RNA editing sites in the Conifer II rps3 gene. Southern blotting analyses of the total DNA and mtDNA, together with the real-time PCR analysis, showed that rps3 exists as a single mitochondrial locus in gymnosperms. It is very likely that the Conifer II rps3 gene has experienced retroprocessing, i.e., the re-integration of its cDNA into the mitochondrial genome, followed by an evolutionary acceleration due to the intron loss. In addition, the phylogenetic analysis of rps3 supports the sister relationship between conifers and Gnetales. In particular, the monophyly of conifer II is strongly supported by the shared loss of two rps3 introns. Our results also indicate that the mitochondrial gene tree would be affected in topology when the “edited” paralogs are analyzed together with their genomic sequences.     

The mitochondrial genome of Chara vulgaris: insights into the mitochondrial DNA architecture of the last common ancestor of green algae and land plants

Mitochondrial DNA (mtDNA) has undergone radical changes during the evolution of green plants, yet little is known about the dynamics of mtDNA evolution in this phylum. Land plant mtDNAs differ from the few green algal mtDNAs that have been analyzed to date by their expanded size, long spacers, and diversity of introns. We have determined the mtDNA sequence of Chara vulgaris (Charophyceae), a green alga belonging to the charophycean order (Charales) that is thought to be the most closely related alga to land plants. This 67,737-bp mtDNA sequence, displaying 68 conserved genes and 27 introns, was compared with those of three angiosperms, the bryophyte Marchantia polymorpha, the charophycean alga Chaetosphaeridium globosum (Coleochaetales), and the green alga Mesostigma viride. Despite important differences in size and intron composition, Chara mtDNA strikingly resembles Marchantia mtDNA; for instance, all except 9 of 68 conserved genes lie within blocks of colinear sequences. Overall, our genome comparisons and phylogenetic analyses provide unequivocal support for a sister-group relationship between the Charales and the land plants. Only four introns in land plant mtDNAs appear to have been inherited vertically from a charalean algar ancestor. We infer that the common ancestor of green algae and land plants harbored a tightly packed, gene-rich, and relatively intron-poor mitochondrial genome. The group II introns in this ancestral genome appear to have spread to new mtDNA sites during the evolution of bryophytes and charalean green algae, accounting for part of the intron diversity found in Chara and land plant mitochondria.     

The evolutionary split of Pinaceae from other conifers: evidence from an intron loss and a multigene phylogeny.

The second intron in the mitochondrial gene nad1 was surveyed using PCR, DNA sequencing, or Southern hybridization in 323 species (313 genera, 212 families) of seed plants. The intron was absent in all 22 species (22 genera, 8 families) of non-Pinaceae conifers studied, in Welwitschia mirabilis, and in seven angiosperms. Whereas absence of the intron in seven angiosperms and Welwitschia is likely due to seven independent losses when evaluated against the recently published multigene phylogenies, the lack of the intron in all non-Pinaceae conifers can be best explained by a single loss. These data suggest that the non-Pinaceae conifers represent a monophyletic group. We also conducted a phylogenetic analysis of seed plants using a combined data set of the partial exon and intron sequences of nad1 generated from this study and published sequences of mitochondrial cox1 and small subunit (SSU) rDNA, chloroplast rbcL, and nuclear 18S rDNA. The results supported the split of conifers into two groups: Pinaceae and non-Pinaceae conifers. The Gnetales were sister to Pinaceae, in agreement with the conclusion from other recent molecular phylogenetic studies that refute the anthophyte hypothesis.     

Molecular phylogenetic analysis of the genus Abies (Pinaceae) based on the nucleotide sequence of chloroplast DNA

A phylogenetic study of firs (Abies Mill.) was conducted using nucleotide sequences of several chloroplast DNA regions with a total length of 5580 bp. The analysis included 37 taxa, which represented the main evolutionary lineages of the genus, and Keteleeria davidiana. According to phylogenetic reconstruction, the Abies species were subdivided into six main groups, generally corresponding to their geographic distribution. The phylogenetic tree had three basal clades. All of these clades contained American species, and only one of them contained Eurasian species. The divergence time calibrations, based on paleobotanical data and the chloroplast DNA mutation rate estimates in Pinaceae, produced similar results. The age of diversification among the basal clades of the present-day Abies was estimated as the end of the Oligocene-beginning of Miocene. The age of the separation of Mediterranean firs from the Asian-North American branch corresponds to the Miocene. The age of diversification within the young groups of Mediterranean, Asian, and “boreal” American firs (A. lasiocarpa, A. balsamea, A. fraseri) was estimated as the Pliocene-Pleistocene. Based on the phylogenetic reconstruction obtained, the most plausible biogeographic scenarios were suggested. It is noted that the existing systematic classification of the genus Abies strongly contradicts with phylogenetic reconstruction and requires revision.     

Mathematical Modeling of the Role of Mitochondrial Fusion and Fission in Mitochondrial DNA Maintenance

Accumulation of mitochondrial DNA (mtDNA) mutations has been implicated in a wide range of human pathologies, including neurodegenerative diseases, sarcopenia, and the aging process itself. In cells, mtDNA molecules are constantly turned over (i.e. replicated and degraded) and are also exchanged among mitochondria during the fusion and fission of these organelles. While the expansion of a mutant mtDNA population is believed to occur by random segregation of these molecules during turnover, the role of mitochondrial fusion-fission in this context is currently not well understood. In this study, an in silico modeling approach is taken to investigate the effects of mitochondrial fusion and fission dynamics on mutant mtDNA accumulation. Here we report model simulations suggesting that when mitochondrial fusion-fission rate is low, the slow mtDNA mixing can lead to an uneven distribution of mutant mtDNA among mitochondria in between two mitochondrial autophagic events leading to more stochasticity in the outcomes from a single random autophagic event. Consequently, slower mitochondrial fusion-fission results in higher variability in the mtDNA mutation burden among cells in a tissue over time, and mtDNA mutations have a higher propensity to clonally expand due to the increased stochasticity. When these mutations affect cellular energetics, nuclear retrograde signalling can upregulate mtDNA replication, which is expected to slow clonal expansion of these mutant mtDNA. However, our simulations suggest that the protective ability of retrograde signalling depends on the efficiency of fusion-fission process. Our results thus shed light on the interplay between mitochondrial fusion-fission and mtDNA turnover and may explain the mechanism underlying the experimentally observed increase in the accumulation of mtDNA mutations when either mitochondrial fusion or fission is inhibited.