Identification and Dissection of a Complex DNA Repair Sensitivity Phenotype in Baker's Yeast

Complex traits typically involve the contribution of multiple gene variants. In this study, we took advantage of a high-density genotyping analysis of the BY (S288c) and RM strains of Saccharomyces cerevisiae and of 123 derived spore progeny to identify the genetic loci that underlie a complex DNA repair sensitivity phenotype. This was accomplished by screening hybrid yeast progeny for sensitivity to a variety of DNA damaging agents. Both the BY and RM strains are resistant to the ultraviolet light–mimetic agent 4-nitroquinoline 1-oxide (4-NQO); however, hybrid progeny from a BY×RM cross displayed varying sensitivities to the drug. We mapped a major quantitative trait locus (QTL), RAD5, and identified the exact polymorphism within this locus responsible for 4-NQO sensitivity. By using a backcrossing strategy along with array-assisted bulk segregant analysis, we identified one other locus, MKT1, and a QTL on Chromosome VII that also link to the hybrid 4-NQO–sensitive phenotype but confer more minor effects. This work suggests an additive model for sensitivity to 4-NQO and provides a strategy for mapping both major and minor QTL that confer background-specific phenotypes. It also provides tools for understanding the effect of genetic background on sensitivity to genotoxic agents

Rapid Evolution of BRCA1 and BRCA2 in Humans and Other Primates

The maintenance of chromosomal integrity is an essential task of every living organism and cellular repair mechanisms exist to guard against insults to DNA. Given the importance of this process, it is expected that DNA repair proteins would be evolutionarily conserved, exhibiting very minimal sequence change over time. However, BRCA1, an essential gene involved in DNA repair, has been reported to be evolving rapidly despite the fact that many protein-altering mutations within this gene convey a significantly elevated risk for breast and ovarian cancers. Results: To obtain a deeper understanding of the evolutionary trajectory of BRCA1, we analyzed complete BRCA1 gene sequences from 23 primate species. We show that specific amino acid sites have experienced repeated selection for amino acid replacement over primate evolution. This selection has been focused specifically on humans and our closest living relatives, chimpanzees (Pan troglodytes) and bonobos (Pan paniscus). After examining BRCA1 polymorphisms in 7 bonobo, 44 chimpanzee, and 44 rhesus macaque (Macaca mulatta) individuals, we find considerable variation within each of these species and evidence for recent selection in chimpanzee populations. Finally, we also sequenced and analyzed BRCA2 from 24 primate species and find that this gene has also evolved under positive selection. Conclusions: While mutations leading to truncated forms of BRCA1 are clearly linked to cancer phenotypes in humans, there is also an underlying selective pressure in favor of amino acid-altering substitutions in this gene. A hypothesis where viruses are the drivers of this natural selection is discussed.National Institutes of Health R01-GM-093086, 8U42OD011197-13National Science Foundation BCS-07115972Burroughs Wellcome FundMolecular Bioscience

Widespread Sequence Variations in VAMP1 across Vertebrates Suggest a Potential Selective Pressure from Botulinum Neurotoxins

Botulinum neurotoxins (BoNT/A-G), the most potent toxins known, act by cleaving three SNARE proteins required for synaptic vesicle exocytosis. Previous studies on BoNTs have generally utilized the major SNARE homologues expressed in brain (VAMP2, syntaxin 1, and SNAP-25). However, BoNTs target peripheral motor neurons and cause death by paralyzing respiratory muscles such as the diaphragm. Here we report that VAMP1, but not VAMP2, is the SNARE homologue predominantly expressed in adult rodent diaphragm motor nerve terminals and in differentiated human motor neurons. In contrast to the highly conserved VAMP2, BoNT-resistant variations in VAMP1 are widespread across vertebrates. In particular, we identified a polymorphism at position 48 of VAMP1 in rats, which renders VAMP1 either resistant (I48) or sensitive (M48) to BoNT/D. Taking advantage of this finding, we showed that rat diaphragms with I48 in VAMP1 are insensitive to BoNT/D compared to rat diaphragms with M48 in VAMP1. This unique intra-species comparison establishes VAMP1 as a physiological toxin target in diaphragm motor nerve terminals, and demonstrates that the resistance of VAMP1 to BoNTs can underlie the insensitivity of a species to members of BoNTs. Consistently, human VAMP1 contains I48, which may explain why humans are insensitive to BoNT/D. Finally, we report that residue 48 of VAMP1 varies frequently between M and I across seventeen closely related primate species, suggesting a potential selective pressure from members of BoNTs for resistance in vertebrates

Incompatibilities Involving Yeast Mismatch Repair Genes: A Role for Genetic Modifiers and Implications for Disease Penetrance and Variation in Genomic Mutation Rates

Genetic background effects underlie the penetrance of most genetically determined phenotypes, including human diseases. To explore how such effects can modify a mutant phenotype in a genetically tractable system, we examined an incompatibility involving the MLH1 and PMS1 mismatch repair genes using a large population sample of geographically and ecologically diverse Saccharomyces cerevisiae strains. The mismatch repair incompatibility segregates into naturally occurring yeast strains, with no strain bearing the deleterious combination. In assays measuring the mutator phenotype conferred by different combinations of MLH1 and PMS1 from these strains, we observed a mutator phenotype only in combinations predicted to be incompatible. Surprisingly, intragenic modifiers could be mapped that specifically altered the strength of the incompatibility over a 20-fold range. Together, these observations provide a powerful model in which to understand the basis of disease penetrance and how such genetic variation, created through mating, could result in new mutations that could be the raw material of adaptive evolution in yeast populations

