The Shu complex prevents mutagenesis and cytotoxicity of single-strand specific alkylation lesions.

Three-methyl cytosine (3meC) are toxic DNA lesions, blocking base pairing. Bacteria and humans express members of the AlkB enzymes family, which directly remove 3meC. However, other organisms, including budding yeast, lack this class of enzymes. It remains an unanswered evolutionary question as to how yeast repairs 3meC, particularly in single-stranded DNA. The yeast Shu complex, a conserved homologous recombination factor, aids in preventing replication-associated mutagenesis from DNA base damaging agents such as methyl methanesulfonate (MMS). We found that MMS-treated Shu complex-deficient cells exhibit a genome-wide increase in A:T and G:C substitutions mutations. The G:C substitutions displayed transcriptional and replicational asymmetries consistent with mutations resulting from 3meC. Ectopic expression of a human AlkB homolog in Shu-deficient yeast rescues MMS-induced growth defects and increased mutagenesis. Thus, our work identifies a novel homologous recombination-based mechanism mediated by the Shu complex for coping with alkylation adducts.

Atypical UV Photoproducts Induce Non-canonical Mutation Classes Associated with Driver Mutations in Melanoma.

Somatic mutations in skin cancers and other ultraviolet (UV)-exposed cells are typified by C>T and CC>TT substitutions at dipyrimidine sequences; however, many oncogenic "driver" mutations in melanoma do not fit this UV signature. Here, we use genome sequencing to characterize mutations in yeast repeatedly irradiated with UV light. Analysis of ~50,000 UV-induced mutations reveals abundant non-canonical mutations, including T>C, T>A, and AC>TT substitutions. These mutations display transcriptional asymmetry that is modulated by nucleotide excision repair (NER), indicating that they are caused by UV photoproducts. Using a sequencing method called UV DNA endonuclease sequencing (UVDE-seq), we confirm the existence of an atypical thymine-adenine photoproduct likely responsible for UV-induced T>A substitutions. Similar non-canonical mutations are present in skin cancers, which also display transcriptional asymmetry and dependence on NER. These include multiple driver mutations, most prominently the recurrent BRAF V600E and V600K substitutions, suggesting that mutations arising from rare, atypical UV photoproducts may play a role in melanomagenesis.

Characterization of systemic genomic instability in budding yeast.

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Mutation signatures specific to DNA alkylating agents in yeast and cancers.

Alkylation is one of the most ubiquitous forms of DNA lesions. However, the motif preferences and substrates for the activity of the major types of alkylating agents defined by their nucleophilic substitution reactions (SN1 and SN2) are still unclear. Utilizing yeast strains engineered for large-scale production of single-stranded DNA (ssDNA), we probed the substrate specificity, mutation spectra and signatures associated with DNA alkylating agents. We determined that SN1-type agents preferably mutagenize double-stranded DNA (dsDNA), and the mutation signature characteristic of the activity of SN1-type agents was conserved across yeast, mice and human cancers. Conversely, SN2-type agents preferably mutagenize ssDNA in yeast. Moreover, the spectra and signatures derived from yeast were detectable in lung cancers, head and neck cancers and tumors from patients exposed to SN2-type alkylating chemicals. The estimates of mutation loads associated with the SN2-type alkylation signature were higher in lung tumors from smokers than never-smokers, pointing toward the mutagenic activity of the SN2-type alkylating carcinogens in cigarettes. In summary, our analysis of mutations in yeast strains treated with alkylating agents, as well as in whole-exome and whole-genome-sequenced tumors identified signatures highly specific to alkylation mutagenesis and indicate the pervasive nature of alkylation-induced mutagenesis in cancers.

Correction: Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

[This corrects the article DOI: 10.1371/journal.pbio.3000464.].

Long transposon-rich centromeres in an oomycete reveal divergence of centromere features in Stramenopila-Alveolata-Rhizaria lineages.

