The mitochondrial and autosomal mutation landscapes of prostate cancer

BACKGROUND: Prostate cancer (PCa) is the most common cancer in men. PCa is strongly age associated; low death rates in surveillance cohorts call into question the widespread use of surgery, which leads to overtreatment and a reduction in quality of life. There is a great need to increase the understanding of tumor characteristics in the context of disease progression.

OBJECTIVE: To perform the first multigenome investigation of PCa through analysis of both autosomal and mitochondrial DNA, and to integrate exome sequencing data, and RNA sequencing and copy-number alteration (CNA) data to investigate how various different tumor characteristics, commonly analyzed separately, are interconnected.

DESIGN, SETTING, AND PARTICIPANTS: Exome sequencing was applied to 64 tumor samples from 55 PCa patients with varying stage and grade. Integrated analysis was performed on a core set of 50 tumors from which exome sequencing, CNA, and RNA sequencing data were available.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Genes, mutated at a significantly higher rate relative to a genomic background, were identified. In addition, mitochondrial and autosomal mutation rates were correlated to CNAs and proliferation, assessed as a cell cycle gene expression signature.

RESULTS AND LIMITATIONS: Genes not previously reported to be significantly mutated in PCa, such as cell division cycle 27 homolog (Saccharomyces cerevisiae) (CDC27), myeloid/lymphoid or mixed-lineage leukemia 3 (MLL3), lysine (K)-specific demethylase 6A (KDM6A), and kinesin family member 5A (KIF5A) were identified. The mutation rate in the mitochondrial genome was 55 times higher than that of the autosomes. Multilevel analysis demonstrated a tight correlation between high reactive-oxygen exposure, chromosomal damage, high proliferation, and in parallel, a transition from multiclonal indolent primary PCa to monoclonal aggressive disease. As we only performed targeted sequence analysis; copy-number neutral rearrangements recently described for PCa were not accounted for.

CONCLUSIONS: The mitochondrial genome displays an elevated mutation rate compared to the autosomal chromosomes. By integrated analysis, we demonstrated that different tumor characteristics are interconnected, providing an increased understanding of PCa etiology.

Male mutation bias and possible long-term effects of human activities.

The ability of a population to adapt to changing environments depends critically on the amount and kind of genetic variability it possesses. Mutations are an important source of new genetic variability and may lead to new adaptations, especially if the population size is large. Mutation rates are extremely variable between and within species, and males usually have higher mutation rates as a result of elevated rates of male germ cell division. This male bias affects the overall mutation rate. We examined the factors that influence male mutation bias, and focused on the effects of classical life-history parameters, such as the average age at reproduction and elevated rates of sperm production in response to sexual selection and sperm competition. We argue that human-induced changes in age at reproduction or in sexual selection will affect male mutation biases and hence overall mutation rates. Depending on the effective population size, these changes are likely to influence the long-term persistence of a population.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Epstein-Barr virus (EBV) nuclear antigen (EBNA)-4 mutation in EBV-associated malignancies in three different populations

Different ethnic groups with a high human leukocyte antigen (HLA)-A11 prevalence have been shown to experience a high rate of Epstein-Barr virus (EBV) infection, EBV-associated malignancies, and Epstein-Barr nuclear antigen (EBNA)-4 mutations. The epitopes 393-408 and 416-424 of EBNA-4 are major antigenic epitopes that elicit an HLA-A11 cytotoxic T lymphocyte (CTL) response to EBV infection. Mutations selectively involving one or more nucleotide residues in these epitopes affect the antigenicity of EBNA-4, because the mutant EBV strains are not recognized by the HLA-A11-restricted CTLs. To investigate these mutations in common EBV-associated malignancies occurring in different populations, we studied the mutation rate of epitopes 393-408 and 416-424 of EBNA-4 in 25 cases of EBV-associated Hodgkin's disease (HD), nine cases of AIDS-related non-Hodgkin's lymphoma, and 37 cases of EBV-associated gastric carcinoma (GC) from the United States, Brazil, and Japan. We found one or more mutations in these two epitopes in 50% (6/12) of United States HD, 15% (2/13) of Brazilian HD, 50% (6/12) United States GC and 28% (7/25) Japanese GC, and 22% (2/9) of United States AIDS-lymphoma. Similar mutations were found in 30% (3/10) of United States reactive, 0% (0/6) of Brazilian reactive, and 25% (2/8) Japanese reactive tissues. The most frequent amino acid substitutions were virtually identical to those seen in previously reported isolates from EBV-associated nasopharyngeal carcinomas and Burkitt's lymphomas occurring in high prevalence HLA-A11 regions. However, only 2/28 (7%) mutations occurred in HLA-A11-positive patients. Our studies suggest that: 1) EBNA-4 mutations are a common phenomenon in EBV-associated HD, GC, and AIDS-lymphoma; 2) the mutation rate does not vary in these geographic areas and ethnic groups; 3) EBNA-4 mutations in EBV-associated United States and Brazilian HD, United States and Japanese GC, and United States AIDS lymphomas are not related to patients' HLA-A11 status.

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

