Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Replication stress affects the fidelity of nucleosome-mediated epigenetic inheritance

The fidelity of epigenetic inheritance or, the precision by which epigenetic information is passed along, is an essential parameter for measuring the effectiveness of the process. How the precision of the process is achieved or modulated, however, remains largely elusive. We have performed quantitative measurement of epigenetic fidelity, using position effect variegation (PEV) in Schizosaccharomyces pombe as readout, to explore whether replication perturbation affects nucleosome-mediated epigenetic inheritance. We show that replication stresses, due to either hydroxyurea treatment or various forms of genetic lesions of the replication machinery, reduce the inheritance accuracy of CENP-A/Cnp1 nucleosome positioning within centromere. Mechanistically, we demonstrate that excessive formation of single-stranded DNA, a common molecular abnormality under these conditions, might have correlation with the reduction in fidelity of centromeric chromatin duplication. Furthermore, we show that replication stress broadly changes chromatin structure at various loci in the genome, such as telomere heterochromatin expanding and mating type locus heterochromatin spreading out of the boundaries. Interestingly, the levels of inheritable expanding at sub-telomeric heterochromatin regions are highly variable among independent cell populations. Finally, we show that HU treatment of the multi-cellular organisms C. elegans and D. melanogaster affects epigenetically programmed development and PEV, illustrating the evolutionary conservation of the phenomenon. Replication stress, in addition to its demonstrated role in genetic instability, promotes variable epigenetic instability throughout the epigenome.

In this study, we found replication stresses reduce the fidelity of nucleosome-mediated epigenetic inheritance. Using Position Effect Variegation (PEV) in centromere as an indicator of chromatin epigenetic stability, we quantified the precision of nucleosomal inheritance and found replication stresses reduce the fidelity of nucleosome-mediated epigenetic inheritance. Further analysis of genome-wide heterochromatin distribution showed that replication stresses affect chromatin structure by expanding of heterochromatin with locus specificity. Mechanistically, we provide evidence suggesting that excessive formation of single-stranded DNA might have correlation with the reduction in fidelity of centromeric chromatin duplication. Finally, we demonstrated replication stress perturb the development process by reducing the fidelity of chromatin organization duplication in fruit fly and worm, illustrating the broadness and the evolutionary conservation of the phenomenon. Together, our results shed light on the importance of replication stresses cause epigenetic instability in addition to genetic stability.

upSET, the Drosophila homologue of SET3, Is Required for Viability and the Proper Balance of Active and Repressive Chromatin Marks

Chromatin plays a critical role in faithful implementation of gene expression programs. Different post-translational modifications (PTMs) of histone proteins reflect the underlying state of gene activity, and many chromatin proteins write, erase, bind, or are repelled by, these histone marks. One such protein is UpSET, the Drosophila homolog of yeast Set3 and mammalian KMT2E (MLL5). Here, we show that UpSET is necessary for the proper balance between active and repressed states. Using CRISPR/Cas-9 editing, we generated S2 cells that are mutant for upSET We found that loss of UpSET is tolerated in S2 cells, but that heterochromatin is misregulated, as evidenced by a strong decrease in H3K9me2 levels assessed by bulk histone PTM quantification. To test whether this finding was consistent in the whole organism, we deleted the upSET coding sequence using CRISPR/Cas-9, which we found to be lethal in both sexes in flies. We were able to rescue this lethality using a tagged upSET transgene, and found that UpSET protein localizes to transcriptional start sites (TSS) of active genes throughout the genome. Misregulated heterochromatin is apparent by suppressed position effect variegation of the wm4 allele in heterozygous upSET-deleted flies. Using nascent-RNA sequencing in the upSET-mutant S2 lines, we show that this result applies to heterochromatin genes generally. Our findings support a critical role for UpSET in maintaining heterochromatin, perhaps by delimiting the active chromatin environment.

[Organization of the Drosophila melanogaster SF1 insulator and its role in transcription regulation in transgenic lines]