Centromeres are chromosomal regions that serve as platforms for kinetochore assembly and spindle attachments, ensuring accurate chromosome segregation during cell division. Despite functional conservation, centromere DNA sequences are diverse and often repetitive, making them challenging to assemble and identify. Here, we describe centromeres in an oomycete Phytophthora sojae by combining long-read sequencing-based genome assembly and chromatin immunoprecipitation for the centromeric histone CENP-A followed by high-throughput sequencing (ChIP-seq). P. sojae centromeres cluster at a single focus at different life stages and during nuclear division. We report an improved genome assembly of the P. sojae reference strain, which enabled identification of 15 enriched CENP-A binding regions as putative centromeres. By focusing on a subset of these regions, we demonstrate that centromeres in P. sojae are regional, spanning 211 to 356 kb. Most of these regions are transposon-rich, poorly transcribed, and lack the histone modification H3K4me2 but are embedded within regions with the heterochromatin marks H3K9me3 and H3K27me3. Strikingly, we discovered a Copia-like transposon (CoLT) that is highly enriched in the CENP-A chromatin. Similar clustered elements are also found in oomycete relatives of P. sojae, and may be applied as a criterion for prediction of oomycete centromeres. This work reveals a divergence of centromere features in oomycetes as compared to other organisms in the Stramenopila-Alveolata-Rhizaria (SAR) supergroup including diatoms and Plasmodium falciparum that have relatively short and simple regional centromeres. Identification of P. sojae centromeres in turn also advances the genome assembly.

Congenital Midline Cervical Cleft: First Report and Genetic Analysis of Two Related Patients.

OBJECTIVES: Congenital midline cervical cleft (CMCC) is a rare congenital anterior neck anatomical anomaly. We present the case of two related patients (grandchild and maternal grandmother) who were both born with a congenital midline cervical cleft along with genetic analysis. METHODS: Clinical examination of both patients and surgical excision of the grandchild was performed. Genetic analysis with exome sequencing (ES) was conducted for both patients. RESULTS: Genetic analysis with exome sequencing (ES) revealed apparently novel single nucleotide variants in 66 genes present in both proband and grandmother. Five of these variants are predicted to cause frameshifting in the coding region of the respective genes and truncated proteins (OVGP1, TYW1B, ZAN, SSPO, FOLR3). Two of these genes (TYW1B and SSPO) have homozygous indel mutations in both patients. CONCLUSIONS: To our knowledge, this is the first case of two related patients with a congenital midline cervical cleft. The results of our genetic analysis reveal potential relevance to CMCC development.

Repair of multiple simultaneous double-strand breaks causes bursts of genome-wide clustered hypermutation.

A single cancer genome can harbor thousands of clustered mutations. Mutation signature analyses have revealed that the origin of clusters are lesions in long tracts of single-stranded (ss) DNA damaged by apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) cytidine deaminases, raising questions about molecular mechanisms that generate long ssDNA vulnerable to hypermutation. Here, we show that ssDNA intermediates formed during the repair of gamma-induced bursts of double-strand breaks (DSBs) in the presence of APOBEC3A in yeast lead to multiple APOBEC-induced clusters similar to cancer. We identified three independent pathways enabling cluster formation associated with repairing bursts of DSBs: 5' to 3' bidirectional resection, unidirectional resection, and break-induced replication (BIR). Analysis of millions of mutations in APOBEC-hypermutated cancer genomes revealed that cancer tolerance to formation of hypermutable ssDNA is similar to yeast and that the predominant pattern of clustered mutagenesis is the same as in resection-defective yeast, suggesting that cluster formation in cancers is driven by a BIR-like mechanism. The phenomenon of genome-wide burst of clustered mutagenesis revealed by our study can play an important role in generating somatic hypermutation in cancers as well as in noncancerous cells.

A Case Study of Genomic Instability in an Industrial Strain of Saccharomyces cerevisiae.

