How important is DNA replication for mutagenesis?

Transcriptional and mutational signatures of the Drosophila ageing germline

Ageing is a complex biological process that is accompanied by changes in gene expression and mutational load. In many species, including humans, older fathers pass on more paternally derived de novo mutations; however, the cellular basis and cell types driving this pattern are still unclear. To explore the root causes of this phenomenon, we performed single-cell RNA sequencing on testes from young and old male Drosophila and genomic sequencing (DNA sequencing) on somatic tissues from the same flies. We found that early germ cells from old and young flies enter spermatogenesis with similar mutational loads but older flies are less able to remove mutations during spermatogenesis. Mutations in old cells may also increase during spermatogenesis. Our data reveal that old and young flies have distinct mutational biases. Many classes of genes show increased postmeiotic expression in the germlines of older flies. Late spermatogenesis-biased genes have higher dN/dS (ratio of non-synonymous to synonymous substitutions) than early spermatogenesis-biased genes, supporting the hypothesis that late spermatogenesis is a source of evolutionary innovation. Surprisingly, genes biased in young germ cells show higher dN/dS than genes biased in old germ cells. Our results provide new insights into the role of the germline in de novo mutation.

Machine Learning Techniques for Classifying the Mutagenic Origins of Point Mutations

There is increasing interest in developing diagnostics that discriminate individual mutagenic mechanisms in a range of applications that include identifying population-specific mutagenesis and resolving distinct mutation signatures in cancer samples. Analyses for these applications assume that mutagenic mechanisms have a distinct relationship with neighboring bases that allows them to be distinguished. Direct support for this assumption is limited to a small number of simple cases, e.g., CpG hypermutability. We have evaluated whether the mechanistic origin of a point mutation can be resolved using only sequence context for a more complicated case. We contrasted single nucleotide variants originating from the multitude of mutagenic processes that normally operate in the mouse germline with those induced by the potent mutagen N-ethyl-N-nitrosourea (ENU). The considerable overlap in the mutation spectra of these two samples make this a challenging problem. Employing a new, robust log-linear modeling method, we demonstrate that neighboring bases contain information regarding point mutation direction that differs between the ENU-induced and spontaneous mutation variant classes. A logistic regression classifier exhibited strong performance at discriminating between the different mutation classes. Concordance between the feature set of the best classifier and information content analyses suggest our results can be generalized to other mutation classification problems. We conclude that machine learning can be used to build a practical classification tool to identify the mutation mechanism for individual genetic variants. Software implementing our approach is freely available under an open-source license.

Genetic mosaicism in haemophilia: A practical review to help evaluate the risk of transmitting the disease

Approximately 70% of patients with haemophilia exhibit a clear inheritance pattern, while for the remaining 30%, patients are the first to be diagnosed in their family and are considered sporadic cases. In such a setting, the determination of carrier status and the risk estimation of disease transmission to another child are major challenges for genetic counselling. Large studies have suggested that genetic testing reveals 70% of sporadic patients' mothers are carriers. In the remaining 30%, in some apparently non-carrier mothers, the pathogenic variant can be detected as low somatic and gonosomal mosaicism. The significance of mosaic pathogenic variants has thus far been underestimated, since conventional Sanger sequencing and other technology are not sufficiently sensitive. The study of various tissue samples and recent extra-sensitive molecular methods have now made it easier to detect both single-nucleotide variants (SNVs) and copy-number variants (CNVs), at a mosaic level in parents, and to predict the probability of disease recurrence. This review seeks to examine various kinds of mosaicism in haemophilia, including the mechanisms by which they arise and the risk of passing these variants on to the next generation. In addition, we focus on the selection of cell tissues and methods to detect these mosaic variants in the haemophilia setting. Taking into account the high rate of mosaicism in mothers of sporadic cases, we propose a diagnostic flow chart that could facilitate better evaluation of the risk of transmitting haemophilia in genetic and prenatal counselling.

Overlooked roles of DNA damage and maternal age in generating human germline mutations

The textbook view that most germline mutations in mammals arise from replication errors is indirectly supported by the fact that there are both more mutations and more cell divisions in the male than in the female germline. When analyzing large de novo mutation datasets in humans, we find multiple lines of evidence that call that view into question. Notably, despite the drastic increase in the ratio of male to female germ cell divisions after the onset of spermatogenesis, even young fathers contribute three times more mutations than young mothers, and this ratio barely increases with parental age. This surprising finding points to a substantial contribution of damage-induced mutations. Indeed, C-to-G transversions and CpG transitions, which together constitute over one-fourth of all base substitution mutations, show genomic distributions and sex-specific age dependencies indicative of double-strand break repair and methylation-associated damage, respectively. Moreover, we find evidence that maternal age at conception influences the mutation rate both because of the accumulation of damage in oocytes and potentially through an influence on the number of postzygotic mutations in the embryo. These findings reveal underappreciated roles of DNA damage and maternal age in the genesis of human germline mutations.

