Single-cell 3D genome structure reveals distinct human pluripotent states.

BACKGROUND: Pluripotent states of embryonic stem cells (ESCs) with distinct transcriptional profiles affect ESC differentiative capacity and therapeutic potential. Although single-cell RNA sequencing has revealed additional subpopulations and specific features of naive and primed human pluripotent stem cells (hPSCs), the underlying mechanisms that regulate their specific transcription and that control their pluripotent states remain elusive. RESULTS: By single-cell analysis of high-resolution, three-dimensional (3D) genomic structure, we herein demonstrate that remodeling of genomic structure is highly associated with the pluripotent states of human ESCs (hESCs). The naive pluripotent state is featured with specialized 3D genomic structures and clear chromatin compartmentalization that is distinct from the primed state. The naive pluripotent state is achieved by remodeling the active euchromatin compartment and reducing chromatin interactions at the nuclear center. This unique genomic organization is linked to enhanced chromatin accessibility on enhancers and elevated expression levels of naive pluripotent genes localized to this region. In contradistinction, the primed state exhibits intermingled genomic organization. Moreover, active euchromatin and primed pluripotent genes are distributed at the nuclear periphery, while repressive heterochromatin is densely concentrated at the nuclear center, reducing chromatin accessibility and the transcription of naive genes. CONCLUSIONS: Our data provide insights into the chromatin structure of ESCs in their naive and primed states, and we identify specific patterns of modifications in transcription and chromatin structure that might explain the genes that are differentially expressed between naive and primed hESCs. Thus, the inversion or relocation of heterochromatin to euchromatin via compartmentalization is related to the regulation of chromatin accessibility, thereby defining pluripotent states and cellular identity.

Heterochromatin Is Not the Only Place for satDNAs: The High Diversity of satDNAs in the Euchromatin of the Beetle Chrysolina americana (Coleoptera, Chrysomelidae).

The satellitome of the beetle Chrysolina americana Linneo, 1758 has been characterized through chromosomal analysis, genomic sequencing, and bioinformatics tools. C-banding reveals the presence of constitutive heterochromatin blocks enriched in A+T content, primarily located in pericentromeric regions. Furthermore, a comprehensive satellitome analysis unveils the extensive diversity of satellite DNA families within the genome of C. americana. Using fluorescence in situ hybridization techniques and the innovative CHRISMAPP approach, we precisely map the localization of satDNA families on assembled chromosomes, providing insights into their organization and distribution patterns. Among the 165 identified satDNA families, only three of them exhibit a remarkable amplification and accumulation, forming large blocks predominantly in pericentromeric regions. In contrast, the remaining, less abundant satDNA families are dispersed throughout euchromatic regions, challenging the traditional association of satDNA with heterochromatin. Overall, our findings underscore the complexity of repetitive DNA elements in the genome of C. americana and emphasize the need for further exploration to elucidate their functional significance and evolutionary implications.

HIV-1 Gag co-localizes with euchromatin histone marks at the nuclear periphery.

IMPORTANCE: The traditional view of retrovirus assembly posits that packaging of gRNA by HIV-1 Gag occurs in the cytoplasm or at the plasma membrane. However, our previous studies showing that HIV-1 Gag enters the nucleus and binds to USvRNA at transcription sites suggest that gRNA selection may occur in the nucleus. In the present study, we observed that HIV-1 Gag trafficked to the nucleus and co-localized with USvRNA within 8 hours of expression. In infected T cells (J-Lat 10.6) reactivated from latency and in a HeLa cell line stably expressing an inducible Rev-dependent HIV-1 construct, we found that Gag preferentially localized with euchromatin histone marks associated with enhancer and promoter regions near the nuclear periphery, which is the favored site HIV-1 integration. These observations support the innovative hypothesis that HIV-1 Gag associates with euchromatin-associated histones to localize to active transcription sites, promoting capture of newly synthesized gRNA for packaging.

The ancestral chromatin landscape of land plants.

Recent studies have shown that correlations between chromatin modifications and transcription vary among eukaryotes. This is the case for marked differences between the chromatin of the moss Physcomitrium patens and the liverwort Marchantia polymorpha. Mosses and liverworts diverged from hornworts, altogether forming the lineage of bryophytes that shared a common ancestor with land plants. We aimed to describe chromatin in hornworts to establish synapomorphies across bryophytes and approach a definition of the ancestral chromatin organization of land plants. We used genomic methods to define the 3D organization of chromatin and map the chromatin landscape of the model hornwort Anthoceros agrestis. We report that nearly half of the hornwort transposons were associated with facultative heterochromatin and euchromatin and formed the center of topologically associated domains delimited by protein coding genes. Transposons were scattered across autosomes, which contrasted with the dense compartments of constitutive heterochromatin surrounding the centromeres in flowering plants. Most of the features observed in hornworts are also present in liverworts or in mosses but are distinct from flowering plants. Hence, the ancestral genome of bryophytes was likely a patchwork of units of euchromatin interspersed within facultative and constitutive heterochromatin. We propose this genome organization was ancestral to land plants.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation.

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