The Saccharomyces cerevisiae strain JAY270/PE2 is a highly efficient biocatalyst used in the production of bioethanol from sugarcane feedstock. This strain is heterothallic and diploid, and its genome is characterized by abundant structural and nucleotide polymorphisms between homologous chromosomes. One of the reasons it is favored by many distilleries is that its cells do not normally aggregate, a trait that facilitates cell recycling during batch-fed fermentations. However, long-term propagation makes the yeast population vulnerable to the effects of genomic instability, which may trigger the appearance of undesirable phenotypes such as cellular aggregation. In pure cultures of JAY270, we identified the recurrent appearance of mutants displaying a mother-daughter cell separation defect resulting in rough colonies in agar media and fast sedimentation in liquid culture. We investigated the genetic basis of the colony morphology phenotype and found that JAY270 is heterozygous for a frameshift mutation in the ACE2 gene (ACE2/ace2-A7), which encodes a transcriptional regulator of mother-daughter cell separation. All spontaneous rough colony JAY270-derived isolates analyzed carried copy-neutral loss-of-heterozygosity (LOH) at the region of chromosome XII where ACE2 is located (ace2-A7/ace2-A7). We specifically measured LOH rates at the ACE2 locus, and at three additional chromosomal regions in JAY270 and in a conventional homozygous diploid laboratory strain. This direct comparison showed that LOH rates at all sites were quite similar between the two strain backgrounds. In this case study of genomic instability in an industrial strain, we showed that the JAY270 genome is dynamic and that structural changes to its chromosomes can lead to new phenotypes. However, our analysis also indicated that the inherent level of genomic instability in this industrial strain is normal relative to a laboratory strain. Our work provides an important frame of reference to contextualize the interpretation of instability processes observed in the complex genomes of industrial yeast strains.

Genome-wide maps of alkylation damage, repair, and mutagenesis in yeast reveal mechanisms of mutational heterogeneity.

DNA base damage is an important contributor to genome instability, but how the formation and repair of these lesions is affected by the genomic landscape and contributes to mutagenesis is unknown. Here, we describe genome-wide maps of DNA base damage, repair, and mutagenesis at single nucleotide resolution in yeast treated with the alkylating agent methyl methanesulfonate (MMS). Analysis of these maps revealed that base excision repair (BER) of alkylation damage is significantly modulated by chromatin, with faster repair in nucleosome-depleted regions, and slower repair and higher mutation density within strongly positioned nucleosomes. Both the translational and rotational settings of lesions within nucleosomes significantly influence BER efficiency; moreover, this effect is asymmetric relative to the nucleosome dyad axis and is regulated by histone modifications. Our data also indicate that MMS-induced mutations at adenine nucleotides are significantly enriched on the nontranscribed strand (NTS) of yeast genes, particularly in BER-deficient strains, due to higher damage formation on the NTS and transcription-coupled repair of the transcribed strand (TS). These findings reveal the influence of chromatin on repair and mutagenesis of base lesions on a genome-wide scale and suggest a novel mechanism for transcription-associated mutation asymmetry, which is frequently observed in human cancers.

APOBEC3B cytidine deaminase targets the non-transcribed strand of tRNA genes in yeast.

Variations in mutation rates across the genome have been demonstrated both in model organisms and in cancers. This phenomenon is largely driven by the damage specificity of diverse mutagens and the differences in DNA repair efficiency in given genomic contexts. Here, we demonstrate that the single-strand DNA-specific cytidine deaminase APOBEC3B (A3B) damages tRNA genes at a 1000-fold higher efficiency than other non-tRNA genomic regions in budding yeast. We found that A3B-induced lesions in tRNA genes were predominantly located on the non-transcribed strand, while no transcriptional strand bias was observed in protein coding genes. Furthermore, tRNA gene mutations were exacerbated in cells where RNaseH expression was completely abolished (Δrnh1Δrnh35). These data suggest a transcription-dependent mechanism for A3B-induced tRNA gene hypermutation. Interestingly, in strains proficient in DNA repair, only 1% of the abasic sites formed upon excision of A3B-deaminated cytosines were not repaired leading to mutations in tRNA genes, while 18% of these lesions failed to be repaired in the remainder of the genome. A3B-induced mutagenesis in tRNA genes was found to be efficiently suppressed by the redundant activities of both base excision repair (BER) and the error-free DNA damage bypass pathway. On the other hand, deficiencies in BER did not have a profound effect on A3B-induced mutations in CAN1, the reporter for protein coding genes. We hypothesize that differences in the mechanisms underlying ssDNA formation at tRNA genes and other genomic loci are the key determinants of the choice of the repair pathways and consequently the efficiency of DNA damage repair in these regions. Overall, our results indicate that tRNA genes are highly susceptible to ssDNA-specific DNA damaging agents. However, increased DNA repair efficacy in tRNA genes can prevent their hypermutation and maintain both genome and proteome homeostasis.

