Leveraging Distant Relatedness to Quantify Human Mutation and Gene-Conversion Rates

The rate at which human genomes mutate is a central biological parameter that has many implications for our ability to understand demographic and evolutionary phenomena. We present a method for inferring mutation and gene-conversion rates by using the number of sequence differences observed in identical-by-descent (IBD) segments together with a reconstructed model of recent population-size history. This approach is robust to, and can quantify, the presence of substantial genotyping error, as validated in coalescent simulations. We applied the method to 498 trio-phased sequenced Dutch individuals and inferred a point mutation rate of 1.66 × 10⁻⁸ per base per generation and a rate of 1.26 × 10⁻⁹ for <20 bp indels. By quantifying how estimates varied as a function of allele frequency, we inferred the probability that a site is involved in non-crossover gene conversion as 5.99 × 10⁻⁶. We found that recombination does not have observable mutagenic effects after gene conversion is accounted for and that local gene-conversion rates reflect recombination rates. We detected a strong enrichment of recent deleterious variation among mismatching variants found within IBD regions and observed summary statistics of local sharing of IBD segments to closely match previously proposed metrics of background selection; however, we found no significant effects of selection on our mutation-rate estimates. We detected no evidence of strong variation of mutation rates in a number of genomic annotations obtained from several recent studies. Our analysis suggests that a mutation-rate estimate higher than that reported by recent pedigree-based studies should be adopted in the context of DNA-based demographic reconstruction.     

Whole Genome Sequencing of Mutation Accumulation Lines Reveals a Low Mutation Rate in the Social Amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10⁻⁹, with a Poisson confidence interval of 4.1×10⁻⁹ − 9.5×10⁻⁹, per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10⁻¹¹, with a Poisson confidence interval ranging from 7.4×10⁻¹³ to 1.6×10⁻¹⁰, is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.     

Probability Distribution of Drug-Resistant Mutants in Unselected Populations of Mycobacterium tuberculosis

The fluctuation test shows that Mycobacterium tuberculosis mutates to resistance to isoniazid, streptomycin, ethambutol and rifampin spontaneously and at random. The average mutation rates for the drugs, in the same order, were calculated to be 2.56 × 10⁻⁸, 2.95 × 10⁻⁸, 10⁻⁷, and 2.25 × 10⁻¹⁰ mutation per bacterium per generation. The relatively high mutation rate to ethambutol resistance and the low mutation rate to rifampin resistance were confirmed by analyzing the increase in the proportion of mutants with time in a growing population of the tubercle bacilli. The highest proportions of mutants to be expected in unselected populations of the tubercle bacilli were calculated from the results of fluctuation tests.     

Phenotypic mutation rates and the abundance of abnormal proteins in yeast

Phenotypic mutations are errors that occur during protein synthesis. These errors lead to amino acid substitutions that give rise to abnormal proteins. Experiments suggest that such errors are quite common. We present a model to study the effect of phenotypic mutation rates on the amount of abnormal proteins in a cell. In our model, genes are regulated to synthesize a certain number of functional proteins. During this process, depending on the phenotypic mutation rate, abnormal proteins are generated. We use data on protein length and abundance in Saccharomyces cerevisiae to parametrize our model. We calculate that for small phenotypic mutation rates most abnormal proteins originate from highly expressed genes that are on average nearly twice as large as the average yeast protein. For phenotypic mutation rates much above 5 × 10⁻⁴, the error-free synthesis of large proteins is nearly impossible and lowly expressed, very large proteins contribute more and more to the amount of abnormal proteins in a cell. This fact leads to a steep increase of the amount of abnormal proteins for phenotypic mutation rates above 5 × 10⁻⁴. Simulations show that this property leads to an upper limit for the phenotypic mutation rate of approximately 2 × 10⁻³ even if the costs for abnormal proteins are extremely low. We also consider the adaptation of individual proteins. Individual genes/proteins can decrease their phenotypic mutation rate by using preferred codons or by increasing their robustness against amino acid substitutions. We discuss the similarities and differences between the two mechanisms and show that they can only slow down but not prevent the rapid increase of the amount of abnormal proteins. Our work allows us to estimate the phenotypic mutation rate based on data on the fraction of abnormal proteins. For S. cerevisiae, we predict that the value for the phenotypic mutation rate is between 2 × 10⁻⁴ and 6 × 10⁻⁴.     

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.     

Empirical estimates of the mutation rate for an alphabaculovirus

