Misalignment-mediated DNA polymerase beta mutations: comparison of microsatellite and frame-shift error rates using a forward mutation assay

Nontraditional Roles of DNA Polymerase Eta Support Genome Duplication and Stability

DNA polymerase eta (Pol η) is a Y-family polymerase and the product of the POLH gene. Autosomal recessive inheritance of POLH mutations is the cause of the xeroderma pigmentosum variant, a cancer predisposition syndrome. This review summarizes mounting evidence for expanded Pol η cellular functions in addition to DNA lesion bypass that are critical for maintaining genome stability. In vitro, Pol η displays efficient DNA synthesis through difficult-to-replicate sequences, catalyzes D-loop extensions, and utilizes RNA-DNA hybrid templates. Human Pol η is constitutively present at the replication fork. In response to replication stress, Pol η is upregulated at the transcriptional and protein levels, and post-translational modifications regulate its localization to chromatin. Numerous studies show that Pol η is required for efficient common fragile site replication and stability. Additionally, Pol η can be recruited to stalled replication forks through protein-protein interactions, suggesting a broader role in replication fork recovery. During somatic hypermutations, Pol η is recruited by mismatch repair proteins and is essential for V_H gene A:T basepair mutagenesis. Within the global context of repeat-dense genomes, the recruitment of Pol η to perform specialized functions during replication could promote genome stability by interrupting pure repeat arrays with base substitutions. Alternatively, not engaging Pol η in genome duplication is costly, as the absence of Pol η leads to incomplete replication and increased chromosomal instability.

Variation in G-quadruplex sequence and topology differentially impacts human DNA polymerase fidelity

G-quadruplexes (G4s), a type of non-B DNA, play important roles in a wide range of molecular processes, including replication, transcription, and translation. Genome integrity relies on efficient and accurate DNA synthesis, and is compromised by various stressors, to which non-B DNA structures such as G4s can be particularly vulnerable. However, the impact of G4 structures on DNA polymerase fidelity is largely unknown. Using an in vitro forward mutation assay, we investigated the fidelity of human DNA polymerases delta (δ4, four-subunit), eta (η), and kappa (κ) during synthesis of G4 motifs representing those in the human genome. The motifs differ in sequence, topology, and stability, features that may affect DNA polymerase errors. Polymerase error rate hierarchy (δ4 < κ < η) is largely maintained during G4 synthesis. Importantly, we observed unique polymerase error signatures during synthesis of VEGF G4 motifs, stable G4s which form parallel topologies. These statistically significant errors occurred within, immediately flanking, and encompassing the G4 motif. For pol δ4, the errors were deletions, insertions and complex errors within the G4 or encompassing the G4 motif and surrounding sequence. For pol η, the errors occurred in 3' sequences flanking the G4 motif. For pol κ, the errors were frameshift mutations within G-tracts of the G4. Because these error signatures were not observed during synthesis of an antiparallel G4 and, to a lesser extent, a hybrid G4, we suggest that G4 topology and/or stability could influence polymerase fidelity. Using in silico analyses, we show that most polymerase errors are predicted to have minimal effects on predicted G4 stability. Our results provide a unique view of G4s not previously elucidated, showing that G4 motif heterogeneity differentially influences polymerase fidelity within the motif and flanking sequences. Thus, our study advances the understanding of how DNA polymerase errors contribute to G4 mutagenesis.

PRIMPOL ready, set, reprime!

DNA replication forks are constantly challenged by DNA lesions induced by endogenous and exogenous sources. DNA damage tolerance mechanisms ensure that DNA replication continues with minimal effects on replication fork elongation either by using specialized DNA polymerases, which have the ability to replicate through the damaged template, or by skipping the damaged DNA, leaving it to be repaired after replication. These mechanisms are evolutionarily conserved in bacteria, yeast, and higher eukaryotes, and are paramount to ensure timely and faithful duplication of the genome. The Primase and DNA-directed Polymerase (PRIMPOL) is a recently discovered enzyme that possesses both primase and polymerase activities. PRIMPOL is emerging as a key player in DNA damage tolerance, particularly in vertebrate and human cells. Here, we review our current understanding of the function of PRIMPOL in DNA damage tolerance by focusing on the structural aspects that define its dual enzymatic activity, as well as on the mechanisms that control its chromatin recruitment and expression levels. We also focus on the latest findings on the mitochondrial and nuclear functions of PRIMPOL and on the impact of loss of these functions on genome stability and cell survival. Defining the function of PRIMPOL in DNA damage tolerance is becoming increasingly important in the context of human disease. In particular, we discuss recent evidence pointing at the PRIMPOL pathway as a novel molecular target to improve cancer cell response to DNA-damaging chemotherapy and as a predictive parameter to stratify patients in personalized cancer therapy.

De novo transcriptome analysis of Rhododendron molle G. Don flowers by Illumina sequencing

Rhododendron molle G. Don occupies an important phylogenetic node in the genus rhododendron with unique yellow flower and medicinal functions. However, only limited genetic resources and their genome information are available for the generation of rhododendron flowers. The next generation sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provide a turning point for rhododendron research. Our goal is to use the genetic information to facilitate the relevant research on flowering and flower color formation in R. molle. In total, 66,026 unigenes were identified, among which 31,298 were annotated in the NCBI non-redundant protein database and 22,410 were annotated in the Swiss-Prot database. Of these annotated unigenes, 9490 and 18,680 unigenes were assigned to clusters of orthologous groups and gene ontology categories, respectively. A total of 7177 genes were mapped to 118 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. In addition, 8266 simple sequence repeats (SSRs) were detected, and these SSRs will undoubtedly benefit rhododendron breeding work. Metabolic pathway analysis revealed that 32 unigenes were predicted to be involved in carotenoid biosynthesis. Our transcriptome revealed 32 engines that encode key enzymes in the carotenoid biosynthesis pathway, including PSY, PDS, LCYB, LCYE, etc. The content of β-carotene was much higher than the other carotenoids throughout the flower development. It was consistent with the key genes expression level in the carotenoid biosynthesis pathway by the Illumina expression profile analysis and the qRT-PCR analysis. Our study identified genes associated with carotenoid biosynthesis in R. molle and provides a valuable resource for understanding the flowering and flower color formation mechanisms in R. molle.

