Striated muscle-specific base editing enables correction of mutations causing dilated cardiomyopathy.

Dilated cardiomyopathy is the second most common cause for heart failure with no cure except a high-risk heart transplantation. Approximately 30% of patients harbor heritable mutations which are amenable to CRISPR-based gene therapy. However, challenges related to delivery of the editing complex and off-target concerns hamper the broad applicability of CRISPR agents in the heart. We employ a combination of the viral vector AAVMYO with superior targeting specificity of heart muscle tissue and CRISPR base editors to repair patient mutations in the cardiac splice factor Rbm20, which cause aggressive dilated cardiomyopathy. Using optimized conditions, we repair >70% of cardiomyocytes in two Rbm20 knock-in mouse models that we have generated to serve as an in vivo platform of our editing strategy. Treatment of juvenile mice restores the localization defect of RBM20 in 75% of cells and splicing of RBM20 targets including TTN. Three months after injection, cardiac dilation and ejection fraction reach wild-type levels. Single-nuclei RNA sequencing uncovers restoration of the transcriptional profile across all major cardiac cell types and whole-genome sequencing reveals no evidence for aberrant off-target editing. Our study highlights the potential of base editors combined with AAVMYO to achieve gene repair for treatment of hereditary cardiac diseases.

Deep phenotyping of two preclinical mouse models and a cohort of RBM20 mutation carriers reveals no sex-dependent disease severity in RBM20 cardiomyopathy.

RBM20 cardiomyopathy is an arrhythmogenic form of dilated cardiomyopathy caused by mutations in the splicing factor RBM20. A recent study found a more severe phenotype in male patients with RBM20 cardiomyopathy patients than in female patients. Here, we aim to determine sex differences in an animal model of RBM20 cardiomyopathy and investigate potential underlying mechanisms. In addition, we aim to determine sex and gender differences in clinical parameters in a novel RBM20 cardiomyopathy patient cohort. We characterized an Rbm20 knockout (KO) mouse model, and show that splicing of key RBM20 targets, cardiac function, and arrhythmia susceptibility do not differ between sexes. Next, we performed deep phenotyping of these mice, and show that male and female Rbm20-KO mice possess transcriptomic and phosphoproteomic differences. Hypothesizing that these differences may influence the heart's ability to compensate for stress, we exposed Rbm20-KO mice to acute catecholaminergic stimulation and again found no functional differences. We also replicate the lack of functional differences in a mouse model with the Rbm20-R636Q mutation. Lastly, we present a patient cohort of 33 RBM20 cardiomyopathy patients and show that these patients do not possess sex and gender differences in disease severity. Current mouse models of RBM20 cardiomyopathy show more pronounced changes in gene expression and phosphorylation of cardiac proteins in male mice, but no sex differences in cardiac morphology and function. Moreover, other than reported before, male RBM20 cardiomyopathy patients do not present with worse cardiac function in a patient cohort from Germany and the Netherlands.NEW & NOTEWORTHY Optimal management of the cardiac disease is increasingly personalized, partly because of differences in outcomes between sexes. RBM20 cardiomyopathy has been described to be more severe in male patients, and this carries the risk that male patients are more scrutinized in the clinic than female patients. Our findings do not support this observation and suggest that treatment should not differ between male and female RBM20 cardiomyopathy patients, but instead should focus on the underlying disease mechanism.

Transcriptional neighborhoods regulate transcript isoform lengths and expression levels.

Sequence features of genes and their flanking regulatory regions are determinants of RNA transcript isoform expression and have been used as context-independent plug-and-play modules in synthetic biology. However, genetic context-including the adjacent transcriptional environment-also influences transcript isoform expression levels and boundaries. We used synthetic yeast strains with stochastically repositioned genes to systematically disentangle the effects of sequence and context. Profiling 120 million full-length transcript molecules across 612 genomic perturbations, we observed sequence-independent alterations to gene expression levels and transcript isoform boundaries that were influenced by neighboring transcription. We identified features of transcriptional context that could predict these alterations and used these features to engineer a synthetic circuit where transcript length was controlled by neighboring transcription. This demonstrates how positional context can be leveraged in synthetic genome engineering.

Patient-derived gene and protein expression signatures of NGLY1 deficiency.

