Double-strand breaks associated with repetitive DNA can reshape the genome

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs) are an important source of genetic variation underlying rapid adaptation and genome evolution. Whereas point mutation rates vary with genomic location and local DNA features, the role of genome architecture in the formation and evolutionary dynamics of CNVs is poorly understood. Previously, we found the GAP1 gene in Saccharomyces cerevisiae undergoes frequent amplification and selection in glutamine-limitation. The gene is flanked by two long terminal repeats (LTRs) and proximate to an origin of DNA replication (autonomously replicating sequence, ARS), which likely promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we evolved engineered strains lacking either the adjacent LTRs, ARS, or all elements in glutamine-limited chemostats. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. Removal of local DNA elements significantly impacts the fitness effect of GAP1 CNVs and the rate of adaptation. In 177 CNV lineages, across all four strains, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, distal ones mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination can mediate gene amplification following de novo retrotransposon events. Our study reveals that template switching during DNA replication is a prevalent source of adaptive CNVs.

Spontaneous and environment induced genomic alterations in yeast model

While genomic alterations are fundamental to biological evolution, enabling adaptation and diversity, they can also result in detrimental outcomes, such as the development of genetic diseases including cancer. The budding yeast Saccharomyces cerevisiae serves as an exemplary model for investigating the mechanisms behind various genomic alterations, including point mutations, chromosomal rearrangements, and whole-chromosome aneuploidy. In this review, we highlight the application of genetic screening systems to assess the mutagenic effects of physical and chemical agents efficiently. Additionally, we discuss the utilization of high-throughput sequencing technologies to uncover comprehensive genomic alterations and rare genetic events. We provide a detailed summary of the features of genomic alterations and discuss the genetic mechanisms driving these changes under both spontaneous and stress-induced conditions. Given the high conservation of DNA replication and repair machinery across different organisms, the insights gained from studies on yeast offer valuable perspectives for understanding the delicate balance between genome plasticity and integrity in other species.

Template switching during DNA replication is a prevalent source of adaptive gene amplification

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Transposable elements regulate thymus development and function

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease activity. Foldback priming at DNA double-strand breaks (DSBs) is one mechanism proposed for the generation of inverted duplications. However, the other pathway steps have not been fully elucidated. Here, we show that a DSB induced near natural inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity, but not Rad1 nuclease, trims the heterologous flaps formed after foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric inversion chromosome. Finally, we show that stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences or by deletion of one centromere.

Nonlethal Furfural Exposure Causes Genomic Alterations and Adaptability Evolution in Saccharomyces cerevisiae

Furfural is a major inhibitor found in lignocellulosic hydrolysate, a promising feedstock for the biofermentation industry. In this study, we aimed to investigate the potential impact of this furan-derived chemical on yeast genome integrity and phenotypic evolution by using genetic screening systems and high-throughput analyses. Our results showed that the rates of aneuploidy, chromosomal rearrangements (including large deletions and duplications), and loss of heterozygosity (LOH) increased by 50-fold, 23-fold, and 4-fold, respectively, when yeast cells were cultured in medium containing a nonlethal dose of furfural (0.6 g/L). We observed significantly different ratios of genetic events between untreated and furfural-exposed cells, indicating that furfural exposure induced a unique pattern of genomic instability. Furfural exposure also increased the proportion of CG-to-TA and CG-to-AT base substitutions among point mutations, which was correlated with DNA oxidative damage. Interestingly, although monosomy of chromosomes often results in the slower growth of yeast under spontaneous conditions, we found that monosomic chromosome IX contributed to the enhanced furfural tolerance. Additionally, terminal LOH events on the right arm of chromosome IV, which led to homozygosity of the SSD1 allele, were associated with furfural resistance. This study sheds light on the mechanisms underlying the influence of furfural on yeast genome integrity and adaptability evolution. IMPORTANCE Industrial microorganisms are often exposed to multiple environmental stressors and inhibitors during their application. This study demonstrates that nonlethal concentrations of furfural in the culture medium can significantly induce genome instability in the yeast Saccharomyces cerevisiae. Notably, furfural-exposed yeast cells displayed frequent chromosome aberrations, indicating the potent teratogenicity of this inhibitor. We identified specific genomic alterations, including monosomic chromosome IX and loss of heterozygosity of the right arm of chromosome IV, that confer furfural tolerance to a diploid S. cerevisiae strain. These findings enhance our understanding of how microorganisms evolve and adapt to stressful environments and offer insights for developing strategies to improve their performance in industrial applications.

Delineation of two multi-invasion-induced rearrangement pathways that differently affect genome stability

Punctuated bursts of structural genomic variations (SVs) have been described in various organisms, but their etiology remains incompletely understood. Homologous recombination (HR) is a template-guided mechanism of repair of DNA double-strand breaks and stalled or collapsed replication forks. We recently identified a DNA break amplification and genome rearrangement pathway originating from the endonucleolytic processing of a multi-invasion (MI) DNA joint molecule formed during HR. Genome-wide approaches confirmed that multi-invasion-induced rearrangement (MIR) frequently leads to several repeat-mediated SVs and aneuploidies. Using molecular and genetic analysis and a novel, highly sensitive proximity ligation-based assay for chromosomal rearrangement quantification, we further delineate two MIR subpathways. MIR1 is a universal pathway occurring in any sequence context, which generates secondary breaks and frequently leads to additional SVs. MIR2 occurs only if recombining donors exhibit substantial homology and results in sequence insertion without additional breaks or SVs. The most detrimental MIR1 pathway occurs late on a subset of persisting DNA joint molecules in a PCNA/Polδ-independent manner, unlike recombinational DNA synthesis. This work provides a refined mechanistic understanding of these HR-based SV formation pathways and shows that complex repeat-mediated SVs can occur without displacement DNA synthesis. Sequence signatures for inferring MIR1 from long-read data are proposed.

Saccharomyces cerevisiae DNA polymerase IV overcomes Rad51 inhibition of DNA polymerase δ in Rad52-mediated direct-repeat recombination

