CDCA7-associated global aberrant DNA hypomethylation translates to localized, tissue-specific transcriptional responses.

Disruption of cell division cycle associated 7 (CDCA7) has been linked to aberrant DNA hypomethylation, but the impact of DNA methylation loss on transcription has not been investigated. Here, we show that CDCA7 is critical for maintaining global DNA methylation levels across multiple tissues in vivo. A pathogenic Cdca7 missense variant leads to the formation of large, aberrantly hypomethylated domains overlapping with the B genomic compartment but without affecting the deposition of H3K9 trimethylation (H3K9me3). CDCA7-associated aberrant DNA hypomethylation translated to localized, tissue-specific transcriptional dysregulation that affected large gene clusters. In the brain, we identify CDCA7 as a transcriptional repressor and epigenetic regulator of clustered protocadherin isoform choice. Increased protocadherin isoform expression frequency is accompanied by DNA methylation loss, gain of H3K4 trimethylation (H3K4me3), and increased binding of the transcriptional regulator CCCTC-binding factor (CTCF). Overall, our in vivo work identifies a key role for CDCA7 in safeguarding tissue-specific expression of gene clusters via the DNA methylation pathway.

PABPN1 loss-of-function causes APA-shift in oculopharyngeal muscular dystrophy.

Alternative polyadenylation (APA) at the 3' UTR of transcripts contributes to the cell transcriptome. APA is suppressed by the nuclear RNA-binding protein PABPN1. Aging-associated reduced PABPN1 levels in skeletal muscles lead to muscle wasting. Muscle weakness in oculopharyngeal muscular dystrophy (OPMD) is caused by short alanine expansion in PABPN1 exon1. The expanded PABPN1 forms nuclear aggregates, an OPMD hallmark. Whether the expanded PABPN1 affects APA and how it contributes to muscle pathology is unresolved. To investigate these questions, we developed a procedure including RNA library preparation and a simple pipeline calculating the APA-shift ratio as a readout for PABPN1 activity. Comparing APA-shift results to previously published PAS utilization and APA-shift results, we validated this procedure. The procedure was then applied on the OPMD cell model and on RNA from OPMD muscles. APA-shift was genome-wide in the mouse OPMD model, primarily affecting muscle transcripts. In OPMD individuals, APA-shift was enriched with muscle transcripts. In an OPMD cell model APA-shift was not significant. APA-shift correlated with reduced expression levels of a subset of PABPN1 isoforms, whereas the expression of the expanded PABPN1 did not correlate with APA-shift. PABPN1 activity is not affected by the expression of expanded PABPN1, but rather by reduced PABPN1 expression levels. In muscles, PABPN1 activity initially affects muscle transcripts. We suggest that muscle weakness in OPMD is caused by PABPN1 loss-of-function leading to APA-shift that primarily affects in muscle transcripts.

MMR vaccination induces trained immunity via functional and metabolic reprogramming of γδ T cells.

The measles, mumps, and rubella (MMR) vaccine protects against all-cause mortality in children, but the immunological mechanisms mediating these effects are poorly known. We systematically investigated whether MMR can induce long-term functional changes in innate immune cells, a process termed trained immunity, that could at least partially mediate this heterologous protection. In a randomized, placebo-controlled trial, 39 healthy adults received either the MMR vaccine or a placebo. Using single-cell RNA-Seq, we found that MMR caused transcriptomic changes in CD14+ monocytes and NK cells, but most profoundly in γδ T cells. Monocyte function was not altered by MMR vaccination. In contrast, the function of γδ T cells was markedly enhanced by MMR vaccination, with higher production of TNF and IFN-γ, as well as upregulation of cellular metabolic pathways. In conclusion, we describe a trained immunity program characterized by modulation of γδ T cell function induced by MMR vaccination.

Antisense oligonucleotide-mediated disruption of HTT caspase-6 cleavage site ameliorates the phenotype of YAC128 Huntington disease mice.