N-Glycanase 1 (NGLY1) deficiency is a rare and complex genetic disorder. Although recent studies have shed light on the molecular underpinnings of NGLY1 deficiency, a systematic characterization of gene and protein expression changes in patient-derived cells has been lacking. Here, we performed RNA-sequencing and mass spectrometry to determine the transcriptomes and proteomes of 66 cell lines representing four different cell types derived from 14 NGLY1 deficient patients and 17 controls. Although NGLY1 protein levels were up to 9.5-fold downregulated in patients compared with parents, residual and likely non-functional NGLY1 protein was detectable in all patient-derived lymphoblastoid cell lines. Consistent with the role of NGLY1 as a regulator of the transcription factor Nrf1, we observed a cell type-independent downregulation of proteasomal genes in NGLY1 deficient cells. In contrast, genes involved in ribosome biogenesis and mRNA processing were upregulated in multiple cell types. In addition, we observed cell type-specific effects. For example, genes and proteins involved in glutathione synthesis, such as the glutamate-cysteine ligase subunits GCLC and GCLM, were downregulated specifically in lymphoblastoid cells. We provide a web application that enables access to all results generated in this study at https://apps.embl.de/ngly1browser. This resource will guide future studies of NGLY1 deficiency in directions that are most relevant to patients.

Loss of N-Glycanase 1 Alters Transcriptional and Translational Regulation in K562 Cell Lines.

N-Glycanase 1 (NGLY1) deficiency is an ultra-rare, complex and devastating neuromuscular disease. Patients display multi-organ symptoms including developmental delays, movement disorders, seizures, constipation and lack of tear production. NGLY1 is a deglycosylating protein involved in the degradation of misfolded proteins retrotranslocated from the endoplasmic reticulum (ER). NGLY1-deficient cells have been reported to exhibit decreased deglycosylation activity and an increased sensitivity to proteasome inhibitors. We show that the loss of NGLY1 causes substantial changes in the RNA and protein landscape of K562 cells and results in downregulation of proteasomal subunits, consistent with its processing of the transcription factor NFE2L1. We employed the CMap database to predict compounds that can modulate NGLY1 activity. Utilizing our robust K562 screening system, we demonstrate that the compound NVP-BEZ235 (Dactosilib) promotes degradation of NGLY1-dependent substrates, concurrent with increased autophagic flux, suggesting that stimulating autophagy may assist in clearing aberrant substrates during NGLY1 deficiency.

Liver-specific deletion of Ngly1 causes abnormal nuclear morphology and lipid metabolism under food stress.

The cytoplasmic peptide:N-glycanase (Ngly1) is a de-N-glycosylating enzyme that cleaves N-glycans from misfolded glycoproteins and is involved in endoplasmic reticulum-associated degradation. The recent discovery of NGLY1-deficiency, which causes severe systemic symptoms, drew attention to the physiological function of Ngly1 in mammals. While several studies have been carried out to reveal the physiological necessity of Ngly1, the semi-lethal nature of Ngly1-deficient animals made it difficult to analyze its function in adults. In this study, we focus on the physiological function of Ngly1 in liver (hepatocyte)-specific Ngly1-deficient mice generated using the cre-loxP system. We found that hepatocyte-specific Ngly1-deficient mice showed abnormal hepatocyte nuclear size/morphology with aging but did not show other notable defects in unstressed conditions. This nuclear phenotype did not appear to be related to the function of the only gene currently reported to rescue Ngly1-deficient murine lethality so far, endo-ß-N-acetylglucosaminidase. We also found that under a high fructose diet induced stress, the hepatocyte-specific Ngly1-deletion resulted in liver transaminases elevation and increased lipid droplet accumulation. We showed that the processing and localization of the transcription factor, nuclear factor erythroid 2-like 1 (Nfe2l1), was impaired in the Ngly1-deficient hepatocytes. Therefore, Nfe2l1, at least partially, contributes to the phenotypes observed in hepatocyte-specific Ngly1-deficient mice. Our results indicate that Ngly1 plays important roles in the adult liver impacting nuclear morphology and lipid metabolism. Hepatocyte-specific Ngly1-deficient mice could thus serve as a valuable animal model for assessing in vivo efficacy of drugs and/or treatment for NGLY1-deficiency.