NRAS and BRAF mutations in primary cutaneous melanoma: A comparison of mutation rates between radial and vertical growth phases in individual tumors

Primary cutaneous melanoma is a cancer arising from melanocytes in the skin. In recent decades the incidence of this malignancy has increased significantly. Mortality rates are high for patients with tumors measuring over a few millimeters in thickness. Response rates to conventional radiation and chemotherapy are very low in patients with metastatic melanoma. New therapies targeting melanoma’s aberrant cell signaling pathways such as the MAP Kinase pathway are being developed. Mutations of NRAS and BRAF genes are quite common in cutaneous melanoma and lead to constitutive activation of the MAP Kinase pathway. This study tests the hypothesis that NRAS and BRAF mutations increase as a tumor progresses from the noninvasive radial growth phase (RGP) to the invasive vertical growth phase (VGP). Laser capture microdissection was used to obtain separate, pure tumor DNA samples from the RGP and VGP of thirty primary cutaneous melanomas. PCR was used to amplify NRAS exon 2 and BRAF exon 15 tumor DNA. The amplified DNA was sequenced and analyzed for mutations. An overall mutation rate of 74% was obtained for the twenty-three melanomas in which there were complete sequence results. With the exception of one melanoma NRAS and BRAF mutations were mutually exclusive. All seven NRAS exon 2 mutations involved codon 61. Three of these melanomas had mutations in both the RGP and VGP. The remaining four tumors were wild type for NRAS exon 2 in the RGP but mutated in the VGP. Of the fifteen BRAF exon 15 mutated melanomas all but one involved codon 600. Twelve of the fifteen BRAF exon 15 mutations were the T1799A type. Nine of the fifteen BRAF mutated tumors had the same mutation in both the RGP and VGP. Five of fifteen melanomas had wild type RGP DNA and BRAF exon 15 mutated VGP DNA. A single melanoma had BRAF exon 15 mutated DNA in the RGP and wild type DNA in the VGP. Overall, these results suggest a trend toward the acquisition of NRAS and BRAF mutations as cutaneous melanomas change from a noninvasive to an invasive, potentially deadly cancer.^

Mutation rates among RNA viruses

The rate of spontaneous mutation is a key parameter in modeling the genetic structure and evolution of populations. The impact of the accumulated load of mutations and the consequences of increasing the mutation rate are important in assessing the genetic health of populations. Mutation frequencies are among the more directly measurable population parameters, although the information needed to convert them into mutation rates is often lacking. A previous analysis of mutation rates in RNA viruses (specifically in riboviruses rather than retroviruses) was constrained by the quality and quantity of available measurements and by the lack of a specific theoretical framework for converting mutation frequencies into mutation rates in this group of organisms. Here, we describe a simple relation between ribovirus mutation frequencies and mutation rates, apply it to the best (albeit far from satisfactory) available data, and observe a central value for the mutation rate per genome per replication of μg ≈ 0.76. (The rate per round of cell infection is twice this value or about 1.5.) This value is so large, and ribovirus genomes are so informationally dense, that even a modest increase extinguishes the population.

Stage specificity, dose response, and doubling dose for mouse minisatellite germ-line mutation induced by acute radiation

