Base composition of mononucleotide runs affects DNA polymerase slippage and removal of frameshift intermediates by mismatch repair in Saccharomyces cerevisiae

Characteristics and Comparative Analysis of Seven Complete Plastomes of Trichoglottis s.l. (Aeridinae, Orchidaceae)

Trichoglottis exhibits a range of rich variations in colors and shapes of flower and is a valuable ornamental orchid genus. The genus Trichoglottis has been expanded by the inclusion of Staurochilus, but this Trichoglottis sensu lato (s.l.) was recovered as a non-monophyletic genus based on molecular sequences from one or a few DNA regions. Here, we present phylogenomic data sets, incorporating complete plastome sequences from seven species (including five species sequenced in this study) of Trichoglottis s.l. (including two species formerly treated as Staurochilus), to compare plastome structure and to reconstruct the phylogenetic relationships of this genus. The seven plastomes possessed the typical quadripartite structure of angiosperms and ranged from 149,402 bp to 149,841 bp with a GC content of 36.6-36.7%. These plastomes contain 120 genes, which comprise 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes, all ndh genes were pseudogenized or lost. A total of 98 (T. philippinensis) to 134 (T. ionosma) SSRs and 33 (T. subviolacea) to 46 (T. ionosma) long repeats were detected. The consistent and robust phylogenetic relationships of Trichoglottis were established using a total of 25 plastid genomes from the Aeridinae subtribe. The genus Trichoglottis s.l. was strongly supported as a monophyletic group, and two species formerly treated as Staurochilus were revealed as successively basal lineages. In addition, five mutational hotspots (trnNGUU-rpl32, trnLUAA, trnSGCU-trnGUCC, rbcL-accD, and trnTGGU-psbD) were identified based on the ranking of PI values. Our research indicates that plastome data is a valuable source for molecular identification and evolutionary studies of Trichoglottis and its related genera.

Characteristics of plastid genomes in the genus Ceratostigma inhabiting arid habitats in China and their phylogenomic implications

BACKGROUND

Ceratostigma, a genus in the Plumbaginaceae, is an ecologically dominant group of shrubs, subshrub and herb mainly distributed in Qinghai-Tibet Plateau and North China. Ceratostigma has been the focal group in several studies, owing to their importance in economic and ecological value and unique breeding styles. Despite this, the genome information is limited and interspecific relationships within the genus Cerotastigma remains unexplored. Here we sequenced, assembled and characterized the 14 plastomes of five species, and conducted phylogenetic analyses of Cerotastigma using plastomes and nuclear ribosomal DNA (nrDNA) data.

RESULTS

Fourteen Cerotastigma plastomes possess typical quadripartite structures with lengths from 164,076 to 168,355 bp that consist of a large single copy, a small single copy and a pair of inverted repeats, and contain 127-128 genes, including 82-83 protein coding genes, 37 transfer RNAs and eight ribosomal RNAs. All plastomes are highly conservative and similar in gene order, simple sequence repeats (SSRs), long repeat repeats and codon usage patterns, but some structural variations in the border of single copy and inverted repeats. Mutation hotspots in coding (Pi values > 0.01: matK, ycf3, rps11, rps3, rpl22 and ndhF) and non-coding regions (Pi values > 0.02: trnH-psbA, rps16-trnQ, ndhF-rpl32 and rpl32-trnL) were identified among plastid genomes that could be served as potential molecular markers for species delimitation and genetic variation studies in Cerotastigma. Gene selective pressure analysis showed that most protein-coding genes have been under purifying selection except two genes. Phylogenetic analyses based on whole plastomes and nrDNA strongly support that the five species formed a monophyletic clade. Moreover, interspecific delimitation was well resolved except C. minus, individuals of which clustered into two main clades corresponding to their geographic distributions. The topology inferred from the nrDNA dataset was not congruent with the tree derived from the analyses of the plastid dataset.

CONCLUSION

These findings represent the first important step in elucidating plastome evolution in this widespread distribution genus Cerotastigma in the Qinghai-Tibet Plateau. The detailed information could provide a valuable resource for understanding the molecular dynamics and phylogenetic relationship in the family Plumbaginaceae. Lineage genetic divergence within C. minus was perhaps promoted by geographic barriers in the Himalaya and Hengduan Mountains region, but introgression or hybridization could not be completely excluded.

Comparative analysis of complete Artemisia subgenus Seriphidium (Asteraceae: Anthemideae) chloroplast genomes: insights into structural divergence and phylogenetic relationships

BACKGROUND

Artemisia subg. Seriphidium, one of the most species-diverse groups within Artemisia, grows mainly in arid or semi-arid regions in temperate climates. Some members have considerable medicinal, ecological, and economic value. Previous studies on this subgenus have been limited by a dearth of genetic information and inadequate sampling, hampering our understanding of their phylogenetics and evolutionary history. We therefore sequenced and compared the chloroplast genomes of this subgenus, and evaluated their phylogenetic relationships.

RESULTS

We newly sequenced 18 chloroplast genomes of 16 subg. Seriphidium species and compared them with one previously published taxon. The chloroplast genomes, at 150,586-151,256 bp in length, comprised 133 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and one pseudogene, with GC content of 37.40-37.46%. Comparative analysis showed that genomic structures and gene order were relatively conserved, with only some variation in IR borders. A total of 2203 repeats (1385 SSRs and 818 LDRs) and 8 highly variable loci (trnK - rps16, trnE - ropB, trnT, ndhC - trnV, ndhF, rpl32 - trnL, ndhG - ndhI and ycf1) were detected in subg. Seriphidium chloroplast genomes. Phylogenetic analysis of the whole chloroplast genomes based on maximum likelihood and Bayesian inference analyses resolved subg. Seriphidium as polyphyletic, and segregated into two main clades, with the monospecific sect. Minchunensa embedded within sect. Seriphidium, suggesting that the whole chloroplast genomes can be used as molecular markers to infer the interspecific relationship of subg. Seriphidium taxa.

CONCLUSION

Our findings reveal inconsistencies between the molecular phylogeny and traditional taxonomy of the subg. Seriphidium and provide new insights into the evolutionary development of this complex taxon. Meanwhile, the whole chloroplast genomes with sufficiently polymorphic can be used as superbarcodes to resolve interspecific relationships in subg. Seriphidium.

The rate of spontaneous mutations in yeast deficient for MutSβ function

Mutations in simple sequence repeat loci underlie many inherited disorders in humans, and are increasingly recognized as important determinants of natural phenotypic variation. In eukaryotes, mutations in these sequences are primarily repaired by the MutSβ mismatch repair complex. To better understand the role of this complex in mismatch repair and the determinants of simple sequence repeat mutation predisposition, we performed mutation accumulation in yeast strains with abrogated MutSβ function. We demonstrate that mutations in simple sequence repeat loci in the absence of mismatch repair are primarily deletions. We also show that mutations accumulate at drastically different rates in short (<8 bp) and longer repeat loci. These data lend support to a model in which the mismatch repair complex is responsible for repair primarily in longer simple sequence repeats.

Histopathologic and Genetic Features of Mismatch Repair-Deficient High-Grade Prostate Cancer

AIMS

Mismatch repair (MMR) deficiency is commonly caused by functional inactivation of MLH1, PMS2, MSH2 or MSH6. The morphologic and molecular correlates of MMR deficiency have been extensively characterized in certain tumor types such as colorectal and endometrial adenocarcinoma. In contrast, the histologic and molecular features of MMR-deficient prostate cancer remain incompletely described. In this study, we evaluated 19 MMR-deficient prostate cancers, including 11 cases without prior systemic treatment.

METHODS AND RESULTS