Biological plasticity rescues target activity in CRISPR knock outs.

Gene knock outs (KOs) are efficiently engineered through CRISPR-Cas9-induced frameshift mutations. While the efficiency of DNA editing is readily verified by DNA sequencing, a systematic understanding of the efficiency of protein elimination has been lacking. Here we devised an experimental strategy combining RNA sequencing and triple-stage mass spectrometry to characterize 193 genetically verified deletions targeting 136 distinct genes generated by CRISPR-induced frameshifts in HAP1 cells. We observed residual protein expression for about one third of the quantified targets, at variable levels from low to original, and identified two causal mechanisms, translation reinitiation leading to N-terminally truncated target proteins or skipping of the edited exon leading to protein isoforms with internal sequence deletions. Detailed analysis of three truncated targets, BRD4, DNMT1 and NGLY1, revealed partial preservation of protein function. Our results imply that systematic characterization of residual protein expression or function in CRISPR-Cas9-generated KO lines is necessary for phenotype interpretation.

Condensin controls cellular RNA levels through the accurate segregation of chromosomes instead of directly regulating transcription.

Condensins are genome organisers that shape chromosomes and promote their accurate transmission. Several studies have also implicated condensins in gene expression, although any mechanisms have remained enigmatic. Here, we report on the role of condensin in gene expression in fission and budding yeasts. In contrast to previous studies, we provide compelling evidence that condensin plays no direct role in the maintenance of the transcriptome, neither during interphase nor during mitosis. We further show that the changes in gene expression in post-mitotic fission yeast cells that result from condensin inactivation are largely a consequence of chromosome missegregation during anaphase, which notably depletes the RNA-exosome from daughter cells. Crucially, preventing karyotype abnormalities in daughter cells restores a normal transcriptome despite condensin inactivation. Thus, chromosome instability, rather than a direct role of condensin in the transcription process, changes gene expression. This knowledge challenges the concept of gene regulation by canonical condensin complexes.

INO80 represses osmostress induced gene expression by resetting promoter proximal nucleosomes.

The conserved INO80 chromatin remodeling complex is involved in regulation of DNA damage repair, replication and transcription. It is commonly recruited to the transcription start region and contributes to the establishment of promoter-proximal nucleosomes. We find a substantial influence of INO80 on nucleosome dynamics and gene expression during stress induced transcription. Transcription induced by osmotic stress leads to genome-wide remodeling of promoter proximal nucleosomes. INO80 function is required for timely return of evicted nucleosomes to the 5΄ end of induced genes. Reduced INO80 function in Arp8-deficient cells leads to correlated prolonged transcription and nucleosome eviction. INO80 and the related complex SWR1 regulate incorporation of the H2A.Z isoform at promoter proximal nucleosomes. However, H2A.Z seems not to influence osmotic stress induced gene regulation. Furthermore, we show that high rates of transcription promote INO80 recruitment to promoter regions, suggesting a connection between active transcription and promoter proximal nucleosome remodeling. In addition, we find that absence of INO80 enhances bidirectional promoter activity at highly induced genes and expression of a number of stress induced transcripts. We suggest that INO80 has a direct repressive role via promoter proximal nucleosome remodeling to limit high levels of transcription in yeast.

Loss of the Yeast SR Protein Npl3 Alters Gene Expression Due to Transcription Readthrough.

Yeast Npl3 is a highly abundant, nuclear-cytoplasmic shuttling, RNA-binding protein, related to metazoan SR proteins. Reported functions of Npl3 include transcription elongation, splicing and RNA 3' end processing. We used UV crosslinking and analysis of cDNA (CRAC) to map precise RNA binding sites, and strand-specific tiling arrays to look at the effects of loss of Npl3 on all transcripts across the genome. We found that Npl3 binds diverse RNA species, both coding and non-coding, at sites indicative of roles in both early pre-mRNA processing and 3' end formation. Tiling arrays and RNAPII mapping data revealed 3' extended RNAPII-transcribed RNAs in the absence of Npl3, suggesting that defects in pre-mRNA packaging events result in termination readthrough. Transcription readthrough was widespread and frequently resulted in down-regulation of neighboring genes. We conclude that the absence of Npl3 results in widespread 3' extension of transcripts with pervasive effects on gene expression.