Saccharomyces cerevisiae DNA polymerase IV (Pol4) like its homolog, human DNA polymerase lambda (Polλ), is involved in Non-Homologous End-Joining and Microhomology-Mediated Repair. Using genetic analysis, we identified an additional role of Pol4 also in homology-directed DNA repair, specifically in Rad52-dependent/Rad51-independent direct-repeat recombination. Our results reveal that the requirement for Pol4 in repeat recombination was suppressed by the absence of Rad51, suggesting that Pol4 counteracts the Rad51 inhibition of Rad52-mediated repeat recombination events. Using purified proteins and model substrates, we reconstituted in vitro reactions emulating DNA synthesis during direct-repeat recombination and show that Rad51 directly inhibits Polδ DNA synthesis. Interestingly, although Pol4 was not capable of performing extensive DNA synthesis by itself, it aided Polδ in overcoming the DNA synthesis inhibition by Rad51. In addition, Pol4 dependency and stimulation of Polδ DNA synthesis in the presence of Rad51 occurred in reactions containing Rad52 and RPA where DNA strand-annealing was necessary. Mechanistically, yeast Pol4 displaces Rad51 from ssDNA independent of DNA synthesis. Together our in vitro and in vivo data suggest that Rad51 suppresses Rad52-dependent/Rad51-independent direct-repeat recombination by binding to the primer-template and that Rad51 removal by Pol4 is critical for strand-annealing dependent DNA synthesis.

High Chromosomal Reorganization and Presence of Microchromosomes in Chactidae Scorpions from the Brazilian Amazon

Scorpions are of particular interest in cytogenomic studies, as they can present a high incidence of chromosomal rearrangements heterozygous in natural populations. In this study, we cytogenetically analyzed four species of Chactidae. In Brotheas, 2n = 40 was observed in Brotheas silvestris, 2n = 48 in Brotheas paraensis, and 2n = 50 (cytotype A) or 2n = 52 (cytotype B) among populations of Brotheas amazonicus. Our results showed a bimodal karyotype in Neochactas parvulus, 2n = 54, with microchromosomes and a concentration of constitutive heterochromatin in macrochromosomes. The 45S rDNA is located in only one pair of the karyotype, with different heteromorphisms of clusters of this rDNA in the cytotype B of B. amazonicus, with NOR-bearing chromosomes involved in multi-chromosomal associations during meiosis I. The U2 snDNA was mapped in the interstitial region of distinct karyotype pairs of three Chactidae species. Our results indicate the possible formation of cryptic species in B. amazonicus; the different 45S rDNA configurations in the genome of this species may result from amplification and degeneration. We suggest that the bimodal karyotype in N. parvulus results from fusion/fission events and that the unequal distribution of repetitive DNAs between macro and microchromosomes contributes to the maintenance of its asymmetry.

Delineation of two multi-invasion-induced rearrangement pathways that differently affect genome stability

Punctuated bursts of structural genomic variations (SVs) have been described in various organisms, but their etiology remains incompletely understood. Homologous recombination (HR) is a template-guided mechanism of repair of DNA double-strand breaks and stalled or collapsed replication forks. We recently identified a DNA break amplification and genome rearrangement pathway originating from the endonucleolytic processing of a multi-invasion (MI) DNA joint molecule formed during HR. Genome-wide sequencing approaches confirmed that multi-invasion-induced rearrangement (MIR) frequently leads to several repeat-mediated SVs and aneuploidies. Using molecular and genetic analysis, and a novel, highly sensitive proximity ligation-based assay for chromosomal rearrangement quantification, we further delineate two MIR sub-pathways. MIR1 is a universal pathway occurring in any sequence context, which generates secondary breaks and frequently leads to additional SVs. MIR2 occurs only if recombining donors exhibit substantial homology, and results in sequence insertion without additional break or SV. The most detrimental MIR1 pathway occurs late on a subset of persisting DNA joint molecules in a PCNA/Polδ-independent manner, unlike recombinational DNA synthesis. This work provides a refined mechanistic understanding of these HR-based SV formation pathways and shows that complex repeat-mediated SVs can occur without displacement DNA synthesis. Sequence signatures for inferring MIR1 from long-read data are proposed.

Shuffling the yeast genome using CRISPR/Cas9-generated DSBs that target the transposable Ty1 elements

Although homologous recombination between transposable elements can drive genomic evolution in yeast by facilitating chromosomal rearrangements, the details of the underlying mechanisms are not fully clarified. In the genome of the yeast Saccharomyces cerevisiae, the most common class of transposon is the retrotransposon Ty1. Here, we explored how Cas9-induced double-strand breaks (DSBs) directed to Ty1 elements produce genomic alterations in this yeast species. Following Cas9 induction, we observed a significant elevation of chromosome rearrangements such as deletions, duplications and translocations. In addition, we found elevated rates of mitotic recombination, resulting in loss of heterozygosity. Using Southern analysis coupled with short- and long-read DNA sequencing, we revealed important features of recombination induced in retrotransposons. Almost all of the chromosomal rearrangements reflect the repair of DSBs at Ty1 elements by non-allelic homologous recombination; clustered Ty elements were hotspots for chromosome rearrangements. In contrast, a large proportion (about three-fourths) of the allelic mitotic recombination events have breakpoints in unique sequences. Our analysis suggests that some of the latter events reflect extensive processing of the broken ends produced in the Ty element that extend into unique sequences resulting in break-induced replication. Finally, we found that haploid and diploid strain have different preferences for the pathways used to repair double-stranded DNA breaks. Our findings demonstrate the importance of DNA lesions in retrotransposons in driving genome evolution.

Double-strand breaks induce inverted duplication chromosome rearrangements by a DNA polymerase δ and Rad51-dependent mechanism

Inverted duplications, also known as foldback inversions, are commonly observed in cancers and are the major class of chromosome rearrangement recovered from yeast cells lacking Mre11 nuclease. Foldback priming at naturally occurring inverted repeats is one mechanism proposed for the generation of inverted duplications. However, the initiating lesion for these events and the mechanism by which they form has not been fully elucidated. Here, we show that a DNA double-strand break (DSB) induced near natural short, inverted repeats drives high frequency inverted duplication in Sae2 and Mre11-deficient cells. We find that DNA polymerase δ proof-reading activity acts non-redundantly with Rad1 nuclease to remove heterologous tails formed during foldback annealing. Additionally, Pol32 is required for the generation of inverted duplications, suggesting that Pol δ catalyzes fill-in synthesis primed from the foldback to create a hairpin-capped chromosome that is subsequently replicated to form a dicentric isochromosome. Stabilization of the dicentric chromosome after breakage involves telomere capture by non-reciprocal translocation mediated by repeat sequences and requires Rad51.

RAD51 Is Implicated in DNA Damage, Chemoresistance and Immune Dysregulation in Solid Tumors