All treatment-naïve cases (11/11, 100%) were Grade Group 4-5 and had predominant cribriform and/or solid growth patterns. Solid components (any amount) and tumor infiltrating lymphocytes were seen in 7/11 (64%) of these cases each. In 68 MMR-proficient Grade Group 5 prostate cancers, predominant cribriform or solid growth patterns, solid components (any amount) and tumor infiltrating lymphocytes were seen at significantly lower frequencies (31/68, 46%; 9/68, 13% and 6/62, 9%, respectively; p<0.001 for all comparisons). Molecular evaluation of 19 cases demonstrated that MMR-deficiency was secondary to functional loss of MSH2/MSH6 and MLH1/PMS2 in 15 cases (79%) and 4 cases (21%), respectively. Definite or likely germline mutations were present in 4 cases (4/19, 21%). TMPRSS2::ERG rearrangements were identified in 2 cases (2/19, 11%). Recurrent cancer-relevant somatic mutations included (but were not limited to) ATM, TP53, FOXA1, RB1, BRCA2 and PTEN.

CONCLUSIONS

MMR deficiency was most commonly secondary to inactivation of MSH2/MSH6 in this study. Importantly, MMR-deficient high-grade prostatic adenocarcinomas had morphologic features that might be useful to identify selected cases for MMR IHC.

Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences

BACKGROUND

Polystachya Hook. is a large pantropical orchid genus (c. 240 species) distributed in Africa, southern Asia and the Americas, with the center of diversity in Africa. Previous studies on species of this genus have not obtained the complete chloroplast genomes, structures and variations. Additionally, the phylogenetic position of the genus in the Orchidaceae is still controversial and uncertain. Therefore, in this study, we sequenced the complete plastomes of six Kenya Polystachya species based on genome skimming, subjected them to comparative genomic analysis, and reconstructed the phylogenetic relationships with other Orchidaceae species.

RESULTS

The results exhibited that the chloroplast genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The plastomes of the six Polystachya species ranged from 145,484 bp to 149,274 bp in length and had an almost similar GC content of 36.9-37.0%. Gene annotation revealed 106-109 single-copy genes. In addition, 19 genes are duplicated in the inverted regions, and 16 genes each possessd one or more introns. Although no large structural variations were observed among the Polystachya plastomes, about 1 kb inversion was found in Polystachya modesta and all 11 ndh genes in the Polystachya plastomes were lost or pseudogenized. Comparative analysis of the overall sequence identity among six complete chloroplast genomes confirmed that for both coding and non-coding regions in Polystachya, SC regions exhibit higher sequence variation than IRs. Furthermore, there were various amplifications in the IR regions among the six Polystachya species. Most of the protein-coding genes of these species had a high degree of codon preference. We screened out SSRs and found seven relatively highly variable loci. Moreover, 13 genes were discovered with significant positive selection. Phylogenetic analysis showed that the six Polystachya species formed a monophyletic clade and were more closely related to the tribe Vandeae. Phylogenetic relationships of the family Orchidaceae inferred from the 85 chloroplast genome sequences were generally consistent with previous studies and robust.

CONCLUSIONS

Our study is the initial report of the complete chloroplast genomes of the six Polystachya species, elucidates the structural characteristics of the chloroplast genome of Polystachya, and filters out highly variable sequences that can contribute to the development of DNA markers for use in the study of genetic variability and evolutionary studies in Polystachya. In addition, the phylogenetic results strongly support that the genus of Polystachya is a part of the tribe Vandeae.

Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

Mutation rate and spectrum in obligately outcrossing Caenorhabditis elegans mutation accumulation lines subjected to RNAi-induced knockdown of the mismatch repair gene msh-2

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutSα homolog, msh-2, in Caenorhabditis elegans male–female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male–female background substantially increased the frequency of nuclear base substitutions (∼23×) and small indels (∼328×) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

AT Homopolymer Strings in Salmonella enterica Subspecies I Contribute to Speciation and Serovar Diversity

Adenine and thymine homopolymer strings of at least 8 nucleotides (AT 8+mers) were characterized in Salmonella enterica subspecies I. The motif differed between other taxonomic classes but not between Salmonella enterica serovars. The motif in plasmids was possibly associated with serovar. Approximately 12.3% of the S. enterica motif loci had mutations. Mutability of AT 8+mers suggests that genomes undergo frequent repair to maintain optimal gene content, and that the motif facilitates self-recognition; in addition, serovar diversity is associated with plasmid content. A theory that genome regeneration accounts for both persistence of predominant Salmonella serovars and serovar diversity provides a new framework for investigating root causes of foodborne illness.

Evolutionary pressures and codon bias in low complexity regions of plasmodia

The biological meaning of low complexity regions in the proteins of Plasmodium species is a topic of discussion in evolutionary biology. There is a debate between selectionists and neutralists, who either attribute or do not attribute an effect of low-complexity regions on the fitness of these parasites, respectively. In this work, we comparatively study 22 Plasmodium species to understand whether their low complexity regions undergo a neutral or, rather, a selective and species-dependent evolution. The focus is on the connection between the codon repertoire of the genetic coding sequences and the occurrence of low complexity regions in the corresponding proteins. The first part of the work concerns the correlation between the length of plasmodial proteins and their propensity at embedding low complexity regions. Relative synonymous codon usage, entropy, and other indicators reveal that the incidence of low complexity regions and their codon bias is species-specific and subject to selective evolutionary pressure. We also observed that protein length, a relaxed selective pressure, and a broad repertoire of codons in proteins, are strongly correlated with the occurrence of low complexity regions. Overall, it seems plausible that the codon bias of low-complexity regions contributes to functional innovation and codon bias enhancement of proteins on which Plasmodium species rest as successful evolutionary parasites.

How does feedback from phage infections influence the evolution of phase variation in Campylobacter?

Campylobacter jejuni (C. jejuni) causes gastroenteritis following the consumption of contaminated poultry meat, resulting in a large health and economic burden worldwide. Phage therapy is a promising technique for eradicating C. jejuni from poultry flocks and chicken carcasses. However, C. jejuni can resist infections by some phages through stochastic, phase-variable ON/OFF switching of the phage receptors mediated by simple sequence repeats (SSR). While selection strength and exposure time influence the evolution of SSR-mediated phase variation (PV), phages offer a more complex evolutionary environment as phage replication depends on having a permissive host organism. Here, we build and explore several continuous culture bacteria-phage computational models, each analysing different phase-variable scenarios calibrated to the experimental SSR rates of C. jejuni loci and replication parameters for the F336 phage. We simulate the evolution of PV rates via the adaptive dynamics framework for varying levels of selective pressures that act on the phage-resistant state. Our results indicate that growth reducing counter-selection on a single PV locus results in the stable maintenance of the phage, while compensatory selection between bacterial states affects the evolutionary stable mutation rates (i.e. very high and very low mutation rates are evolutionarily disadvantageous), whereas, in the absence of either selective pressure the evolution of PV rates results in mutation rates below the basal values. Contrastingly, a biologically-relevant model with two phase-variable loci resulted in phage extinction and locking of the bacteria into a phage-resistant state suggesting that another counter-selective pressure is required, instance, the use of a distinct phage whose receptor is an F336-phage-resistant state. We conclude that a delicate balance between counter-selection and phage-attack can result in both the evolution of phase-variable phage receptors and persistence of PV-receptor-specific phage.

Globally rising rates of antibiotic resistance have renewed interest in phage therapy. Bacteriophages (phages) act on bacteria to select for resistance mechanisms such as loss of phage receptors by phase variation (PV). Phase-variable genes mediate rapid adaption by stochastic switching of gene expression. Campylobacter jejuni is a common commensal of birds but also causes serious gastrointestinal infections in humans. Optimisation of phage therapy against C. jejuni requires an in-depth understanding of how PV has evolved and mediates phage resistance. Here, we use a detailed continuous culture model for nutrient-limited bacteria-phage interactions, with PV rates calibrated to match the experimental observations for C.jejuni and phage F336. Evolution within a model accounting for two phase-variable loci closely matches the experimental results when growth reducing counter-selection is imposed on all phage-resistant states, but, not when restricted to the particular states associated with resistance to immune effectors. Our results emphasize that delicate balancing of selective pressures, imposed by single and multiple distinct phages, are necessary for effective use of phage therapy against C. jejuni.