The SF1 insulator was found to contain a polyadenylation signal, which corresponded to the functional polyadenylation signal in embryonic S2 cells and the transgenic lines of Drosophila and bi-directional promoter that functioned in S2 cells. The studies performed did not confirm the ability of the SF1 insulator to protect expression of reporter gene white from the chromosome position effect in transgenic lines.

Mutations in CG8878, a novel putative protein kinase, enhance P element dependent silencing (PDS) and position effect variegation (PEV) in Drosophila melanogaster

Genes in multicellular organisms are expressed as part of a developmental program that is largely dependent on self-perpetuating higher-order chromatin states. The mechanism of establishing and maintaining these epigenetic events is well studied in Drosophila. The first known example of an epigenetic effect was that of (PEV) in Drosophila, which has been shown to be due to gene silencing via heterochromatin formation. We are investigating a process similar to Position Effect Variegation (PEV) using a mini-w transgene, called Pci, inserted in the upstream regulatory region of ci. The mini-white+ transgene in Pci is expressed throughout the adult eye; however, when other P or KP elements are present, a variegated eye phenotype results indicating random w+ silencing during development. This P element dependent silencing (PDS) can be modified by the haplo-suppressors/triplo-enhancers, Su(var)205 and Su(var)3-7, indicating that these heterochromatic modifiers also act dose dependently in PDS. Here we use a spontaneous derivative mutation of Pci called PciE1 (E1) that variegates like PDS in the absence of P elements, presumably due to an adjacent gypsy element insertion, to screen for second-site modifier mutations that enhance variable silencing of white+ in E1. We isolated 7 mutations in CG8878, an essential gene, that enhance the E1 variegated phenotype. CG8878, a previously uncharacterized gene, potentially encodes a serine/threonine kinase whose closest Drosophila paralogue, ballchen (nhk-1), phosphorylates histones. These mutant alleles enhance both PDS at E1 and Position Effect Variegation (PEV) at w(m4), indicating a previously unknown common silencing mechanism between the two.

Relative transgene expression frequencies in homozygous versus hemizygous transgenic mice

We have used a simple binomial model of stochastic transgene inactivation at the level of the chromosome or transgene, rather than the cellular level, for the analysis of two mouse transgenic lines that show variegated patterns of expression. This predicts the percentages of cells that express one, both or neither alleles of the transgene in homozygotes from the observed percentages of cells, which express the transgene in hemizygotes. It adequately explained the relationship between the numbers of cells expressing the transgene in hemizygous and homozygous mosaic 21OH/LacZ mouse adrenals and mosaic BLG/7 mouse mammary glands. The binomial model also predicted that a small proportion of cells in mosaic mammary glands of BLG/7 homozygotes would express both BLG/7 alleles but published data indicated that all cells expressing the transgene showed monoallelic expression. Although it didn't fit all of the BLG/7 data as precisely as a more complex model, which used several ad hoc assumptions to explain these results, the simple binomial model was able to explain the relationship in observed transgene expression frequencies between hemizygous and homozygous mosaic tissues for both 21OH/LacZ and BLG/7 mice. It may prove to be a useful general model for analysing other transgenic animals showing mosaic transgene expression.

[Correlation on a cellular level of gene transcriptional silencing and heterochromatin compartment dragging in case of PEV-producing eu-heterochromatin rearrangement in Drosophila melanogaster]

Eu-heterochromatic rearrangements transfer genes into the heterochromatin and cause their variegated inactivation (PEV). Genes affected by PEV often demonstrate association with heterochromatic nuclear compartment (a distinct area composed of heterochromatin sequences like satellite DNA and enriched in specific chromatin proteins e.g. HP1). Here, we investigate the nuclear localization and the expression levels of the genes subjected to PEV caused by chromosome inversion, In(2)A4. We demonstrate that the degree of PEV-caused gene inactivation depends on a developmental stage, and the maximum of repression corresponds to the gene expression activation period. In the case of In(2)A4 rearrangement we detect the dragging of affected euchromatic region into heterochromatic nuclear compartment and the increase in HP1 occupancy in this region. We developed a protocol of simultaneous RNA-DNA-protein staining to demonstrate firstly in a single cell a strong correlation between transcriptional activity of affected gene and its distance from chromosome 2 satellite DNA.