BACKGROUND

In normal cells, homologous recombination (HR) is tightly regulated and plays an important role in the maintenance of genomic integrity and stability through precise repair of DNA damage. RAD51 is a recombinase that mediates homologous base pairing and strand exchange during DNA repair by HR. Our previous data in multiple myeloma and esophageal adenocarcinoma (EAC) show that dysregulated HR mediates genomic instability. Purpose of this study was to investigate role of HR in genomic instability, chemoresistance and immune dysregulation in solid tumors including colon and breast cancers.

METHODS

The GEO dataset were used to investigate correlation of RAD51 expression with patient survival and expression of various immune markers in EAC, breast and colorectal cancers. RAD51 was inhibited in cancer cell lines using shRNAs and a small molecule inhibitor. HR activity was evaluated using a plasmid-based assay, DNA breaks assessed by evaluating expression of γ-H2AX (a marker of DNA breaks) and p-RPA32 (a marker of DNA end resection) using Western blotting. Genomic instability was monitored by investigating micronuclei (a marker of genomic instability). Impact of RAD51 inhibitor and/or a DNA-damaging agent was assessed on viability and apoptosis in EAC, breast and colon cancer cell lines in vitro and in a subcutaneous tumor model of EAC. Impact of RAD51 inhibitor on expression profile was monitored by RNA sequencing.

RESULTS

Elevated RAD51 expression correlated with poor survival of EAC, breast and colon cancer patients. RAD51 knockdown in cancer cell lines inhibited DNA end resection and strand exchange activity (key steps in the initiation of HR) as well as spontaneous DNA breaks, whereas its overexpression increased DNA breaks and genomic instability. Treatment of EAC, colon and breast cancer cell lines with a small molecule inhibitor of RAD51 inhibited DNA breaking agent-induced DNA breaks and genomic instability. RAD51 inhibitor potentiated cytotoxicity of DNA breaking agent in all cancer cell types tested in vitro as well as in a subcutaneous model of EAC. Evaluation by RNA sequencing demonstrated that DNA repair and cell cycle related pathways were induced by DNA breaking agent whereas their induction either prevented or reversed by RAD51 inhibitor. In addition, immune-related pathways such as PD-1 and Interferon Signaling were also induced by DNA breaking agent whereas their induction prevented by RAD51 inhibitor. Consistent with these observations, elevated RAD51 expression also correlated with that of genes involved in inflammation and other immune surveillance.

CONCLUSIONS

Elevated expression of RAD51 and associated HR activity is involved in spontaneous and DNA damaging agent-induced DNA breaks and genomic instability thus contributing to chemoresistance, immune dysregulation and poor prognosis in cancer. Therefore, inhibitors of RAD51 have great potential as therapeutic agents for EAC, colon, breast and probably other solid tumors.

The contribution of DNA repair pathways to genome editing and evolution in filamentous pathogens

DNA double-strand breaks require repair or risk corrupting the language of life. To ensure genome integrity and viability, multiple DNA double-strand break repair pathways function in eukaryotes. Two such repair pathways, canonical non-homologous end joining and homologous recombination, have been extensively studied, while other pathways such as microhomology-mediated end joint and single-strand annealing, once thought to serve as back-ups, now appear to play a fundamental role in DNA repair. Here, we review the molecular details and hierarchy of these four DNA repair pathways, and where possible, a comparison for what is known between animal and fungal models. We address the factors contributing to break repair pathway choice, and aim to explore our understanding and knowledge gaps regarding mechanisms and regulation in filamentous pathogens. We additionally discuss how DNA double-strand break repair pathways influence genome engineering results, including unexpected mutation outcomes. Finally, we review the concept of biased genome evolution in filamentous pathogens, and provide a model, termed Biased Variation, that links DNA double-strand break repair pathways with properties of genome evolution. Despite our extensive knowledge for this universal process, there remain many unanswered questions, for which the answers may improve genome engineering and our understanding of genome evolution.

Ribodysgenesis: sudden genome instability in the yeast Saccharomyces cerevisiae arising from RNase H2 cleavage at genomic-embedded ribonucleotides

Ribonucleotides can be incorporated into DNA during replication by the replicative DNA polymerases. These aberrant DNA subunits are efficiently recognized and removed by Ribonucleotide Excision Repair, which is initiated by the heterotrimeric enzyme RNase H2. While RNase H2 is essential in higher eukaryotes, the yeast Saccharomyces cerevisiae can survive without RNase H2 enzyme, although the genome undergoes mutation, recombination and other genome instability events at an increased rate. Although RNase H2 can be considered as a protector of the genome from the deleterious events that can ensue from recognition and removal of embedded ribonucleotides, under conditions of high ribonucleotide incorporation and retention in the genome in a RNase H2-negative strain, sudden introduction of active RNase H2 causes massive DNA breaks and genome instability in a condition which we term 'ribodysgenesis'. The DNA breaks and genome instability arise solely from RNase H2 cleavage directed to the ribonucleotide-containing genome. Survivors of ribodysgenesis have massive loss of heterozygosity events stemming from recombinogenic lesions on the ribonucleotide-containing DNA, with increases of over 1000X from wild-type. DNA breaks are produced over one to two divisions and subsequently cells adapt to RNase H2 and ribonucleotides in the genome and grow with normal levels of genome instability.

Bursts of Genomic Instability Potentiate Phenotypic and Genomic Diversification in Saccharomyces cerevisiae

How microbial cells leverage their phenotypic potential to survive in a changing environment is a complex biological problem, with important implications for pathogenesis and species evolution. Stochastic phenotype switching, a particularly fascinating adaptive approach observed in numerous species across the tree of life, introduces phenotypic diversity into a population through mechanisms which have remained difficult to define. Here we describe our investigations into the mechanistic basis of colony morphology phenotype switching which occurs in populations of a pathogenic isolate of Saccharomyces cerevisiae, YJM311. We observed that clonal populations of YJM311 cells produce variant colonies that display altered morphologies and, using whole genome sequence analysis, discovered that these variant clones harbored an exceptional collection of karyotypes newly altered by de novo structural genomic variations (SVs). Overall, our analyses indicate that copy number alterations, more often than changes in allelic identity, provide the causative basis of this phenotypic variation. Individual variants carried between 1 and 16 de novo copy number variations, most of which were whole chromosomal aneuploidies. Notably, we found that the inherent stability of the diploid YJM311 genome is comparable to that of domesticated laboratory strains, indicating that the collections of SVs harbored by variant clones did not arise by a chronic chromosomal instability (CIN) mechanism. Rather, our data indicate that these variant clones acquired such complex karyotypic configurations simultaneously, during stochastic and transient episodes of punctuated systemic genomic instability (PSGI). Surprisingly, we found that the majority of these highly altered variant karyotypes were propagated with perfect fidelity in long-term passaging experiments, demonstrating that high aneuploidy burdens can often be conducive with prolonged genomic integrity. Together, our results demonstrate that colony morphology switching in YJM311 is driven by a stochastic process in which genome stability and plasticity are integrally coupled to phenotypic heterogeneity. Consequently, this system simultaneously introduces both phenotypic and genomic variation into a population of cells, which can, in turn perpetuate population diversity for many generations thereafter.