Exact mapping of Illumina blind spots in the Mycobacterium tuberculosis genome reveals platform-wide and workflow-specific biases

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

The fitness cost of mismatch repair mutators in Saccharomyces cerevisiae: partitioning the mutational load

Mutational load is the depression in a population's mean fitness that results from the continual influx of deleterious mutations. Here, we directly estimate the mutational load in a population of haploid Saccharomyces cerevisiae that are deficient for mismatch repair. We partition the load in haploids into two components. To estimate the load due to nonlethal mutations, we measure the competitive fitness of hundreds of randomly selected clones from both mismatch-repair-deficient and -proficient populations. Computation of the mean clone fitness for the mismatch-repair-deficient strain permits an estimation of the nonlethal load, and the histogram of fitness provides an interesting visualization of a loaded population. In a separate experiment, in order to estimate the load due to lethal mutations (i.e. the lethal mutation rate), we manipulate thousands of individual pairs of mother and daughter cells and track their fates. These two approaches yield point estimates for the two contributors to load, and the addition of these estimates is nearly equal to the separately measured short-term competitive fitness deficit for the mismatch-repair-deficient strain. This correspondence suggests that there is no need to invoke direct fitness effects to explain the fitness difference between mismatch-repair-deficient and -proficient strains. Assays in diploids are consistent with deleterious mutations in diploids tending towards recessivity. These results enhance our understanding of mutational load, a central population genetics concept, and we discuss their implications for the evolution of mutation rates.

Cooperation between non-essential DNA polymerases contributes to genome stability in Saccharomyces cerevisiae

DNA polymerases influence genome stability through their involvement in DNA replication, response to DNA damage, and DNA repair processes. Saccharomyces cerevisiae possess four non-essential DNA polymerases, Pol λ, Pol η, Pol ζ, and Rev1, which have varying roles in genome stability. In order to assess the contribution of the non-essential DNA polymerases in genome stability, we analyzed the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant in microhomology mediated repair, due to recent studies linking some of these DNA polymerases to this repair pathway. Our results suggest that the length and quality of microhomology influence both the overall efficiency of repair and the involvement of DNA polymerases. Furthermore, the non-essential DNA polymerases demonstrate overlapping and redundant functions when repairing double-strand breaks using short microhomologies containing mismatches. Then, we examined genome-wide mutation accumulation in the pol4Δ rev1Δ rev3Δ rad30Δ quadruple mutant compared to wild type cells. We found a significant decrease in the overall rate of mutation accumulation in the quadruple mutant cells compared to wildtype, but an increase in frameshift mutations and a shift towards transversion base-substitution with a preference for G:C to T:A or C:G. Thus, the non-essential DNA polymerases have an impact on the nature of the mutational spectrum. The sequence and functional homology shared between human and S. cerevisiae non-essential DNA polymerases suggest these DNA polymerases may have a similar role in human cells.

Guidelines for DNA recombination and repair studies: Cellular assays of DNA repair pathways

Understanding the plasticity of genomes has been greatly aided by assays for recombination, repair and mutagenesis. These assays have been developed in microbial systems that provide the advantages of genetic and molecular reporters that can readily be manipulated. Cellular assays comprise genetic, molecular, and cytological reporters. The assays are powerful tools but each comes with its particular advantages and limitations. Here the most commonly used assays are reviewed, discussed, and presented as the guidelines for future studies.

Involvement of C-terminal truncation mutation of kinesin-5 in resistance to kinesin-5 inhibitor

Cultured cells easily develop resistance to kinesin-5 inhibitors (K5Is) often by overexpressing a related motor protein, kinesin-12/KIF15, or by acquiring mutations in the N-terminal motor domain of kinesin-5/KIF11 itself. We aimed to identify novel mechanisms responsible for resistance to S-trityl L-cysteine (STLC), one of the K5Is, using human osteosarcoma cell lines. Among six lines examined, U-2OS and HOS survived chronic STLC treatment and gave rise to resistant cells with IC₅₀s at least 10-fold higher than those of the respective parental lines. Depletion of KIF15 largely eliminated the acquired K5I resistance in both cases, consistent with the proposed notion that KIF15 is indispensable for it. In contrast to the KIF11-independent property of the cells derived from HOS, those derived from U-2OS still required KIF11 for their growth and, intriguingly, expressed a C-terminal truncated variant of KIF11 resulting from a frame shift mutation (S1017fs). All of the isolated clones harbored the same mutation, suggesting its clonal expansion in the cell population due to the growth advantage during chronic STLC treatment. Transgenic expression of KIF11^S1017fs in the parental U-2OS cells, as well as in HeLa cells, conferred a moderate but reproducible STLC resistance, probably owing to STLC-resistant localization of the mutant KIF11 on mitotic spindle. Our observations indicate that both KIF15 and the C-terminal-truncated KIF11 contributes to the STLC resistance of the U-2OS derived cells.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

BAF180: Its Roles in DNA Repair and Consequences in Cancer

In 2011, Varela et al. reported that the PBRM1 gene is mutated in approximately 40% of clear cell renal cell carcinoma cases. Since then, the number of studies relating PBRM1 mutations to cancers has substantially increased. BAF180 has now been linked to more than 30 types of cancers, including ccRCC, cholangiocarcinomas, esophageal squamous cell carcinoma, bladder cancer, and breast cancer. The mutations associated with BAF180 are most often truncations, which result in a loss of protein expression. This loss has been shown to adversely affect the expression of genes, likely because BAF180 is the chromatin recognition subunit of PBAF. In addition, BAF180 functions in numerous DNA repair mechanisms. Its roles in mediating DNA repair are likely the mechanism by which BAF180 acts a tumor suppressor protein. As research on this protein gains more interest, scientists will begin to piece together the complicated puzzle of the BAF180 protein and why its loss often results in cancer.

DNA mutation motifs in the genes associated with inherited diseases

Mutations in human genes can be responsible for inherited genetic disorders and cancer. Mutations can arise due to environmental factors or spontaneously. It has been shown that certain DNA sequences are more prone to mutate. These sites are termed hotspots and exhibit a higher mutation frequency than expected by chance. In contrast, DNA sequences with lower mutation frequencies than expected by chance are termed coldspots. Mutation hotspots are usually derived from a mutation spectrum, which reflects particular population where an effect of a common ancestor plays a role. To detect coldspots/hotspots unaffected by population bias, we analysed the presence of germline mutations obtained from HGMD database in the 5-nucleotide segments repeatedly occurring in genes associated with common inherited disorders, in particular, the PAH, LDLR, CFTR, F8, and F9 genes. Statistically significant sequences (mutational motifs) rarely associated with mutations (coldspots) and frequently associated with mutations (hotspots) exhibited characteristic sequence patterns, e.g. coldspots contained purine tract while hotspots showed alternating purine-pyrimidine bases, often with the presence of CpG dinucleotide. Using molecular dynamics simulations and free energy calculations, we analysed the global bending properties of two selected coldspots and two hotspots with a G/T mismatch. We observed that the coldspots were inherently more flexible than the hotspots. We assume that this property might be critical for effective mismatch repair as DNA with a mutation recognized by MutSα protein is noticeably bent.

Mechanisms of glycosylase induced genomic instability