Somatic microsatellite variability as a predictive marker for colorectal cancer and liver cancer progression

Microsatellites (MSTs) are short tandem repeated genetic motifs that comprise ~3% of the genome. MST instability (MSI), defined as acquired/lost primary alleles at a small subset of microsatellite loci (e.g. Bethesda markers), is a clinically relevant marker for colorectal cancer. However, these markers are not applicable to other types of cancers, specifically, for liver cancer which has a high mortality rate. Here we show that somatic MST variability (SMV), defined as the presence of additional, non-primary (aka minor) alleles at MST loci, is a complementary measure of MSI, and a genetic marker for colorectal and liver cancer. Re-analysis of Illumina sequenced exomes from The Cancer Genome Atlas indicates that SMV may distinguish a subpopulation of African American patients with colorectal cancer, which represents ~33% of the population in this study. Further, for liver cancer, a higher rate of SMV may be indicative of an earlier age of onset. The work presented here suggests that classical MSI should be expanded to include SMV, going beyond alterations of the primary alleles at a small number of microsatellite loci. This measure of SMV may represent a potential new diagnostic for a variety of cancers and may provide new information for colorectal cancer patients.

The complete mitochondrial genome sequence of the green microalga Lobosphaera (Parietochloris) incisa reveals a new type of palindromic repetitive repeat

Background

Lobosphaera incisa, formerly known as Myrmecia incisa and then Parietochloris incisa, is an oleaginous unicellular green alga belonging to the class Trebouxiophyceae (Chlorophyta). It is the richest known plant source of arachidonic acid, an ω-6 poly-unsaturated fatty acid valued by the pharmaceutical and baby-food industries. It is therefore an organism of high biotechnological interest, and we recently reported the sequence of its chloroplast genome.

Results

We now report the complete sequence of the mitochondrial genome of L. incisa from high-throughput Illumina short-read sequencing. The circular chromosome of 69,997 bp is predicted to encode a total of 64 genes, some harboring specific self-splicing group I and group II introns. Overall, the gene content is highly similar to that of the mitochondrial genomes of other Trebouxiophyceae, with 34 protein-coding, 3 rRNA, and 27 tRNA genes. Genes are distributed in two clusters located on different DNA strands, a bipartite arrangement that suggests expression from two divergent promoters yielding polycistronic primary transcripts. The L. incisa mitochondrial genome contains families of intergenic dispersed DNA repeat sequences that are not shared with other known mitochondrial genomes of Trebouxiophyceae. The most peculiar feature of the genome is a repetitive palindromic repeat, the LIMP (L. Incisa Mitochondrial Palindrome), found 19 times in the genome. It is formed by repetitions of an AACCA pentanucleotide, followed by an invariant 7-nt loop and a complementary repeat of the TGGTT motif. Analysis of the genome sequencing reads indicates that the LIMP can be a substrate for large-scale genomic rearrangements. We speculate that LIMPs can act as origins of replication. Deep sequencing of the L. incisa transcriptome also suggests that the LIMPs with long stems are sites of transcript processing. The genome also contains five copies of a related palindromic repeat, the HyLIMP, with a 10-nt motif related to that of the LIMP.

Conclusions

The mitochondrial genome of L. incisa encodes a unique type of repetitive palindromic repeat sequence, the LIMP, which can mediate genome rearrangements and play a role in mitochondrial gene expression. Experimental studies are needed to confirm and further characterize the functional role(s) of the LIMP.

Accurate typing of short tandem repeats from genome-wide sequencing data and its applications

Short tandem repeats (STRs) are implicated in dozens of human genetic diseases and contribute significantly to genome variation and instability. Yet profiling STRs from short-read sequencing data is challenging because of their high sequencing error rates. Here, we developed STR-FM, short tandem repeat profiling using flank-based mapping, a computational pipeline that can detect the full spectrum of STR alleles from short-read data, can adapt to emerging read-mapping algorithms, and can be applied to heterogeneous genetic samples (e.g., tumors, viruses, and genomes of organelles). We used STR-FM to study STR error rates and patterns in publicly available human and in-house generated ultradeep plasmid sequencing data sets. We discovered that STRs sequenced with a PCR-free protocol have up to ninefold fewer errors than those sequenced with a PCR-containing protocol. We constructed an error correction model for genotyping STRs that can distinguish heterozygous alleles containing STRs with consecutive repeat numbers. Applying our model and pipeline to Illumina sequencing data with 100-bp reads, we could confidently genotype several disease-related long trinucleotide STRs. Utilizing this pipeline, for the first time we determined the genome-wide STR germline mutation rate from a deeply sequenced human pedigree. Additionally, we built a tool that recommends minimal sequencing depth for accurate STR genotyping, depending on repeat length and sequencing read length. The required read depth increases with STR length and is lower for a PCR-free protocol. This suite of tools addresses the pressing challenges surrounding STR genotyping, and thus is of wide interest to researchers investigating disease-related STRs and STR evolution.

Genetic evidence that both dNTP-stabilized and strand slippage mechanisms may dictate DNA polymerase errors within mononucleotide microsatellites