Statistical Methods for Identifying Sequence Motifs Affecting Point Mutations

Mutation processes differ between types of point mutation, genomic locations, cells, and biological species. For some point mutations, specific neighboring bases are known to be mechanistically influential. Beyond these cases, numerous questions remain unresolved, including: what are the sequence motifs that affect point mutations? How large are the motifs? Are they strand symmetric? And, do they vary between samples? We present new log-linear models that allow explicit examination of these questions, along with sequence logo style visualization to enable identifying specific motifs. We demonstrate the performance of these methods by analyzing mutation processes in human germline and malignant melanoma. We recapitulate the known CpG effect, and identify novel motifs, including a highly significant motif associated with A[Formula: see text]G mutations. We show that major effects of neighbors on germline mutation lie within [Formula: see text] of the mutating base. Models are also presented for contrasting the entire mutation spectra (the distribution of the different point mutations). We show the spectra vary significantly between autosomes and X-chromosome, with a difference in T[Formula: see text]C transition dominating. Analyses of malignant melanoma confirmed reported characteristic features of this cancer, including statistically significant strand asymmetry, and markedly different neighboring influences. The methods we present are made freely available as a Python library https://bitbucket.org/pycogent3/mutationmotif.

Breaking bad: The mutagenic effect of DNA repair

Species survival depends on the faithful replication of genetic information, which is continually monitored and maintained by DNA repair pathways that correct replication errors and the thousands of lesions that arise daily from the inherent chemical lability of DNA and the effects of genotoxic agents. Nonetheless, neutrally evolving DNA (not under purifying selection) accumulates base substitutions with time (the neutral mutation rate). Thus, repair processes are not 100% efficient. The neutral mutation rate varies both between and within chromosomes. For example it is 10-50 fold higher at CpGs than at non-CpG positions. Interestingly, the neutral mutation rate at non-CpG sites is positively correlated with CpG content. Although the basis of this correlation was not immediately apparent, some bioinformatic results were consistent with the induction of non-CpG mutations by DNA repair at flanking CpG sites. Recent studies with a model system showed that in vivo repair of preformed lesions (mismatches, abasic sites, single stranded nicks) can in fact induce mutations in flanking DNA. Mismatch repair (MMR) is an essential component for repair-induced mutations, which can occur as distant as 5 kb from the introduced lesions. Most, but not all, mutations involved the C of TpCpN (G of NpGpA) which is the target sequence of the C-preferring single-stranded DNA specific APOBEC deaminases. APOBEC-mediated mutations are not limited to our model system: Recent studies by others showed that some tumors harbor mutations with the same signature, as can intermediates in RNA-guided endonuclease-mediated genome editing. APOBEC deaminases participate in normal physiological functions such as generating mutations that inactivate viruses or endogenous retrotransposons, or that enhance immunoglobulin diversity in B cells. The recruitment of normally physiological error-prone processes during DNA repair would have important implications for disease, aging and evolution. This perspective briefly reviews both the bioinformatic and biochemical literature relevant to repair-induced mutagenesis and discusses future directions required to understand the mechanistic basis of this process.

Chromothripsis: potential origin in gametogenesis and preimplantation cell divisions. A review

OBJECTIVE

To review the discovery of chromothripsis and analyze its impact on human reproduction.

DESIGN

Database and literature analysis.

SETTING

University hospital.

PATIENT(S)

Carriers of massive and complex chromosomal rearrangements.

INTERVENTION(S)

Cytogenetic analysis and molecular testing (fluorescence in situ hybridization, microarray, whole-genome sequencing).

MAIN OUTCOME MEASURE(S)

Chromothripsis occurrence in human gametes and preimplantation embryos, with regard to the potential causative mechanisms described in literature.

RESULT(S)