Germ-line mutation induction at mouse minisatellite loci by acute irradiation with x-rays was studied at premeiotic and postmeiotic stages of spermatogenesis. An elevated paternal mutation rate was found after irradiation of premeiotic spermatogonia and stem cells, whereas the frequency of minisatellite mutation after postmeiotic irradiation of spermatids was similar to that in control litters. In contrast, paternal irradiation did not affect the maternal mutation rate. A linear dose–response curve for paternal mutation induced at premeiotic stages was found, with a doubling dose of 0.33 Gy, a value close to those obtained in mice after acute spermatogonia irradiation using other systems for mutation detection. High frequencies of spontaneous and induced mutations at minisatellite loci allow mutation induction to be evaluated at low doses of exposure in very small population samples, which currently makes minisatellite DNA the most powerful tool for monitoring radiation-induced germ-line mutation.

Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius

Microbes whose genomes are encoded by DNA and for which adequate information is available display similar genomic mutation rates (average 0.0034 mutations per chromosome replication, range 0.0025 to 0.0046). However, this value currently is based on only a few well characterized microbes reproducing within a narrow range of environmental conditions. In particular, no genomic mutation rate has been determined either for a microbe whose natural growth conditions may extensively damage DNA or for any member of the archaea, a prokaryotic lineage deeply diverged from both bacteria and eukaryotes. Both of these conditions are met by the extreme thermoacidophile Sulfolobus acidocaldarius. We determined the genomic mutation rate for this species when growing at pH 3.5 and 75°C based on the rate of forward mutation at the pyrE gene and the nucleotide changes identified in 101 independent mutants. The observed value of about 0.0018 extends the range of DNA-based microbes with rates close to the standard rate simultaneously to an archaeon and to an extremophile whose cytoplasmic pH and normal growth temperature greatly accelerate the spontaneous decomposition of DNA. The mutations include base pair substitutions (BPSs) and additions and deletions of various sizes, but the S. acidocaldarius spectrum differs from those of other DNA-based organisms in being relatively poor in BPSs. The paucity of BPSs cannot yet be explained by known properties of DNA replication or repair enzymes of Sulfolobus spp. It suggests, however, that molecular evolution per genome replication may proceed more slowly in S. acidocaldarius than in other DNA-based organisms examined to date.

Enhanced somatic mutation rates induced in stem cells of mice by low chronic exposure to ethylnitrosourea.

We have found that the somatic mutation rate at the Dlb-1 locus increases exponentially during low daily exposure to ethylnitrosourea over 4 months. This effect, enhanced mutagenesis, was not observed at a lacI transgene in the same tissue, although the two loci respond very similarly to acute doses. Since both mutations are neutral, the mutant frequency was expected to increase linearly with time in response to a constant mutagenic exposure, as it did for lacI. Enhanced mutagenesis does not result from an overall sensitization of the animals, since mice that had first been treated with a low daily dose for 90 days and then challenged with a large acute dose were not sensitized to the acute dose. Nor was the increased mutant frequency due to selection, since animals that were treated for 90 days and then left untreated for up to 60 days showed little change from the 90-day frequency. The effect is substantial: about 8 times as many Dlb-1 mutants were induced between 90 and 120 days as in the first 30 days. This resulted in a reverse dose rate effect such that 90 mg/kg induced more mutants when delivered at 1 mg/kg per day than at 3 mg/kg per day. We postulate that enhanced mutagenesis arises from increased stem cell proliferation and the preferential repair of transcribed genes. Enhanced mutagenesis may be important for risk evaluation, as the results show that chronic exposures can be more mutagenic than acute ones and raise the possibility of synergism between chemicals at low doses.

Modeling and optimization of populations subject to time-dependent mutation.

It has become clear that many organisms possess the ability to regulate their mutation rate in response to environmental conditions. So the question of finding an optimal mutation rate must be replaced by that of finding an optimal mutation schedule. We show that this task cannot be accomplished with standard population-dynamic models. We then develop a "hybrid" model for populations experiencing time-dependent mutation that treats population growth as deterministic but the time of first appearance of new variants as stochastic. We show that the hybrid model agrees well with a Monte Carlo simulation. From this model, we derive a deterministic approximation, a "threshold" model, that is similar to standard population dynamic models but differs in the initial rate of generation of new mutants. We use these techniques to model antibody affinity maturation by somatic hypermutation. We had previously shown that the optimal mutation schedule for the deterministic threshold model is phasic, with periods of mutation between intervals of mutation-free growth. To establish the validity of this schedule, we now show that the phasic schedule that optimizes the deterministic threshold model significantly improves upon the best constant-rate schedule for the hybrid and Monte Carlo models.