The Impact of Environmental and Endogenous Damage on Somatic Mutation Load in Human Skin Fibroblasts.

Accumulation of somatic changes, due to environmental and endogenous lesions, in the human genome is associated with aging and cancer. Understanding the impacts of these processes on mutagenesis is fundamental to understanding the etiology, and improving the prognosis and prevention of cancers and other genetic diseases. Previous methods relying on either the generation of induced pluripotent stem cells, or sequencing of single-cell genomes were inherently error-prone and did not allow independent validation of the mutations. In the current study we eliminated these potential sources of error by high coverage genome sequencing of single-cell derived clonal fibroblast lineages, obtained after minimal propagation in culture, prepared from skin biopsies of two healthy adult humans. We report here accurate measurement of genome-wide magnitude and spectra of mutations accrued in skin fibroblasts of healthy adult humans. We found that every cell contains at least one chromosomal rearrangement and 600­13,000 base substitutions. The spectra and correlation of base substitutions with epigenomic features resemble many cancers. Moreover, because biopsies were taken from body parts differing by sun exposure, we can delineate the precise contributions of environmental and endogenous factors to the accrual of genetic changes within the same individual. We show here that UV-induced and endogenous DNA damage can have a comparable impact on the somatic mutation loads in skin fibroblasts. Trial Registration: ClinicalTrials.gov NCT01087307.

APOBEC3A and APOBEC3B Preferentially Deaminate the Lagging Strand Template during DNA Replication.

APOBEC family cytidine deaminases have recently been implicated as powerful mutators of cancer genomes. How APOBECs, which are ssDNA-specific enzymes, gain access to chromosomal DNA is unclear. To ascertain the chromosomal ssDNA substrates of the APOBECs, we expressed APOBEC3A and APOBEC3B, the two most probable APOBECs mediating cancer mutagenesis, in a yeast model system. We demonstrate, using mutation reporters and whole genome sequencing, that APOBEC3A- and APOBEC3B-induced mutagenesis primarily results from the deamination of the lagging strand template during DNA replication. Moreover, our results indicate that both genetic deficiencies in replication fork-stabilizing proteins and chemical induction of replication stress greatly augment the mutagenesis of APOBEC3A and APOBEC3B. Taken together, these results strongly indicate that ssDNA formed during DNA lagging strand synthesis is a major substrate for APOBECs and may be the principal substrate in human cancers experiencing replication stress.

An APOBEC3A hypermutation signature is distinguishable from the signature of background mutagenesis by APOBEC3B in human cancers.

Elucidation of mutagenic processes shaping cancer genomes is a fundamental problem whose solution promises insights into new treatment, diagnostic and prevention strategies. Single-strand DNA-specific APOBEC cytidine deaminase(s) are major source(s) of mutation in several cancer types. Previous indirect evidence implicated APOBEC3B as the more likely major mutator deaminase, whereas the role of APOBEC3A is not established. Using yeast models enabling the controlled generation of long single-strand genomic DNA substrates, we show that the mutation signatures of APOBEC3A and APOBEC3B are statistically distinguishable. We then apply three complementary approaches to identify cancer samples with mutation signatures resembling either APOBEC. Strikingly, APOBEC3A-like samples have over tenfold more APOBEC-signature mutations than APOBEC3B-like samples. We propose that APOBEC3A-mediated mutagenesis is much more frequent because APOBEC3A itself is highly proficient at generating DNA breaks, whose repair can trigger the formation of single-strand hypermutation substrates.