Mononucleotide microsatellites are tandem repeats of a single base pair, abundant within coding exons and frequent sites of mutation in the human genome. Because the repeated unit is one base pair, multiple mechanisms of insertion/deletion (indel) mutagenesis are possible, including strand-slippage, dNTP-stabilized, and misincorportion-misalignment. Here, we examine the effects of polymerase identity (mammalian Pols α, β, κ, and η), template sequence, dNTP pool size, and reaction temperature on indel errors during in vitro synthesis of mononucleotide microsatellites. We utilized the ratio of insertion to deletion errors as a genetic indicator of mechanism. Strikingly, we observed a statistically significant bias toward deletion errors within mononucleotide repeats for the majority of the 28 DNA template and polymerase combinations examined, with notable exceptions based on sequence and polymerase identity. Using mutator forms of Pol β did not substantially alter the error specificity, suggesting that mispairing-misalignment mechanism is not a primary mechanism. Based on our results for mammalian DNA polymerases representing three structurally distinct families, we suggest that dNTP-stabilized mutagenesis may be an alternative mechanism for mononucleotide microsatellite indel mutation. The change from a predominantly dNTP-stabilized mechanism to a strand-slippage mechanism with increasing microsatellite length may account for the differential rates of tandem repeat mutation that are observed genome-wide.

Human PrimPol is a highly error-prone polymerase regulated by single-stranded DNA binding proteins

PrimPol is a recently identified polymerase involved in eukaryotic DNA damage tolerance, employed in both re-priming and translesion synthesis mechanisms to bypass nuclear and mitochondrial DNA lesions. In this report, we investigate how the enzymatic activities of human PrimPol are regulated. We show that, unlike other TLS polymerases, PrimPol is not stimulated by PCNA and does not interact with it in vivo. We identify that PrimPol interacts with both of the major single-strand binding proteins, RPA and mtSSB in vivo. Using NMR spectroscopy, we characterize the domains responsible for the PrimPol-RPA interaction, revealing that PrimPol binds directly to the N-terminal domain of RPA70. In contrast to the established role of SSBs in stimulating replicative polymerases, we find that SSBs significantly limit the primase and polymerase activities of PrimPol. To identify the requirement for this regulation, we employed two forward mutation assays to characterize PrimPol's replication fidelity. We find that PrimPol is a mutagenic polymerase, with a unique error specificity that is highly biased towards insertion-deletion errors. Given the error-prone disposition of PrimPol, we propose a mechanism whereby SSBs greatly restrict the contribution of this enzyme to DNA replication at stalled forks, thus reducing the mutagenic potential of PrimPol during genome replication.

Exome-wide somatic microsatellite variation is altered in cells with DNA repair deficiencies

Microsatellites (MST), tandem repeats of 1–6 nucleotide motifs, are mutational hot-spots with a bias for insertions and deletions (INDELs) rather than single nucleotide polymorphisms (SNPs). The majority of MST instability studies are limited to a small number of loci, the Bethesda markers, which are only informative for a subset of colorectal cancers. In this paper we evaluate non-haplotype alleles present within next-gen sequencing data to evaluate somatic MST variation (SMV) within DNA repair proficient and DNA repair defective cell lines. We confirm that alleles present within next-gen data that do not contribute to the haplotype can be reliably quantified and utilized to evaluate the SMV without requiring comparisons of matched samples. We observed that SMV patterns found in DNA repair proficient cell lines without DNA repair defects, MCF10A, HEK293 and PD20 RV:D2, had consistent patterns among samples. Further, we were able to confirm that changes in SMV patterns in cell lines lacking functional BRCA2, FANCD2 and mismatch repair were consistent with the different pathways perturbed. Using this new exome sequencing analysis approach we show that DNA instability can be identified in a sample and that patterns of instability vary depending on the impaired DNA repair mechanism, and that genes harboring minor alleles are strongly associated with cancer pathways. The MST Minor Allele Caller used for this study is available at https://github.com/zalmanv/MST_minor_allele_caller.

Microsatellite interruptions stabilize primate genomes and exist as population-specific single nucleotide polymorphisms within individual human genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.

Microsatellites are short tandem repeat DNA sequences located throughout the human genome that display a high degree of inter-individual variation. This characteristic makes microsatellites an attractive tool for population genetics and forensics research. Some microsatellites affect gene expression, and mutations within such microsatellites can cause disease. Interruption mutations disrupt the perfect repeated array and are frequently associated with altered disease risk, but they have not been thoroughly studied in human genomes. We identified interrupted mono-, di-, tri- and tetranucleotide MSs (iMS) within individual genomes from African, European, Asian and American population groups. We show that many iMSs, including some within disease-associated genes, are unique to a single population group. By measuring the conservation of microsatellites between human and chimpanzee genomes, we demonstrate that interruptions decrease the probability of microsatellite mutations throughout the genome. We demonstrate that iMSs arise in the human genome by single base changes within the DNA, and provide biochemical data suggesting that these stabilizing changes may be created by error-prone DNA polymerases. Our genome-wide study supports the model in which iMSs act to stabilize individual genomes, and suggests that population-specific differences in microsatellite architecture may be an avenue by which genetic ancestry impacts individual disease risk.

De novo assembly and annotation of Salvia splendens transcriptome using the Illumina platform

BACKGROUND

As an important perennial herbaceous flower, Salvia splendens possesses high ornamental value. Understanding its branching processes may help scientists select the best plant type. Although Salvia splendens is a frequently-used horticultural flower, only limited transcriptomic or genomic research is available in public databases. In the present study, we, for the first time, constructed a comprehensive dataset for Salvia splendens through de novo high-throughput transcriptome sequencing.

METHODOLOGY/PRINCIPAL FINDINGS

We performed de novo transcriptome sequencing on two different branching type plants (Strain 35 and Cailinghong) using the Illumina paired-end sequencing technology. For Strain 35, a total of 16,488,829 reads were generated and assembled into 38,498 unigenes, with a mean length of approximately 779 bp. For Cailinghong, 16,464,713 reads were generated and assembled into 34,302 unigenes, with a mean length of approximately 812 bp. Moreover, a total of 49,310 unigenes for Salvia splendens were identified, among them 33,925 (68.80%) were annotated in the non-redundant NCBI database, 25,371 (51.45%) were annotated in the Swiss-Prot database, while 24,888 (50.47%) and 9,896 (20.07%) unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. Using the Kyoto Encyclopedia of Genes and Genomes pathway database, we identified 134 differently expressed unigenes between Strain 35 and Cailinghong, and then these unigenes were mapped to 79 pathways. In addition, we detected 2,453 simple sequence repeats (SSRs).