Databases were searched for the literature published up to March 2014. Chromothripsis is characterized by the shattering of one (or a few) chromosome segments followed by a haphazard reassembly of the fragments generated, arising through a single initial catastrophic event. Several mechanisms involving abortive apoptosis, telomere erosion, mitotic errors, micronuclei formation, and p53 inactivation might cause chromothripsis. The remarkable point is that all these plausible mechanisms have been identified in the field of human reproduction as causal factors for reproductive failures and the genesis of chromosomal abnormalities. Specific features of gametogenesis and early embryonic development such as the weakness of cell cycle and mitosis checkpoints and the rapid kinetics of division in germ cells and early cleavage embryos may contribute to the emergence of chromothripsis.

CONCLUSION(S)

The discovery of this new class of massive chromosomal rearrangement has deeply modified our understanding on the genesis of complex genomic rearrangements. Data presented in this review support the assumption that chromothripsis could operate in human germlines and during early embryonic development. Chromothripsis might arise more frequently than previously thought in both gametogenesis and early human embryogenesis.

Repair of naturally occurring mismatches can induce mutations in flanking DNA

'Normal' genomic DNA contains hundreds of mismatches that are generated daily by the spontaneous deamination of C (U/G) and methyl-C (T/G). Thus, a mutagenic effect of their repair could constitute a serious genetic burden. We show here that while mismatches introduced into human cells on an SV40-based episome were invariably repaired, this process induced mutations in flanking DNA at a significantly higher rate than no mismatch controls. Most mutations involved the C of TpC, the substrate of some single strand-specific APOBEC cytidine deaminases, similar to the mutations that can typify the 'mutator phenotype' of numerous tumors. siRNA knockdowns and chromatin immunoprecipitation showed that TpC preferring APOBECs mediate the mutagenesis, and siRNA knockdowns showed that both the base excision and mismatch repair pathways are involved. That naturally occurring mispairs can be converted to mutators, represents an heretofore unsuspected source of genetic changes that could underlie disease, aging, and evolutionary change.DOI: http://dx.doi.org/10.7554/eLife.02001.001.

Mosaicism and haemophilia

Mosaicism may affect the haemophilia phenotype. Well-known instances include chromosomal mosaicism due to aneuploidy and pseudo-mosaicism due to varying patterns of X-chromosome inactivation. Chromosomal mosaicism in a chimera is a potential source of phenotypic variation. Gene mosaicism is commonplace. Its pattern and effect depend on the stage of development at which a mutation occurs. Proven or possible genetic mosaicism is an important consideration when predicting the likelihood of transmission of haemophilia to a future generation.

Evidence that replication-associated mutation alone does not explain between-chromosome differences in substitution rates

Since Haldane first noticed an excess of paternally derived mutations, it has been considered that most mutations derive from errors during germ line replication. Miyata et al. (1987) proposed that differences in the rate of neutral evolution on X, Y, and autosome can be employed to measure the extent of this male bias. This commonly applied method assumes replication to be the sole source of between-chromosome variation in substitution rates. We propose a simple test of this assumption: If true, estimates of the male bias should be independent of which two chromosomal classes are compared. Prior evidence from rodents suggested that this might not be true, but conclusions were limited by a lack of rat Y-linked sequence. We therefore sequenced two rat Y-linked bacterial artificial chromosomes and determined evolutionary rate by comparison with mouse. For estimation of rates we consider both introns and synonymous rates. Surprisingly, for both data sets the prediction of congruent estimates of alpha is strongly rejected. Indeed, some comparisons suggest a female bias with autosomes evolving faster than Y-linked sequence. We conclude that the method of Miyata et al. (1987) has the potential to provide incorrect estimates. Correcting the method requires understanding of the other causes of substitution that might differ between chromosomal classes. One possible cause is recombination-associated substitution bias for which we find some evidence. We note that if, as some suggest, this association is dominantly owing to male recombination, the high estimates of alpha seen in birds is to be expected as Z chromosomes recombine in males.

Mosaics and haemophilia

Some mosaic conditions may affect the haemophilia phenotype. Well-known instances include chromosomal mosaicism because of aneuploidy and pseudo-mosaicism because of varying patterns of X-chromosome inactivation. Chromosomal mosaicism in a chimera is a potential source of phenotypic variation. Gene mosaicism is commonplace. Its pattern and effect depend on the stage of development at which a mutation occurs. Proven or possible genetic mosaicism is an important consideration when predicting the likelihood of transmission of haemophilia to a future generation. A mosaic is an individual who has two or more cell lines, genetically different with regard to chromosomes or genes. As techniques improve and studies accumulate, mosaics are being found to be more common than hitherto believed. Some mosaic conditions may affect the phenotype of haemophilia in males and of the carrier state in females. Cells may be mosaic with regard to chromosomes, as in some instances of aneuploidy, and in chimeras, and in females owing to the pattern of X-chromosome inactivation. Cells may be mosaic with regard to new gene mutations. The pattern of mosaicism depends upon the stage in embryogenesis or in germ-cell formation in which the mutation arose.