Evolutionary Reconstructions of the Transferrin Receptor of Caniforms Supports Canine Parvovirus Being a Re-emerged and Not a Novel Pathogen in Dogs

Parvoviruses exploit transferrin receptor type-1 (TfR) for cellular entry in carnivores, and specific interactions are key to control of host range. We show that several key mutations acquired by TfR during the evolution of Caniforms (dogs and related species) modified the interactions with parvovirus capsids by reducing the level of binding. These data, along with signatures of positive selection in the TFRC gene, are consistent with an evolutionary arms race between the TfR of the Caniform clade and parvoviruses. As well as the modifications of amino acid sequence which modify binding, we found that a glycosylation site mutation in the TfR of dogs which provided resistance to the carnivore parvoviruses which were in circulation prior to about 1975 predates the speciation of coyotes and dogs. Because the closely-related black-backed jackal has a TfR similar to their common ancestor and lacks the glycosylation site, reconstructing this mutation into the jackal TfR shows the potency of that site in blocking binding and infection and explains the resistance of dogs until recent times. This alters our understanding of this well-known example of viral emergence by indicating that canine parvovirus emergence likely resulted from the re-adaptation of a parvovirus to the resistant receptor of a former host

Ancient and Recent Adaptive Evolution of Primate Non-Homologous End Joining Genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes

Detection of Heterozygous Mutations in the Genome of Mismatch Repair Defective Diploid Yeast Using a Bayesian Approach

DNA replication errors that escape polymerase proofreading and mismatch repair (MMR) can lead to base substitution and frameshift mutations. Such mutations can disrupt gene function, reduce fitness, and promote diseases such as cancer and are also the raw material of molecular evolution. To analyze with limited bias genomic features associated with DNA polymerase errors, we performed a genome-wide analysis of mutations that accumulate in MMR-deficient diploid lines of Saccharomyces cerevisiae. These lines were derived from a common ancestor and were grown for 160 generations, with bottlenecks reducing the population to one cell every 20 generations. We sequenced to between 8- and 20-fold coverage one wild-type and three mutator lines using Illumina Solexa 36-bp reads. Using an experimentally aware Bayesian genotype caller developed to pool experimental data across sequencing runs for all strains, we detected 28 heterozygous single-nucleotide polymorphisms (SNPs) and 48 single-nt insertion/deletions (indels) from the data set. This method was evaluated on simulated data sets and found to have a very low false-positive rate (∼6 × 10−5) and a false-negative rate of 0.08 within the unique mapping regions of the genome that contained at least sevenfold coverage. The heterozygous mutations identified by the Bayesian genotype caller were confirmed by Sanger sequencing. All of the mutations were unique to a given line, except for a single-nt deletion mutation which occurred independently in two lines. All 48 indels, composed of 46 deletions and two insertions, occurred in homopolymer (HP) tracts [i.e., 47 poly(A) or (T) tracts, 1 poly(G) or (C) tract] between 5 and 13 bp long. Our findings are of interest because HP tracts are present at high levels in the yeast genome (>77,400 for 5- to 20-nt HP tracts), and frameshift mutations in these regions are likely to disrupt gene function. In addition, they demonstrate that the mutation pattern seen previously in mismatch repair defective strains using a limited number of reporters holds true for the entire genome

A putative SUMO interacting motif in the B30.2/SPRY domain of rhesus macaque TRIM5α important for NF-κB/AP-1 signaling and HIV-1 restriction

TRIM5α from the rhesus macaque (TRIM5αRh) is a restriction factor that shows strong activity against HIV-1. TRIM5αRh binds specifically to HIV-1 capsid (CA) through its B30.2/PRYSPRY domain shortly after entry of the virus into the cytoplasm. Recently, three putative SUMO interacting motifs (SIMs) have been identified in the PRYSPRY domain of human and macaque TRIM5α. However, structural modeling of this domain suggested that two of them were buried in the hydrophobic core of the protein, implying that interaction with SUMO was implausible, while the third one was not relevant to restriction. In light of these results, we re-analyzed the TRIM5αRh PRYSPRY sequence and identified an additional putative SIM (435VIIC438) which we named SIM4. This motif is exposed at the surface of the PRYSPRY domain, allowing potential interactions with SUMO or SUMOylated proteins. Introducing a double mutation in SIM4 (V435K, I436K) did not alter stability, unlike mutations in SIM1. SIM4-mutated TRIM5αRh failed to bind HIV-1CA and lost the ability to restrict this virus. Accordingly, SIM4 undergoes significant variation among primates and substituting this motif with naturally occurring SIM4 variants affected HIV-1 restriction by TRIM5αRh, suggesting a direct role in capsid recognition. Interestingly, SIM4-mutated TRIM5αRh also failed to activate NF-κB and AP-1-mediated transcription. Although there is no direct evidence that SIM4 is involved in direct interaction with SUMO or a SUMOylated protein, mutating this motif strongly reduced co-localization of TRIM5αRh with SUMO-1 and with PML, a SUMOylated nuclear protein. In conclusion, this new putative SIM is crucial for both direct interaction with incoming capsids and for NF-κB/AP-1 signaling. We speculate that the latter function is mediated by interactions of SIM4 with a SUMOylated protein involved in the NF-κB/AP-1 signaling pathways

Host species examined in this analysis.

<p><i>TFRC</i> sequences were determined from cDNA prepared from mRNA isolated from the samples indicated.</p