Chromatin Dynamics and the RNA Exosome Function in Concert to Regulate Transcriptional Homeostasis.

The histone variant H2A.Z is a hallmark of nucleosomes flanking promoters of protein-coding genes and is often found in nucleosomes that carry lysine 56-acetylated histone H3 (H3-K56Ac), a mark that promotes replication-independent nucleosome turnover. Here, we find that H3-K56Ac promotes RNA polymerase II occupancy at many protein-coding and noncoding loci, yet neither H3-K56Ac nor H2A.Z has a significant impact on steady-state mRNA levels in yeast. Instead, broad effects of H3-K56Ac or H2A.Z on RNA levels are revealed only in the absence of the nuclear RNA exosome. H2A.Z is also necessary for the expression of divergent, promoter-proximal noncoding RNAs (ncRNAs) in mouse embryonic stem cells. Finally, we show that H2A.Z functions with H3-K56Ac to facilitate formation of chromosome interaction domains (CIDs). Our study suggests that H2A.Z and H3-K56Ac work in concert with the RNA exosome to control mRNA and ncRNA expression, perhaps in part by regulating higher-order chromatin structures.

Expression of nuclear and mitochondrial genes encoding ATP synthase is synchronized by disassembly of a multisynthetase complex.

In eukaryotic cells, oxidative phosphorylation involves multisubunit complexes of mixed genetic origin. Assembling these complexes requires an organelle-independent synchronizing system for the proper expression of nuclear and mitochondrial genes. Here we show that proper expression of the F1FO ATP synthase (complex V) depends on a cytosolic complex (AME) made of two aminoacyl-tRNA synthetases (cERS and cMRS) attached to an anchor protein, Arc1p. When yeast cells adapt to respiration the Snf1/4 glucose-sensing pathway inhibits ARC1 expression triggering simultaneous release of cERS and cMRS. Free cMRS and cERS relocate to the nucleus and mitochondria, respectively, to synchronize nuclear transcription and mitochondrial translation of ATP synthase genes. Strains releasing asynchronously the two aminoacyl-tRNA synthetases display aberrant expression of nuclear and mitochondrial genes encoding subunits of complex V resulting in severe defects of the oxidative phosphorylation mechanism. This work shows that the AME complex coordinates expression of enzymes that require intergenomic control.

High-density tiling microarray analysis of the full transcriptional activity of yeast.

Understanding the relationship between DNA sequence variation and phenotypic variation in complex or quantitative traits is one of the major challenges in modern biology. We are witnessing a deluge of DNA sequence information and association studies of genetic polymorphisms with phenotypes of interest in families and populations. In addition, it has become clear that large portions of eukaryotic genomes beyond protein-coding genes are transcribed, generating numerous noncoding RNA (ncRNA) molecules whose functions remain mostly unknown.DNA oligonucleotide microarrays constitute a powerful technology for studying the expression of genes in different organisms. The Saccharomyces cerevisiae tiling array presents a significant advance over previous array-based platforms. It has a high density of overlapping probes that start on average every 8 bp along each strand of the genome, enabling precise definition of transcript structure. Furthermore, the array includes probes specific for the polymorphic positions of another, distantly related yeast strain, allowing accurate measurement of allele-specific expression in a hybrid of the two strains. This technology thus allows high-resolution, quantitative, strand- and allele-specific measurements of transcription from a full eukaryotic genome. In this chapter, we describe the methods for extracting RNA, synthesizing first-strand cDNA, fragmenting, and labeling of samples for hybridization to the tiling array. Combining genome-wide information on variation in DNA sequence with variation in transcript structure and levels promises to increase our understanding of the genotype-to-phenotype relationship.

Alternative polyadenylation diversifies post-transcriptional regulation by selective RNA-protein interactions.