CONCLUSIONS

We obtained a comprehensive transcriptomic information from this work and provided a valuable resource of transcript sequences of Salvia splendens in public databases. Moreover, some candidate genes potentially involved in branching were identified. Furthermore, numerous obtained SSRs might contribute to marker-assisted selection. These data could be further utilized in functional genomics studies on Salvia splendens.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

Distinct mutational behaviors differentiate short tandem repeats from microsatellites in the human genome

A tandem repeat's (TR) propensity to mutate increases with repeat number, and can become very pronounced beyond a critical boundary, transforming it into a microsatellite (MS). However, a clear understanding of the mutational behavior of different TR classes and motifs and related mechanisms is lacking, as is a consensus on the existence of a boundary separating short TRs (STRs) from MSs. This hinders our understanding of MSs' mutational properties and their effective use as genetic markers. Using indel calls for 179 individuals from 1000 Genomes Pilot-1 Project, we determined polymorphism incidence for four major TR classes, and formalized its varying relationship with repeat number using segmented regression. We observed a biphasic regime with a transition from a faster to a slower exponential growth at 9, 5, 4, and 4 repeats for mono-, di-, tri-, and tetranucleotide TRs, respectively. We used an in vitro mutagenesis assay to evaluate the contribution of strand slippage errors to mutability. STRs and MSs differ in their absolute polymorphism levels, but more importantly in their rates of mutability growth. Although strand slippage is a major factor driving mononucleotide polymorphism incidence, dinucleotide polymorphism incidence is greater than that expected due to strand slippage alone, indicating that additional cellular factors might be driving dinucleotide mutability in the human genome. Leveraging on hundreds of human genomes, we present the first comprehensive, genome-wide analysis of TR mutational behavior, encompassing several motif sizes and compositions.

Mature microsatellites: mechanisms underlying dinucleotide microsatellite mutational biases in human cells

Dinucleotide microsatellites are dynamic DNA sequences that affect genome stability. Here, we focused on mature microsatellites, defined as pure repeats of lengths above the threshold and unlikely to mutate below it in a single mutational event. We investigated the prevalence and mutational behavior of these sequences by using human genome sequence data, human cells in culture, and purified DNA polymerases. Mature dinucleotides (≥10 units) are present within exonic sequences of >350 genes, resulting in vulnerability to cellular genetic integrity. Mature dinucleotide mutagenesis was examined experimentally using ex vivo and in vitro approaches. We observe an expansion bias for dinucleotide microsatellites up to 20 units in length in somatic human cells, in agreement with previous computational analyses of germ-line biases. Using purified DNA polymerases and human cell lines deficient for mismatch repair (MMR), we show that the expansion bias is caused by functional MMR and is not due to DNA polymerase error biases. Specifically, we observe that the MutSα and MutLα complexes protect against expansion mutations. Our data support a model wherein different MMR complexes shift the balance of mutations toward deletion or expansion. Finally, we show that replication fork progression is stalled within long dinucleotides, suggesting that mutational mechanisms within long repeats may be distinct from shorter lengths, depending on the biochemistry of fork resolution. Our work combines computational and experimental approaches to explain the complex mutational behavior of dinucleotide microsatellites in humans.

Tumor-specific microsatellite instability: do distinct mechanisms underlie the MSI-L and EMAST phenotypes?

Microsatellite DNA sequences display allele length alterations or microsatellite instability (MSI) in tumor tissues, and MSI is used diagnostically for tumor detection and classification. We discuss the known types of tumor-specific MSI patterns and the relevant mechanisms underlying each pattern. Mutation rates of individual microsatellites vary greatly, and the intrinsic DNA features of motif size, sequence, and length contribute to this variation. MSI is used for detecting mismatch repair (MMR)-deficient tumors, which display an MSI-high phenotype due to genome-wide microsatellite destabilization. Because several pathways maintain microsatellite stability, tumors that have undergone other events associated with moderate genome instability may display diagnostic MSI only at specific di- or tetranucleotide markers. We summarize evidence for such alternative MSI forms (A-MSI) in sporadic cancers, also referred to as MSI-low and EMAST. While the existence of A-MSI is not disputed, there is disagreement about the origin and pathologic significance of this phenomenon. Although ambiguities due to PCR methods may be a source, evidence exists for other mechanisms to explain tumor-specific A-MSI. Some portion of A-MSI tumors may result from random mutational events arising during neoplastic cell evolution. However, this mechanism fails to explain the specificity of A-MSI for di- and tetranucleotide instability. We present evidence supporting the alternative argument that some A-MSI tumors arise by a distinct genetic pathway, and give examples of DNA metabolic pathways that, when altered, may be responsible for instability at specific microsatellite motifs. Finally, we suggest that A-MSI in tumors could be molecular signatures of environmental influences and DNA damage. Importantly, A-MSI occurs in several pre-neoplastic inflammatory states, including inflammatory bowel diseases, consistent with a role of oxidative stress in A-MSI. Understanding the biochemical basis of A-MSI tumor phenotypes will advance the development of new diagnostic tools and positively impact the clinical management of individual cancers.

DNA polymerase kappa microsatellite synthesis: two distinct mechanisms of slippage-mediated errors