Neo-sex chromosomes in the black muntjac recapitulate incipient evolution of mammalian sex chromosomes

BACKGROUND

The regular mammalian X and Y chromosomes diverged from each other at least 166 to 148 million years ago, leaving few traces of their early evolution, including degeneration of the Y chromosome and evolution of dosage compensation.

RESULTS

We studied the intriguing case of black muntjac, in which a recent X-autosome fusion and a subsequent large autosomal inversion within just the past 0.5 million years have led to inheritance patterns identical to the traditional X-Y (neo-sex chromosomes). We compared patterns of genome evolution in 35-kilobase noncoding regions and 23 gene pairs on the homologous neo-sex chromosomes. We found that neo-Y alleles have accumulated more mutations, comprising a wide variety of mutation types, which indicates cessation of recombination and is consistent with an ongoing neo-Y degeneration process. Putative deleterious mutations were observed in coding regions of eight investigated genes as well as cis-regulatory regions of two housekeeping genes. In vivo assays characterized a neo-Y insertion in the promoter of the CLTC gene that causes a significant reduction in allelic expression. A neo-Y-linked deletion in the 3'-untranslated region of gene SNX22 abolished a microRNA target site. Finally, expression analyses revealed complex patterns of expression divergence between neo-Y and neo-X alleles.

CONCLUSION

The nascent neo-sex chromosome system of black muntjacs is a valuable model in which to study the evolution of sex chromosomes in mammals. Our results illustrate the degeneration scenarios in various genomic regions. Of particular importance, we report--for the first time--that regulatory mutations were probably able to accelerate the degeneration process of Y and contribute to further evolution of dosage compensation.

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.

Rates of genome evolution and branching order from whole genome analysis

Accurate estimation of any phylogeny is important as a framework for evolutionary analysis of form and function at all levels of organization from sequence to whole organism. Using alignments of nonrepetitive components of opossum, human, mouse, rat, and dog genomes we evaluated two alternative tree topologies for eutherian evolution. We show with very high confidence that there is a basal split between rodents (as represented by the mouse and rat) and a branch joining primates (as represented by humans) and carnivores (as represented by dogs), consistent with some but not the most widely accepted mammalian phylogenies. The result was robust to substitution model choice with equivalent inference returned from a spectrum of models ranging from a general time reversible model, a model that treated nucleotides as either purines and pyrimidines, and variants of these that incorporated rate heterogeneity among sites. By determining this particular branching order we are able to show that the rate of molecular evolution is almost identical in rodent and carnivore lineages and that sequences evolve approximately 11%-14% faster in these lineages than in the primate lineage. In addition by applying the chicken as outgroup the analyses suggested that the rate of evolution in all eutherian lineages is approximately 30% slower than in the opossum lineage. This pattern of relative rates is inconsistent with the hypothesis that generation time is an important determinant of substitution rates and, by implication, mutation rates. Possible factors causing rate differences between the lineages include differences in DNA repair and replication enzymology, and shifts in nucleotide pools. Our analysis demonstrates the importance of using multiple sequences from across the genome to estimate phylogeny and relative evolutionary rate in order to reduce the influence of distorting local effects evident even in relatively long sequences.

Characteristics, causes and evolutionary consequences of male-biased mutation

Mutation has traditionally been considered a random process, but this paradigm is challenged by recent evidence of divergence rate heterogeneity in different genomic regions. One facet of mutation rate variation is the propensity for genetic change to correlate with the number of germ cell divisions, reflecting the replication-dependent origin of many mutations. Haldane was the first to connect this association of replication and mutation to the difference in the number of cell divisions in oogenesis (low) and spermatogenesis (usually high), and the resulting sex difference in the rate of mutation. The concept of male-biased mutation has been thoroughly analysed in recent years using an evolutionary approach, in which sequence divergence of autosomes and/or sex chromosomes are compared to allow inference about the relative contribution of mothers and fathers in the accumulation of mutations. For instance, assuming that a neutral sequence is analysed, that rate heterogeneity owing to other factors is cancelled out by the investigation of many loci and that the effect of ancestral polymorphism is properly taken into account, the male-to-female mutation rate ratio, alpham, can be solved from the observed difference in rate of X and Y chromosome divergence. The male mutation bias is positively correlated with the relative excess of cell divisions in the male compared to the female germ line, as evidenced by a generation time effect: in mammals, alpham is estimated at approximately 4-6 in primates, approximately 3 in carnivores and approximately 2 in small rodents. Another life-history correlate is sexual selection: when there is intense sperm competition among males, increased sperm production will be associated with a larger number of mitotic cell divisions in spermatogenesis and hence an increase in alpham. Male-biased mutation has implications for important aspects of evolutionary biology such as mate choice in relation to mutation load, sexual selection and the maintenance of genetic diversity despite strong directional selection, the tendency for a disproportionate large role of the X (Z) chromosome in post-zygotic isolation, and the evolution of sex.

