Estimation of spontaneous mutation rates

Spontaneous or randomly occurring mutations play a key role in cancer progression. Estimation of the mutation rate of cancer cells can provide useful information about the disease. To ascertain these mutation rates, we need mathematical models that describe the distribution of mutant cells. In this investigation, we develop a discrete time stochastic model for a mutational birth process. We assume that mutations occur concurrently with mitosis so that when a nonmutant parent cell splits into two progeny, one of these daughter cells could carry a mutation. We propose an estimator for the mutation rate and investigate its statistical properties via theory and simulations. A salient feature of this estimator is the ease with which it can be computed. The methods developed herein are applied to a human colorectal cancer cell line and compared to existing continuous time models.     

Statistical and algorithmic methods for fluctuation analysis with SALVADOR as an implementation

This paper aims at removing certain long-standing impediments to more effective and widespread use of fluctuation analysis. The paper presents a method of constructing confidence intervals for mutation rates using data from fluctuation experiments. The method was inspired by a rediscovery of a little-known, not fully developed method of Lea and Coulson; substantial modifications have been made both to enhance computational efficiency and to widen the scope of the original method's applicability. A computer package named SALVADOR is presented that can be used for Monte Carlo simulation, for point and interval estimation of mutation rates, and for exploration of various hypotheses spawned by the directed mutation controversy. In addition to the maximum likelihood method, methods of considerable historical interest are also examined and included in SALVADOR to help the reader compare and assess some of the most popular methods for estimating mutation rates.     

Rapid Method for Assessing Oxygen Consumption Rate of Cells from Transient-state Measurements of Pericellular Dissolved Oxygen Concentration

Recently we described a method for estimating the oxygen consumption rate (OCR) of cells in static culture from equilibrium measurements of dissolved oxygen concentration (dO₂), using an oxygen-sensing microplate and the steady-state solution to Fick's Law (Guarino et al. 2004). Here we describe a complementary method for estimating OCR from the transient-state rate of change of measured dO₂. Although the system is open to the atmosphere and subject to a significant lag in sensor response, the rate of change of the measured dO₂ immediately after seeding correlated directly with both cell number and steady-state OCR. This transient-state method is linear with cell number to a much higher density than is possible with the steady-state method because it derives from measurements made before diffusion limitations can be established. For a given sensor thickness, the same correlation line between the transient and non-diffusion-limited steady-state estimates of OCR was found to apply for various preparations of rat hepatocytes. The correlation slope varied predictably with sensor thickness. Thus, despite the non-idealities of this system, the initial rate measurement offers a rapid method to obtain an estimate of absolute OCR. To demonstrate the utility of this method, we purposefully treated rat hepatocytes in ways expected to change OCR. Cells deprived of oxygen by storage under several centimeters of medium showed decreases in both OCR and viability with time. Likewise, the OCR of hepatocytes exposed to the oxidative phosphorylation inhibitor rotenone decreased, whereas those exposed to the uncoupler dinoseb increased.     

The adequacy of aging techniques in vertebrates for rapid estimation of population mortality rates from age distributions

As a key parameter in population dynamics, mortality rates are frequently estimated using mark–recapture data, which requires extensive, long‐term data sets. As a potential rapid alternative, we can measure variables correlated to age, allowing the compilation of population age distributions, from which mortality rates can be derived. However, most studies employing such techniques have ignored their inherent inaccuracy and have thereby failed to provide reliable mortality estimates. In this study, we present a general statistical model linking birth rate, mortality rate, and population age distributions. We next assessed the reliability and data needs (i.e., sample size) for estimating mortality rate of eight different aging techniques. The results revealed that for half of the aging techniques, correlations with age varied considerably, translating into highly variable accuracies when used to estimate mortality rate from age distributions. Telomere length is generally not sufficiently correlated to age to provide reliable mortality rate estimates. DNA methylation, signal‐joint T‐cell recombination excision circle (sjTREC), and racemization are generally more promising techniques to ultimately estimate mortality rate, if a sufficiently high sample size is available. Otolith ring counts, otolithometry, and age‐length keys in fish, and skeletochronology in reptiles, mammals, and amphibians, outperformed all other aging techniques and generated relatively accurate mortality rate estimation with a sample size that can be feasibly obtained. Provided the method chosen is minimizing and estimating the error in age estimation, it is possible to accurately estimate mortality rates from age distributions. The method therewith has the potential to estimate a critical, population dynamic parameter to inform conservation efforts within a limited time frame as opposed to mark–recapture analyses.     