Microsatellite tandem repeats are frequent sites of strand slippage mutagenesis in the human genome. Microsatellite mutations often occur as insertion/deletion of a repeat motif (unit-based indels), and increase in frequency with increasing repeat length after a threshold is reached. We recently demonstrated that DNA polymerase κ (Pol κ) produces fewer unit-based indel errors within dinucleotide microsatellites than does polymerase δ. Here, we examined human Pol κ's error profile within microsatellite alleles of varying sequence composition and length, using an in vitro HSV-tk gap-filling assay. We observed that Pol κ displays relatively accurate synthesis for unit-based indels, using di- and tetranucleotide repeat templates longer than the threshold length. We observed an abrupt increase in the unit-based indel frequency when the total microsatellite length exceeds 28 nucleotides, suggesting that extended Pol κ protein-DNA interactions enhance fidelity of the enzyme when synthesizing these microsatellite alleles. In contrast, Pol κ is error-prone within the HSV-tk coding sequence, producing frequent single-base errors in a manner that is highly biased with regard to sequence context. Single-nucleotide errors are also created by Pol κ within di- and tetranucleotide repeats, independently of the microsatellite allele length and at a frequency per nucleotide similar to the frequency of single base errors within the coding sequence. These single-base errors represent the mutational signature of Pol κ, and we propose them a mechanism independent of homology-stabilized slippage. Pol κ's dual fidelity nature provides a unique research tool to explore the distinct mechanisms of slippage-mediated mutagenesis.

Characterisation of tetra-nucleotide microsatellite loci in the ostrich (Struthio camelus)

1. Eleven polymorphic tetra-nucleotide microsatellite loci were identified in the ostrich (Struthio camelus) using a selective enrichment protocol. 2. The average number of alleles observed was 6·6 with an average heterozygosity of 0·4. 3. The population was found to be in Hardy-Weinberg equilibrium and two of the loci had a greater than 95% probability of having null alleles. 4. These microsatellite loci will add to the existing pool of markers available for the ostrich and help to facilitate analysis of population structure and pedigree determination.

A plasmid-based lacZα gene assay for DNA polymerase fidelity measurement

A significantly improved DNA polymerase fidelity assay, based on a gapped plasmid containing the lacZα reporter gene in a single-stranded region, is described. Nicking at two sites flanking lacZα, and removing the excised strand by thermocycling in the presence of complementary competitor DNA, is used to generate the gap. Simple methods are presented for preparing the single-stranded competitor. The gapped plasmid can be purified, in high amounts and in a very pure state, using benzoylated-naphthoylated DEAE-cellulose, resulting in a low background mutation frequency (~1 × 10(-4)). Two key parameters, the number of detectable sites and the expression frequency, necessary for measuring polymerase error rates have been determined. DNA polymerase fidelity is measured by gap filling in vitro, followed by transformation into Escherichia coli and scoring of blue/white colonies and converting the ratio to error rate. Several DNA polymerases have been used to fully validate this straightforward and highly sensitive system.

Transcriptome sequencing and de novo analysis for Ma bamboo (Dendrocalamus latiflorus Munro) using the Illumina platform

BACKGROUND

Bamboo occupies an important phylogenetic node in the grass family with remarkable sizes, woodiness and a striking life history. However, limited genetic research has focused on bamboo partially because of the lack of genomic resources. The advent of high-throughput sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provides a turning point for bamboo research. In the present study, we performed de novo transcriptome sequencing for the first time to produce a comprehensive dataset for the Ma bamboo (Dendrocalamus latiflorus Munro).

RESULTS

The Ma bamboo transcriptome was sequenced using the Illumina paired-end sequencing technology. We produced 15,138,726 reads and assembled them into 103,354 scaffolds. A total of 68,229 unigenes were identified, among which 46,087 were annotated in the NCBI non-redundant protein database and 28,165 were annotated in the Swiss-Prot database. Of these annotated unigenes, 11,921 and 10,147 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. We could map 45,649 unigenes onto 292 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. The annotated unigenes were compared against Moso bamboo, rice and millet. Unigenes that did not match any of those three sequence datasets are considered to be Ma bamboo unique. We predicted 105 unigenes encoding eight key enzymes involved in lignin biosynthesis. In addition, 621 simple sequence repeats (SSRs) were detected.

CONCLUSION

Our data provide the most comprehensive transcriptomic resource currently available for D. latiflorus Munro. Candidate genes potentially involved in growth and development were identified, and those predicted to be unique to Ma bamboo are expected to give a better insight on Ma bamboo gene diversity. Numerous SSRs characterized contributed to marker development. These data constitute a new valuable resource for genomic studies on D. latiflorus Munro and, more generally, bamboo.

Generation and analysis of expressed sequence tags from a normalized cDNA library of young leaf from Ma bamboo (Dendrocalamus latiflorus Munro)

Ma bamboo (Dendrocalamus latiflorus Munro) belongs to Dendrocalamus genus, Bambusease tribe, Bambusoideae subfamily, Poaceae family. It is a representative species of clumping bamboo, and a principal commercial species for various construction purposes using mature culms and for human consumption using young shoots. A normalized cDNA library was constructed from young leaves of Ma bamboo and 9,574 high-quality ESTs were generated, from which 5,317 unigenes including 1,502 contigs and 3,815 singletons were assembled. The unigenes were assigned into different gene ontology (GO) categories and summarized into 13 broad biologically functional groups according to similar functional characteristics or cellular roles by BLAST search against public databases. Eight hundred and ninety-one unigenes were assigned by KO identifiers and mapped to six KEGG biochemical pathways. The transcripts involved in biosynthesis of secondary metabolites such as cytochrome 450, flavonol synthase/flavanone 3-hydroxylase, and dihydroflavonol-4-reductase were well represented by 14 unigenes in the unigene set. The candidate genes involved in phytohormone metabolism, signal transduction and encoding cell wall-associated receptor kinases were also identified. Sixty-seven unigenes related to plant resistance (R) genes, including RPP genes, RGAs and RDL/RF genes, were discovered. These results will provide genome-wide knowledge about the molecular physiology of Ma bamboo young leaves and tools for advanced studies of molecular mechanism underlying leaf growth and development.

Beyond translesion synthesis: polymerase κ fidelity as a potential determinant of microsatellite stability