A pangenome reference of 36 Chinese populations.

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference1. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

In skeletal muscle and neural crest cells, SMCHD1 regulates biological pathways relevant for Bosma syndrome and facioscapulohumeral dystrophy phenotype.

Many genetic syndromes are linked to mutations in genes encoding factors that guide chromatin organization. Among them, several distinct rare genetic diseases are linked to mutations in SMCHD1 that encodes the structural maintenance of chromosomes flexible hinge domain containing 1 chromatin-associated factor. In humans, its function as well as the impact of its mutations remains poorly defined. To fill this gap, we determined the episignature associated with heterozygous SMCHD1 variants in primary cells and cell lineages derived from induced pluripotent stem cells for Bosma arhinia and microphthalmia syndrome (BAMS) and type 2 facioscapulohumeral dystrophy (FSHD2). In human tissues, SMCHD1 regulates the distribution of methylated CpGs, H3K27 trimethylation and CTCF at repressed chromatin but also at euchromatin. Based on the exploration of tissues affected either in FSHD or in BAMS, i.e. skeletal muscle fibers and neural crest stem cells, respectively, our results emphasize multiple functions for SMCHD1, in chromatin compaction, chromatin insulation and gene regulation with variable targets or phenotypical outcomes. We concluded that in rare genetic diseases, SMCHD1 variants impact gene expression in two ways: (i) by changing the chromatin context at a number of euchromatin loci or (ii) by directly regulating some loci encoding master transcription factors required for cell fate determination and tissue differentiation.

Long interspersed nuclear element 1 and B1/Alu repeats blueprint genome compartmentalization.

The organization of the genome into euchromatin and heterochromatin has been known for almost 100 years [1]. More than 50% of mammalian genomes contain repetitive sequences [2,3]. Recently, a functional link between the genome and its folding has been identified [4,5]. Homotypic clustering of long interspersed nuclear element 1 (LINE1 or L1) and B1/Alu retrotransposons forms grossly exclusive nuclear domains that characterize and predict heterochromatin and euchromatin, respectively. The spatial segregation of L1 and B1/Alu-rich compartments is conserved in mammalian cells and can be rebuilt during the cell cycle and established de novo in early embryogenesis. Inhibition of L1 RNA drastically weakened homotypic repeat contacts and compartmental segregation, indicating that L1 plays a more significant role than just being a compartmental marker. This simple and inclusive genetic coding model of L1 and B1/Alu in shaping the macroscopic structure of the genome provides a plausible explanation for the remarkable conservation and robustness of its folding in mammalian cells. It also proposes a conserved core structure on which subsequent dynamic regulation takes place.

Euchromatin factors HULC and Set1C affect heterochromatin organization and mating-type switching in fission yeast Schizosaccharomyces pombe.

Mating-type (P or M) of fission yeast Schizosaccharomyces pombe is determined by the transcriptionally active mat1 cassette and is switched by gene conversion using a donor, either mat2 or mat3, located in an adjacent heterochromatin region (mating-type switching; MTS). In the switching process, heterochromatic donors of genetic information are selected based on the P or M cell type and on the action of two recombination enhancers, SRE2 promoting the use of mat2-P and SRE3 promoting the use of mat3-M, leading to replacement of the content of the expressed mat1 cassette. Recently, we found that the histone H3K4 methyltransferase complex Set1C participates in donor selection, raising the question of how a complex best known for its effects in euchromatin controls recombination in heterochromatin. Here, we report that the histone H2BK119 ubiquitin ligase complex HULC functions with Set1C in MTS, as mutants in the shf1, brl1, brl2 and rad6 genes showed defects similar to Set1C mutants and belonged to the same epistasis group as set1Δ. Moreover, using H3K4R and H2BK119R histone mutants and a Set1-Y897A catalytic mutant, we found that ubiquitylation of histone H2BK119 by HULC and methylation of histone H3K4 by Set1C are functionally coupled in MTS. Cell-type biases in MTS in these mutants suggested that HULC and Set1C inhibit the use of the SRE3 recombination enhancer in M cells, thus favoring SRE2 and mat2-P. Consistent with this, imbalanced switching in the mutants was traced to compromised association of the directionality factor Swi6 with the recombination enhancers in M cells. Based on their known effects at other chromosomal locations, we speculate that HULC and Set1C control nucleosome mobility and strand invasion near the SRE elements. In addition, we uncovered distinct effects of HULC and Set1C on histone H3K9 methylation and gene silencing, consistent with additional functions in the heterochromatic domain.