Limitations of mouse models for sickle cell disease conferred by their human globin transgene configurations

We characterized the human β-like globin transgenes in two mouse models of sickle cell disease (SCD) and tested a genome-editing strategy to induce red blood cell fetal hemoglobin (HbF; α₂γ₂). Berkeley SCD mice contain four to 22 randomly arranged, fragmented copies of three human transgenes (HBA1, HBG2-HBG1-HBD-HBB^S and a mini-locus control region) integrated into a single site of mouse chromosome 1. Cas9 disruption of the BCL11A repressor binding motif in the γ-globin gene (HBG1 and HBG2; HBG) promoters of Berkeley mouse hematopoietic stem cells (HSCs) caused extensive death from multiple double-strand DNA breaks. Long-range sequencing of Townes SCD mice verified that the endogenous Hbb genes were replaced by single-copy segments of human HBG1 and HBB^S including proximal but not some distal gene-regulatory elements. Townes mouse HSCs were viable after Cas9 disruption of the HBG1 BCL11A binding motif but failed to induce HbF to therapeutic levels, contrasting with human HSCs. Our findings provide practical information on the genomic structures of two common mouse SCD models, illustrate their limitations for analyzing therapies to induce HbF and confirm the importance of distal DNA elements in human globin regulation.

This article has an associated First Person interview with the first author of the paper.

Editor's choice: This study describes the genomic structures of two common sickle cell disease mouse models, illustrates their limitations for analyzing some genetic therapies and confirms the importance of distal DNA elements in human globin gene regulation.

Genomic structural plasticity of rodent-associated Bartonella in nature

Rodent‐associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1–1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter‐ and intraspecies competition and adaptation.

Genome-wide sequencing analysis of Sgs1, Exo1, Rad51, and Srs2 in DNA repair by homologous recombination

Homologous recombination is essential to maintain genome stability in response to DNA damage. Here, we have used genome-wide sequencing to quantitatively analyze at nucleotide resolution the dynamics of DNA end resection, re-synthesis, and gene conversion at a double-strand break. Resection initiates asymmetrically in an MRX-independent manner before proceeding steadily in both directions. Sgs1, Exo1, Rad51, and Srs2 differently regulate the rate and symmetry of early and late resection. Exo1 also ensures the coexistence of resection and re-synthesis, while Srs2 guarantees a constant and symmetrical DNA re-polymerization. Gene conversion is MMR independent, spans only a minor fraction of the resected region, and its unidirectionality depends on Srs2. Finally, these repair factors prevent the development of alterations remote from the DNA lesion, such as subtelomeric instability, duplication of genomic regions, and over-replication of Ty elements. Altogether, this approach allows a quantitative analysis and a direct genome-wide visualization of DNA repair by homologous recombination.

Genomic characterization of a wild diploid isolate of Saccharomyces cerevisiae reveals an extensive and dynamic landscape of structural variation

The budding yeast Saccharomyces cerevisiae has been extensively characterized for many decades and is a critical resource for the study of numerous facets of eukaryotic biology. Recently, whole genome sequence analysis of over 1000 natural isolates of S. cerevisiae has provided critical insights into the evolutionary landscape of this species by revealing a population structure comprised of numerous genomically diverse lineages. These survey-level analyses have been largely devoid of structural genomic information, mainly because short read sequencing is not suitable for detailed characterization of genomic architecture. Consequently, we still lack a complete perspective of the genomic variation the exists within the species. Single molecule long read sequencing technologies, such as Oxford Nanopore and PacBio, provide sequencing-based approaches with which to rigorously define the structure of a genome, and have empowered yeast geneticists to explore this poorly described realm of eukaryotic genomics. Here, we present the comprehensive genomic structural analysis of a wild diploid isolate of S. cerevisiae, YJM311. We used long read sequence analysis to construct a haplotype-phased, telomere-to-telomere length assembly of the YJM311 genome and characterized the structural variations (SVs) therein. We discovered that the genome of YJM311 contains significant intragenomic structural variation, some of which imparts notable consequences to the genomic stability and developmental biology of the strain. Collectively, we outline a new methodology for creating accurate haplotype-phased genome assemblies and highlight how such genomic analyses can define the structural architectures of S. cerevisiae isolates. It is our hope that continued structural characterization of S. cerevisiae genomes, such as we have reported here for YJM311, will comprehensively advance our understanding of eukaryotic genome structure-function relationships, structural genomic diversity, and evolution.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

Roles of homologous recombination in response to ionizing radiation-induced DNA damage

PURPOSE

Ionizing radiation induces a vast array of DNA lesions including base damage, and single- and double-strand breaks (SSB, DSB). DSBs are among the most cytotoxic lesions, and mis-repair causes small- and large-scale genome alterations that can contribute to carcinogenesis. Indeed, ionizing radiation is a 'complete' carcinogen. DSBs arise immediately after irradiation, termed 'frank DSBs,' as well as several hours later in a replication-dependent manner, termed 'secondary' or 'replication-dependent DSBs. DSBs resulting from replication fork collapse are single-ended and thus pose a distinct problem from two-ended, frank DSBs. DSBs are repaired by error-prone nonhomologous end-joining (NHEJ), or generally error-free homologous recombination (HR), each with sub-pathways. Clarifying how these pathways operate in normal and tumor cells is critical to increasing tumor control and minimizing side effects during radiotherapy.

CONCLUSIONS