DNA repair in mammalian embryos

Mammalian cells have developed complex mechanisms to identify DNA damage and activate the required response to maintain genome integrity. Those mechanisms include DNA damage detection, DNA repair, cell cycle arrest and apoptosis which operate together to protect the conceptus from DNA damage originating either in parental gametes or in the embryo's somatic cells. DNA repair in the newly fertilized preimplantation embryo is believed to rely entirely on the oocyte's machinery (mRNAs and proteins deposited and stored prior to ovulation). DNA repair genes have been shown to be expressed in the early stages of mammalian development. The survival of the embryo necessitates that the oocyte be sufficiently equipped with maternal stored products and that embryonic gene expression commences at the correct time. A Medline based literature search was performed using the keywords 'DNA repair' and 'embryo development' or 'gametogenesis' (publication dates between 1995 and 2006). Mammalian studies which investigated gene expression were selected. Further articles were acquired from the citations in the articles obtained from the preliminary Medline search. This paper reviews mammalian DNA repair from gametogenesis to preimplantation embryos to late gestational stages.

Evolution on the X chromosome: unusual patterns and processes

Although the X chromosome is usually similar to the autosomes in size and cytogenetic appearance, theoretical models predict that its hemizygosity in males may cause unusual patterns of evolution. The sequencing of several genomes has indeed revealed differences between the X chromosome and the autosomes in the rates of gene divergence, patterns of gene expression and rates of gene movement between chromosomes. A better understanding of these patterns should provide valuable information on the evolution of genes located on the X chromosome. It could also suggest solutions to more general problems in molecular evolution, such as detecting selection and estimating mutational effects on fitness.

Molecular clocks: four decades of evolution

During the past four decades, the molecular-clock hypothesis has provided an invaluable tool for building evolutionary timescales, and has served as a null model for testing evolutionary and mutation rates in different species. Molecular clocks have also influenced the development of theories of molecular evolution. As DNA-sequencing technologies have progressed, the use of molecular clocks has increased, with a profound effect on our understanding of the temporal diversification of species and genomes.

Developmental regulation and in vitro culture effects on expression of DNA repair and cell cycle checkpoint control genes in rhesus monkey oocytes and embryos

DNA repair is essential for maintaining genomic integrity, and may be required in the early embryo to correct damage inherited via the gametes, damage that arises during DNA replication, or damage that arises in response to exposure to genotoxic agents. The capacity of preimplantation stage mammalian embryos to repair damaged DNA has not been well characterized, particularly in primate embryos. In this study, we examined the expression of 48 mRNAs related to sensing different kinds of DNA damage, repairing that DNA damage, and controlling the cell cycle to provide an opportunity for DNA repair. The expression data reveal dynamic temporal changes, indicating a changing ability of the rhesus embryo to detect and repair different kinds of DNA damage. Low expression or overexpression of specific DNA repair genes may limit the ability of the embryo to respond to DNA damage at certain stages. Additionally, our data reveal that in vitro culture may lead to dysregulation of many such genes and a potentially impaired ability to repair DNA damage, thus affecting cellular viability and long-term embryo viability via effects on genome integrity. This effect of in vitro culture on nonhuman primate embryos may be relevant to assessing the potential advantages and disadvantages of prolonged in vitro culture of human embryos.

Effects of scrotal heating on sperm surface protein PH-20 expression in sheep