Birth and population prevalence of Duchenne muscular dystrophy in the Netherlands

Summary Mutations causing Duchenne muscular dystrophy (DMD) have a short survival. Therefore, birth and population prevalence are maintained by new mutations. The present inventory was made to estimate the birth and population prevalence rates of DMD in the Netherlands. Seven methods of case identification were used. Data on 496 definite, probable or possible DMD patients born since 1961, or alive on January 1, 1983, were obtained. Several methods gave an estimated ascertainment of more than 95%. The prevalence rate at birth of DMD was estimated at 23.7×10^–5 (14215) male live births (MLB) yearly. The prevalence rate in the male population on January 1, 1983 was 5.4×10^–5 (118496). About 1% of the males in this study may have autosomal recessive Duchenne-like muscular dystrophy. Until now there has been no convincing evidence for geographic differences in DMD prevalence at birth. A list of frequency studies of Duchenne muscular dystrophy is included. The DMD mutation rate calculated by the indirect method is 7.9×10^–5 genes per generation. However, this may well be an over-estimate, as this method does not account for germline mosaicism. Using a modified sex ratio method the proportion of sporadic DMD among all cases was estimated to be 0.106 (range 0–0.332). High frequency of germline mosaicism in DMD is a likely cause for the apparent lack of sporadic cases as found in previous studies, if mutation rates in male and female gametes are equal. Therefore, methods for estimating the proportion of new mutants in DMD should take germline mosaicism into account. The modified sex ratio method allows incorporation of data on germline mosaicism if available.     

THE METAPHASE ARREST TECHNIQUE

The accuracy and precision of the stathmokinetic experiment for the determination of cell birth rate are discussed from practical and theoretical viewpoints. Topics covered include the determination of the optimal dosage of the metaphase arresting agent, the times at which observations should be taken and which cells should be counted, the correct method of estimating cell birth rate and the dependence of any derived estimate of the turnover time on the age distribution of cycling cells. Failure to consider any of these points can lead to considerable inaccuracy. It is further shown that the stathmokinetic method is a rather imprecise one for measuring the turnover time and the duration of mitosis. Data are present comparing the results of using the technique with those obtained from autoradiographic studies.     

一种确定rm的简单方法的验证与分析

本文对一种测定蚜虫和螨的内禀增长率(rm)的简单方法进行了验证与分析。这种方法不要求详细的生殖力表资料,计算rm的公式如下:rm=0.74 ln(Md)/d。式中,d表示生殖前期,即从出生到第一次生殖之间的天数;Md表示在2d期间每个雌虫繁殖的平均雌性后代数。这种方法是Wyatt和White提出的。我们的研究表明:当Md小于1,这种方法不宜采用;当Md大于1和2d期间的繁殖贡献接近或大于70％时,可以用这种方法计算rm;不过,所有蚜虫和螨种群在2d期间对rm的繁殖贡献是否能够达到70％或70％以上还有待研究。研究繁殖分布特征对于这种简单方法的实际应用有重要意义。     

The Coalescent Process with Selfing

A method for estimating the selfing rate using DNA sequence data was recently proposed by Milligan. Unfortunately, a number of errors make interpretation of his results problematic. In the present paper we first show how the usual coalescent process can be adapted to models that include selfing, and then use this result to find moment estimators as well as the likelihood surface for the selfing rate, s, and the scaled mutation rate, θ. We conclude that, regardless of the method used, large sample sizes are necessary to estimate s with any degree of certainty, and that the estimate is always highly sensitive to recent changes in the true value.     