The choice between NHEJ and HR is regulated during the cell cycle and by other factors. DSB repair pathways are major contributors to cell survival after ionizing radiation, including tumor-resistance to radiotherapy. Several nucleases are important for HR-mediated repair of replication-dependent DSBs and thus replication fork restart. These include three structure-specific nucleases, the 3' MUS81 nuclease, and two 5' nucleases, EEPD1 and Metnase, as well as three end-resection nucleases, MRE11, EXO1, and DNA2. The three structure-specific nucleases evolved at very different times, suggesting incremental acceleration of replication fork restart to limit toxic HR intermediates and genome instability as genomes increased in size during evolution, including the gain of large numbers of HR-prone repetitive elements. Ionizing radiation also induces delayed effects, observed days to weeks after exposure, including delayed cell death and delayed HR. In this review we highlight the roles of HR in cellular responses to ionizing radiation, and discuss the importance of HR as an exploitable target for cancer radiotherapy.

Recombination-mediated genome rearrangements

Homologous recombination (HR) is a universal DNA double-strand break (DSB) repair pathway that uses an intact DNA molecule as a template. Signature HR reactions are homology search and DNA strand invasion catalyzed by the prototypical RecA-ssDNA filament (Rad51 and Dmc1 in eukaryotes), which produces heteroduplex DNA-containing joint molecules (JMs). These reactions uniquely infringe on the DNA strands association established at replication, on the basis of substantial sequence similarity. For that reason, and despite the high fidelity of its templated nature, DSB repair by HR authorizes the alteration of genome structure, guided by repetitive DNA elements. The resulting structural variations (SVs) can involve vast genomic regions, potentially affecting multiple coding sequences and regulatory elements at once, with possible pathological consequences. Here, we discuss recent advances in our understanding of genetic and molecular vulnerabilities of HR leading to SVs, and of the various fidelity-enforcing factors acting across scales on the balancing act of this complex pathway. An emphasis is put on extra-chomosomal DNAs, both product of, and substrate for HR-mediated chromosomal rearrangements.

The dark side of homology-directed repair

DNA double strand breaks (DSB) are cytotoxic lesions that can lead to genome rearrangements and genomic instability, which are hallmarks of cancer. The two main DSB repair pathways are non-homologous end joining and homologous recombination (HR). While HR is generally highly accurate, it has the potential for rearrangements that occur directly or through intermediates generated during the repair process. Whole genome sequencing of cancers has revealed numerous types of structural rearrangement signatures that are often indicative of repair mediated by sequence homology. However, it can be challenging to delineate repair mechanisms from sequence analysis of rearrangement end products from cancer genomes, or even model systems, because the same rearrangements can be generated by different pathways. Here, we review homology-directed repair pathways and their consequences. Exploring those pathways can lead to a greater understanding of rearrangements that occur in cancer cells.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Multi-omic analysis elucidates the genetic basis of hydrocephalus

We conducted PrediXcan analysis of hydrocephalus risk in ten neurological tissues and whole blood. Decreased expression of MAEL in the brain was significantly associated (Bonferroni-adjusted p < 0.05) with hydrocephalus. PrediXcan analysis of brain imaging and genomics data in the independent UK Biobank (N = 8,428) revealed that MAEL expression in the frontal cortex is associated with white matter and total brain volumes. Among the top differentially expressed genes in brain, we observed a significant enrichment for gene-level associations with these structural phenotypes, suggesting an effect on disease risk through regulation of brain structure and integrity. We found additional support for these genes through analysis of the choroid plexus transcriptome of a murine model of hydrocephalus. Finally, differential protein expression analysis in patient cerebrospinal fluid recapitulated disease-associated expression changes in neurological tissues, but not in whole blood. Our findings provide convergent evidence highlighting the importance of tissue-specific pathways and mechanisms in the pathophysiology of hydrocephalus.

Learning Yeast Genetics from Miro Radman

Miroslav Radman's far-sighted ideas have penetrated many aspects of our study of the repair of broken eukaryotic chromosomes. For over 35 years my lab has studied different aspects of the repair of chromosomal breaks in the budding yeast, Saccharomyces cerevisiae. From the start, we have made what we thought were novel observations that turned out to have been predicted by Miro's extraordinary work in the bacterium Escherichia coli and then later in the radiation-resistant Dienococcus radiodurans. In some cases, we have been able to extend some of his ideas a bit further.

Origin, Regulation, and Fitness Effect of Chromosomal Rearrangements in the Yeast Saccharomyces cerevisiae

Chromosomal rearrangements comprise unbalanced structural variations resulting in gain or loss of DNA copy numbers, as well as balanced events including translocation and inversion that are copy number neutral, both of which contribute to phenotypic evolution in organisms. The exquisite genetic assay and gene editing tools available for the model organism Saccharomyces cerevisiae facilitate deep exploration of the mechanisms underlying chromosomal rearrangements. We discuss here the pathways and influential factors of chromosomal rearrangements in S. cerevisiae. Several methods have been developed to generate on-demand chromosomal rearrangements and map the breakpoints of rearrangement events. Finally, we highlight the contributions of chromosomal rearrangements to drive phenotypic evolution in various S. cerevisiae strains. Given the evolutionary conservation of DNA replication and recombination in organisms, the knowledge gathered in the small genome of yeast can be extended to the genomes of higher eukaryotes.

Genome-Wide Analysis of Mitotic Recombination in Budding Yeast

DNA break lesions pose a serious threat to the integrity of the genome. Eukaryotic cells can repair these lesions using the homologous recombination pathway that guides the repair reaction by using a homologous DNA template. The budding yeast Saccharomyces cerevisiae is an excellent model system with which to study this repair mechanism and the resulting patterns of genomic change resulting from it. In this chapter, we describe an approach that utilizes whole-genome sequencing data to support the analysis of tracts of loss-of-heterozygosity (LOH) that can arise from mitotic recombination in the context of the entire diploid yeast genome. The workflow and the discussion in this chapter are intended to enable classically trained molecular biologists and geneticists with limited experience in computational methods to conceptually understand and execute the steps of genome-wide LOH analysis as well as to adapt and apply them to their own specific studies and experimental models.

A Versatile Protocol to Generate Translocations in Yeast Genomes Using CRISPR/Cas9