In Huntington disease, cellular toxicity is particularly caused by toxic protein fragments generated from the mutant huntingtin (HTT) protein. By modifying the HTT protein, we aim to reduce proteolytic cleavage and ameliorate the consequences of mutant HTT without lowering total HTT levels. To that end, we use an antisense oligonucleotide (AON) that targets HTT pre-mRNA and induces partial skipping of exon 12, which contains the critical caspase-6 cleavage site. Here, we show that AON-treatment can partially restore the phenotype of YAC128 mice, a mouse model expressing the full-length human HTT gene including 128 CAG-repeats. Wild-type and YAC128 mice were treated intracerebroventricularly with AON12.1, scrambled AON or vehicle starting at 6 months of age and followed up to 12 months of age, when MRI was performed and mice were sacrificed. AON12.1 treatment induced around 40% exon skip and protein modification. The phenotype on body weight and activity, but not rotarod, was restored by AON treatment. Genes differentially expressed in YAC128 striatum changed toward wild-type levels and striatal volume was preserved upon AON12.1 treatment. However, scrambled AON also showed a restorative effect on gene expression and appeared to generally increase brain volume.

DMD antisense oligonucleotide mediated exon skipping efficiency correlates with flanking intron retention time and target position within the exon.

Mutations in the DMD gene are causative for Duchenne muscular dystrophy (DMD). Antisense oligonucleotide (AON) mediated exon skipping to restore disrupted dystrophin reading frame is a therapeutic approach that allows production of a shorter but functional protein. As DMD causing mutations can affect most of the 79 exons encoding dystrophin, a wide variety of AONs are needed to treat the patient population. Design of AONs is largely guided by trial-and-error, and it is yet unclear what defines the skippability of an exon. Here, we use a library of phosphorodiamidate morpholino oligomer (PMOs) AONs of similar physical properties to test the skippability of a large number of DMD exons. The DMD transcript is non-sequentially spliced, meaning that certain introns are retained longer in the transcript than downstream introns. We tested whether the relative intron retention time has a significant effect on AON efficiency, and found that targeting an out-of-frame exon flanked at its 5'-end by an intron that is retained in the transcript longer ('slow' intron) leads to overall higher exon skipping efficiency than when the 5'-end flanking intron is 'fast'. Regardless of splicing speed of flanking introns, we find that positioning an AON closer to the 5'-end of the target exon leads to higher exon skipping efficiency opposed to targeting an exons 3'-end. The data enclosed herein can be of use to guide future target selection and preferential AON binding sites for both DMD and other disease amenable by exon skipping therapies.

STRAIGHT-IN enables high-throughput targeting of large DNA payloads in human pluripotent stem cells.

Inserting large DNA payloads (>10 kb) into specific genomic sites of mammalian cells remains challenging. Applications ranging from synthetic biology to evaluating the pathogenicity of disease-associated variants for precision medicine initiatives would greatly benefit from tools that facilitate this process. Here, we merge the strengths of different classes of site-specific recombinases and combine these with CRISPR-Cas9-mediated homologous recombination to develop a strategy for stringent site-specific replacement of genomic fragments at least 50 kb in size in human induced pluripotent stem cells (hiPSCs). We demonstrate the versatility of STRAIGHT-IN (serine and tyrosine recombinase-assisted integration of genes for high-throughput investigation) by (1) inserting various combinations of fluorescent reporters into hiPSCs to assess the excitation-contraction coupling cascade in derivative cardiomyocytes and (2) simultaneously targeting multiple variants associated with inherited cardiac arrhythmic disorders into a pool of hiPSCs. STRAIGHT-IN offers a precise approach to generate genetically matched panels of hiPSC lines efficiently and cost effectively.

Iron accumulation induces oxidative stress, while depressing inflammatory polarization in human iPSC-derived microglia.

Iron accumulation in microglia has been observed in Alzheimer's disease and other neurodegenerative disorders and is thought to contribute to disease progression through various mechanisms, including neuroinflammation. To study this interaction, we treated human induced pluripotent stem cell-derived microglia (iPSC-MG) with iron, in combination with inflammatory stimuli such as interferon gamma (IFN-γ) and amyloid ß. Both IFN-γ and iron treatment increased labile iron levels, but only iron treatment led to a consistent increase of ferritin levels, reflecting long-term iron storage. Therefore, in iPSC-MG, ferritin appeared to be regulated by iron revels rather than inflammation. Further investigation showed that while IFN-γ induced pro-inflammatory activation, iron treatment dampened both classic pro- and anti-inflammatory activation on a transcriptomic level. Notably, iron-loaded microglia showed strong upregulation of cellular stress response pathways, the NRF2 pathway, and other oxidative stress pathways. Functionally, iPSC-MG exhibited altered phagocytosis and impaired mitochondrial metabolism following iron treatment. Collectively, these data suggest that in MG, in contrast to current hypotheses, iron treatment does not result in pro-inflammatory activation, but rather dampens it and induces oxidative stress.