Human alkyladenine DNA glycosylase (AAG) initiates base excision repair (BER) to guard against mutations by excising alkylated and deaminated purines. Counterintuitively, increased expression of AAG has been implicated in increased rates of spontaneous mutation in microsatellite repeats. This microsatellite mutator phenotype is consistent with a model in which AAG excises bulged (unpaired) bases, altering repeat length. To directly test the role of base excision in AAG-induced mutagenesis, we conducted mutation accumulation experiments in yeast overexpressing different variants of AAG and detected mutations via high-depth genome resequencing. We also developed a new software tool, hp_caller, to perform accurate genotyping at homopolymeric repeat loci. Overexpression of wild-type AAG elevated indel mutations in homopolymeric sequences distributed throughout the genome. However, catalytically inactive variants (E125Q/E125A) caused equal or greater increases in frameshift mutations. These results disprove the hypothesis that base excision is the key step in mutagenesis by overexpressed wild-type AAG. Instead, our results provide additional support for the previously published model wherein overexpressed AAG interferes with the mismatch repair (MMR) pathway. In addition to the above results, we observed a dramatic mutator phenotype for N169S AAG, which has increased rates of excision of undamaged purines. This mutant caused a 10-fold increase in point mutations at G:C base pairs and a 50-fold increase in frameshifts in A:T homopolymers. These results demonstrate that it is necessary to consider the relative activities and abundance of many DNA replication and repair proteins when considering mutator phenotypes, as they are relevant to the development of cancer and its resistance to treatment.

Identification and characterization of a -1 reading frameshift in the heavy chain constant region of an IgG1 recombinant monoclonal antibody produced in CHO cells

Frameshifts lead to complete alteration of the intended amino acid sequences, and therefore may affect the biological activities of protein therapeutics and pose potential immunogenicity risks. We report here the identification and characterization of a novel -1 frameshift variant in a recombinant IgG1 therapeutic monoclonal antibody (mAb) produced in Chinese hamster ovary cells during the cell line selection studies. The variant was initially observed as an atypical post-monomer fragment peak in size exclusion chromatography. Characterization of the fragment peak using intact and reduced liquid chromatography-mass spectrometry (LC-MS) analyses determined that the fragment consisted of a normal light chain disulfide-linked to an aberrant 26 kDa fragment that could not be assigned to any HC fragment even after considering common modifications. Further analysis using LC-MS/MS peptide mapping revealed that the aberrant fragment contained the expected HC amino acid sequence (1-232) followed by a 20-mer novel sequence corresponding to expression of heavy chain DNA sequence in the -1 reading frame. Examination of the DNA sequence around the frameshift initiation site revealed that a mononucleotide repeat GGGGGG located in the IgG1 HC constant region was most likely the structural root cause of the frameshift. Rapid identification of the frameshift allowed us to avoid use of a problematic cell line containing the frameshift as the production cell line. The frameshift reported here may be observed in other mAb products and the hypothesis-driven analytical approaches employed here may be valuable for rapid identification and characterization of frameshift variants in other recombinant proteins.

The peculiar landscape of repetitive sequences in the olive (Olea europaea L.) genome

Analyzing genome structure in different species allows to gain an insight into the evolution of plant genome size. Olive (Olea europaea L.) has a medium-sized haploid genome of 1.4 Gb, whose structure is largely uncharacterized, despite the growing importance of this tree as oil crop. Next-generation sequencing technologies and different computational procedures have been used to study the composition of the olive genome and its repetitive fraction. A total of 2.03 and 2.3 genome equivalents of Illumina and 454 reads from genomic DNA, respectively, were assembled following different procedures, which produced more than 200,000 differently redundant contigs, with mean length higher than 1,000 nt. Mapping Illumina reads onto the assembled sequences was used to estimate their redundancy. The genome data set was subdivided into highly and medium redundant and nonredundant contigs. By combining identification and mapping of repeated sequences, it was established that tandem repeats represent a very large portion of the olive genome (∼31% of the whole genome), consisting of six main families of different length, two of which were first discovered in these experiments. The other large redundant class in the olive genome is represented by transposable elements (especially long terminal repeat-retrotransposons). On the whole, the results of our analyses show the peculiar landscape of the olive genome, related to the massive amplification of tandem repeats, more than that reported for any other sequenced plant genome.

Assembly of a phased diploid Candida albicans genome facilitates allele-specific measurements and provides a simple model for repeat and indel structure

Background

Candida albicans is a ubiquitous opportunistic fungal pathogen that afflicts immunocompromised human hosts. With rare and transient exceptions the yeast is diploid, yet despite its clinical relevance the respective sequences of its two homologous chromosomes have not been completely resolved.

Results

We construct a phased diploid genome assembly by deep sequencing a standard laboratory wild-type strain and a panel of strains homozygous for particular chromosomes. The assembly has 700-fold coverage on average, allowing extensive revision and expansion of the number of known SNPs and indels. This phased genome significantly enhances the sensitivity and specificity of allele-specific expression measurements by enabling pooling and cross-validation of signal across multiple polymorphic sites. Additionally, the diploid assembly reveals pervasive and unexpected patterns in allelic differences between homologous chromosomes. Firstly, we see striking clustering of indels, concentrated primarily in the repeat sequences in promoters. Secondly, both indels and their repeat-sequence substrate are enriched near replication origins. Finally, we reveal an intimate link between repeat sequences and indels, which argues that repeat length is under selective pressure for most eukaryotes. This connection is described by a concise one-parameter model that explains repeat-sequence abundance in C. albicans as a function of the indel rate, and provides a general framework to interpret repeat abundance in species ranging from bacteria to humans.

Conclusions

The phased genome assembly and insights into repeat plasticity will be valuable for better understanding allele-specific phenomena and genome evolution.

Increased and imbalanced dNTP pools symmetrically promote both leading and lagging strand replication infidelity

The fidelity of DNA replication requires an appropriate balance of dNTPs, yet the nascent leading and lagging strands of the nuclear genome are primarily synthesized by replicases that differ in subunit composition, protein partnerships and biochemical properties, including fidelity. These facts pose the question of whether imbalanced dNTP pools differentially influence leading and lagging strand replication fidelity. Here we test this possibility by examining strand-specific replication infidelity driven by a mutation in yeast ribonucleotide reductase, rnr1-Y285A, that leads to elevated dTTP and dCTP concentrations. The results for the CAN1 mutational reporter gene present in opposite orientations in the genome reveal that the rates, and surprisingly even the sequence contexts, of replication errors are remarkably similar for leading and lagging strand synthesis. Moreover, while many mismatches driven by the dNTP pool imbalance are efficiently corrected by mismatch repair, others are repaired less efficiently, especially those in sequence contexts suggesting reduced proofreading due to increased mismatch extension driven by the high dTTP and dCTP concentrations. Thus the two DNA strands of the nuclear genome are at similar risk of mutations resulting from this dNTP pool imbalance, and this risk is not completely suppressed even when both major replication error correction mechanisms are genetically intact.

Comparative analysis of error-prone replication mononucleotide repeats across baculovirus genomes

Genome replication by the baculovirus DNA polymerase often generates errors in mononucleotide repeat (MNR) sequences due to replication slippage. This results in the inactivation of genes that affects different stages of the cell infection cycle. Here we mapped these MNRs in the 59 baculovirus genomes. We found that the MNR frequencies of baculovirus genomes are different and not correlated with the genome sizes. Although the average A/T content of baculoviruses is 58.67%, the A/T MNR frequency is significantly higher than that of the G/C MNRs. Furthermore, the A7/T7 MNRs are the most frequent of those we studied. Finally, MNR frequencies in different classes of baculovirus genes, such as immediate early genes, show differences between baculovirus genomes, suggesting that the distribution and frequency of different MNRs are unique to each baculovirus species or strain. Therefore, the results of this study can help select appropriate baculoviruses for the development of biological insecticides.

DNA repair mechanisms and the bypass of DNA damage in Saccharomyces cerevisiae