Sperm surface protein PH-20 expression was studied during spermatogenesis in pubertal and adult sheep, using molecular and histological methods. The effects of 24 hr of insulation raising scrotal temperatures to 39 degrees C on PH-20 expression in ejaculated sheep sperm were also determined. A 282 nt cDNA fragment of ovine PH-20 was identified in total RNA extracts of sheep testes, which exhibited 76% identity at the nucleotide level with the equivalent region of the human sequence. Ovine PH-20 mRNA and immunoreactivity were identified only in adult ram testis and not in peri-pubertal ram testis tubules lacking round spermatids, nor in adult sheep brain, pituitary, heart, spleen, lung, liver, kidney, epididymis, or ovary. Ovine PH-20 protein was distributed predominantly on the postacrosomal membrane and was also present on the anterior membrane of the sperm head in fresh, unheated sheep semen. Scrotal heating caused a significant, transient decrease in the percentage of PH-20 immunoreactive sperm, but did not change the pattern of PH-20 staining on the sperm head. The results strongly suggest that ovine PH-20 is postmeiotically expressed in haploid germ cells in sheep testis and is arrayed on the membrane of the mature ovine spermatozoon. Scrotal heating appears to have few effects on PH-20 expression and distribution on ejaculated sperm.

Modeling the impact of DNA methylation on the evolution of BRCA1 in mammals

The modified base 5-methylcytosine ((m)C) plays an important functional role in the biology of mammals as an epigenetic modification and appears to exert a striking impact on the molecular evolution of mammal genomes. The collective epigenetic functions of (m)C revolve around its effect on gene transcription, while the influence of this modified base on the evolution of mammal genomes derives from the greatly elevated spontaneous mutation rate of (m)C to T. In mammals, (m)C occurs at the dinucleotides CpG, CpA, and CpT. As a step toward a comprehensive statistical examination of the role of (m)C in mammal molecular evolution, we have developed novel Markov models of codon substitution that incorporate dinucleotide-level terms relevant to (m)C mutation. We apply these models to two data sets of aligned BRCA1 exon 11 sequences from bats and primates. In all cases, terms specific to mutations that affect the dinucleotides CpG, CpA, and CpT significantly improved model fit. For the CpG-specific terms, both transition and transversion substitution rates were elevated. These rates differed between the data sets. Bats exhibited a lower relative rate of substitutions at CpG-containing codons. Transition substitutions were significantly less than 1 at CpA-containing codons but greater than 1 at CpT-containing codons. The inclusion of interaction terms in the codon models to represent possible confounding with the effect of natural selection were supported for codons that contained CpG and CpT, but not CpA. From the results, we infer that mutation of (m)C is a probable factor that affects BRCA1 codons containing the dinucleotide CpG, a possible factor for CpA-containing codons, and an unlikely factor that affects CpT-containing codons. The confounding of estimated terms with the effect of natural selection indicate this confounding must be addressed for comparisons between different coding and noncoding regions.

The population genetics of multistage carcinogenesis

Despite the many successes of cancer research, we lack the framework necessary to predict the ratio of familial (inherited) to sporadic (non-inherited) cancers. An evolutionary model of multistage carcinogenesis provides this framework by demonstrating that the number of tumour suppressor loci (TSLs) preventing cancer in a given tissue is expected to depend upon the tissue's vulnerability to pre-reproductive somatic mutation. Since this vulnerability increases with tissue size, single gene control of human cancer may be restricted to retinoblastoma, a cancer of the tiny embryonic retina. The model is used to estimate the frequency of mutant alleles causing inherited cancers, based on the population genetics of the mutation-selection balance between new mutations arising and selection that eliminates them. For each specific cancer, this balance is determined by the effectiveness with which pre-reproductive cancer is suppressed in the non-mutant genotype characteristic of that population. Effectiveness depends on an interaction between the number of TSLs suppressing the cancer and factors determining the tissue-wide somatic mutation rate, such as tissue size and number of pre-reproductive cell divisions. The model predicts that the commonest pre-reproductive cancers will have the lowest proportion of familial cases, and that cancers associated with the most TSLs will have the highest post-reproductive incidence but no elevated pre-reproductive risk (a pattern seen in human epithelial cancers).

Male-biased transmission of deleterious mutations to the progeny in Arabidopsis thaliana

The extent and cause of male-biased mutation rates, the higher number of mutations in sperm than in eggs, is currently an active and controversial subject. Recent evidence indicates that this male (sperm) bias not only occurs in animals but also in plants. The higher mutation rate in plant sperm was inferred from rates of evolution of neutral DNA regions, and the results were confined to the mitochondria and chloroplasts of gymnosperms. However, the relative transmission rates of deleterious mutations, which have substantial evolutionary consequences, have rarely been studied. Here, an investigation is described by using the hermaphroditic self-compatible flowering plant Arabidopsis thaliana, in which we artificially increased the rate of mutation in pollen (i.e., sperm donor) and maternal (i.e., egg donor) parents, by using two kinds of UV irradiation in parallel and separate experiments, and assessed the deleterious effects on fitness of the F(2) generation. The results show that more deleterious induced mutations are transmitted to the progeny by a sperm than by an egg. These findings provide the first experimental evidence that more deleterious mutations are inherited from sperm than from an egg in any organism. Possible causes underlying this male bias are discussed.