Cavitation Enhancing Nanodroplets Mediate Efficient DNA Fragmentation in a Bench Top Ultrasonic Water Bath.

A perfluorocarbon nanodroplet formulation is shown to be an effective cavitation enhancement agent, enabling rapid and consistent fragmentation of genomic DNA in a standard ultrasonic water bath. This nanodroplet-enhanced method produces genomic DNA libraries and next-generation sequencing results indistinguishable from DNA samples fragmented in dedicated commercial acoustic sonication equipment, and with higher throughput. This technique thus enables widespread access to fast bench-top genomic DNA fragmentation.

Tracking replication enzymology in vivo by genome-wide mapping of ribonucleotide incorporation.

Ribonucleotides are frequently incorporated into DNA during replication in eukaryotes. Here we map genome-wide distribution of these ribonucleotides as markers of replication enzymology in budding yeast, using a new 5' DNA end-mapping method, hydrolytic end sequencing (HydEn-seq). HydEn-seq of DNA from ribonucleotide excision repair-deficient strains reveals replicase- and strand-specific patterns of ribonucleotides in the nuclear genome. These patterns support the roles of DNA polymerases α and δ in lagging-strand replication and of DNA polymerase ɛ in leading-strand replication. They identify replication origins, termination zones and variations in ribonucleotide incorporation frequency across the genome that exceed three orders of magnitude. HydEn-seq also reveals strand-specific 5' DNA ends at mitochondrial replication origins, thus suggesting unidirectional replication of a circular genome. Given the conservation of enzymes that incorporate and process ribonucleotides in DNA, HydEn-seq can be used to track replication enzymology in other organisms.

Heterogeneous polymerase fidelity and mismatch repair bias genome variation and composition.

Mutational heterogeneity must be taken into account when reconstructing evolutionary histories, calibrating molecular clocks, and predicting links between genes and disease. Selective pressures and various DNA transactions have been invoked to explain the heterogeneous distribution of genetic variation between species, within populations, and in tissue-specific tumors. To examine relationships between such heterogeneity and variations in leading- and lagging-strand replication fidelity and mismatch repair, we accumulated 40,000 spontaneous mutations in eight diploid yeast strains in the absence of selective pressure. We found that replicase error rates vary by fork direction, coding state, nucleosome proximity, and sequence context. Further, error rates and DNA mismatch repair efficiency both vary by mismatch type, responsible polymerase, replication time, and replication origin proximity. Mutation patterns implicate replication infidelity as one driver of variation in somatic and germline evolution, suggest mechanisms of mutual modulation of genome stability and composition, and predict future observations in specific cancers.

Variant calling in low-coverage whole genome sequencing of a Native American population sample.

BACKGROUND: The reduction in the cost of sequencing a human genome has led to the use of genotype sampling strategies in order to impute and infer the presence of sequence variants that can then be tested for associations with traits of interest. Low-coverage Whole Genome Sequencing (WGS) is a sampling strategy that overcomes some of the deficiencies seen in fixed content SNP array studies. Linkage-disequilibrium (LD) aware variant callers, such as the program Thunder, may provide a calling rate and accuracy that makes a low-coverage sequencing strategy viable. RESULTS: We examined the performance of an LD-aware variant calling strategy in a population of 708 low-coverage whole genome sequences from a community sample of Native Americans. We assessed variant calling through a comparison of the sequencing results to genotypes measured in 641 of the same subjects using a fixed content first generation exome array. The comparison was made using the variant calling routines GATK Unified Genotyper program and the LD-aware variant caller Thunder. Thunder was found to improve concordance in a coverage dependent fashion, while correctly calling nearly all of the common variants as well as a high percentage of the rare variants present in the sample. CONCLUSIONS: Low-coverage WGS is a strategy that appears to collect genetic information intermediate in scope between fixed content genotyping arrays and deep-coverage WGS. Our data suggests that low-coverage WGS is a viable strategy with a greater chance of discovering novel variants and associations than fixed content arrays for large sample association analyses.