Adenine base editing of the DUX4 polyadenylation signal for targeted genetic therapy in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is caused by chromatin relaxation of the D4Z4 repeat resulting in misexpression of the D4Z4-encoded DUX4 gene in skeletal muscle. One of the key genetic requirements for the stable production of full-length DUX4 mRNA in skeletal muscle is a functional polyadenylation signal (ATTAAA) in exon three of DUX4 that is used in somatic cells. Base editors hold great promise to treat DNA lesions underlying genetic diseases through their ability to carry out specific and rapid nucleotide mutagenesis even in postmitotic cells such as skeletal muscle. In this study, we present a simple and straightforward strategy for mutagenesis of the somatic DUX4 polyadenylation signal by adenine base editing in immortalized myoblasts derived from independent FSHD-affected individuals. We show that mutating this critical cis-regulatory element results in downregulation of DUX4 mRNA and its direct transcriptional target genes. Our findings identify the somatic DUX4 polyadenylation signal as a therapeutic target and represent the first step toward clinical application of the CRISPR-Cas9 base editing platform for FSHD gene therapy.

Integrating Whole-Genome Sequencing in Clinical Genetics: A Novel Disruptive Structural Rearrangement Identified in the Dystrophin Gene (DMD).

While in most patients the identification of genetic alterations causing dystrophinopathies is a relatively straightforward task, a significant number require genomic and transcriptomic approaches that go beyond a routine diagnostic set-up. In this work, we present a Becker Muscular Dystrophy patient with elevated creatinine kinase levels, progressive muscle weakness, mild intellectual disability and a muscle biopsy showing dystrophic features and irregular dystrophin labelling. Routine molecular techniques (Southern-blot analysis, multiplex PCR, MLPA and genomic DNA sequencing) failed to detect a defect in the DMD gene. Muscle DMD transcript analysis (RT-PCR and cDNA-MLPA) showed the absence of exons 75 to 79, seen to be present at the genomic level. These results prompted the application of low-coverage linked-read whole-genome sequencing (WGS), revealing a possible rearrangement involving DMD intron 74 and a region located upstream of the PRDX4 gene. Breakpoint PCR and Sanger sequencing confirmed the presence of a ~8 Mb genomic inversion. Aberrant DMD transcripts were subsequently identified, some of which contained segments from the region upstream of PRDX4. Besides expanding the mutational spectrum of the disorder, this study reinforces the importance of transcript analysis in the diagnosis of dystrophinopathies and shows how WGS has a legitimate role in clinical laboratory genetics.

Detailed genetic and functional analysis of the hDMDdel52/mdx mouse model.

Duchenne muscular dystrophy (DMD) is a severe, progressive neuromuscular disorder caused by reading frame disrupting mutations in the DMD gene leading to absence of functional dystrophin. Antisense oligonucleotide (AON)-mediated exon skipping is a therapeutic approach aimed at restoring the reading frame at the pre-mRNA level, allowing the production of internally truncated partly functional dystrophin proteins. AONs work in a sequence specific manner, which warrants generating humanized mouse models for preclinical tests. To address this, we previously generated the hDMDdel52/mdx mouse model using transcription activator like effector nuclease (TALEN) technology. This model contains mutated murine and human DMD genes, and therefore lacks mouse and human dystrophin resulting in a dystrophic phenotype. It allows preclinical evaluation of AONs inducing the skipping of human DMD exons 51 and 53 and resulting in restoration of dystrophin synthesis. Here, we have further characterized this model genetically and functionally. We discovered that the hDMD and hDMDdel52 transgene is present twice per locus, in a tail-to-tail-orientation. Long-read sequencing revealed a partial deletion of exon 52 (first 25 bp), and a 2.3 kb inversion in intron 51 in both copies. These new findings on the genomic make-up of the hDMD and hDMDdel52 transgene do not affect exon 51 and/or 53 skipping, but do underline the need for extensive genetic analysis of mice generated with genome editing techniques to elucidate additional genetic changes that might have occurred. The hDMDdel52/mdx mice were also evaluated functionally using kinematic gait analysis. This revealed a clear and highly significant difference in overall gait between hDMDdel52/mdx mice and C57BL6/J controls. The motor deficit detected in the model confirms its suitability for preclinical testing of exon skipping AONs for human DMD at both the functional and molecular level.