Four new mitochondrial genomes and the increased stability of evolutionary trees of mammals from improved taxon sampling

We have sequenced four new mitochondrial genomes to improve the stability of the tree for placental mammals; they are two insectivores (a gymnure, Echinosorex gymnurus and Formosan shrew Soriculus fumidus); a Formosan lesser horseshoe bat (Rhinolophus monoceros); and the New Zealand fur seal (Arctocephalus forsteri). A revision to the hedgehog sequence (Erinaceus europaeus) is also reported. All five are from the Laurasiatheria grouping of eutherian mammals. On this new data set there is a strong tendency for the hedgehog and its relative, the gymnure, to join with the other Laurasiatherian insectivores (mole and shrews). To quantify the stability of trees from this data we define, based on nuclear sequences, a major four-way split in Laurasiatherians. This ([Xenarthra, Afrotheria], [Laurasiatheria, Supraprimates]) split is also found from mitochondrial genomes using either protein-coding or RNA (rRNA and tRNA) data sets. The high similarity of the mitochondrial and nuclear-derived trees allows a quantitative estimate of the stability of trees from independent data sets, as detected from a triplet Markov analysis. There are significant changes in the mutational processes within placental mammals that are ignored by current tree programs. On the basis of our quantitative results, we expect the evolutionary tree for mammals to be resolved quickly, and this will allow other problems to be solved.

The observed human sperm mutation frequency cannot explain the achondroplasia paternal age effect

The lifelong spermatogonial stem cell divisions unique to male germ cell production are thought to contribute to a higher mutation frequency in males. The fact that certain de novo human genetic conditions (e.g., achondroplasia) increase in incidence with the age of the father is consistent with this idea. Although it is assumed that the paternal age effect is the result of an increasing frequency of mutant sperm as a man grows older, no direct molecular measurement of the germ-line mutation frequency has been made to confirm this hypothesis. Using sperm DNA from donors of different ages, we determined the frequency of the nucleotide substitution in the fibroblast growth factor receptor 3 (FGFR3) gene that causes achondroplasia. Surprisingly, the magnitude of the increase in mutation frequency with age appears insufficient to explain why older fathers have a greater chance of having a child with this condition. A number of alternatives may explain this discrepancy, including selection for sperm that carry the mutation or an age-dependent increase in premutagenic lesions that remain unrepaired in sperm and are inefficiently detected by the PCR assay.

Male-driven evolution of mitochondrial and chloroplastidial DNA sequences in plants

Although there is substantial evidence that, in animals, male-inherited neutral DNA evolves at a higher rate than female-inherited DNA, the relative evolutionary rate of male- versus female-inherited DNA has not been investigated in plants. We compared the substitution rates at neutral sites of maternally and paternally inherited organellar DNA in gymnosperms. The analysis provided substantial support for the presence of a higher evolutionary rate in both the mitochondrial and chloroplastidial DNA when the organelle was inherited paternally than when inherited maternally. These results suggest that, compared with eggs, sperm tend to carry a greater number of mutations in mitochondrial and chloroplastidial DNA. The existence of a male mutation bias in plants is remarkable because, unlike animals, the germ-lines are not separated from the somatic cells throughout an individual's lifetime. The data therefore suggest that even a brief period of male and female germ-line separation can cause gender-specific mutation rates. These results are the first to show that, at least in some species, germ-lines influence the number of mutations carried in the gametes. Possible causes of male mutation bias in plants are discussed.

Mutation rates in mammalian genomes

Knowledge of the rate of point mutation is of fundamental importance, because mutations are a vital source of genetic novelty and a significant cause of human diseases. Currently, mutation rate is thought to vary many fold among genes within a genome and among lineages in mammals. We have conducted a computational analysis of 5,669 genes (17,208 sequences) from species representing major groups of placental mammals to characterize the extent of mutation rate differences among genes in a genome and among diverse mammalian lineages. We find that mutation rate is approximately constant per year and largely similar among genes. Similarity of mutation rates among lineages with vastly different generation lengths and physiological attributes points to a much greater contribution of replication-independent mutational processes to the overall mutation rate. Our results suggest that the average mammalian genome mutation rate is 2.2 x 10(-9) per base pair per year, which provides further opportunities for estimating species and population divergence times by using molecular clocks.