Genomic engineering methods represent powerful tools to examine chromosomal modifications and to subsequently study their impacts on cellular phenotypes. However, quantifying the fitness impact of translocations, independently from base substitutions or the insertion of genetic markers, remains a challenge. Here we report a rapid and straightforward protocol for engineering either targeted reciprocal translocations at the base pair level of resolution between two chromosomes or multiple simultaneous rearrangements in the yeast genome, without inserting any marker sequence in the chromosomes. Our CRISPR/Cas9-based method consists of inducing either (1) two double-strand breaks (DSBs) in two different chromosomes with two distinct guide RNAs (gRNAs) while providing specifically designed homologous donor DNA forcing the trans-repair of chromosomal extremities to generate a targeted reciprocal translocation or (2) multiple DSBs with a single gRNA targeting dispersed repeated sequences and leaving endogenous uncut copies of the repeat to be used as donor DNA, thereby generating multiple translocations, often associated with large segmental duplications (Fleiss, et al. PLoS Genet 15:e1008332, 2019).

Genomic Convergence in the Adaptation to Extreme Environments

Convergent evolution is especially common in plants that have independently adapted to the same extreme environments (i.e., extremophile plants). The recent burst of omics data has alleviated many limitations that have hampered molecular convergence studies of non-model extremophile plants. In this review, we summarize cases of genomic convergence in these taxa to examine the extent and type of genomic convergence during the process of adaptation to extreme environments. Despite being well studied by candidate gene approaches, convergent evolution at individual sites is rare and often has a high false-positive rate when assessed in whole genomes. By contrast, genomic convergence at higher genetic levels has been detected during adaptation to the same extreme environments. Examples include the convergence of biological pathways and changes in gene expression, gene copy number, amino acid usage, and GC content. Higher convergence levels play important roles in the adaptive evolution of extremophiles and may be more frequent and involve more genes. In several cases, multiple types of convergence events have been found to co-occur. However, empirical and theoretical studies of this higher level convergent evolution are still limited. In conclusion, both the development of powerful approaches and the detection of convergence at various genetic levels are needed to further reveal the genetic mechanisms of plant adaptation to extreme environments.

What do we know about the biology of the emerging fungal pathogen of humans Candida auris?

Candida auris is a worrisome fungal pathogen of humans which emerged merely about a decade ago. Ever since then the scientific community worked hard to understand clinically relevant traits, such as virulence factors, antifungal resistance mechanisms, and its ability to adhere to human skin and medical devices. Whole-genome sequencing of clinical isolates and epidemiological studies outlining the path of nosocomial outbreaks have been the focus of research into this pathogenic and multidrug-resistant yeast since its first description in 2009. More recently, work was started by several laboratories to explore the biology of C. auris. Here, we review the insights of studies characterizing the mechanisms underpinning antifungal drug resistance, biofilm formation, morphogenetic switching, cell aggregation, virulence, and pathogenicity of C. auris. We conclude that, although some progress has been made, there is still a long journey ahead of us, before we fully understand this novel pathogen. Critically important is the development of molecular tools for C. auris to make this fungus genetically tractable and traceable. This will allow an in-depth molecular dissection of the life cycle of C. auris, of its characteristics while interacting with the human host, and the mechanisms it employs to avoid being killed by antifungals and the immune system.

Long-read only assembly of Drechmeria coniospora genomes reveals widespread chromosome plasticity and illustrates the limitations of current nanopore methods

BACKGROUND

Long-read sequencing is increasingly being used to determine eukaryotic genomes. We used nanopore technology to generate chromosome-level assemblies for 3 different strains of Drechmeria coniospora, a nematophagous fungus used extensively in the study of innate immunity in Caenorhabditis elegans.

RESULTS

One natural geographical isolate demonstrated high stability over decades, whereas a second isolate not only had a profoundly altered genome structure but exhibited extensive instability. We conducted an in-depth analysis of sequence errors within the 3 genomes and established that even with state-of-the-art tools, nanopore methods alone are insufficient to generate eukaryotic genome sequences of sufficient accuracy to merit inclusion in public databases.

CONCLUSIONS

Although nanopore long-read sequencing is not accurate enough to produce publishable eukaryotic genomes, in our case, it has revealed new information about genome plasticity in D. coniospora and provided a backbone that will permit future detailed study to characterize gene evolution in this important model fungal pathogen.

Widespread and tissue-specific expression of endogenous retroelements in human somatic tissues

BACKGROUND

Endogenous retroelements (EREs) constitute about 42% of the human genome and have been implicated in common human diseases such as autoimmunity and cancer. The dominant paradigm holds that EREs are expressed in embryonic stem cells (ESCs) and germline cells but are repressed in differentiated somatic cells. Despite evidence that some EREs can be expressed at the RNA and protein levels in specific contexts, a system-level evaluation of their expression in human tissues is lacking.

METHODS

Using RNA sequencing data, we analyzed ERE expression in 32 human tissues and cell types, including medullary thymic epithelial cells (mTECs). A tissue specificity index was computed to identify tissue-restricted ERE families. We also analyzed the transcriptome of mTECs in wild-type and autoimmune regulator (AIRE)-deficient mice. Finally, we developed a proteogenomic workflow combining RNA sequencing and mass spectrometry (MS) in order to evaluate whether EREs might be translated and generate MHC I-associated peptides (MAP) in B-lymphoblastoid cell lines (B-LCL) from 16 individuals.

RESULTS

We report that all human tissues express EREs, but the breadth and magnitude of ERE expression are very heterogeneous from one tissue to another. ERE expression was particularly high in two MHC I-deficient tissues (ESCs and testis) and one MHC I-expressing tissue, mTECs. In mutant mice, we report that the exceptional expression of EREs in mTECs was AIRE-independent. MS analyses identified 103 non-redundant ERE-derived MAPs (ereMAPs) in B-LCLs. These ereMAPs preferentially derived from sense translation of intronic EREs. Notably, detailed analyses of their amino acid composition revealed that ERE-derived MAPs presented homology to viral MAPs.

CONCLUSIONS

This study shows that ERE expression in somatic tissues is more pervasive and heterogeneous than anticipated. The high and diversified expression of EREs in mTECs and their ability to generate MAPs suggest that EREs may play an important role in the establishment of self-tolerance. The viral-like properties of ERE-derived MAPs suggest that those not expressed in mTECs can be highly immunogenic.

Diverse transposable element landscapes in pathogenic and nonpathogenic yeast models: the value of a comparative perspective