Recent research has uncovered extensive variability in the boundaries of transcript isoforms, yet the functional consequences of this variation remain largely unexplored. Here, we systematically discriminate between the molecular phenotypes of overlapping coding and non-coding transcriptional events from each genic locus using a novel genome-wide, nucleotide-resolution technique to quantify the half-lives of 3' transcript isoforms in yeast. Our results reveal widespread differences in stability among isoforms for hundreds of genes in a single condition, and that variation of even a single nucleotide in the 3' untranslated region (UTR) can affect transcript stability. While previous instances of negative associations between 3' UTR length and transcript stability have been reported, here, we find that shorter isoforms are not necessarily more stable. We demonstrate the role of RNA-protein interactions in conditioning isoform-specific stability, showing that PUF3 binds and destabilizes specific polyadenylation isoforms. Our findings indicate that although the functional elements of a gene are encoded in DNA sequence, the selective incorporation of these elements into RNA through transcript boundary variation allows a single gene to have diverse functional consequences.

Role of histone modifications and early termination in pervasive transcription and antisense-mediated gene silencing in yeast.

Most genomes, including yeast Saccharomyces cerevisiae, are pervasively transcribed producing numerous non-coding RNAs, many of which are unstable and eliminated by nuclear or cytoplasmic surveillance pathways. We previously showed that accumulation of PHO84 antisense RNA (asRNA), in cells lacking the nuclear exosome component Rrp6, is paralleled by repression of sense transcription in a process dependent on the Hda1 histone deacetylase (HDAC) and the H3K4 histone methyl transferase Set1. Here we investigate this process genome-wide and measure the whole transcriptome of various histone modification mutants in a Δrrp6 strain using tiling arrays. We confirm widespread occurrence of potentially antisense-dependent gene regulation and identify three functionally distinct classes of genes that accumulate asRNAs in the absence of Rrp6. These classes differ in whether the genes are silenced by the asRNA and whether the silencing is HDACs and histone methyl transferase-dependent. Among the distinguishing features of asRNAs with regulatory potential, we identify weak early termination by Nrd1/Nab3/Sen1, extension of the asRNA into the open reading frame promoter and dependence of the silencing capacity on Set1 and the HDACs Hda1 and Rpd3 particularly at promoters undergoing extensive chromatin remodelling. Finally, depending on the efficiency of Nrd1/Nab3/Sen1 early termination, asRNA levels are modulated and their capability of silencing is changed.

Set3 HDAC mediates effects of overlapping noncoding transcription on gene induction kinetics.

The Set3 histone deacetylase complex (Set3C) binds histone H3 dimethylated at lysine 4 (H3K4me2) to mediate deacetylation of histones in 5'-transcribed regions. To discern how Set3C affects gene expression, genome-wide transcription was analyzed in yeast undergoing a series of carbon source shifts. Deleting SET3 primarily caused changes during transition periods, as genes were induced or repressed. Surprisingly, a majority of Set3-affected genes are overlapped by noncoding RNA (ncRNA) transcription. Many Set3-repressed genes have H3K4me2 instead of me3 over promoter regions, due to either reduced H3K4me3 or ncRNA transcription from distal or antisense promoters. Set3C also represses internal cryptic promoters, but in different regions of genes than the Set2/Rpd3S pathway. Finally, Set3C stimulates some genes by repressing an overlapping antagonistic antisense transcript. These results show that overlapping noncoding transcription can fine-tune gene expression, not via the ncRNA but by depositing H3K4me2 to recruit the Set3C deacetylase.

Experimental relocation of the mitochondrial ATP9 gene to the nucleus reveals forces underlying mitochondrial genome evolution.

Only a few genes remain in the mitochondrial genome retained by every eukaryotic organism that carry out essential functions and are implicated in severe diseases. Experimentally relocating these few genes to the nucleus therefore has both therapeutic and evolutionary implications. Numerous unproductive attempts have been made to do so, with a total of only 5 successes across all organisms. We have taken a novel approach to relocating mitochondrial genes that utilizes naturally nuclear versions from other organisms. We demonstrate this approach on subunit 9/c of ATP synthase, successfully relocating this gene for the first time in any organism by expressing the ATP9 genes from Podospora anserina in Saccharomyces cerevisiae. This study substantiates the role of protein structure in mitochondrial gene transfer: expression of chimeric constructs reveals that the P. anserina proteins can be correctly imported into mitochondria due to reduced hydrophobicity of the first transmembrane segment. Nuclear expression of ATP9, while permitting almost fully functional oxidative phosphorylation, perturbs many cellular properties, including cellular morphology, and activates the heat shock response. Altogether, our study establishes a novel strategy for allotopic expression of mitochondrial genes, demonstrates the complex adaptations required to relocate ATP9, and indicates a reason that this gene was only transferred to the nucleus during the evolution of multicellular organisms.