Microsatellite DNA synthesis represents a significant component of human genome replication that must occur faithfully. However, yeast replicative DNA polymerases do not possess high fidelity for microsatellite synthesis. We hypothesized that the structural features of Y-family polymerases that facilitate accurate translesion synthesis may promote accurate microsatellite synthesis. We compared human polymerases κ (Pol κ) and η (Pol η) fidelities to that of replicative human polymerase δ holoenzyme (Pol δ4), using the in vitro HSV-tk assay. Relative polymerase accuracy for insertion/deletion (indel) errors within 2-3 unit repeats internal to the HSV-tk gene concurred with the literature: Pol δ4 >> Pol κ or Pol η. In contrast, relative polymerase accuracy for unit-based indel errors within [GT](10) and [TC](11) microsatellites was: Pol κ ≥ Pol δ4 > Pol η. The magnitude of difference was greatest between Pols κ and δ4 with the [GT] template. Biochemically, Pol κ displayed less synthesis termination within the [GT] allele than did Pol δ4. In dual polymerase reactions, Pol κ competed with either a stalled or moving Pol δ4, thereby reducing termination. Our results challenge the ideology that pol κ is error prone, and suggest that DNA polymerases with complementary biochemical properties can function cooperatively at repetitive sequences.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

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.

DNA pol λ's extraordinary ability to stabilize misaligned DNA

DNA polymerases have the venerable task of maintaining genome stability during DNA replication and repair. Errors, nonetheless, occur with error propensities that are polymerase specific. For example, DNA polymerase λ (pol λ) generates single-base deletions through template-strand slippage within short repetitive DNA regions much more readily than does the closely related polymerase β (pol β). Here we present in silico evidence to help interpret pol λ's greater tendency for deletion errors than pol β by its more favorable protein/DNA electrostatic interactions immediately around the extrahelical nucleotide on the template strand. Our molecular dynamics and free energy analyses suggest that pol λ provides greater stabilization to misaligned DNA than aligned DNA. Our study of several pol λ mutants of Lys544 (Ala, Phe, Glu) probes the interactions between the extrahelical nucleotide and the adjacent Lys544 to show that the charge of the 544 residue controls stabilization of the DNA misalignment. In addition, we identify other thumb residues (Arg538, Lys521, Arg517, and Arg514) that play coordinating roles in stabilizing pol λ's interactions with misaligned DNA. Interestingly, their aggregate stabilization effect is more important than that of any one component residue, in contrast to aligned DNA systems, as we determined from mutations of these key residues and energetic analyses. No such comparable network of stabilizing misaligned DNA exists in pol β. Evolutionary needs for DNA repair on substrates with minimal base-pairing, such as those encountered by pol λ in the non-homologous end-joining pathway, may have been solved by a greater tolerance to deletion errors. Other base-flipping proteins share similar binding properties and motions for extrahelical nucleotides.

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.

Plasmid-based lacZalpha assay for DNA polymerase fidelity: application to archaeal family-B DNA polymerase

The preparation of a gapped pUC18 derivative, containing the lacZα reporter gene in the single-stranded region, is described. Gapping is achieved by flanking the lacZα gene with sites for two related nicking endonucleases, enabling the excision of either the coding or non-coding strand. However, the excised strand remains annealed to the plasmid through non-covalent Watson–Crick base-pairing; its removal, therefore, requires a heat–cool cycle in the presence of an exactly complementary competitor DNA. The gapped plasmids can be used to assess DNA polymerase fidelity using in vitro replication, followed by transformation into Escherichia coli and scoring the blue/white colony ratio. Results found with plasmids are similar to the well established method based on gapped M13, in terms of background (∼0.08% in both cases) and the mutation frequencies observed with a number of DNA polymerases, providing validation for this straightforward and technically uncomplicated approach. Several error prone variants of the archaeal family-B DNA polymerase from Pyrococcus furiosus have been investigated, illuminating the potential of the method.

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.

Loop II of DNA polymerase beta is important for discrimination during substrate binding

Loop II of DNA polymerase beta (pol beta) consists of 14 amino acid residues and is highly flexible and solvent exposed. Previous research from our laboratory has shown that this loop is important for polymerase activity and fidelity. In the study presented here, we demonstrate that a shortened five amino acid residue loop compromises the fidelity of pol beta. This five-residue loop, termed ENEYP, induces one base frameshift errors and A-C transversions within a specific sequence context. We demonstrate that ENEYP misincorporates dGTP opposite template A at higher efficiencies than wild-type pol beta. The kinetic basis for misincorporation is a defect in discrimination of the correct from incorrect dNTP substrate at the level of ground-state binding. Our results are consistent with the idea that loop II of pol beta functions to maintain accurate DNA synthesis by a direct or indirect influence on the nucleotide binding pocket.

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.

DNA triple helices: biological consequences and therapeutic potential

DNA structure is a critical element in determining its function. The DNA molecule is capable of adopting a variety of non-canonical structures, including three-stranded (i.e. triplex) structures, which will be the focus of this review. The ability to selectively modulate the activity of genes is a long-standing goal in molecular medicine. DNA triplex structures, either intermolecular triplexes formed by binding of an exogenously applied oligonucleotide to a target duplex sequence, or naturally occurring intramolecular triplexes (H-DNA) formed at endogenous mirror repeat sequences, present exploitable features that permit site-specific alteration of the genome. These structures can induce transcriptional repression and site-specific mutagenesis or recombination. Triplex-forming oligonucleotides (TFOs) can bind to duplex DNA in a sequence-specific fashion with high affinity, and can be used to direct DNA-modifying agents to selected sequences. H-DNA plays important roles in vivo and is inherently mutagenic and recombinogenic, such that elements of the H-DNA structure may be pharmacologically exploitable. In this review we discuss the biological consequences and therapeutic potential of triple helical DNA structures. We anticipate that the information provided will stimulate further investigations aimed toward improving DNA triplex-related gene targeting strategies for biotechnological and potential clinical applications.

DNA polymerase kappa produces interrupted mutations and displays polar pausing within mononucleotide microsatellite sequences