Phase Separation and Correlated Motions in Motorized Genome.

The human genome is arranged in the cell nucleus nonrandomly, and phase separation has been proposed as an important driving force for genome organization. However, the cell nucleus is an active system, and the contribution of nonequilibrium activities to phase separation and genome structure and dynamics remains to be explored. We simulated the genome using an energy function parametrized with chromosome conformation capture (Hi-C) data with the presence of active, nondirectional forces that break the detailed balance. We found that active forces that may arise from transcription and chromatin remodeling can dramatically impact the spatial localization of heterochromatin. When applied to euchromatin, active forces can drive heterochromatin to the nuclear envelope and compete with passive interactions among heterochromatin that tend to pull them in opposite directions. Furthermore, active forces induce long-range spatial correlations among genomic loci beyond single chromosome territories. We further showed that the impact of active forces could be understood from the effective temperature defined as the fluctuation-dissipation ratio. Our study suggests that nonequilibrium activities can significantly impact genome structure and dynamics, producing unexpected collective phenomena.

Human 5-lipoxygenase regulates transcription by association to euchromatin.

Human 5-lipoxygenase (5-LO) is the key enzyme of leukotriene biosynthesis, mostly expressed in leukocytes and thus a crucial component of the innate immune system. In this study, we show that 5-LO, besides its canonical function as an arachidonic acid metabolizing enzyme, is a regulator of gene expression associated with euchromatin. By Crispr-Cas9-mediated 5-LO knockout (KO) in MonoMac6 (MM6) cells and subsequent RNA-Seq analysis, we identified 5-LO regulated genes which could be clustered to immune/defense response, cell adhesion, transcription and growth/developmental processes. Analysis of differentially expressed genes identified cyclooxygenase-2 (COX2, PTGS2) and kynureninase (KYNU) as strongly regulated 5-LO target genes. 5-LO knockout affected MM6 cell adhesion and tryptophan metabolism via inhibition of the degradation of the immunoregulator kynurenine. By subsequent FAIRE-Seq and 5-LO ChIP-Seq analyses, we found an association of 5-LO with euchromatin, with prominent 5-LO binding to promoter regions in actively transcribed genes. By enrichment analysis of the ChIP-Seq results, we identified potential 5-LO interaction partners. Furthermore, 5-LO ChIP-Seq peaks resemble patterns of H3K27ac histone marks, suggesting that 5-LO recruitment mainly takes place at acetylated histones. In summary, we demonstrate a noncanonical function of 5-LO as transcriptional regulator in monocytic cells.

The chromatin remodeller ATRX facilitates diverse nuclear processes, in a stochastic manner, in both heterochromatin and euchromatin.

The chromatin remodeller ATRX interacts with the histone chaperone DAXX to deposit the histone variant H3.3 at sites of nucleosome turnover. ATRX is known to bind repetitive, heterochromatic regions of the genome including telomeres, ribosomal DNA and pericentric repeats, many of which are putative G-quadruplex forming sequences (PQS). At these sites ATRX plays an ancillary role in a wide range of nuclear processes facilitating replication, chromatin modification and transcription. Here, using an improved protocol for chromatin immunoprecipitation, we show that ATRX also binds active regulatory elements in euchromatin. Mutations in ATRX lead to perturbation of gene expression associated with a reduction in chromatin accessibility, histone modification, transcription factor binding and deposition of H3.3 at the sequences to which it normally binds. In erythroid cells where downregulation of α-globin expression is a hallmark of ATR-X syndrome, perturbation of chromatin accessibility and gene expression occurs in only a subset of cells. The stochastic nature of this process suggests that ATRX acts as a general facilitator of cell specific transcriptional and epigenetic programmes, both in heterochromatin and euchromatin.

Topoisomerase VI participates in an insulator-like function that prevents H3K9me2 spreading.