Increased mutation frequency of feline immunodeficiency virus lacking functional deoxyuridine-triphosphatase.

Feline immunodeficiency virus (FIV) encodes the enzyme deoxyuridine-triphosphatase (DU; EC 3.6.1.23) between the coding regions for reverse transcriptase and integrase in the pol gene. Here, we report the in vivo infection of cats with a DU- variant of the PPR strain of FIV and compare its growth properties and tissue distribution with those of wild-type FIV-PPR. The results reveal several important points: (i) DU- FIV is able to infect the cat, with kinetics similar to that observed with wild-type FIV; (ii) both wild-type and DU- FIV-infected specific-pathogen free cats mount a strong humoral antibody response which is able to limit the virus burden in both groups of animals; (iii) the virus burden is reduced in the DU- FIV-infected cats, particularly in tissues such as spleen and salivary gland; and (iv) the mutation frequency in DU- FIVs integrated in the DNA of primary macrophages after 9 months of infection is approximately 5-fold greater than the frequency observed in DU- FIV DNA integrated in T lymphocytes. Mutation rate with wild-type FIV remains the same in both cell types in vivo. The dominant mutations seen in macrophages with DU- FIV are G-->A base changes, consistent with an increased misincorporation of deoxyuridine into viral DNA of DU- FIVs during reverse transcription. Because this enzyme is absent from human immunodeficiency virus type 1 and other primate lentiviruses, virus replication in cell environments with low DU activity may lead to increased mutation and contribute to the rapid expansion of the viral repertoire.

Mutation at the Human D1S80 Minisatellite Locus

Little is known about the general biology of minisatellites. The purpose of this study is to examine repeat mutations from the D1S80 minisatellite locus by sequence analysis to elucidate the mutational process at this locus. This is a highly polymorphic minisatellite locus, located in the subtelomeric region of chromosome 1. We have analyzed 90,000 human germline transmission events and found seven (7) mutations at this locus. The D1S80 alleles of the parentage trio, the child, mother, and the alleged father were sequenced and the origin of the mutation was determined. Using American Association of Blood Banks (AABB) guidelines, we found a male mutation rate of 1 . 0 4 × 1 0− 4 and a female mutation rate of 5 . 1 8 × 1 0− 5 with an overall mutation rate of approximately 7 . 7 7 × 1 0− 5. Also, in this study, we found that the identified mutations are in close proximity to the center of the repeat array rather than at the ends of the repeat array. Several studies have examined the mutational mechanisms of the minisatellites according to infinite allele model (IAM) and the one-step stepwise mutation model (SMM). In this study, we found that this locus fits into the one-step mutation model (SMM) mechanism in six out of seven instances similar to STR loci.

Variable mutation rates as an adaptive strategy in replicator populations

For evolving populations of replicators, there is much evidence that the effect of mutations on fitness depends on the degree of adaptation to the selective pressures at play. In optimized populations, most mutations have deleterious effects, such that low mutation rates are favoured. In contrast to this, in populations thriving in changing environments a larger fraction of mutations have beneficial effects, providing the diversity necessary to adapt to new conditions. What is more, non-adapted populations occasionally benefit from an increase in the mutation rate. Therefore, there is no optimal universal value of the mutation rate and species attempt to adjust it to their momentary adaptive needs. In this work we have used stationary populations of RNA molecules evolving in silico to investigate the relationship between the degree of adaptation of an optimized population and the value of the mutation rate promoting maximal adaptation in a short time to a new selective pressure. Our results show that this value can significantly differ from the optimal value at mutation-selection equilibrium, being strongly influenced by the structure of the population when the adaptive process begins. In the short-term, highly optimized populations containing little variability respond better to environmental changes upon an increase of the mutation rate, whereas populations with a lower degree of optimization but higher variability benefit from reducing the mutation rate to adapt rapidly. These findings show a good agreement with the behaviour exhibited by actual organisms that replicate their genomes under broadly different mutation rates. © 2010 Stich et al.