Extensive phenotypic variation among allelic T-DNA inserts in Arabidopsis thaliana

T-DNA insertion mutants are a tool used widely in Arabidopsis thaliana to disrupt gene function. We phenotyped multiple homozygous T-DNA A. thaliana mutants at each of two loci (AT1G11060 and AT4G00210). We measured life history traits, including germination, size at reproduction and fruit production. Allelic T-DNA lines differed for most traits at AT1G11060 but not at AT4G00210. However, insertions in exons differed from other insertion positions in AT4G00210 but not in AT1G11060. We found evidence for additional insertions in approximately half of the lines, but found few phenotypic consequences. In general, our results suggest that a cautious interpretation of T-DNA phenotypes is warranted.

[Nuclear envelope attachment sites of interphase chromosomes: barrier elements but not insulators]

A neDNA fragment that protects a transgene from position effect variegation when flanking it was tested for insulator properties. The fragment did not act as an insulator. A similarity of neDNA and various barrier elements was examined, and Drosophila melanogaster and Mus musculus chromosome DNA regions homologous to neDNA were analyzed. Additional conserved sites and DNA duplex destabilization sites were found in the neDNA sequence, and DNA conformational specifics were assumed for the chromosomal neighborhood of neDNA sites of the D. melanogaster and M. musculus genomes.

Chromosomal position effects are linked to sir2-mediated variation in transcriptional burst size

Gene expression noise varies with genomic position and is a driving force in the evolution of chromosome organization. Nevertheless, position effects remain poorly characterized. Here, we present a systematic analysis of chromosomal position effects by characterizing single-cell gene expression from euchromatic positions spanning the length of a eukaryotic chromosome. We demonstrate that position affects gene expression by modulating the size of transcriptional bursts, rather than their frequency, and that the histone deacetylase Sir2 plays a role in this process across the chromosome.

Superhelical duplex destabilization and the recombination position effect

The susceptibility to recombination of a plasmid inserted into a chromosome varies with its genomic position. This recombination position effect is known to correlate with the average G+C content of the flanking sequences. Here we propose that this effect could be mediated by changes in the susceptibility to superhelical duplex destabilization that would occur. We use standard nonparametric statistical tests, regression analysis and principal component analysis to identify statistically significant differences in the destabilization profiles calculated for the plasmid in different contexts, and correlate the results with their measured recombination rates. We show that the flanking sequences significantly affect the free energy of denaturation at specific sites interior to the plasmid. These changes correlate well with experimentally measured variations of the recombination rates within the plasmid. This correlation of recombination rate with superhelical destabilization properties of the inserted plasmid DNA is stronger than that with average G+C content of the flanking sequences. This model suggests a possible mechanism by which flanking sequence base composition, which is not itself a context-dependent attribute, can affect recombination rates at positions within the plasmid.

Sub-telomeric core X and Y' elements in S. cerevisiae suppress extreme variations in gene silencing

Telomere Position Effect (TPE) is governed by strong repression signals emitted by telomeres via the Sir2/3/4 Histone Deacetylase complex. These signals are then relayed by weak proto-silencers residing in the subtelomeric core X and Y' elements. Subtelomeres also contain Sub-Telomeric Anti-silencing Regions (STARs). In this study we have prepared telomeres built of different combinations of core X, Y' and STARs and have analyzed them in strains lacking Histone-Acetyltransferase genes as well as in cdc6-1 and Δrif1 strains. We show that core X and Y' dramatically reduce both positive and negative variations in TPE, that are caused by these mutations. We also show that the deletion of Histone-Acetyltransferase genes reduce the silencing activity of an ACS proto-silencer, but also reduce the anti-silencing activity of a STAR. We postulate that core X and Y' act as epigenetic “cushioning” cis-elements.