Recombination and the Evolution of Diploidy

With two copies of every gene, a diploid organism is able to mask recessive deleterious mutations. In this paper we present the analysis of a two-locus model designed to determine when the masking of deleterious alleles favors the evolution of a dominant diploid phase in organisms that alternate between haploid and diploid phases ("alternation of generations"). It is hypothesized that diploidy will be favored whenever masking occurs ("the masking hypothesis"). Using analytical methods, we confirm that this masking hypothesis is essentially correct under free recombination: as long as the heterozygous expression of deleterious alleles is sufficiently masked by the wild-type allele, diploidy is favored over haploidy. When the rate of recombination is lower, however, diploidy is much less likely to be favored over haploidy. In fact, according to our model, the evolution of diploidy is impossible without significant levels of recombination even when masking is fairly strong.     

Confidence-based Somatic Mutation Evaluation and Prioritization

Next generation sequencing (NGS) has enabled high throughput discovery of somatic mutations. Detection depends on experimental design, lab platforms, parameters and analysis algorithms. However, NGS-based somatic mutation detection is prone to erroneous calls, with reported validation rates near 54% and congruence between algorithms less than 50%. Here, we developed an algorithm to assign a single statistic, a false discovery rate (FDR), to each somatic mutation identified by NGS. This FDR confidence value accurately discriminates true mutations from erroneous calls. Using sequencing data generated from triplicate exome profiling of C57BL/6 mice and B16-F10 melanoma cells, we used the existing algorithms GATK, SAMtools and SomaticSNiPer to identify somatic mutations. For each identified mutation, our algorithm assigned an FDR. We selected 139 mutations for validation, including 50 somatic mutations assigned a low FDR (high confidence) and 44 mutations assigned a high FDR (low confidence). All of the high confidence somatic mutations validated (50 of 50), none of the 44 low confidence somatic mutations validated, and 15 of 45 mutations with an intermediate FDR validated. Furthermore, the assignment of a single FDR to individual mutations enables statistical comparisons of lab and computation methodologies, including ROC curves and AUC metrics. Using the HiSeq 2000, single end 50 nt reads from replicates generate the highest confidence somatic mutation call set.     

Population Studies in Microorganisms I. Evolution of Diploidy in SACCHAROMYCES CEREVISIAE

The relative adaptation of isogenic haploid and diploid strains of yeast was investigated in different sets of physiological conditions. When all nutrients were present in excess, no difference in the reproductive rates of isogenic haploid and diploid strains of yeast was detected in both optimal and non-optimal growth conditions. Competition between haploid and diploid strains of yeast was observed when growth was limited by the concentration of a single nutrilite. Under certain conditions when fitness (reproductive rate) is determined by transport of an essential nutrilite that exists in very low concentrations, diploid cells were selected against. These environmental conditions are similar to those found in offshore marine environments where nutrients are often present in extremely low concentrations. The fitness of diploid cells was estimated to be.93 +/-.02 (haploid fitness = 1). The reduced fitness of diploid cells in this environment can be explained by the reduced surface area/volume ratio possessed by diploid cells in comparison to haploid cells. The fitnesses of haploid and diploid cells in these environments are closely correlated with geometric variations in these strains. These results are consistent with the hypothesis that diploid cells are simply double haploids, and diploidy per se does not confer any direct adaptive advantage. The mechanism of the evolution of diploidy as a dominant phase in the life cycle of higher plants and animals remains obscure.     

Heightened Exposure to Parasites Favors the Evolution of Immunity in Brood Parasitic Cowbirds