Parental haplotype-specific single-cell transcriptomics reveal incomplete epigenetic reprogramming in human female germ cells.

In contrast to mouse, human female germ cells develop asynchronously. Germ cells transition to meiosis, erase genomic imprints, and reactivate the X chromosome. It is unknown if these events all appear asynchronously, and how they relate to each other. Here we combine exome sequencing of human fetal and maternal tissues with single-cell RNA-sequencing of five donors. We reconstruct full parental haplotypes and quantify changes in parental allele-specific expression, genome-wide. First we distinguish primordial germ cells (PGC), pre-meiotic, and meiotic transcriptional stages. Next we demonstrate that germ cells from various stages monoallelically express imprinted genes and confirm this by methylation patterns. Finally, we show that roughly 30% of the PGCs are still reactivating their inactive X chromosome and that this is related to transcriptional stage rather than fetal age. Altogether, we uncover the complexity and cell-to-cell heterogeneity of transcriptional and epigenetic remodeling in female human germ cells.

Immune stimuli shape the small non-coding transcriptome of extracellular vesicles released by dendritic cells.

The release and uptake of nano-sized extracellular vesicles (EV) is a highly conserved means of intercellular communication. The molecular composition of EV, and thereby their signaling function to target cells, is regulated by cellular activation and differentiation stimuli. EV are regarded as snapshots of cells and are, therefore, in the limelight as biomarkers for disease. Although research on EV-associated RNA has predominantly focused on microRNAs, the transcriptome of EV consists of multiple classes of small non-coding RNAs with potential gene-regulatory functions. It is not known whether environmental cues imposed on cells induce specific changes in a broad range of EV-associated RNA classes. Here, we investigated whether immune-activating or -suppressing stimuli imposed on primary dendritic cells affected the release of various small non-coding RNAs via EV. The small RNA transcriptomes of highly pure EV populations free from ribonucleoprotein particles were analyzed by RNA sequencing and RT-qPCR. Immune stimulus-specific changes were found in the miRNA, snoRNA, and Y-RNA content of EV from dendritic cells, whereas tRNA and snRNA levels were much less affected. Only part of the changes in EV-RNA content reflected changes in cellular RNA, which urges caution in interpreting EV as snapshots of cells. By comprehensive analysis of RNA obtained from highly purified EV, we demonstrate that multiple RNA classes contribute to genetic messages conveyed via EV. The identification of multiple RNA classes that display cell stimulation-dependent association with EV is the prelude to unraveling the function and biomarker potential of these EV-RNAs.

Full-length mRNA sequencing uncovers a widespread coupling between transcription initiation and mRNA processing.

BACKGROUND: The multifaceted control of gene expression requires tight coordination of regulatory mechanisms at transcriptional and post-transcriptional level. Here, we studied the interdependence of transcription initiation, splicing and polyadenylation events on single mRNA molecules by full-length mRNA sequencing. RESULTS: In MCF-7 breast cancer cells, we find 2700 genes with interdependent alternative transcription initiation, splicing and polyadenylation events, both in proximal and distant parts of mRNA molecules, including examples of coupling between transcription start sites and polyadenylation sites. The analysis of three human primary tissues (brain, heart and liver) reveals similar patterns of interdependency between transcription initiation and mRNA processing events. We predict thousands of novel open reading frames from full-length mRNA sequences and obtained evidence for their translation by shotgun proteomics. The mapping database rescues 358 previously unassigned peptides and improves the assignment of others. By recognizing sample-specific amino-acid changes and novel splicing patterns, full-length mRNA sequencing improves proteogenomics analysis of MCF-7 cells. CONCLUSIONS: Our findings demonstrate that our understanding of transcriptome complexity is far from complete and provides a basis to reveal largely unresolved mechanisms that coordinate transcription initiation and mRNA processing.