DNA repair mechanisms are critical for maintaining the integrity of genomic DNA, and their loss is associated with cancer predisposition syndromes. Studies in Saccharomyces cerevisiae have played a central role in elucidating the highly conserved mechanisms that promote eukaryotic genome stability. This review will focus on repair mechanisms that involve excision of a single strand from duplex DNA with the intact, complementary strand serving as a template to fill the resulting gap. These mechanisms are of two general types: those that remove damage from DNA and those that repair errors made during DNA synthesis. The major DNA-damage repair pathways are base excision repair and nucleotide excision repair, which, in the most simple terms, are distinguished by the extent of single-strand DNA removed together with the lesion. Mistakes made by DNA polymerases are corrected by the mismatch repair pathway, which also corrects mismatches generated when single strands of non-identical duplexes are exchanged during homologous recombination. In addition to the true repair pathways, the postreplication repair pathway allows lesions or structural aberrations that block replicative DNA polymerases to be tolerated. There are two bypass mechanisms: an error-free mechanism that involves a switch to an undamaged template for synthesis past the lesion and an error-prone mechanism that utilizes specialized translesion synthesis DNA polymerases to directly synthesize DNA across the lesion. A high level of functional redundancy exists among the pathways that deal with lesions, which minimizes the detrimental effects of endogenous and exogenous DNA damage.

Mutation rates, spectra, and genome-wide distribution of spontaneous mutations in mismatch repair deficient yeast

DNA mismatch repair is a highly conserved DNA repair pathway. In humans, germline mutations in hMSH2 or hMLH1, key components of mismatch repair, have been associated with Lynch syndrome, a leading cause of inherited cancer mortality. Current estimates of the mutation rate and the mutational spectra in mismatch repair defective cells are primarily limited to a small number of individual reporter loci. Here we use the yeast Saccharomyces cerevisiae to generate a genome-wide view of the rates, spectra, and distribution of mutation in the absence of mismatch repair. We performed mutation accumulation assays and next generation sequencing on 19 strains, including 16 msh2 missense variants implicated in Lynch cancer syndrome. The mutation rate for DNA mismatch repair null strains was approximately 1 mutation per genome per generation, 225-fold greater than the wild-type rate. The mutations were distributed randomly throughout the genome, independent of replication timing. The mutation spectra included insertions/deletions at homopolymeric runs (87.7%) and at larger microsatellites (5.9%), as well as transitions (4.5%) and transversions (1.9%). Additionally, repeat regions with proximal repeats are more likely to be mutated. A bias toward deletions at homopolymers and insertions at (AT)_n microsatellites suggests a different mechanism for mismatch generation at these sites. Interestingly, 5% of the single base pair substitutions might represent double-slippage events that occurred at the junction of immediately adjacent repeats, resulting in a shift in the repeat boundary. These data suggest a closer scrutiny of tumor suppressors with homopolymeric runs with proximal repeats as the potential drivers of oncogenesis in mismatch repair defective cells.

Cell-dependent production of polyhedra and virion occlusion of Autographa californica multiple nucleopolyhedrovirus fp25k mutants in vitro and in vivo

Members of the family Baculoviridae are insect-specific dsDNA viruses that have been used for biological control of insect pests in agriculture and forestry, as well as in research and pharmaceutical protein expression in insect cells and larvae. Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is the type species of the family Baculoviridae. During infection of AcMNPV in permissive cells, fp25k mutants are positively selected, leading to the formation of the few polyhedra (FP) phenotype with reduced yield of polyhedra and reduced virion occlusion efficiency, which leads to decreased oral infectivity for insects. Here we report that polyhedra of AcMNPV fp25k mutants produced from different insect cell lines and insects have differences in larval per os infectivity, and that these variations are due to different virion occlusion efficiencies in these cell lines and insects. Polyhedra of AcMNPV fp25k mutants produced from Sf cells (Sf21 and Sf9, derived from Spodoptera frugiperda) and S. frugiperda larvae had poorer virion occlusion efficiency than those from Hi5 cells (derived from Trichoplusia ni) and T. ni larvae, based on immunoblots, DNA isolation and larval oral infection analysis. AcMNPV fp25k mutants formed clusters of FP and many polyhedra (MP) in the fat body cells of both T. ni and S. frugiperda larvae. Transmission electron microscopy revealed that the nature of virion occlusion of AcMNPV fp25k mutants was dependent on the different cells of the T. ni fat body tissue. Taken together, these results indicate that the FP phenotype and virion occlusion efficiency of fp25k mutants are influenced by the host insect cells.

Beyond junk-variable tandem repeats as facilitators of rapid evolution of regulatory and coding sequences

Copy Number Variations (CNVs) and Single Nucleotide Polymorphisms (SNPs) have been the major focus of most large-scale comparative genomics studies to date. Here, we discuss a third, largely ignored, type of genetic variation, namely changes in tandem repeat number. Historically, tandem repeats have been designated as non functional “junk” DNA, mostly as a result of their highly unstable nature. With the exception of tandem repeats involved in human neurodegenerative diseases, repeat variation was often believed to be neutral with no phenotypic consequences. Recent studies, however, have shown that as many as 10% to 20% of coding and regulatory sequences in eukaryotes contain an unstable repeat tract. Contrary to initial suggestions, tandem repeat variation can have useful phenotypic consequences. Examples include rapid variation in microbial cell surface, tuning of internal molecular clocks in flies and the dynamic morphological plasticity in mammals. As such, tandem repeats can be useful functional elements that facilitate evolvability and rapid adaptation.

Phase variable genes of Campylobacter jejuni exhibit high mutation rates and specific mutational patterns but mutability is not the major determinant of population structure during host colonization

Phase variation of surface structures occurs in diverse bacterial species due to stochastic, high frequency, reversible mutations. Multiple genes of Campylobacter jejuni are subject to phase variable gene expression due to mutations in polyC/G tracts. A modal length of nine repeats was detected for polyC/G tracts within C. jejuni genomes. Switching rates for these tracts were measured using chromosomally-located reporter constructs and high rates were observed for cj1139 (G8) and cj0031 (G9). Alteration of the cj1139 tract from G8 to G11 increased mutability 10-fold and changed the mutational pattern from predominantly insertions to mainly deletions. Using a multiplex PCR, major changes were detected in ‘on/off’ status for some phase variable genes during passage of C. jejuni in chickens. Utilization of observed switching rates in a stochastic, theoretical model of phase variation demonstrated links between mutability and genetic diversity but could not replicate observed population diversity. We propose that modal repeat numbers have evolved in C. jejuni genomes due to molecular drivers associated with the mutational patterns of these polyC/G repeats, rather than by selection for particular switching rates, and that factors other than mutational drift are responsible for generating genetic diversity during host colonization by this bacterial pathogen.

Escherichia coli frameshift mutation rate depends on the chromosomal context but not on the GATC content near the mutation site

Different studies have suggested that mutation rate varies at different positions in the genome. In this work we analyzed if the chromosomal context and/or the presence of GATC sites can affect the frameshift mutation rate in the Escherichia coli genome. We show that in a mismatch repair deficient background, a condition where the mutation rate reflects the fidelity of the DNA polymerization process, the frameshift mutation rate could vary up to four times among different chromosomal contexts. Furthermore, the mismatch repair efficiency could vary up to eight times when compared at different chromosomal locations, indicating that detection and/or repair of frameshift events also depends on the chromosomal context. Also, GATC sequences have been proved to be essential for the correct functioning of the E. coli mismatch repair system. Using bacteriophage heteroduplexes molecules it has been shown that GATC influence the mismatch repair efficiency in a distance- and number-dependent manner, being almost nonfunctional when GATC sequences are located at 1 kb or more from the mutation site. Interestingly, we found that in E. coli genomic DNA the mismatch repair system can efficiently function even if the nearest GATC sequence is located more than 2 kb away from the mutation site. The results presented in this work show that even though frameshift mutations can be efficiently generated and/or repaired anywhere in the genome, these processes can be modulated by the chromosomal context that surrounds the mutation site.