Mutation rates are of key importance for understanding evolutionary processes and predicting their outcomes. Empirical mutation rate estimates are available for a number of RNA viruses, but few are available for DNA viruses, which tend to have larger genomes. Whilst some viruses have very high mutation rates, lower mutation rates are expected for viruses with large genomes to ensure genome integrity. Alphabaculoviruses are insect viruses with large genomes and often have high levels of polymorphism, suggesting high mutation rates despite evidence of proofreading activity by the replication machinery. Here, we report an empirical estimate of the mutation rate per base per strand copying (s/n/r) of Autographa californica multiple nucleopolyhedrovirus (AcMNPV). To avoid biases due to selection, we analyzed mutations that occurred in a stable, non-functional genomic insert after five serial passages in Spodoptera exigua larvae. Our results highlight that viral demography and the stringency of mutation calling affect mutation rate estimates, and that using a population genetic simulation model to make inferences can mitigate the impact of these processes on estimates of mutation rate. We estimated a mutation rate of μ = 1×10⁻⁷ s/n/r when applying the most stringent criteria for mutation calling, and estimates of up to μ = 5×10⁻⁷ s/n/r when relaxing these criteria. The rates at which different classes of mutations accumulate provide good evidence for neutrality of mutations occurring within the inserted region. We therefore present a robust approach for mutation rate estimation for viruses with stable genomes, and strong evidence of a much lower alphabaculovirus mutation rate than supposed based on the high levels of polymorphism observed.     

Rates of Microbial Metabolism in Deep Coastal Plain Aquifers

Rates of microbial metabolism in deep anaerobic aquifers of the Atlantic coastal plain of South Carolina were investigated by both microbiological and geochemical techniques. Rates of [2-¹⁴C]acetate and [U-¹⁴C]glucose oxidation as well as geochemical evidence indicated that metabolic rates were faster in the sandy sediments composing the aquifers than in the clayey sediments of the confining layers. In the sandy aquifer sediments, estimates of the rates of CO₂ production (millimoles of CO₂ per liter per year) based on the oxidation of [2-¹⁴C] acetate were 9.4 × 10⁻³ to 2.4 × 10⁻¹ for the Black Creek aquifer, 1.1 × 10⁻² for the Middendorf aquifer, and <7 × 10⁻⁵ for the Cape Fear aquifer. These estimates were at least 2 orders of magnitude lower than previously published estimates that were based on the accumulation of CO₂ in laboratory incubations of similar deep subsurface sediments. In contrast, geochemical modeling of groundwater chemistry changes along aquifer flowpaths gave rate estimates that ranged from 10⁻⁴ to 10⁻⁶ mmol of CO₂ per liter per year. The age of these sediments (ca. 80 million years) and their organic carbon content suggest that average rates of CO₂ production could have been no more than 10⁻⁴ mmol per liter per year. Thus, laboratory incubations may greatly overestimate the in situ rates of microbial metabolism in deep subsurface environments. This has important implications for the use of laboratory incubations in attempts to estimate biorestoration capacities of deep aquifers. The rate estimates from geochemical modeling indicate that deep aquifers are among the most oligotrophic aquatic environments in which there is ongoing microbial metabolism.     

Estimation of the Spontaneous Mutation Rate per Nucleotide Site in a Drosophila melanogaster Full-Sib Family

We employed deep genome sequencing of two parents and 12 of their offspring to estimate the mutation rate per site per generation in a full-sib family of Drosophila melanogaster recently sampled from a natural population. Sites that were homozygous for the same allele in the parents and heterozygous in one or more offspring were categorized as candidate mutations and subjected to detailed analysis. In 1.23 × 10⁹ callable sites from 12 individuals, we confirmed six single nucleotide mutations. We estimated the false negative rate in the experiment by generating synthetic mutations using the empirical distributions of numbers of nonreference bases at heterozygous sites in the offspring. The proportion of synthetic mutations at callable sites that we failed to detect was <1%, implying that the false negative rate was extremely low. Our estimate of the point mutation rate is 2.8 × 10⁻⁹ (95% confidence interval = 1.0 × 10⁻⁹ − 6.1 × 10⁻⁹) per site per generation, which is at the low end of the range of previous estimates, and suggests an effective population size for the species of ∼1.4 × 10⁶. At one site, point mutations were present in two individuals, indicating that there had been a premeiotic mutation cluster, although surprisingly one individual had a G→A transition and the other a G→T transversion, possibly associated with error-prone mismatch repair. We also detected three short deletion mutations and no insertions, giving a deletion mutation rate of 1.2 × 10⁻⁹ (95% confidence interval = 0.7 × 10⁻⁹ − 11 × 10⁻⁹).     

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10⁻⁷ mutations/site/year (95% confidence interval 3.75 × 10⁻⁸–1.12 × 10⁻⁶). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods.     