The organization of the genome into transcriptionally active and inactive chromatin domains requires well-delineated chromatin boundaries and insulator functions in order to maintain the identity of adjacent genomic loci with antagonistic chromatin marks and functionality. In plants that lack known chromatin insulators, the mechanisms that prevent heterochromatin spreading into euchromatin remain to be identified. Here, we show that DNA Topoisomerase VI participates in a chromatin boundary function that safeguards the expression of genes in euchromatin islands within silenced heterochromatin regions. While some transposable elements are reactivated in mutants of the Topoisomerase VI complex, genes insulated in euchromatin islands within heterochromatic regions of the Arabidopsis thaliana genome are specifically down-regulated. H3K9me2 levels consistently increase at euchromatin island loci and decrease at some transposable element loci. We further show that Topoisomerase VI physically interacts with S-adenosylmethionine synthase methionine adenosyl transferase 3 (MAT3), which is required for H3K9me2. A Topoisomerase VI defect affects MAT3 occupancy on heterochromatic elements and its exclusion from euchromatic islands, thereby providing a possible mechanistic explanation to the essential role of Topoisomerase VI in the delimitation of chromatin domains.

Comparative 3D genome architecture in vertebrates.

BACKGROUND: The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. RESULTS: This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. CONCLUSIONS: Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Nanodosimetric Calculations of Radiation-Induced DNA Damage in a New Nucleus Geometrical Model Based on the Isochore Theory.

Double-strand breaks (DSBs) in nuclear DNA represents radiation-induced damage that has been identified as particularly deleterious. Calculating this damage using Monte Carlo track structure modeling could be a suitable indicator to better assess and anticipate the side-effects of radiation therapy. However, as already demonstrated in previous work, the geometrical description of the nucleus and the DNA content used in the simulation significantly influence damage calculations. Therefore, in order to obtain accurate results, this geometry must be as realistic as possible. In this study, a new geometrical model of an endothelial cell nucleus and DNA distribution according to the isochore theory are presented and used in a Monte Carlo simulation chain based on the Geant4-DNA toolkit. In this theory, heterochromatin and euchromatin compaction are distributed along the genome according to five different families (L1, L2, H1, H2, and H3). Each of these families is associated with a different hetero/euchromatin rate related to its compaction level. In order to compare the results with those obtained using a previous nuclear geometry, simulations were performed for protons with linear energy transfers (LETs) of 4.29 keV/µm, 19.51 keV/µm, and 43.25 keV/µm. The organization of the chromatin fibers at different compaction levels linked to isochore families increased the DSB yield by 6-10%, and it allowed the most affected part of the genome to be identified. These new results indicate that the genome core is more radiosensitive than the genome desert, with a 3-8% increase in damage depending on the LET. This work highlights the importance of using realistic distributions of chromatin compaction levels to calculate radio-induced damage using Monte Carlo simulation methods.

H3K56 deacetylation and H2A.Z deposition are required for aberrant heterochromatin spreading.

Crucial mechanisms are required to restrict self-propagating heterochromatin spreading within defined boundaries and prevent euchromatic gene silencing. In the filamentous fungus Neurospora crassa, the JmjC domain protein DNA METHYLATION MODULATOR-1 (DMM-1) prevents aberrant spreading of heterochromatin, but the molecular details remain unknown. Here, we revealed that DMM-1 is highly enriched in a well-defined 5-kb heterochromatin domain upstream of the cat-3 gene, hereby called 5H-cat-3 domain, to constrain aberrant heterochromatin spreading. Interestingly, aberrant spreading of the 5H-cat-3 domain observed in the dmm-1KO strain is accompanied by robust deposition of histone variant H2A.Z, and deletion of H2A.Z abolishes aberrant spreading of the 5H-cat-3 domain into adjacent euchromatin. Furthermore, lysine 56 of histone H3 is deacetylated at the expanded heterochromatin regions, and mimicking H3K56 acetylation with an H3K56Q mutation effectively blocks H2A.Z-mediated aberrant spreading of the 5H-cat-3 domain. Importantly, genome-wide analyses demonstrated the general roles of H3K56 deacetylation and H2A.Z deposition in aberrant spreading of heterochromatin. Together, our results illustrate a previously unappreciated regulatory process that mediates aberrant heterochromatin spreading.

Comparison of structural variants in the whole genome sequences of two Medicago truncatula ecotypes: Jemalong A17 and R108.