Widespread occurrence of power-law distributions in inter-repeat distances shaped by genome dynamics

Repetitive DNA sequences derived from transposable elements (TE) are distributed in a non-random way, co-clustering with other classes of repeat elements, genes and other genomic components. In a previous work we reported power-law-like size distributions (linearity in log-log scale) in the spatial arrangement of Alu and LINE1 elements in the human genome. Here we investigate the large-scale features of the spatial arrangement of all principal classes of TEs in 14 genomes from phylogenetically distant organisms by studying the size distribution of inter-repeat distances. Power-law-like size distributions are found to be widespread, extending up to several orders of magnitude. In order to understand the emergence of this distributional pattern, we introduce an evolutionary scenario, which includes (i) Insertions of DNA segments (e.g., more recent repeats) into the considered sequence and (ii) Eliminations of members of the studied TE family. In the proposed model we also incorporate the potential for transposition events (characteristic of the DNA transposons' life-cycle) and segmental duplications. Simulations reproduce the main features of the observed size distributions. Furthermore, we investigate the effects of various genomic features on the presence and extent of power-law size distributions including TE class and age, mode of parental TE transmission, GC content, deletion and recombination rates in the studied genomic region, etc. Our observations corroborate the hypothesis that insertions of genomic material and eliminations of repeats are at the basis of power-laws in inter-repeat distances. The existence of these power-laws could facilitate the formation of the recently proposed "fractal globule" for the confined chromatin organization.

Mutation hot spots in yeast caused by long-range clustering of homopolymeric sequences

Evolutionary theory assumes that mutations occur randomly in the genome; however, studies performed in a variety of organisms indicate the existence of context-dependent mutation biases. Sources of mutagenesis variation across large genomic contexts (e.g., hundreds of bases) have not been identified. Here, we use high-coverage whole-genome sequencing of a conditional mismatch repair mutant line of diploid yeast to identify mutations that accumulated after 160 generations of growth. The vast majority of the mutations accumulated as insertion/deletions (in/dels) in homopolymeric [poly(dA:dT)] and repetitive DNA tracts. Surprisingly, the likelihood of an in/del mutation in a given poly(dA:dT) tract is increased by the presence of nearby poly(dA:dT) tracts in up to a 1,000 bp region centered on the given tract. Our work suggests that specific mutation hot spots can contribute disproportionately to the genetic variation that is introduced into populations and provides long-range genomic sequence context that contributes to mutagenesis.

Evolutionary pressures on simple sequence repeats in prokaryotic coding regions

Simple sequence repeats (SSRs) are indel mutational hotspots in genomes. In prokaryotes, SSR loci can cause phase variation, a microbial survival strategy that relies on stochastic, reversible on–off switching of gene activity. By analyzing multiple strains of 42 fully sequenced prokaryotic species, we measure the relative variability and density distribution of SSRs in coding regions. We demonstrate that repeat type strongly influences indel mutation rates, and that the most mutable types are most strongly avoided across genomes. We thoroughly characterize SSR density and variability as a function of N→C position along protein sequences. Using codon-shuffling algorithms that preserve amino acid sequence, we assess evolutionary pressures on SSRs. We find that coding sequences suppress repeats in the middle of proteins, and enrich repeats near termini, yielding U-shaped SSR density curves. We show that for many species this characteristic shape can be attributed to purely biophysical constraints of protein structure. In multiple cases, however, particularly in certain pathogenic bacteria, we observe over enrichment of SSRs near protein N-termini significantly beyond expectation based on structural constraints. This increases the probability that frameshifts result in non-functional proteins, revealing that these species may evolutionarily tune SSR positions in coding regions to facilitate phase variation.

Transient suppression of MLH1 allows effective single-nucleotide substitution by single-stranded DNA oligonucleotides

Short synthetic single-stranded oligodeoxyribonucleotides (ssODNs) can be used to introduce subtle modifications into the genome of mouse embryonic stem cells (ESCs). We have previously shown that effective application of ssODN-mediated gene targeting in ESC requires (transient) suppression of DNA mismatch repair (MMR). However, whereas transient down-regulation of the mismatch recognition protein MSH2 allowed substitution of 3 or 4 nucleotides, 1 or 2 nucleotide substitutions were still suppressed. We now demonstrate that single- or dinucleotide substitution can effectively be achieved by transient down-regulation of the downstream MMR protein MLH1. By exploiting highly specific real-time PCR, we demonstrate the feasibility of substituting a single basepair in a non-selectable gene. However, disabling the MMR machinery may lead to inadvertent mutations. To obtain insight into the mutation rate associated with transient MMR suppression, we have compared the impact of transient and constitutive MMR deficiency on the repair of frameshift intermediates at mono- and dinucleotide repeats. Repair at these repeats relied on the substrate specificity and functional redundancy of the MSH2/MSH6 and MSH2/MSH3 MMR complexes. MLH1 knockdown increased the level of spontaneous mutagenesis, but modified ESCs remained germ line competent. Thus, transient MLH1 suppression provides a valuable extension of the MSH2 knockdown strategy, allowing rapid generation of mice carrying single basepair alterations in their genome.

The in vitro fidelity of yeast DNA polymerase δ and polymerase ε holoenzymes during dinucleotide microsatellite DNA synthesis

Elucidating the sources of genetic variation within microsatellite alleles has important implications for understanding the etiology of human diseases. Mismatch repair is a well described pathway for the suppression of microsatellite instability. However, the cellular polymerases responsible for generating microsatellite errors have not been fully described. We address this gap in knowledge by measuring the fidelity of recombinant yeast polymerase δ (Pol δ) and ɛ (Pol ɛ) holoenzymes during synthesis of a [GT/CA] microsatellite. The in vitro HSV-tk forward assay was used to measure DNA polymerase errors generated during gap-filling of complementary GT(10) and CA(10)-containing substrates and ∼90 nucleotides of HSV-tk coding sequence surrounding the microsatellites. The observed mutant frequencies within the microsatellites were 4 to 30-fold higher than the observed mutant frequencies within the coding sequence. More specifically, the rate of Pol δ and Pol ɛ misalignment-based insertion/deletion errors within the microsatellites was ∼1000-fold higher than the rate of insertion/deletion errors within the HSV-tk gene. Although the most common microsatellite error was the deletion of a single repeat unit, ∼ 20% of errors were deletions of two or more units for both polymerases. The differences in fidelity for wild type enzymes and their exonuclease-deficient derivatives were ∼2-fold for unit-based microsatellite insertion/deletion errors. Interestingly, the exonucleases preferentially removed potentially stabilizing interruption errors within the microsatellites. Since Pol δ and Pol ɛ perform not only the bulk of DNA replication in eukaryotic cells but also are implicated in performing DNA synthesis associated with repair and recombination, these results indicate that microsatellite errors may be introduced into the genome during multiple DNA metabolic pathways.

Simple sequence repeat variation in the Daphnia pulex genome

Background

Simple sequence repeats (SSRs) are highly variable features of all genomes. Their rapid evolution makes them useful for tracing the evolutionary history of populations and investigating patterns of selection and mutation across gnomes. The recently sequenced Daphnia pulex genome provides us with a valuable data set to study the mode and tempo of SSR evolution, without the inherent biases that accompany marker selection.

Results

Here we catalogue SSR loci in the Daphnia pulex genome with repeated motif sizes of 1-100 nucleotides with a minimum of 3 perfect repeats. We then used whole genome shotgun reads to determine the average heterozygosity of each SSR type and the relationship that it has to repeat number, motif size, motif sequence, and distribution of SSR loci. We find that SSR heterozygosity is motif specific, and positively correlated with repeat number as well as motif size. For non-repeat unit polymorphisms, we identify a motif-dependent end-nucleotide polymorphism bias that may contribute to the patterns of abundance for specific homopolymers, dimers, and trimers. Our observations confirm the high frequency of multiple unit variation (multistep) at large microsatellite loci, and further show that the occurrence of multiple unit variation is dependent on both repeat number and motif size. Using the Daphnia pulex genetic map, we show a positive correlation between dimer and trimer frequency and recombination.

Conclusions

This genome-wide analysis of SSR variation in Daphnia pulex indicates that several aspects of SSR variation are motif dependent and suggests that a combination of unit length variation and end repeat biased base substitution contribute to the unique spectrum of SSR repeat loci.

Scaling properties and fractality in the distribution of coding segments in eukaryotic genomes revealed through a block entropy approach

Statistical methods, including block entropy based approaches, have already been used in the study of long-range features of genomic sequences seen as symbol series, either considering the full alphabet of the four nucleotides or the binary purine or pyrimidine character set. Here we explore the alternation of short protein-coding segments with long noncoding spacers in entire chromosomes, focusing on the scaling properties of block entropy. In previous studies, it has been shown that the sizes of noncoding spacers follow power-law-like distributions in most chromosomes of eukaryotic organisms from distant taxa. We have developed a simple evolutionary model based on well-known molecular events (segmental duplications followed by elimination of most of the duplicated genes) which reproduces the observed linearity in log-log plots. The scaling properties of block entropy H(n) have been studied in several works. Their findings suggest that linearity in semilogarithmic scale characterizes symbol sequences which exhibit fractal properties and long-range order, while this linearity has been shown in the case of the logistic map at the Feigenbaum accumulation point. The present work starts with the observation that the block entropy of the Cantor-like binary symbol series scales in a similar way. Then, we perform the same analysis for the full set of human chromosomes and for several chromosomes of other eukaryotes. A similar but less extended linearity in semilogarithmic scale, indicating fractality, is observed, while randomly formed surrogate sequences clearly lack this type of scaling. Genomic sequences always present entropy values much lower than their random surrogates. Symbol sequences produced by the aforementioned evolutionary model follow the scaling found in genomic sequences, thus corroborating the conjecture that "segmental duplication-gene elimination" dynamics may have contributed to the observed long rangeness in the coding or noncoding alternation in genomes.

The baker's yeast diploid genome is remarkably stable in vegetative growth and meiosis

Accurate estimates of mutation rates provide critical information to analyze genome evolution and organism fitness. We used whole-genome DNA sequencing, pulse-field gel electrophoresis, and comparative genome hybridization to determine mutation rates in diploid vegetative and meiotic mutation accumulation lines of Saccharomyces cerevisiae. The vegetative lines underwent only mitotic divisions while the meiotic lines underwent a meiotic cycle every ∼20 vegetative divisions. Similar base substitution rates were estimated for both lines. Given our experimental design, these measures indicated that the meiotic mutation rate is within the range of being equal to zero to being 55-fold higher than the vegetative rate. Mutations detected in vegetative lines were all heterozygous while those in meiotic lines were homozygous. A quantitative analysis of intra-tetrad mating events in the meiotic lines showed that inter-spore mating is primarily responsible for rapidly fixing mutations to homozygosity as well as for removing mutations. We did not observe 1-2 nt insertion/deletion (in-del) mutations in any of the sequenced lines and only one structural variant in a non-telomeric location was found. However, a large number of structural variations in subtelomeric sequences were seen in both vegetative and meiotic lines that did not affect viability. Our results indicate that the diploid yeast nuclear genome is remarkably stable during the vegetative and meiotic cell cycles and support the hypothesis that peripheral regions of chromosomes are more dynamic than gene-rich central sections where structural rearrangements could be deleterious. This work also provides an improved estimate for the mutational load carried by diploid organisms.

Avoidance of long mononucleotide repeats in codon pair usage

Protein is an essential component for life, and its synthesis is mediated by codons in any organisms on earth. While some codons encode the same amino acid, their usage is often highly biased. There are many factors that can cause the bias, but a potential effect of mononucleotide repeats, which are known to be highly mutable, on codon usage and codon pair preference is largely unknown. In this study we performed a genomic survey on the relationship between mononucleotide repeats and codon pair bias in 53 bacteria, 68 archaea, and 13 eukaryotes. By distinguishing the codon pair bias from the codon usage bias, four general patterns were revealed: strong avoidance of five or six mononucleotide repeats in codon pairs; lower observed/expected (o/e) ratio for codon pairs with C or G repeats (C/G pairs) than that with A or T repeats (A/T pairs); a negative correlation between genomic GC contents and the o/e ratios, particularly for C/G pairs; and avoidance of C/G pairs in highly conserved genes. These results support natural selection against long mononucleotide repeats, which could induce frameshift mutations in coding sequences. The fact that these patterns are found in all kingdoms of life suggests that this is a general phenomenon in living organisms. Thus, long mononucleotide repeats may play an important role in base composition and genetic stability of a gene and gene functions.

What is a microsatellite: a computational and experimental definition based upon repeat mutational behavior at A/T and GT/AC repeats

Microsatellites are abundant in eukaryotic genomes and have high rates of strand slippage-induced repeat number alterations. They are popular genetic markers, and their mutations are associated with numerous neurological diseases. However, the minimal number of repeats required to constitute a microsatellite has been debated, and a definition of a microsatellite that considers its mutational behavior has been lacking. To define a microsatellite, we investigated slippage dynamics for a range of repeat sizes, utilizing two approaches. Computationally, we assessed length polymorphism at repeat loci in ten ENCODE regions resequenced in four human populations, assuming that the occurrence of polymorphism reflects strand slippage rates. Experimentally, we determined the in vitro DNA polymerase-mediated strand slippage error rates as a function of repeat number. In both approaches, we compared strand slippage rates at tandem repeats with the background slippage rates. We observed two distinct modes of mutational behavior. At small repeat numbers, slippage rates were low and indistinguishable from background measurements. A marked transition in mutability was observed as the repeat array lengthened, such that slippage rates at large repeat numbers were significantly higher than the background rates. For both mononucleotide and dinucleotide microsatellites studied, the transition length corresponded to a similar number of nucleotides (approximately 10). Thus, microsatellite threshold is determined not by the presence/absence of strand slippage at repeats but by an abrupt alteration in slippage rates relative to background. These findings have implications for understanding microsatellite mutagenesis, standardization of genome-wide microsatellite analyses, and predicting polymorphism levels of individual microsatellite loci.

Detection of heterozygous mutations in the genome of mismatch repair defective diploid yeast using a Bayesian approach

DNA replication errors that escape polymerase proofreading and mismatch repair (MMR) can lead to base substitution and frameshift mutations. Such mutations can disrupt gene function, reduce fitness, and promote diseases such as cancer and are also the raw material of molecular evolution. To analyze with limited bias genomic features associated with DNA polymerase errors, we performed a genome-wide analysis of mutations that accumulate in MMR-deficient diploid lines of Saccharomyces cerevisiae. These lines were derived from a common ancestor and were grown for 160 generations, with bottlenecks reducing the population to one cell every 20 generations. We sequenced to between 8- and 20-fold coverage one wild-type and three mutator lines using Illumina Solexa 36-bp reads. Using an experimentally aware Bayesian genotype caller developed to pool experimental data across sequencing runs for all strains, we detected 28 heterozygous single-nucleotide polymorphisms (SNPs) and 48 single-nt insertion/deletions (indels) from the data set. This method was evaluated on simulated data sets and found to have a very low false-positive rate (∼6 × 10(-5)) and a false-negative rate of 0.08 within the unique mapping regions of the genome that contained at least sevenfold coverage. The heterozygous mutations identified by the Bayesian genotype caller were confirmed by Sanger sequencing. All of the mutations were unique to a given line, except for a single-nt deletion mutation which occurred independently in two lines. All 48 indels, composed of 46 deletions and two insertions, occurred in homopolymer (HP) tracts [i.e., 47 poly(A) or (T) tracts, 1 poly(G) or (C) tract] between 5 and 13 bp long. Our findings are of interest because HP tracts are present at high levels in the yeast genome (>77,400 for 5- to 20-nt HP tracts), and frameshift mutations in these regions are likely to disrupt gene function. In addition, they demonstrate that the mutation pattern seen previously in mismatch repair defective strains using a limited number of reporters holds true for the entire genome.

The nucleotide composition of microsatellites impacts both replication fidelity and mismatch repair in human colorectal cells

Microsatellite instability is a key mechanism of colon carcinogenesis. We have previously studied mutations within a (CA)13 microsatellite using an enhanced green fluorescent protein (EGFP)-based reporter assay that allows the distinction of replication errors and mismatch repair (MMR) activity. Here we utilize this assay to compare mutations of mono- and dinucleotide repeats in human colorectal cells. HCT116 and HCT116+chr3 cells were stably transfected with EGFP-based plasmids harboring A10, G10, G16, (CA)13 and (CA)26 repeats. EGFP-positive mutant fractions were quantitated by flow cytometry, mutation rates were calculated and the mutant spectrum was analyzed by cycle sequencing. EGFP fluorescence pattern changed with the microsatellite's nucleotide sequence and cell type and clonal variations were observed in mononucleotide repeats. Replication errors (as calculated in HCT116) at A10 repeats were 5–10-fold higher than in G10, G16 were 30-fold higher than G10 and (CA)26 were 10-fold higher than (CA)13. The mutation rates in hMLH1-proficient HCT116+chr3 were 30–230-fold lower than in HCT116. MMR was more efficient in G16 than in A10 clones leading to a higher stability of poly-G tracts. Mutation spectra revealed predominantly 1-unit deletions in A10, (CA)13 and G10 and 2-unit deletions or 1-unit insertion in (CA)26. These findings indicate that both replication fidelity and MMR are affected by the microsatellite's nucleotide composition.

Every microsatellite is different: Intrinsic DNA features dictate mutagenesis of common microsatellites present in the human genome

Microsatellite sequences are ubiquitous in the human genome and are important regulators of genome function. Here, we examine the mutational mechanisms governing the stability of highly abundant mono-, di-, and tetranucleotide microsatellites. Microsatellite mutation rate estimates from pedigree analyses and experimental models range from a low of approximately 10(-6) to a high of approximately 10(-2) mutations per locus per generation. The vast majority of observed mutational variation can be attributed to features intrinsic to the allele itself, including motif size, length, and sequence composition. A greater than linear relationship between motif length and mutagenesis has been observed in several model systems. Motif sequence differences contribute up to 10-fold to the variation observed in human cell mutation rates. The major mechanism of microsatellite mutagenesis is strand slippage during DNA synthesis. DNA polymerases produce errors within microsatellites at a frequency that is 10- to 100-fold higher than the frequency of frameshifts in coding sequences. Motif sequence significantly affects both polymerase error rate and specificity, resulting in strand biases within complementary microsatellites. Importantly, polymerase errors within microsatellites include base substitutions, deletions, and complex mutations, all of which produced interrupted alleles from pure microsatellites. Postreplication mismatch repair efficiency is affected by microsatellite motif size and sequence, also contributing to the observed variation in microsatellite mutagenesis. Inhibition of DNA synthesis within common microsatellites is highly sequence-dependent, and is positively correlated with the production of errors. DNA secondary structure within common microsatellites can account for some DNA polymerase pause sites, and may be an important factor influencing mutational specificity.

Human postmeiotic segregation 2 exhibits biased repair at tetranucleotide microsatellite sequences

The mismatch repair (MMR) system plays a major role in removing DNA polymerization errors, and loss of this pathway results in hereditary cancers characterized by microsatellite instability. We investigated microsatellite stability during DNA replication within human postmeiotic segregation 2 (hPMS2)-deficient and proficient human lymphoblastoid cell lines. Using a shuttle vector assay, we measured mutation rates at reporter cassettes containing defined mononucleotide, dinucleotide, and tetranucleotide microsatellite sequences. A mutator phenotype was observed in the hPMS2-deficient cell line. The mutation rate of vectors containing [G/C](10) or [GT/CA](10) alleles was elevated 20-fold to 40-fold in hPMS2-deficient cells, relative to an hPMS2-expressing cell line. We observed a 6-fold and 12-fold relative increase in mutation rate of [TTTC/AAAG](9) and [TTCC/AAGG](9) sequences, respectively, in hPMS2-deficient cells. Mutational specificity analyses suggested that repair by hPMS2 is biased. In the absence of hPMS2, a greater number of microsatellite expansion versus deletion mutations was observed, and expansion rates of the tetranucleotide alleles were similar. In the presence of hPMS2, we observed a 29-fold decrease in the [TTCC/AAGG](9) expansion rate but only a 6-fold decrease for the [TTTC/AAAG](9) allele. Our data indicate that hPMS2 is more protective of tetranucleotide expansions than deletions and that hPMS2 displays a sequence bias, wherein [TTCC/AAGG] sequences are stabilized to a greater extent than [TTTC/AAAG]. Our results allow for greater accuracy during identification of MMR defects by providing a mutational signature characteristic of hPMS2 defect. This study also provides clues to possible mechanisms of repair by hPMS2 in the context of the MMR system.

The rate and character of spontaneous mutation in Thermus thermophilus

Selection of spontaneous, loss-of-function mutations at two chromosomal loci (pyrF and pyrE) enabled the first molecular-level analysis of replication fidelity in the extremely thermophilic bacterium Thermus thermophilus. Two different methods yielded similar mutation rates, and mutational spectra determined by sequencing of independent mutants revealed a variety of replication errors distributed throughout the target genes. The genomic mutation rate estimated from these targets, 0.00097 +/- 0.00052 per replication, was lower than corresponding estimates from mesophilic microorganisms, primarily because of a low rate of base substitution. However, both the rate and spectrum of spontaneous mutations in T. thermophilus resembled those of the thermoacidophilic archaeon Sulfolobus acidocaldarius, despite important molecular differences between these two thermophiles and their genomes.

The rate and spectrum of microsatellite mutation in Caenorhabditis elegans and Daphnia pulex

The effective use of microsatellite loci as tools for microevolutionary analysis requires knowledge of the factors influencing the rate and pattern of mutation, much of which is derived from indirect inference from population samples. Interspecific variation in microsatellite stability also provides a glimpse into aspects of phylogenetic constancy of mutational processes. Using long-term series of mutation-accumulation lines, we have obtained direct estimates of the spectrum of microsatellite mutations in two model systems: the nematode Caenorhabditis elegans and the microcrustacean Daphnia pulex. Although the scaling of the mutation rate with the number of tandem repeats is highly consistent across distantly related species, including yeast and human, the per-cell-division mutation rate appears to be elevated in multicellular species. Contrary to the expectations under the stepwise mutation model, most microsatellite mutations in C. elegans and D. pulex involve changes of multiple repeat units, with expansions being much more common than contractions.

A genome-wide view of the spectrum of spontaneous mutations in yeast

The mutation process ultimately defines the genetic features of all populations and, hence, has a bearing on a wide range of issues involving evolutionary genetics, inheritance, and genetic disorders, including the predisposition to cancer. Nevertheless, formidable technical barriers have constrained our understanding of the rate at which mutations arise and the molecular spectrum of their effects. Here, we report on the use of complete-genome sequencing in the characterization of spontaneously arising mutations in the yeast Saccharomyces cerevisiae. Our results confirm some findings previously obtained by indirect methods but also yield numerous unexpected findings, in particular a very high rate of point mutation and skewed distribution of base-substitution types in the mitochondrion, a very high rate of segmental duplication and deletion in the nuclear genome, and substantial deviations in the mutational profile among various model organisms.