Microsatellites are ubiquitously present in eukaryotic genomes and are implicated as positive factors in evolution. At the nucleotide level, microsatellites undergo slippage events that alter allele length and base changes that interrupt the repetitive tract. We examined DNA polymerase errors within a [T]₁₁ microsatellite using an in vitro assay that preferentially detects mutations other than unit changes. We observed that human DNA polymerase kappa (Pol κ) inserts dGMP and dCMP within the [T]₁₁ mononucleotide repeat, producing an interrupted 12-bp allele. Polymerase β produced such interruptions at a lower frequency. These data demonstrate that DNA polymerases are capable of directly producing base interruptions within microsatellites. At the molecular level, expanded microsatellites have been implicated in DNA replication fork stalling. Using an in vitro primer extension assay, we observed sequence-specific synthesis termination by DNA polymerases within mononucleotides. Quantitatively, intense, polar pausing was observed for both pol κ and polymerase α-primase within a [T]₁₁ allele. A mechanism is proposed in which pausing results from DNA bending within the duplex stem of the nascent DNA. Our data support the concept of a microsatellite life-cycle, and are consistent with the models in which DNA sequence or secondary structures contributes to non-uniform rates of replication fork progression.

The genome-wide determinants of human and chimpanzee microsatellite evolution

Mutation rates of microsatellites vary greatly among loci. The causes of this heterogeneity remain largely enigmatic yet are crucial for understanding numerous human neurological diseases and genetic instability in cancer. In this first genome-wide study, the relative contributions of intrinsic features and regional genomic factors to the variation in mutability among orthologous human-chimpanzee microsatellites are investigated with resampling and regression techniques. As a result, we uncover the intricacies of microsatellite mutagenesis as follows. First, intrinsic features (repeat number, length, and motif size), which all influence the probability and rate of slippage, are the strongest predictors of mutability. Second, mutability increases nonuniformly with length, suggesting that processes additional to slippage, such as faulty repair, contribute to mutations. Third, mutability varies among microsatellites with different motif composition likely due to dissimilarities in secondary DNA structure formed by their slippage intermediates. Fourth, mutability of mononucleotide microsatellites is impacted by their location on sex chromosomes vs. autosomes and inside vs. outside of Alu repeats, the former confirming the importance of replication and the latter suggesting a role for gene conversion. Fifth, transcription status and location in a particular isochore do not influence microsatellite mutability. Sixth, compared with intrinsic features, regional genomic factors have only minor effects. Finally, our regression models explain approximately 90% of variation in microsatellite mutability and can generate useful predictions for the studies of human diseases, forensics, and conservation genetics.

Loop II of DNA polymerase beta is important for polymerization activity and fidelity

The accurate replication and transmission of genetic information is critical in the life of an organism. During its entire lifespan, the genetic information is constantly under attack from endogenous and exogenous sources of damage. To ensure that the content of its genetic information is faithfully preserved for synthesis and transmission, eukaryotic cells have developed a complex system of genomic quality control. Key players in this process are DNA polymerases, the enzymes responsible for synthesizing the DNA, because errors introduced into the genome by polymerase can result in mutations. We use DNA polymerase beta (pol β) as a model system to investigate mechanisms of preserving fidelity during nucleotide incorporation. In the study described here, we characterized the role that loop II of pol β plays in maintaining the activity and fidelity of pol β. We report here that the absence or shortening of loop II compromises the catalytic activity of pol β. Our data also show that loop variants of a specific length have a lower fidelity when compared to the wild-type polymerase. Taken together, our results indicate that loop II is important for the catalytic activity and fidelity of pol β.

Cytological diploidization and rapid genome changes of the newly synthesized allotetraploids Cucumis x hytivus

We used a newly synthesized allotetraploid between C. sativus (2n = 2x = 14, n gametic chromosome number, x haploid chromosome number) and C. hystrix (2n = 2x = 24) to study the genomic events in its early generations. Results from cytological characterization of the F(1) and the allotetraploid progenies showed that the rate of bivalents in meiotic metaphase I of the F(1) was greatly improved by chromosome doubling, and further improved during the selfing process of allopolyploid resulting into relatively diploid-like meiosis. Extensive genomic changes were detected by amplified fragment length polymorphism analysis. The changes mainly involved loss of parental restriction fragments and gaining of novel fragments. The total detectable changes were from 11.1 to 32.1%, and the frequency of losing parental fragments was much higher than that of gaining novel fragments. Some of the changes were initiated as early as in the F(1) hybrid, whereas others occurred after chromosome doubling (polyploid formation). No significant differences were detected in the reciprocal F(1) hybrids and S(0) generations. But the data showed that the frequency of sequence losing in C. sativus was about two times higher than in the C. hystrix. Our results demonstrated that the sequence elimination was the major event of genomic changes, and it might provide the physical basis for the diploid-like meiotic behavior in the diploidization of the newly formed allopolyploids. Moreover, the results suggest that the sequence elimination was not caused by cytoplasmic factors, and might relate to genomic recombination and to the numbers of parental chromosome.

Sequence-based methods for identifying epidemiologically linked herpes simplex virus type 2 strains

Traditional methods for confirming the identity of herpes simplex virus (HSV) isolates use restriction fragment length polymorphism (RFLP). However, RFLP is less amenable to high-throughput analyses of many samples, and the extent to which small differences in RFLP patterns distinguish between different viral strains remains unclear. Viral HSV type 2 (HSV-2) DNA isolates from 14 persons experiencing a primary HSV-2 infection and from their sexual partners were analyzed by RFLP and heteroduplex mobility assays. We also compared the HSV-2 sequences from seven regions, including noncoding regions between UL19 and UL20, UL24 and UL25, UL37 and UL38, and UL41 and UL42 and coding segments of the gC, gB, and gG genes. Although the resulting RFLP patterns of the couples were almost identical, minor banding differences existed between the source and susceptible partners in five couples. Heteroduplex mobility assays were unable to distinguish between unrelated strains. Overall, 22 sites of sequence variation were found in 1,482 bp of analyzed sequence. The DNA sequences differentiated between all unrelated infections, and epidemiologically related isolates had identical sequences in all but two pairs. Our results suggest that a multilocus assay based on several DNA sequences has the potential to be an informative tool for identifying epidemiologically related HSV-2 strains.