Searching for avidity by chemical ligation of combinatorially self-assembled DNA-encoded ligand libraries.

DNA encoded ligands are self-assembled into bivalent complexes and chemically ligated to link their identities. To demonstrate their potential as a combinatorial screening platform for avidity interactions, the optimal bivalent aptamer design (examplar ligands) for human alpha-thrombin is determined in a single round of selection and the DNA scaffold replaced with minimal impact on the final design.

Non-sequential and multi-step splicing of the dystrophin transcript.

The dystrophin protein encoding DMD gene is the longest human gene. The 2.2 Mb long human dystrophin transcript takes 16 hours to be transcribed and is co-transcriptionally spliced. It contains long introns (24 over 10kb long, 5 over 100kb long) and the heterogeneity in intron size makes it an ideal transcript to study different aspects of the human splicing process. Splicing is a complex process and much is unknown regarding the splicing of long introns in human genes. Here, we used ultra-deep transcript sequencing to characterize splicing of the dystrophin transcripts in 3 different human skeletal muscle cell lines, and explored the order of intron removal and multi-step splicing. Coverage and read pair analyses showed that around 40% of the introns were not always removed sequentially. Additionally, for the first time, we report that non-consecutive intron removal resulted in 3 or more joined exons which are flanked by unspliced introns and we defined these joined exons as an exon block. Lastly, computational and experimental data revealed that, for the majority of dystrophin introns, multistep splicing events are used to splice out a single intron. Overall, our data show for the first time in a human transcript, that multi-step intron removal is a general feature of mRNA splicing.

KeyGenes, a Tool to Probe Tissue Differentiation Using a Human Fetal Transcriptional Atlas.

Differentiated derivatives of human pluripotent stem cells in culture are generally phenotypically immature compared to their adult counterparts. Their identity is often difficult to determine with certainty because little is known about their human fetal equivalents in vivo. Cellular identity and signaling pathways directing differentiation are usually determined by extrapolating information from either human adult tissue or model organisms, assuming conservation with humans. To resolve this, we generated a collection of human fetal transcriptional profiles at different developmental stages. Moreover, we developed an algorithm, KeyGenes, which uses this dataset to quantify the extent to which next-generation sequencing or microarray data resemble specific cell or tissue types in the human fetus. Using KeyGenes combined with the human fetal atlas, we identified multiple cell and tissue samples unambiguously on a limited set of features. We thus provide a flexible and expandable platform to monitor and evaluate the efficiency of differentiation in vitro.

Huntington's disease biomarker progression profile identified by transcriptome sequencing in peripheral blood.

With several therapeutic approaches in development for Huntington's disease, there is a need for easily accessible biomarkers to monitor disease progression and therapy response. We performed next-generation sequencing-based transcriptome analysis of total RNA from peripheral blood of 91 mutation carriers (27 presymptomatic and, 64 symptomatic) and 33 controls. Transcriptome analysis by DeepSAGE identified 167 genes significantly associated with clinical total motor score in Huntington's disease patients. Relative to previous studies, this yielded novel genes and confirmed previously identified genes, such as H2AFY, an overlap in results that has proven difficult in the past. Pathway analysis showed enrichment of genes of the immune system and target genes of miRNAs, which are downregulated in Huntington's disease models. Using a highly parallelized microfluidics array chip (Fluidigm), we validated 12 of the top 20 significant genes in our discovery cohort and 7 in a second independent cohort. The five genes (PROK2, ZNF238, AQP9, CYSTM1 and ANXA3) that were validated independently in both cohorts present a candidate biomarker panel for stage determination and therapeutic readout in Huntington's disease. Finally we suggest a first empiric formula predicting total motor score from the expression levels of our biomarker panel. Our data support the view that peripheral blood is a useful source to identify biomarkers for Huntington's disease and monitor disease progression in future clinical trials.