Genome-wide transcriptome analysis in yeast using high-density tiling arrays.

In the last decade, it became clear that transcription goes far beyond that of protein-coding genes. Most RNA molecules are transcribed from intergenic regions or introns and exhibit much variability in size, expression level, secondary structure, and evolutionary conservation. While for several types of non-coding RNAs some cellular functions have been reported, like for micro-RNAs and small nucleolar RNAs, for most others no indications of function or regulation have so far been found. Therefore, the RNA population inside a cell is diverse and cryptic and, thus, demands powerful methods to study its composition, abundance, and structure. DNA oligonucleotide microarrays have proven to be of great utility to study transcription of genes in various organisms. Recently, due to advancement in microarray technology, tiling microarrays that extend transcription measurement to genomic regions beyond protein-coding genes were designed for several species. The Saccharomyces cerevisiae yeast tiling array contains overlapping probes across the full genomic sequence, with consecutive probes starting every 8 bp on average on each strand, enabling strand-specific measurement of transcription from a full eukaryotic genome. Here, we describe the methods used to extract yeast RNA, convert it into first-strand cDNA, fragment, and label it for hybridization to the tiling array. This protocol will enable researchers not only to study which genes are expressed and to what levels, but also to identify non-coding RNAs and to study the structure of transcripts including their untranslated regions, alternative start, stop, and processing sites. This information will allow understanding their roles inside cells.

Accumulation of noncoding RNA due to an RNase P defect in Saccharomyces cerevisiae.

Ribonuclease P (RNase P) is an essential endoribonuclease that catalyzes the cleavage of the 5' leader of pre-tRNAs. In addition, a growing number of non-tRNA substrates have been identified in various organisms. RNase P varies in composition, as bacterial RNase P contains a catalytic RNA core and one protein subunit, while eukaryotic nuclear RNase P retains the catalytic RNA but has at least nine protein subunits. The additional eukaryotic protein subunits most likely provide additional functionality to RNase P, with one possibility being additional RNA recognition capabilities. To investigate the possible range of additional RNase P substrates in vivo, a strand-specific, high-density microarray was used to analyze what RNA accumulates with a mutation in the catalytic RNA subunit of nuclear RNase P in Saccharomyces cerevisiae. A wide variety of noncoding RNAs were shown to accumulate, suggesting that nuclear RNase P participates in the turnover of normally unstable nuclear RNAs. In some cases, the accumulated noncoding RNAs were shown to be antisense to transcripts that commensurately decreased in abundance. Pre-mRNAs containing introns also accumulated broadly, consistent with either compromised splicing or failure to efficiently turn over pre-mRNAs that do not enter the splicing pathway. Taken together with the high complexity of the nuclear RNase P holoenzyme and its relatively nonspecific capacity to bind and cleave mixed sequence RNAs, these data suggest that nuclear RNase P facilitates turnover of nuclear RNAs in addition to its role in pre-tRNA biogenesis.

Antisense expression increases gene expression variability and locus interdependency.

Genome-wide transcription profiling has revealed extensive expression of non-coding RNAs antisense to genes, yet their functions, if any, remain to be understood. In this study, we perform a systematic analysis of sense-antisense expression in response to genetic and environmental changes in yeast. We find that antisense expression is associated with genes of larger expression variability. This is characterized by more 'switching off' at low levels of expression for genes with antisense compared to genes without, yet similar expression at maximal induction. By disrupting antisense transcription, we demonstrate that antisense expression confers an on-off switch on gene regulation for the SUR7 gene. Consistent with this, genes that must respond in a switch-like manner, such as stress-response and environment-specific genes, are enriched for antisense expression. In addition, our data provide evidence that antisense expression initiated from bidirectional promoters enables the spreading of regulatory signals from one locus to neighbouring genes. These results indicate a general regulatory effect of antisense expression on sense genes and emphasize the importance of antisense-initiating regions downstream of genes in models of gene regulation.