Variable Spontaneous Mutation and Loss of Heterozygosity among Heterozygous Genomes in Yeast

Mutation and recombination are the primary sources of genetic variation. To better understand the evolution of genetic variation, it is crucial to comprehensively investigate the processes involving mutation accumulation and recombination. In this study, we performed mutation accumulation experiments on four heterozygous diploid yeast species in the Saccharomycodaceae family to determine spontaneous mutation rates, mutation spectra, and losses of heterozygosity (LOH). We observed substantial variation in mutation rates and mutation spectra. We also observed high LOH rates (1.65–11.07×10⁻⁶ events per heterozygous site per cell division). Biases in spontaneous mutation and LOH together with selection ultimately shape the variable genome-wide nucleotide landscape in yeast species.     

Calibrating the Human Mutation Rate via Ancestral Recombination Density in Diploid Genomes

The human mutation rate is an essential parameter for studying the evolution of our species, interpreting present-day genetic variation, and understanding the incidence of genetic disease. Nevertheless, our current estimates of the rate are uncertain. Most notably, recent approaches based on counting de novo mutations in family pedigrees have yielded significantly smaller values than classical methods based on sequence divergence. Here, we propose a new method that uses the fine-scale human recombination map to calibrate the rate of accumulation of mutations. By comparing local heterozygosity levels in diploid genomes to the genetic distance scale over which these levels change, we are able to estimate a long-term mutation rate averaged over hundreds or thousands of generations. We infer a rate of 1.61 ± 0.13 × 10⁻⁸ mutations per base per generation, which falls in between phylogenetic and pedigree-based estimates, and we suggest possible mechanisms to reconcile our estimate with previous studies. Our results support intermediate-age divergences among human populations and between humans and other great apes.     

Impact of Genetic Background on Allele Selection in a Highly Mutable Candida albicans Gene,PNG2

In many microbes rapid mutation of highly mutable contingency genes continually replenishes a pool of variant alleles from which the most suitable are selected, assisting in rapid adaptation and evasion of the immune response. In some contingency genes mutability is achieved through DNA repeats within the coding region. The fungal human pathogen Candida albicans has 2600 repeat-containing ORFs. For those investigated (ALS genes, HYR1, HYR2, CEK1, RLM1) many protein variants with differing amino acid repeat regions exist, as expected for contingency genes. However, specific alleles dominate in different clades, which is unexpected if allele variation is used for short-term adaptation. Generation of new alleles of repeat-containing C. albicans ORFs has never been observed directly. Here we present evidence for restrictions on the emergence of new alleles in a highly mutable C. albicans repeat-containing ORF, PNG2, encoding a putative secreted or cell surface glycoamidase. In laboratory cultures new PNG2 alleles arose at a rate of 2.8×10⁻⁵ (confidence interval 3.3×10⁻⁶−9. 9×10⁻⁵) per cell per division, comparable to rates measured for contingency genes. Among 80 clinical isolates 17 alleles of different length and 23 allele combinations were distinguishable; sequence differences between repeat regions of identical size suggest the existence of 36 protein variants. Specific allele combinations predominated in different genetic backgrounds, as defined by DNA fingerprinting and multilocus sequence typing. Given the PNG2 mutation rate, this is unexpected, unless in different genetic backgrounds selection favors different alleles. Specific alleles or allele combinations were not preferentially associated with C. albicans isolates from particular body sites or geographical regions. Our results suggest that the mutability of PNG2 is not used for short-term adaptation or evasion of the immune system. Nevertheless the large number of alleles observed indicates that mutability of PNG2 may assist C. albicans strains from different genetic backgrounds optimize their interaction with the host in the long term.     

Amino Acids Transitioning of 2009 H1N1pdm in Taiwan from 2009 to 2011