Estimating mutation rate and generation time from longitudinal samples of DNA sequences

We present in this paper a simple method for estimating the mutation rate per site per year which also yields an estimate of the length of a generation when mutation rate per site per generation is known. The estimator, which takes advantage of DNA polymorphisms in longitudinal samples, is unbiased under a number of population models, including population structure and variable population size over time. We apply the new method to a longitudinal sample of DNA sequences of the env gene of human immunodeficiency virus type 1 (HIV-1) from a single patient and obtain 1.62 x 10(-2) as the mutation rate per site per year for HIV-1. Using an independent data set to estimate the mutation rate per generation, we obtain 1.8 days as the length of a generation of HIV-1, which agrees well with recent estimates based on viral load data. Our estimate of generation time differs considerably from a recent estimate by Rodrigo et al. when the same mutation rate per site per generation is used. Some factors that may contribute to the difference among different estimators are discussed.     

The influence of selection on the evolutionary distance estimated from the base changes observed between homologous nucleotide sequences

In most studies of molecular evolution, the nucleotide base at a site is assumed to change with the apparent rate under functional constraint, and the comparison of base changes between homologous genes is thought to yield the evolutionary distance corresponding to the site-average change rate multiplied by the divergence time. However, this view is not sufficiently successful in estimating the divergence time of species, but mostly results in the construction of tree topology without a time-scale. In the present paper, this problem is investigated theoretically by considering that observed base changes are the results of comparing the survivals through selection of mutated bases. In the case of weak selection, the time course of base changes due to mutation and selection can be obtained analytically, leading to a theoretical equation showing how the selection has influence on the evolutionary distance estimated from the enumeration of base changes. This result provides a new method for estimating the divergence time more accurately from the observed base changes by evaluating both the strength of selection and the mutation rate. The validity of this method is verified by analysing the base changes observed at the third codon positions of amino acid residues with four-fold codon degeneracy in the protein genes of mammalian mitochondria; i.e. the ratios of estimated divergence times are fairly well consistent with a series of fossil records of mammals. Throughout this analysis, it is also suggested that the mutation rates in mitochondrial genomes are almost the same in different lineages of mammals and that the lineage-specific base-change rates indicated previously are due to the selection probably arising from the preference of transfer RNAs to codons.     

Mutation Rate Variability across Human Y-Chromosome Haplogroups

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.     

鲤鱼微卫星突变速率和模式(英文)

利用150个微卫星分子标记在F1代家系的基因型分析过程中,共有27600个等位基因从亲本向子代传递,其中在5个微卫星座位上检测到6个突变的等位基因。对突变的等位基因数目进行统计分析后得出:鲤鱼平均每个世代每个微卫星座位的突变速率为2.53×10-4。在发现突变的5个位点中,经测序发现,突变序列中插入1个以上的重复单元就导致了突变的发生。这些突变表明,鲤鱼的微卫星突变没有遵循严格的渐变突变模型(stepwise mutation model,SMM)。该文关于鲤鱼微卫星突变速率和模式的研究将会对统计鲤鱼有效群体的统计提供有效参数。     

Estimation of mutation rate from rare protein variants.

A method for estimating the mutation rate for protein loci from the number of rare alleles in the population is presented. It seems to have a number of advantages compared with Kimura and Ohta's method. Applying this method to Neel's data from American Indians in South America and to Nozawa's data from Japanese macaques, the mutation rate for electrophoretically detectable alleles is estimated to be (2 approximately 3) x 10(-6) per locus per generation. This estimate may not include many severely or substantially deleterious mutations.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Sequence diversity under the multispecies coalescent with Yule process and constant population size