Immunologists and evolutionary biologists are interested in how the immune system evolves to fit an ecological niche. We studied the relationship between exposure to parasites and strength of immunity by investigating the response of two species of New World cowbirds (genus Molothrus, Icteridae), obligate brood parasites with contrasting life history strategies, to experimental arboviral infection. The South American shiny cowbird (M. bonariensis) is an extreme host-generalist that lays its eggs in the nests of >225 different avian species. The Central American bronzed cowbird (M. aeneus) is a relative host-specialist that lays its eggs preferentially in the nests of approximately 12 orioles in a single sister genus. West Nile virus provided a strong challenge and delineated immune differences between these species. The extreme host-generalist shiny cowbird, like the North American host-generalist, the brown-headed cowbird, showed significantly lower viremia to three arboviruses than related icterid species that were not brood parasites. The bronzed cowbird showed intermediate viremia. These findings support the interpretation that repeated exposure to a high diversity of parasites favors the evolution of enhanced immunity in brood parasitic cowbirds and makes them useful models for future studies of innate immunity.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.     

Somatic Point Mutation Calling in Low Cellularity Tumors

Somatic mutation calling from next-generation sequencing data remains a challenge due to the difficulties of distinguishing true somatic events from artifacts arising from PCR, sequencing errors or mis-mapping. Tumor cellularity or purity, sub-clonality and copy number changes also confound the identification of true somatic events against a background of germline variants. We have developed a heuristic strategy and software (http://www.qcmg.org/bioinformatics/qsnp/) for somatic mutation calling in samples with low tumor content and we show the superior sensitivity and precision of our approach using a previously sequenced cell line, a series of tumor/normal admixtures, and 3,253 putative somatic SNVs verified on an orthogonal platform.     

Parental Age Affects Somatic Mutation Rates in the Progeny of Flowering Plants

In humans, it is well known that the parental reproductive age has a strong influence on mutations transmitted to their progeny. Meiotic nondisjunction is known to increase in older mothers, and base substitutions tend to go up with paternal reproductive age. Hence, it is clear that the germinal mutation rates are a function of both maternal and paternal ages in humans. In contrast, it is unknown whether the parental reproductive age has an effect on somatic mutation rates in the progeny, because these are rare and difficult to detect. To address this question, we took advantage of the plant model system Arabidopsis (Arabidopsis thaliana), where mutation detector lines allow for an easy quantitation of somatic mutations, to test the effect of parental age on somatic mutation rates in the progeny. Although we found no significant effect of parental age on base substitutions, we found that frameshift mutations and transposition events increased in the progeny of older parents, an effect that is stronger through the maternal line. In contrast, intrachromosomal recombination events in the progeny decrease with the age of the parents in a parent-of-origin-dependent manner. Our results clearly show that parental reproductive age affects somatic mutation rates in the progeny and, thus, that some form of age-dependent information, which affects the frequency of double-strand breaks and possibly other processes involved in maintaining genome integrity, is transmitted through the gametes.In humans, it has long been recognized that the reproductive age of the parents has an influence on the health of their progeny. An older reproductive age of the mother is known to increase the fraction of aneuploid gamete formation (Hurles, 2012). For instance, the risk for a trisomy increases from 2% to 3% for mothers in their 20s to more than 30% for mothers in their 40s (Hassold and Hunt, 2009). The age of the father also has an effect on the frequency of spontaneous congenital disorders and common complex diseases, such as autism and some cancers (Goriely and Wilkie, 2012). Indeed, sperm from 36- to 57-year-old men have more double-strand breaks (DSBs) than those of 20- to 35-year-old individuals (Singh et al., 2003). Similarly, the efficiency of DSB repair was reported to decrease with age in vegetative tissues of the plant model system Arabidopsis (Arabidopsis thaliana; Boyko et al., 2006).Owing to the continuous divisions of spermatogonial stem cells, the male germline of humans is thought to be more mutagenic than the female germline. Indeed, it was shown that the paternal germline is more mutagenic than the maternal one with respect to base substitutions (Kong et al., 2012) and replication slippage errors at microsatellites (Sun et al., 2012). It is also known that carriers of germline mutations in mismatch repair (MMR) genes in humans are prone to get colorectal cancer and that the risk depends on the parent-of-origin of the mutation (van Vliet et al., 2011). The molecular basis of these parental effects is not entirely clear but is likely to involve higher rates of nondisjunction during female meiosis, higher mutation rates during spermatogenesis, and probably additional effects of aging.In contrast to the effect of parental age on germline mutations, not much is known about potential effects of parental reproductive age on somatic mutation rates in the offspring. However, it has been shown in animal studies that radiation of males can lead to somatic mutations in their progeny—and subsequent generations—that cannot be attributed to mutations in the paternal germline (for review, see Little et al., 2013). Moreover, several recent studies have illustrated the existence of complex parental and transgenerational effects in humans, although their molecular basis is not clear (Grossniklaus et al., 2013). These effects can be of either genetic nature (but the effect is seen even in offspring that did not inherit the genetic variant from their parents; for review, see Nadeau, 2009) or epigenetic nature (where environmental influences can possibly exert effects on subsequent generations; for review, see Pembrey et al., 2006; Pembrey, 2010; Curley et al., 2011). It is currently not known whether such parental effects affect the somatic mutation rates in the offspring or whether the effects are modulated by parental age.Taking advantage of the plant model system Arabidopsis, in which various somatic mutation rates can readily be assessed (Bashir et al., 2014), we investigated the effects of parental reproductive age on somatic mutation rates in the progeny. We report that there is a pronounced effect of parental age on somatic mutation rates in their offspring in a parent-of-origin-dependent fashion. Thus, some form of parental information, which is inherited through the gametes to the next generation, seems to alter the somatic mutation rates in the progeny and changes with parental reproductive age.     

Human Microsatellites: Mutation and Evolution

Microsatellites (MSs) are short tandem DNA repeats with the repetitive motif of two to six nucleotides, forming tracts up to hundreds of nucleotides long. Notwithstanding the active use of MSs in genetic studies of various biological problems, the reasons for their wide occurrence in the genome, their possible functions, and mutational behavior are still unclear. The mutation rate in MS repeats is on average several orders of magnitude higher than in the remaining DNA, which allows for direct estimation of evolutionary transformation rate in nucleotide sequences of the genome. Mutation process in MSs is highly heterogeneous, with distinct differences between species; furthermore, within a species it differs among loci with different repeat size, among alleles of one locus, and among individuals of different sex and age. Most MS mutations are caused by DNA slippage during replication but the probability of this event depends on the locus. In this review, a number of models of MS evolution are discussed, which account for the relationship between mutation rate and allele size, different mutation direction in alleles of different size, and the appearance of point mutations within repeat tracts restricting allele size. The MS evolution is considered mainly in the context of selective neutrality, although there is evidence showing functional significance of some variants of tandem repeats and thus their possible selective value.     

Mutation Rate Evolution in Replicator Dynamics

The mutation rate of an organism is itself evolvable. In stable environments, if faithful replication is costless, theory predicts that mutation rates will evolve to zero. However, positive mutation rates can evolve in novel or fluctuating environments, as analytical and empirical studies have shown. Previous work on this question has focused on environments that fluctuate independently of the evolving population. Here we consider fluctuations that arise from frequency-dependent selection in the evolving population itself. We investigate how the dynamics of competing traits can induce selective pressure on the rates of mutation between these traits. To address this question, we introduce a theoretical framework combining replicator dynamics and adaptive dynamics. We suppose that changes in mutation rates are rare, compared to changes in the traits under direct selection, so that the expected evolutionary trajectories of mutation rates can be obtained from analysis of pairwise competition between strains of different rates. Depending on the nature of frequency-dependent trait dynamics, we demonstrate three possible outcomes of this competition. First, if trait frequencies are at a mutation–selection equilibrium, lower mutation rates can displace higher ones. Second, if trait dynamics converge to a heteroclinic cycle—arising, for example, from “rock-paper-scissors” interactions—mutator strains succeed against non-mutators. Third, in cases where selection alone maintains all traits at positive frequencies, zero and nonzero mutation rates can coexist indefinitely. Our second result suggests that relatively high mutation rates may be observed for traits subject to cyclical frequency-dependent dynamics.     

RADIA: RNA and DNA Integrated Analysis for Somatic Mutation Detection

The detection of somatic single nucleotide variants is a crucial component to the characterization of the cancer genome. Mutation calling algorithms thus far have focused on comparing the normal and tumor genomes from the same individual. In recent years, it has become routine for projects like The Cancer Genome Atlas (TCGA) to also sequence the tumor RNA. Here we present RADIA (RNA and DNA Integrated Analysis), a novel computational method combining the patient-matched normal and tumor DNA with the tumor RNA to detect somatic mutations. The inclusion of the RNA increases the power to detect somatic mutations, especially at low DNA allelic frequencies. By integrating an individual’s DNA and RNA, we are able to detect mutations that would otherwise be missed by traditional algorithms that examine only the DNA. We demonstrate high sensitivity (84%) and very high precision (98% and 99%) for RADIA in patient data from endometrial carcinoma and lung adenocarcinoma from TCGA. Mutations with both high DNA and RNA read support have the highest validation rate of over 99%. We also introduce a simulation package that spikes in artificial mutations to patient data, rather than simulating sequencing data from a reference genome. We evaluate sensitivity on the simulation data and demonstrate our ability to rescue back mutations at low DNA allelic frequencies by including the RNA. Finally, we highlight mutations in important cancer genes that were rescued due to the incorporation of the RNA.     

Imperfect Genes, Fisherian Mutation and the Evolution of Sex

In this paper we present a mathematical model of mutation and selection that allows for the coexistence of multiple alleles at a locus with very small selective differences between alleles. The model also allows for the determination of fitness by multiple loci. Models of this sort are biologically plausible. However, some previous attempts to construct similar models have assumed that all mutations produce a decrease in fitness, and this has led to a tendency for the average fitness of population members to decline when population numbers are finite. In our model we incorporate some of the ideas of R. A. FISHER, so that both deleterious and beneficial mutations are possible. As a result, average fitness tends to approach a stationary distribution. We have used computer simulation methods to apply the Fisherian mutation model to the problem of the evolution of sex and recombination. The results suggest that sex and recombination can provide very large benefits in terms of average fitness. The results also suggest that obligately sexual species will win ecological competitions with species that produce a substantial fraction of their offspring asexually, so long as the number of sites under selection within the genomes of the competing species is not too small and the population sizes are not too large. Our model focuses on fertility selection in an hermaphroditic plant. However, the results are likely to generalize to a wide variety of other situations as well.     

Fluctuating Mutation Bias and the Evolution of Base Composition in Drosophila

The idea that the pattern of point mutation in Drosophila has remained constant during the evolution of the genus has recently been challenged. A study of the nucleotide composition focused on the Drosophila saltans group has evidenced unsuspected nucleotide composition differences among lineages. Compositional differences are associated with an accelerated rate of amino acid replacement in functionally less constrained regions. Here we reassess this issue from a different perspective. Adopting a maximum-likelihood estimation approach, we focus on the different predictions that mutation and selection make about the nonsynonymous-to-synonymous rate ratio. We investigate two gene regions, alcohol dehydrogenase (Adh) and xanthine dehydrogenase (Xdh), using a balanced data set that comprises representatives from the melangaster, obscura, saltans, and willistoni groups. We also consider representatives of the Hawaiian picture-winged group. These Hawaiian species are known to have experienced repeated bottlenecks and are included as a reference for comparison. Our results confirm patterns previously detected. The branch ancestral to the fast-evolving willistoni/saltans lineage, where most of the change in GC content has occurred, exhibits an excess of synonymous substitutions. The shift in mutation bias has affected the extent of the rate variation among sites in Xdh. Received: 4 May 1999 / Accepted: 26 July 1999