A swine-origin influenza A was detected in April 2009 and soon became the 2009 H1N1 pandemic strain (H1N1pdm). The current study revealed the genetic diversity of H1N1pdm, based on 77 and 70 isolates which we collected, respectively, during the 2009/2010 and 2010/2011 influenza seasons in Taiwan. We focused on tracking the amino acid transitioning of hemagglutinin (HA) and neuraminidase (NA) genes in the early diversification of the virus and compared them with H1N1pdm strains reported worldwide. We identified newly emerged mutation markers based on A/California/04/2009, described how these markers shifted from the first H1N1pdm season to the one that immediately followed, and discussed how these observations may relate to antigenicity, receptor-binding, and drug susceptibility. It was found that the amino acid mutation rates of H1N1pdm were elevated, from 9.29×10⁻³ substitutions per site in the first season to 1.46×10⁻² in the second season in HA, and from 5.23×10⁻³ to 1.10×10⁻² in NA. Many mutation markers were newly detected in the second season, including 11 in HA and 8 in NA, and some were found having statistical correlation to disease severity. There were five noticeable HA mutations made to antigenic sites. No significant titer changes, however, were detected based on hemagglutination inhibition tests. Only one isolate with H275Y mutation known to reduce susceptibility to NA inhibitors was detected. As limited Taiwanese H1N1pdm viruses were isolated after our sampling period, we gathered 8,876 HA and 6,017 NA H1N1pdm sequences up to April 2012 from NCBI to follow up the dynamics of mentioned HA mutations. While some mutations described in this study seemed to either settle in or die out in the 2011–2012 season, a number of them still showed signs of transitioning, prompting the importance of continuous monitoring of this virus for more seasons to come.     

Kinetics of Leptin Binding to the Q223R Leptin Receptor

Studies in human populations and mouse models of disease have linked the common leptin receptor Q223R mutation to obesity, multiple forms of cancer, adverse drug reactions, and susceptibility to enteric and respiratory infections. Contradictory results cast doubt on the phenotypic consequences of this variant. We set out to determine whether the Q223R substitution affects leptin binding kinetics using surface plasmon resonance (SPR), a technique that allows sensitive real-time monitoring of protein-protein interactions. We measured the binding and dissociation rate constants for leptin to the extracellular domain of WT and Q223R murine leptin receptors expressed as Fc-fusion proteins and found that the mutant receptor does not significantly differ in kinetics of leptin binding from the WT leptin receptor. (WT: k_a 1.76×10⁶±0.193×10⁶ M⁻¹ s⁻¹, k_d 1.21×10⁻⁴±0.707×10⁻⁴ s⁻¹, K_D 6.47×10⁻¹¹±3.30×10⁻¹¹ M; Q223R: k_a 1.75×10⁶±0.0245×10⁶ M⁻¹ s⁻¹, k_d 1.47×10⁻⁴±0.0505×10⁻⁴ s⁻¹, K_D 8.43×10⁻¹¹±0.407×10⁻¹¹ M). Our results support earlier findings that differences in affinity and kinetics of leptin binding are unlikely to explain mechanistically the phenotypes that have been linked to this common genetic variant. Future studies will seek to elucidate the mechanism by which this mutation influences susceptibility to metabolic, infectious, and malignant pathologies.     

Upper-limit mutation rate estimation for a plant RNA virus

It is generally accepted that mutation rates of RNA viruses are inherently high due to the lack of proofreading mechanisms. However, direct estimates of mutation rate are surprisingly scarce, in particular for plant viruses. Here, based on the analysis of in vivo mutation frequencies in tobacco etch virus, we calculate an upper-bound mutation rate estimation of 3×10⁻⁵ per site and per round of replication; a value which turns out to be undistinguishable from the methodological error. Nonetheless, the value is barely on the lower side of the range accepted for RNA viruses, although in good agreement with the only direct estimate obtained for other plant viruses. These observations suggest that, perhaps, differences in the selective pressures operating during plant virus evolution may have driven their mutation rates towards values lower than those characteristic of other RNA viruses infecting bacteria or animals.     

Rates of transposition in Escherichia coli

The evolutionary role of transposable elements (TEs) is still highly controversial. Two key parameters, the transposition rate (u and w, for replicative and non-replicative transposition) and the excision rate (e) are fundamental to understanding their evolution and maintenance in populations. We have estimated u, w and e for six families of TEs (including eight members: IS1, IS2, IS3, IS4, IS5, IS30, IS150 and IS186) in Escherichia coli, using a mutation accumulation (MA) experiment. In this experiment, mutations accumulate essentially at the rate at which they appear, during a period of 80 500 (1610 generations × 50 lines) generations, and spontaneous transposition events can be detected. This differs from other experiments in which insertions accumulated under strong selective pressure or over a limited genomic target. We therefore provide new estimates for the spontaneous rates of transposition and excision in E. coli. We observed 25 transposition and three excision events in 50 MA lines, leading to overall rate estimates of u ∼ 1.15 × 10^–5, w ∼ 4 × 10⁻⁸ and e ∼ 1.08 × 10⁻⁶ (per element, per generation). Furthermore, extensive variation between elements was found, consistent with previous knowledge of the mechanisms and regulation of transposition for the different elements.     