Microsatellite evolution: Markov transition functions for a suite of models

This paper takes from the collection of models considered by Whittaker et al. [2003. Likelihood-based estimation of microsatellite mutation rates. Genetics 164, 781-787] derived from direct observation of microsatellite mutation in parent-child pairs and provides analytical expressions for the probability distributions for the change in number of repeats over any given number of generations. The mathematical framework for this analysis is the theory of Markov processes. We find these expressions using two approaches, approximating by circulant matrices and solving a partial differential equation satisfied by the generating function. The impact of the differing choice of models is examined using likelihood estimates for time to most recent common ancestor. The analysis presented here may play a role in elucidating the connections between these two approaches and shows promise in reconciling differences between estimates for mutation rates based on Whittaker's approach and methods based on phylogenetic analyses.

Review of capture-recapture methods applicable to noninvasive genetic sampling

The use of noninvasive genetic sampling to identify individual animals for capture-recapture studies has become widespread in the past decade. Strong emphasis has been placed on the field protocols and genetic analyses with fruitful results. Little attention has been paid to the capture-recapture application for this specific type of data beyond stating the effects of assumption violations. Here, we review the broad class of capture-recapture methods that are available for use with DNA-based capture-recapture data, noting the array of biologically interesting parameters such as survival, emigration rates, state transition rates and the finite rate of population change that can be estimated from such data. We highlight recent developments in capture-recapture theory specifically designed for noninvasive genetic sampling data.

Patterns of sequence loss and cytosine methylation within a population of newly resynthesized Brassica napus allopolyploids

Allopolyploid formation requires the adaptation of two nuclear genomes within a single cytoplasm, which may involve programmed genetic and epigenetic changes during the initial generations following genome fusion. To study the dynamics of genome change, we synthesized 49 isogenic Brassica napus allopolyploids and surveyed them with 76 restriction fragment length polymorphism (RFLP) probes and 30 simple sequence repeat (SSR) primer pairs. Here, we report on the types and distribution of genetic and epigenetic changes within the S(1) genotypes. We found that insertion/deletion (indel) events were rare, but not random. Of the 57,710 (54,383 RFLP and 3,327 SSR) parental fragments expected among the amphidiploids, we observed 56,676 or 99.9%. Three loci derived from Brassica rapa had indels, and one indel occurred repeatedly across 29% (14/49) of the lines. Loss of one parental fragment was due to the 400-bp reduction of a guanine-adenine dinucleotide repeat-rich sequence. In contrast to the 4% (3/76) RFLP probes that detected indels, 48% (35/73) detected changes in the CpG methylation status between parental genomes and the S1 lines. Some loci were far more likely than others to undergo epigenetic change, but the number of methylation changes within each synthetic polyploid was remarkably similar to others. Clear de novo methylation occurred at a much higher frequency than de novo demethylation within allopolyploid sequences derived from B. rapa. Our results suggest that there is little genetic change in the S(0) generation of resynthesized B. napus polyploids. In contrast, DNA methylation was altered extensively in a pattern that indicates tight regulation of epigenetic changes.

Positive correlation between DNA polymerase alpha-primase pausing and mutagenesis within polypyrimidine/polypurine microsatellite sequences

Microsatellite DNA sequences are ubiquitous in the human genome, and mutation rates of these repetitive sequences vary with respect to DNA sequence as well as length. We have analyzed polymerase-DNA interactions as a function of microsatellite sequence, using polypyrimidine/polypurine di- and tetranucleotide alleles representative of those found in the human genome. Using an in vitro primer extension assay and the mammalian DNA polymerase alpha-primase complex, we have observed a polymerase termination profile for each microsatellite that is unique to that allele. Interestingly, a periodic termination profile with an interval size (9-11 nucleotides) unrelated to microsatellite unit length was observed for the [TC](20) and [TTCC](9) templates. In contrast, a unit-punctuated polymerase termination profile was found for the longer polypurine templates. We detected strong polymerase pauses within the [TC](20) allele at low reaction pH which were eliminated by the addition of deaza-dGTP, consistent with these specific pauses being a consequence of triplex DNA formation during DNA synthesis. Quantitatively, a strand bias was observed in the primer extension assay, in that polymerase synthesis termination is more intense when the polypurine sequence serves as the template, relative to its complementary polypyrimidine sequence. The HSV-tk forward mutation assay was utilized to determine the corresponding polymerase alpha-primase error frequencies and specificities at the microsatellite alleles. A higher microsatellite polymerase error frequency (50x10(-4) to 60x10(-4)) was measured when polypurine sequences serve as templates for DNA synthesis, relative to the polypyrimidine template (18x10(-4)). Thus, a positive correlation exists between polymerase alpha-primase pausing and mutagenesis within microsatellite DNA alleles.

An Analysis of Microsatellite Loci in Arabidopsis thaliana: Mutational Dynamics and Application

Microsatellite loci are among the most commonly used molecular markers. These loci typically exhibit variation for allele frequency distribution within a species. However, the factors contributing to this variation are not well understood. To expand on the current knowledge of microsatellite evolution, 20 microsatellite loci were examined for 126 accessions of the flowering plant, Arabidopsis thaliana. Substantial variability in mutation pattern among loci was found, most of which cannot be explained by the assumptions of the traditional stepwise mutation model or infinite alleles model. Here it is shown that the degree of locus diversity is strongly correlated with the number of contiguous repeats, more so than with the total number of repeats. These findings support a strong role for repeat disruptions in stabilizing microsatellite loci by reducing the substrate for polymerase slippage and recombination. Results of cluster analyses are also presented, demonstrating the potential of microsatellite loci for resolving relationships among accessions of A. thaliana.