Determining the quality and complexity of next-generation sequencing data without a reference genome.

We describe an open-source kPAL package that facilitates an alignment-free assessment of the quality and comparability of sequencing datasets by analyzing k-mer frequencies. We show that kPAL can detect technical artefacts such as high duplication rates, library chimeras, contamination and differences in library preparation protocols. kPAL also successfully captures the complexity and diversity of microbiomes and provides a powerful means to study changes in microbial communities. Together, these features make kPAL an attractive and broadly applicable tool to determine the quality and comparability of sequence libraries even in the absence of a reference sequence. kPAL is freely available at https://github.com/LUMC/kPAL webcite.

New insights into domestication of carrot from root transcriptome analyses.

BACKGROUND: Understanding the molecular basis of domestication can provide insights into the processes of rapid evolution and crop improvement. Here we demonstrated the processes of carrot domestication and identified genes under selection based on transcriptome analyses. RESULTS: The root transcriptomes of widely differing cultivated and wild carrots were sequenced. A method accounting for sequencing errors was introduced to optimize SNP (single nucleotide polymorphism) discovery. 11,369 SNPs were identified. Of these, 622 (out of 1000 tested SNPs) were validated and used to genotype a large set of cultivated carrot, wild carrot and other wild Daucus carota subspecies, primarily of European origin. Phylogenetic analysis indicated that eastern carrot may originate from Western Asia and western carrot may be selected from eastern carrot. Different wild D. carota subspecies may have contributed to the domestication of cultivated carrot. Genetic diversity was significantly reduced in western cultivars, probably through bottlenecks and selection. However, a high proportion of genetic diversity (more than 85% of the genetic diversity in wild populations) is currently retained in western cultivars. Model simulation indicated high and asymmetric gene flow from wild to cultivated carrots, spontaneously and/or by introgression breeding. Nevertheless, high genetic differentiation exists between cultivated and wild carrots (Fst = 0.295) showing the strong effects of selection. Expression patterns differed radically for some genes between cultivated and wild carrot roots which may be related to changes in root traits. The up-regulation of water-channel-protein gene expression in cultivars might be involved in changing water content and transport in roots. The activated expression of carotenoid-binding-protein genes in cultivars could be related to the high carotenoid accumulation in roots. The silencing of allergen-protein-like genes in cultivated carrot roots suggested strong human selection to reduce allergy. These results suggest that regulatory changes of gene expressions may have played a predominant role in domestication. CONCLUSIONS: Western carrots may originate from eastern carrots. The reduction in genetic diversity in western cultivars due to domestication bottleneck/selection may have been offset by introgression from wild carrot. Differential gene expression patterns between cultivated and wild carrot roots may be a signature of strong selection for favorable cultivation traits.

DeepSAGE reveals genetic variants associated with alternative polyadenylation and expression of coding and non-coding transcripts.

Many disease-associated variants affect gene expression levels (expression quantitative trait loci, eQTLs) and expression profiling using next generation sequencing (NGS) technology is a powerful way to detect these eQTLs. We analyzed 94 total blood samples from healthy volunteers with DeepSAGE to gain specific insight into how genetic variants affect the expression of genes and lengths of 3'-untranslated regions (3'-UTRs). We detected previously unknown cis-eQTL effects for GWAS hits in disease- and physiology-associated traits. Apart from cis-eQTLs that are typically easily identifiable using microarrays or RNA-sequencing, DeepSAGE also revealed many cis-eQTLs for antisense and other non-coding transcripts, often in genomic regions containing retrotransposon-derived elements. We also identified and confirmed SNPs that affect the usage of alternative polyadenylation sites, thereby potentially influencing the stability of messenger RNAs (mRNA). We then combined the power of RNA-sequencing with DeepSAGE by performing a meta-analysis of three datasets, leading to the identification of many more cis-eQTLs. Our results indicate that DeepSAGE data is useful for eQTL mapping of known and unknown transcripts, and for identifying SNPs that affect alternative polyadenylation. Because of the inherent differences between DeepSAGE and RNA-sequencing, our complementary, integrative approach leads to greater insight into the molecular consequences of many disease-associated variants.