BACKGROUND: Structural variants (SVs) constitute a large proportion of the genomic variation that results in phenotypic variation in plants. However, they are still a largely unexplored feature in most plant genomes. Here, we present the whole-genome landscape of SVs between two model legume Medicago truncatula ecotypes-Jemalong A17 and R108- that have been extensively used in various legume biology studies. RESULTS: To catalogue SVs, we first resolved the previously published R108 genome assembly (R108 v1.0) to chromosome-scale using 124 × Hi-C data, resulting in a high-quality genome assembly. The inter-chromosomal reciprocal translocations between chromosomes 4 and 8 were confirmed by performing syntenic analysis between the two genomes. Combined with the Hi-C data, it appears that these translocation events had a significant effect on chromatin organization. Using both whole-genome and short-read alignments, we identified the genomic landscape of SVs between the two genomes, some of which may account for several phenotypic differences, including their differential responses to aluminum toxicity and iron deficiency, and the development of different anthocyanin leaf markings. We also found extensive SVs within the nodule-specific cysteine-rich gene family which encodes antimicrobial peptides essential for terminal bacteroid differentiation during nitrogen-fixing symbiosis. CONCLUSIONS: Our results provide a near-complete R108 genome assembly and the first genomic landscape of SVs obtained by comparing two M. truncatula ecotypes. This may provide valuable genomic resources for the functional and molecular research of legume biology in the future.

Chromatin Velocity reveals epigenetic dynamics by single-cell profiling of heterochromatin and euchromatin.

Recent efforts have succeeded in surveying open chromatin at the single-cell level, but high-throughput, single-cell assessment of heterochromatin and its underlying genomic determinants remains challenging. We engineered a hybrid transposase including the chromodomain (CD) of the heterochromatin protein-1α (HP-1α), which is involved in heterochromatin assembly and maintenance through its binding to trimethylation of the lysine 9 on histone 3 (H3K9me3), and developed a single-cell method, single-cell genome and epigenome by transposases sequencing (scGET-seq), that, unlike single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), comprehensively probes both open and closed chromatin and concomitantly records the underlying genomic sequences. We tested scGET-seq in cancer-derived organoids and human-derived xenograft (PDX) models and identified genetic events and plasticity-driven mechanisms contributing to cancer drug resistance. Next, building upon the differential enrichment of closed and open chromatin, we devised a method, Chromatin Velocity, that identifies the trajectories of epigenetic modifications at the single-cell level. Chromatin Velocity uncovered paths of epigenetic reorganization during stem cell reprogramming and identified key transcription factors driving these developmental processes. scGET-seq reveals the dynamics of genomic and epigenetic landscapes underlying any cellular processes.

Autism risk gene POGZ promotes chromatin accessibility and expression of clustered synaptic genes.

Deleterious genetic variants in POGZ, which encodes the chromatin regulator Pogo Transposable Element with ZNF Domain protein, are strongly associated with autism spectrum disorder (ASD). Although it is a high-confidence ASD risk gene, the neurodevelopmental functions of POGZ remain unclear. Here we reveal the genomic binding of POGZ in the developing forebrain at euchromatic loci and gene regulatory elements (REs). We profile chromatin accessibility and gene expression in Pogz-/- mice and show that POGZ promotes the active chromatin state and transcription of clustered synaptic genes. We further demonstrate that POGZ forms a nuclear complex and co-occupies loci with ADNP, another high-confidence ASD risk gene, and provide evidence that POGZ regulates other neurodevelopmental disorder risk genes as well. Our results reveal a neurodevelopmental function of an ASD risk gene and identify molecular targets that may elucidate its function in ASD.

Targeted Inter-Homologs Recombination in Arabidopsis Euchromatin and Heterochromatin.

Homologous recombination (HR) typically occurs during meiosis between homologs, at a few unplanned locations along the chromosomes. In this study, we tested whether targeted recombination between homologous chromosomes can be achieved via Clustered Regulatory Interspaced Short Palindromic Repeat associated protein Cas9 (CRISPR-Cas9)-induced DNA double-strand break (DSB) repair in Arabidopsis thaliana. Our experimental system includes targets for DSB induction in euchromatic and heterochromatic genomic regions of hybrid F1 plants, in one or both parental chromosomes, using phenotypic and molecular markers to measure Non-Homologous End Joining and HR repair. We present a series of evidence showing that targeted DSBs can be repaired via HR using a homologous chromosome as the template in various chromatin contexts including in pericentric regions. Targeted crossover was rare, but gene conversion events were the most frequent outcome of HR and were found in both "hot and cold" regions. The length of the conversion tracts was variable, ranging from 5 to 7505 bp. In addition, a typical feature of these tracks was that they often were interrupted. Our findings pave the way for the use of targeted gene-conversion for precise breeding.