Genomics and other large-scale analyses have drawn increasing attention to the potential impacts of transposable elements (TEs) on their host genomes. However, it remains challenging to transition from identifying potential roles to clearly demonstrating the level of impact TEs have on genome evolution and possible functions that they contribute to their host organisms. I summarize TE content and distribution in four well-characterized yeast model systems in this review: the pathogens Candida albicans and Cryptococcus neoformans, and the nonpathogenic species Saccharomyces cerevisiae and Schizosaccharomyces pombe. I compare and contrast their TE landscapes to their lifecycles, genomic features, as well as the presence and nature of RNA interference pathways in each species to highlight the valuable diversity represented by these models for functional studies of TEs. I then review the regulation and impacts of the Ty1 and Ty3 retrotransposons from Saccharomyces cerevisiae and Tf1 and Tf2 retrotransposons from Schizosaccharomyces pombe to emphasize parallels and distinctions between these well-studied elements. I propose that further characterization of TEs in the pathogenic yeasts would enable this set of four yeast species to become an excellent set of models for comparative functional studies to address outstanding questions about TE-host relationships.

Genomic Instability in Fungal Plant Pathogens

Fungi and fungal-like organisms (oomycetes) that cause diseases in plants have impacted human communities for centuries and probably from the dawn of agriculture. In modern agriculture, there is a constant race between new strategies to manage fungal plant pathogens and their ability to adapt. An important component in this race is fungal genetic diversity. Mechanisms such as sexual and parasexual recombination that contribute to the creation of novel allele combinations in fungal plant pathogens are briefly discussed in the first part of this review. Advances in genomics have enabled the investigation of chromosomal aberrations of agriculturally important fungal isolates at the nucleotide level. Some of these cases are summarized in the second part of this review; it is claimed that the effect of chromosomal aberrations on pathogenicity should be studied mechanistically. More data on the effect of gene copy number variations on phenotypes that are relevant to agriculture are especially needed. Genome rearrangements through translocations have shaped the genome of fungal plant pathogens by creating lineage-specific chromosome territories encoding for genes participating in plant diseases. Pathogenicity chromosomes are unique cases of such lineage-specific genetic elements, interestingly these chromosomes can be transferred horizontally and thus transforming a non-pathogenic strain to a pathogenic one. The third part of this review describes our attempts to reveal mutators in fungal plant pathogens by identifying fungi that lack important DNA repair genes or respond to DNA damage in an unconventional way. We found that a group of fungal plant pathogens lack conserved genes that are needed for an important Holliday junction resolution pathway. In addition, in Fusarium oxysporum, the rate-limiting step in dNTP production is not induced under DNA replication stress. This is very different from organisms from bacteria to humans. It remains to be seen if these mechanisms promote genetic instability in fungal plant pathogens.

Centromere scission drives chromosome shuffling and reproductive isolation

A fundamental characteristic of eukaryotic organisms is the generation of genetic variation via sexual reproduction. Conversely, significant large-scale genome structure variations could hamper sexual reproduction, causing reproductive isolation and promoting speciation. The underlying processes behind large-scale genome rearrangements are not well understood and include chromosome translocations involving centromeres. Recent genomic studies in the Cryptococcus species complex revealed that chromosome translocations generated via centromere recombination have reshaped the genomes of different species. In this study, multiple DNA double-strand breaks (DSBs) were generated via the CRISPR/Cas9 system at centromere-specific retrotransposons in the human fungal pathogen Cryptococcus neoformans The resulting DSBs were repaired in a complex manner, leading to the formation of multiple interchromosomal rearrangements and new telomeres, similar to chromothripsis-like events. The newly generated strains harboring chromosome translocations exhibited normal vegetative growth but failed to undergo successful sexual reproduction with the parental wild-type strain. One of these strains failed to produce any spores, while another produced ∼3% viable progeny. The germinated progeny exhibited aneuploidy for multiple chromosomes and showed improved fertility with both parents. All chromosome translocation events were accompanied without any detectable change in gene sequences and thus suggest that chromosomal translocations alone may play an underappreciated role in the onset of reproductive isolation and speciation.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Computational search for UV radiation resistance strategies in Deinococcus swuensis isolated from Paramo ecosystems

Ultraviolet radiation (UVR) is widely known as deleterious for many organisms since it can cause damage to biomolecules either directly or indirectly via the formation of reactive oxygen species. The goal of this study was to analyze the capacity of high-mountain Espeletia hartwegiana plant phyllosphere microorganisms to survive UVR and to identify genes related to resistance strategies. A strain of Deinococcus swuensis showed a high survival rate of up to 60% after UVR treatment at 800J/m2 and was used for differential expression analysis using RNA-seq after exposing cells to 400J/m2 of UVR (with >95% survival rate). Differentially expressed genes were identified using the R-Bioconductor package NOISeq and compared with other reported resistance strategies reported for this genus. Genes identified as being overexpressed included transcriptional regulators and genes involved in protection against damage by UVR. Non-coding (nc)RNAs were also differentially expressed, some of which have not been previously implicated. This study characterized the immediate radiation response of D. swuensis and indicates the involvement of ncRNAs in the adaptation to extreme environmental conditions.

Partner Choice in Spontaneous Mitotic Recombination in Wild Type and Homologous Recombination Mutants of Candida albicans

Candida albicans, the most common fungal pathogen, is a diploid with a genome that is rich in repeats and has high levels of heterozygosity. To study the role of different recombination pathways on direct-repeat recombination, we replaced either allele of the RAD52 gene (Chr6) with the URA-blaster cassette (hisG-URA3-hisG), measured rates of URA3 loss as resistance to 5-fluoroorotic acid (5FOA^R) and used CHEF Southern hybridization and SNP-RFLP analysis to identify recombination mechanisms and their frequency in wildtype and recombination mutants. FOA^R rates varied little across different strain backgrounds. In contrast, the type and frequency of mechanisms underlying direct repeat recombination varied greatly. For example, wildtype, rad59 and lig4 strains all displayed a bias for URA3 loss via pop-out/deletion vs. inter-homolog recombination and this bias was reduced in rad51 mutants. In addition, in rad51-derived 5FOA^R strains direct repeat recombination was associated with ectopic translocation (5%), chromosome loss/truncation (14%) and inter-homolog recombination (6%). In the absence of RAD52, URA3 loss was mostly due to chromosome loss and truncation (80–90%), and the bias of retained allele frequency points to the presence of a recessive lethal allele on Chr6B. However, a few single-strand annealing (SSA)-like events were identified and these were independent of either Rad59 or Lig4. Finally, the specific sizes of Chr6 truncations suggest that the inserted URA-blaster could represent a fragile site.

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.