Mutations within a furin consensus sequence block proteolytic release of ectodysplasin-A and cause X-linked hypohidrotic ectodermal dysplasia

X-linked hypohidrotic ectodermal dysplasia (XLHED) is a heritable disorder of the ED-1 gene disrupting the morphogenesis of ectodermal structures. The ED-1 gene product, ectodysplasin-A (EDA), is a tumor necrosis factor (TNF) family member and is synthesized as a membrane-anchored precursor protein with the TNF core motif located in the C-terminal domain. The stalk region of EDA contains the sequence -Arg-Val-Arg-Arg156-Asn-Lys-Arg159-, representing overlapping consensus cleavage sites (Arg-X-Lys/Arg-Arg( downward arrow)) for the proprotein convertase furin. Missense mutations in four of the five basic residues within this sequence account for approximately 20% of all known XLHED cases, with mutations occurring most frequently at Arg156, which is shared by the two consensus furin sites. These analyses suggest that cleavage at the furin site(s) in the stalk region is required for the EDA-mediated cell-to-cell signaling that regulates the morphogenesis of ectodermal appendages. Here we show that the 50-kDa EDA parent molecule is cleaved at -Arg156Asn-Lys-Arg(159 downward arrow)- to release the soluble C-terminal fragment containing the TNF core domain. This cleavage appears to be catalyzed by furin, as release of the TNF domain was blocked either by expression of the furin inhibitor alpha1-PDX or by expression of EDA in furin-deficient LoVo cells. These results demonstrate that mutation of a functional furin cleavage site in a developmental signaling molecule is a basis for human disease (XLHED) and raise the possibility that furin cleavage may regulate the ability of EDA to act as a juxtacrine or paracrine factor.

Measurement of mutational flow implies both a high new-mutation rate for Huntington disease and substantial underascertainment of late-onset cases

We describe a new approach for analysis of the epidemiology of progressive genetic disorders that quantifies the rate of progression of the disease in the population by measuring the mutational flow. The framework is applied to Huntington disease (HD), a dominant neurological disorder caused by the expansion of a CAG-trinucleotide sequence to >35 repeats. The disease is 100% penetrant in individuals with > or = 42 repeats. Measurement of the flow from disease alleles provides a minimum estimate of the flow in the whole population and implies that the new mutation rate for HD in each generation is > or = 10% of currently known cases (95% confidence limits 6%-14%). Analysis of the pattern of flow demonstrates systematic underascertainment for repeat lengths <44. Ascertainment falls to <50% for individuals with 40 repeats and to <5% for individuals with 36-38 repeats. Clinicians should not assume that HD is rare outside known pedigrees or that most cases have onset at age <50 years.

A case for evolutionary genomics and the comprehensive examination of sequence biodiversity

Comparative analysis is one of the most powerful methods available for understanding the diverse and complex systems found in biology, but it is often limited by a lack of comprehensive taxonomic sampling. Despite the recent development of powerful genome technologies capable of producing sequence data in large quantities (witness the recently completed first draft of the human genome), there has been relatively little change in how evolutionary studies are conducted. The application of genomic methods to evolutionary biology is a challenge, in part because gene segments from different organisms are manipulated separately, requiring individual purification, cloning, and sequencing. We suggest that a feasible approach to collecting genome-scale data sets for evolutionary biology (i.e., evolutionary genomics) may consist of combination of DNA samples prior to cloning and sequencing, followed by computational reconstruction of the original sequences. This approach will allow the full benefit of automated protocols developed by genome projects to be realized; taxon sampling levels can easily increase to thousands for targeted genomes and genomic regions. Sequence diversity at this level will dramatically improve the quality and accuracy of phylogenetic inference, as well as the accuracy and resolution of comparative evolutionary studies. In particular, it will be possible to make accurate estimates of normal evolution in the context of constant structural and functional constraints (i.e., site-specific substitution probabilities), along with accurate estimates of changes in evolutionary patterns, including pairwise coevolution between sites, adaptive bursts, and changes in selective constraints. These estimates can then be used to understand and predict the effects of protein structure and function on sequence evolution and to predict unknown details of protein structure, function, and functional divergence. In order to demonstrate the practicality of these ideas and the potential benefit for functional genomic analysis, we describe a pilot project we are conducting to simultaneously sequence large numbers of vertebrate mitochondrial genomes.