The study of sequence diversity under phylogenetic models is now classic. Theoretical studies of diversity under the Kingman coalescent appeared shortly after the introduction of the coalescent. In this paper we revisit this topic under the multispecies coalescent, an extension of the single population model to multiple populations. We derive exact formulas for the sequence dissimilarity of two sequences drawn at random under a basic multispecies setup. The multispecies model uses three parameters—the species tree birth rate under the pure birth process (Yule), the species effective population size and the mutation rate. We also discuss the effects of relaxing some of the model assumptions.     

High mutation rate of TPE repeats: a microsatellite in the putative transposase of the hobo element in Drosophila melanogaster

The hobo transposable element contains a polymorphic microsatellite sequence located in its coding region, the TPE repeats. Previous surveys of natural populations of Drosophila melanogaster have detected at least seven different hobo transposons. These natural populations are geographically structured with regard to TPE polymorphism, and a scenario has been proposed for the invasion process. Natural populations have recently been completely invaded by hobo elements with three TPE repeats. New elements then appeared by mutation, triggering a new stage of invasion by other elements. Since TPE polymorphism appeared over a short period of time, we focused on estimating the mutation rate of these TPE repeats. We used transgenic lines harboring three TPE and/or five TPE hobo elements that had been evolving for at least 16 generations to search for a new TPE repeat polymorphism. We detected three mutants, with four, seven, and eight TPE repeats, respectively. The estimated mutation rate of the TPE repeats is therefore higher than that of neutral microsatellites in D. melanogaster (4.2 x 10-4 versus 6.5 x 10-6). The role of the transposition mechanism and the particular structure of the TPE repeats of the hobo element in this increase in the mutation rate are discussed.     

A method for estimating cell volume of amoebae based on measurements of cell length of motile forms: physiological and ecological applications

A novel method of estimating the cell volume of non-testate, free-living amoebae is presented. Measurements of the length (L) of motile forms of some amoebae, obtained by applying some basic rules, can be used to calculate diameter (D) using a relatively simple formula, derived from regression analysis (N=39): D=0.6L. Important physiological and ecological variables, including carbon content and respiratory rate, can then be easily calculated using this estimated cell volume.     

2, 3, 5-triphenyl tetrazolium chloride (TTC) reduction as exponential growth phase marker for mammalian cells in culture and for myeloma hybridization experiments

Triphenyl tetrazolium chloride in vitro reduction by cells produces a red formazan pellet which can be extracted and measured. We have shown that such reduction is associated with animal cell growth, and particularly with the specific growth rate, so the measurement of Triphenyl tetrazolium chloride reduction is proposed as a physiological marker of the exponential growth of cultured cells. Further application of this technique is shown using this Redox reaction for estimating plasmacytoma fusion potential for hybridoma cell line production.Abbreviations TTC 2,3,5-Triphenyl Tetrazolium Chloride     

The influence of premeiotic clusters of mutation on indirect estimations of mutation rate

Based on the hypothesis that rare alleles are in mutation and random loss equilibrium, mutation rate can be indirectly estimated by measuring the number of rare variants and the average existing time of a mutant allele. This method can be applied to estimate the mutation rate in humans. However, this estimation of mutation rate is affected by the presence of premeiotic clusters of mutation. Mutation clusters change both the number of initial mutants and the average existing time of a mutant allele. As a result, the formula indirectly estimating mutation rate should be modified. The influence of premeiotic clusters is more obvious when the population size is small or the average cluster size is big. For example, if the population size is 3,000 and average cluster size is two, instead of one, the mutation rate is increased by about 9.4%.     

胚胎移植技术对特种突变小鼠的净化研究

目的通过净化特种突变小鼠来建立SPF级种子库。方法利用胚胎移植技术对四种突变无毛小鼠和一种白内障小鼠的桑葚胚进行子宫内移植。结果四种突变无毛小鼠和白内障小鼠均产出仔鼠,获得移植成功,其中白内障小鼠的产仔率最高为25.5%。结论建立了SPF级特种突变小鼠种群,为进一步研究与开发利用提供了可靠保证。