The Comparative Genomics of Human Respiratory Syncytial Virus Subgroups A and B: Genetic Variability and Molecular Evolutionary Dynamics

Genomic variation and related evolutionary dynamics of human respiratory syncytial virus (RSV), a common causative agent of severe lower respiratory tract infections, may affect its transmission behavior. RSV evolutionary patterns are likely to be influenced by a precarious interplay between selection favoring variants with higher replicative fitness and variants that evade host immune responses. Studying RSV genetic variation can reveal both the genes and the individual codons within these genes that are most crucial for RSV survival. In this study, we conducted genetic diversity and evolutionary rate analyses on 36 RSV subgroup B (RSV-B) whole-genome sequences. The attachment protein, G, was the most variable protein; accordingly, the G gene had a higher substitution rate than other RSV-B genes. Overall, less genetic variability was found among the available RSV-B genome sequences than among RSV-A genome sequences in a comparable sample. The mean substitution rates of the two subgroups were, however, similar (for subgroup A, 6.47 × 10⁻⁴ substitutions/site/year [95% credible interval {CI 95%}, 5.56 × 10⁻⁴ to 7.38 × 10⁻⁴]; for subgroup B, 7.76 × 10⁻⁴ substitutions/site/year [CI 95%, 6.89 × 10⁻⁴ to 8.58 × 10⁻⁴]), with the time to their most recent common ancestors (TMRCAs) being much lower for RSV-B (19 years) than for RSV-A (46.8 years). The more recent RSV-B TMRCA is apparently the result of a genetic bottleneck that, over longer time scales, is still compatible with neutral population dynamics. Whereas the immunogenic G protein seems to require high substitution rates to ensure immune evasion, strong purifying selection in conserved proteins such as the fusion protein and nucleocapsid protein is likely essential to preserve RSV viability.     

Variation in RNA Virus Mutation Rates across Host Cells

It is well established that RNA viruses exhibit higher rates of spontaneous mutation than DNA viruses and microorganisms. However, their mutation rates vary amply, from 10⁻⁶ to 10⁻⁴ substitutions per nucleotide per round of copying (s/n/r) and the causes of this variability remain poorly understood. In addition to differences in intrinsic fidelity or error correction capability, viral mutation rates may be dependent on host factors. Here, we assessed the effect of the cellular environment on the rate of spontaneous mutation of the vesicular stomatitis virus (VSV), which has a broad host range and cell tropism. Luria-Delbrück fluctuation tests and sequencing showed that VSV mutated similarly in baby hamster kidney, murine embryonic fibroblasts, colon cancer, and neuroblastoma cells (approx. 10⁻⁵ s/n/r). Cell immortalization through p53 inactivation and oxygen levels (1–21%) did not have a significant impact on viral replication fidelity. This shows that previously published mutation rates can be considered reliable despite being based on a narrow and artificial set of laboratory conditions. Interestingly, we also found that VSV mutated approximately four times more slowly in various insect cells compared with mammalian cells. This may contribute to explaining the relatively slow evolution of VSV and other arthropod-borne viruses in nature.     

Elemental Composition (C, N, P) and Cell Volume of Exponentially Growing and Nutrient-Limited Bacterioplankton

Marine bacterioplankton were isolated and grown in batch cultures until their growth became limited by organic carbon (C), nitrogen (N), or phosphorus (P). Samples were taken from the cultures at both the exponential and stationary phases. The elemental composition of individual bacterial cells was analyzed by X-ray microanalysis with an electron microscope. The cell size was also measured. The elemental content was highest in exponentially growing cells (149 ± 8 fg of C cell⁻¹, 35 ± 2 fg of N cell⁻¹, and 12 ± 1 fg of P cell⁻¹; average of all isolates ± standard error). The lowest C content was found in C-limited cells (39 ± 3 fg of C cell⁻¹), the lowest N content in C- and P-limited cells (12 ± 1 and 12 ± 2 fg of N cell⁻¹, respectively), and the lowest P content in P-limited cells (2.3 ± 0.6 fg of P cell⁻¹). The atomic C:N ratios varied among treatments between 3.8 ± 0.1 and 9.5 ± 1.0 (average ± standard error), the C:P ratios between 35 ± 2 and 178 ± 28, and the N:P ratios between 6.7 ± 0.3 and 18 ± 3. The carbon-volume ratios showed large variation among isolates due to different types of nutrient limitation (from 51± 4 to 241 ± 38 fg of C μm⁻¹; average of individual isolates and treatments ± standard error). The results show that different growth conditions and differences in the bacterial community may explain some of the variability of previously reported elemental and carbon-volume ratios.