This study uses yeast expressing a human cytidine deaminase to reveal simultaneous stretches of long single-strand DNA and multiple vast mutation clusters in a single eukaryotic cell repairing multiple double-strand breaks. This is reminiscent of the phenomenon of “kataegis” or hypermutation observed in cancer genomes, suggesting that a similar mechanism is involved.

Genome architecture and stability in the Saccharomyces cerevisiae knockout collection

Despite major progress in defining the functional roles of genes, a complete understanding of their influences is far from being realized, even in relatively simple organisms. A major milestone in this direction arose via the completion of the yeast Saccharomyces cerevisiae gene-knockout collection (YKOC), which has enabled high-throughput reverse genetics, phenotypic screenings and analyses of synthetic-genetic interactions1-3. Ensuing experimental work has also highlighted some inconsistencies and mistakes in the YKOC, or genome instability events that rebalance the effects of specific knockouts4-6, but a complete overview of these is lacking. The identification and analysis of genes that are required for maintaining genomic stability have traditionally relied on reporter assays and on the study of deletions of individual genes, but whole-genome-sequencing technologies now enable-in principle-the direct observation of genome instability globally and at scale. To exploit this opportunity, we sequenced the whole genomes of nearly all of the 4,732 strains comprising the homozygous diploid YKOC. Here, by extracting information on copy-number variation of tandem and interspersed repetitive DNA elements, we describe-for almost every single non-essential gene-the genomic alterations that are induced by its loss. Analysis of this dataset reveals genes that affect the maintenance of various genomic elements, highlights cross-talks between nuclear and mitochondrial genome stability, and shows how strains have genetically adapted to life in the absence of individual non-essential genes.

Global Analysis of Furfural-Induced Genomic Instability Using a Yeast Model

Furfural is an important renewable precursor for multiple commercial chemicals and fuels; a main inhibitor existing in cellulosic hydrolysate, which is used for bioethanol fermentation; and a potential carcinogen, as well. Using a genetic system in Saccharomyces cerevisiae that allows detection of crossover events, we observed that the frequency of mitotic recombination was elevated by 1.5- to 40-fold when cells were treated with 0.1 g/liter to 20 g/liter furfural. Analysis of the gene conversion tracts associated with crossover events suggested that most furfural-induced recombination resulted from repair of DNA double-strand breaks (DSBs) that occurred in the G₁ phase. Furfural was incapable of breaking DNA directly in vitro but could trigger DSBs in vivo related to reactive oxygen species accumulation. By whole-genome single nucleotide polymorphism (SNP) microarray and sequencing, furfural-induced genomic alterations that range from single base substitutions, loss of heterozygosity, and chromosomal rearrangements to aneuploidy were explored. At the whole-genome level, furfural-induced events were evenly distributed across 16 chromosomes but were enriched in high-GC-content regions. Point mutations, particularly the C-to-T/G-to-A transitions, were significantly elevated in furfural-treated cells compared to wild-type cells. This study provided multiple novel insights into the global effects of furfural on genomic stability.IMPORTANCE Whether and how furfural affects genome integrity have not been clarified. Using a Saccharomyces cerevisiae model, we found that furfural exposure leads to in vivo DSBs and elevation in mitotic recombination by orders of magnitude. Gross chromosomal rearrangements and aneuploidy events also occurred at a higher frequency in furfural-treated cells. In a genome-wide analysis, we show that the patterns of mitotic recombination and point mutations differed dramatically in furfural-treated cells and wild-type cells.

Reshuffling yeast chromosomes with CRISPR/Cas9

Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions.

Chromosomes are highly dynamic objects that often undergo large structural variations such as reciprocal translocations. Such rearrangements can have dramatic functional consequences, as they can disrupt genes, change their regulation or create novel fusion genes at their breakpoints. For instance, 90–95% of patients diagnosed with chronic myeloid leukemia carry the Philadelphia chromosome characterized by a reciprocal translocation between chromosomes 9 and 22. In addition, translocations reorganize the genetic information along chromosomes, which in turn can modify the 3D architecture of the genome and potentially affect its functioning. Quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. Here, we report a novel CRISPR/Cas9-based technology allowing to generate with high efficiency and at a base-pair precision either uniquely targeted or multiple reciprocal translocations in yeast, without leaving any marker or scar in the genome. Engineering targeted reciprocal translocations allowed us for the first time to untangle the phenotypic impacts of large chromosomal rearrangements from that of point mutations. In addition, the generation of multiple translocations led to a large reorganization of the genetic information along the chromosomes, often including unanticipated large segmental duplications. We showed that reshuffling the genome resulted in the emergence of fitness advantage in stressful environmental conditions, even in strains where no gene was disrupted or amplified by the translocations.

Allele-specific genome editing using CRISPR-Cas9 is associated with loss of heterozygosity in diploid yeast

Targeted DNA double-strand breaks (DSBs) with CRISPR–Cas9 have revolutionized genetic modification by enabling efficient genome editing in a broad range of eukaryotic systems. Accurate gene editing is possible with near-perfect efficiency in haploid or (predominantly) homozygous genomes. However, genomes exhibiting polyploidy and/or high degrees of heterozygosity are less amenable to genetic modification. Here, we report an up to 99-fold lower gene editing efficiency when editing individual heterozygous loci in the yeast genome. Moreover, Cas9-mediated introduction of a DSB resulted in large scale loss of heterozygosity affecting DNA regions up to 360 kb and up to 1700 heterozygous nucleotides, due to replacement of sequences on the targeted chromosome by corresponding sequences from its non-targeted homolog. The observed patterns of loss of heterozygosity were consistent with homology directed repair. The extent and frequency of loss of heterozygosity represent a novel mutagenic side-effect of Cas9-mediated genome editing, which would have to be taken into account in eukaryotic gene editing. In addition to contributing to the limited genetic amenability of heterozygous yeasts, Cas9-mediated loss of heterozygosity could be particularly deleterious for human gene therapy, as loss of heterozygous functional copies of anti-proliferative and pro-apoptotic genes is a known path to cancer.

Mating-System Evolution: All Roads Lead to Selfing

Fungi reproduce via both heterothallic outcrossing and homothallic selfing modes, and transitions between the two are common throughout the tree of life. A new study reports that the transition from heterothallism to homothallism is common and has repeatedly punctuated the evolutionary trajectory across a major lineage of the fungal kingdom.