A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage 1-3 . As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days 4,5 , the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2 , a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

Herpesviruses mimic zygotic genome activation to promote viral replication

DUX4 is a germline transcription factor and a master regulator of zygotic genome activation. During early embryogenesis, DUX4 is crucial for maternal to zygotic transition at the 2-8-cell stage in order to overcome silencing of genes and enable transcription from the zygotic genome. In adult somatic cells, DUX4 expression is silenced and its activation in adult muscle cells causes the genetic disorder Facioscapulohumeral Muscular Dystrophy (FSHD). Here we show that herpesviruses from alpha-, beta- and gamma-herpesvirus subfamilies as well as papillomaviruses actively induce DUX4 expression to promote viral transcription and replication. We demonstrate that HSV-1 immediate early proteins directly induce expression of DUX4 and its target genes including endogenous retroelements, which mimics zygotic genome activation. We further show that DUX4 directly binds to the viral genome and promotes viral transcription. DUX4 is functionally required for herpesvirus infection, since genetic depletion of DUX4 by CRISPR/Cas9 abrogates viral replication. Our results show that herpesviruses induce DUX4 expression and its downstream germline-specific genes and retroelements, thus mimicking an early embryonic-like transcriptional program that prevents epigenetic silencing of the viral genome and facilitates herpesviral gene expression.

The SPOC proteins DIDO3 and PHF3 co-regulate gene expression and neuronal differentiation

Transcription is regulated by a multitude of activators and repressors, which bind to the RNA polymerase II (Pol II) machinery and modulate its progression. Death-inducer obliterator 3 (DIDO3) and PHD finger protein 3 (PHF3) are paralogue proteins that regulate transcription elongation by docking onto phosphorylated serine-2 in the C-terminal domain (CTD) of Pol II through their SPOC domains. Here, we show that DIDO3 and PHF3 form a complex that bridges the Pol II elongation machinery with chromatin and RNA processing factors and tethers Pol II in a phase-separated microenvironment. Their SPOC domains and C-terminal intrinsically disordered regions are critical for transcription regulation. PHF3 and DIDO exert cooperative and antagonistic effects on the expression of neuronal genes and are both essential for neuronal differentiation. In the absence of PHF3, DIDO3 is upregulated as a compensatory mechanism. In addition to shared gene targets, DIDO specifically regulates genes required for lipid metabolism. Collectively, our work reveals multiple layers of gene expression regulation by the DIDO3 and PHF3 paralogues, which have specific, co-regulatory and redundant functions in transcription.

Transcriptional reprogramming by mutated IRF4 in lymphoma

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

The SPOC domain is a phosphoserine binding module that bridges transcription machinery with co- and post-transcriptional regulators

The heptad repeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulators as well as RNA capping, splicing and 3'end processing factors. The SPOC domain of PHF3 was recently identified as a CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP (also known as SPEN) and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP SPOC in complex with CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with writers and readers of m6A, the most abundant RNA modification. RBM15 positively regulates m6A levels and mRNA stability in a SPOC-dependent manner, while SHARP SPOC is essential for its localization to inactive X-chromosomes. Our findings suggest that the SPOC domain is a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

SARS-CoV-2 infection dynamics revealed by wastewater sequencing analysis and deconvolution

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater. While applying our tool on a wide range of different datasets (from different types of wastewater sample locations and sequenced with different methods), we show that PiGx SARS-CoV-2 can be used to identify new mutations and detect any emerging new lineages in a highly automated and scalable way. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen.

Cleidocervical muscle: a mini literature survey of a human muscle variation

Cleidocervical muscles (CCM) or levator claviculae muscles in humans can be found as supernumerary unilaterally or bilaterally on the neck attached proximally to the clavicle and distally to the transverse process of cervical vertebrae at various levels. Altogether 20 case reports from year 1994 till present including 25 subjects related to CCM were found and analysed where parameters such as cervical insertion level, clavicular insertion at the middle vs. lateral third, unilateral vs. bilateral presence of the muscle, study type, reported gender of the subjects were extracted. Our literature survey shows that the prevalence of CCM in male and female was equally presented in radiological studies whereas almost 3-fold higher prevalence of males was found in cadaver reports. Since body donor system worldwide is male dominant, a 1:1 proportion of male and female in radiological studies could show more reality-based distribution of this muscle. Nevertheless, the presentation of this muscle was found in over 90% of the case reports unilaterally with higher left sided dominance. Even though the attachment points of CCM varied from case to case, the proximal attachment was found slightly more frequent on the middle third of the clavicle whereas the distal insertion was present more often on the superior cervical vertebrae than the lower ones. With prevalence of CCM in the population around 2.0-2.5%, the clinical, radiological and surgical relevance of this variation has to be highlighted to avoid potential misleading diagnostics in the neck.

Identifying tumor cells at the single-cell level using machine learning

Tumors are complex tissues of cancerous cells surrounded by a heterogeneous cellular microenvironment with which they interact. Single-cell sequencing enables molecular characterization of single cells within the tumor. However, cell annotation—the assignment of cell type or cell state to each sequenced cell—is a challenge, especially identifying tumor cells within single-cell or spatial sequencing experiments. Here, we propose ikarus, a machine learning pipeline aimed at distinguishing tumor cells from normal cells at the single-cell level. We test ikarus on multiple single-cell datasets, showing that it achieves high sensitivity and specificity in multiple experimental contexts.

Rapid Inflammasome Activation Is Attenuated in Post-Myocardial Infarction Monocytes

Inflammasomes are crucial gatekeepers of the immune response, but their maladaptive activation associates with inflammatory pathologies. Besides canonical activation, monocytes can trigger non-transcriptional or rapid inflammasome activation that has not been well defined in the context of acute myocardial infarction (AMI). Rapid transcription-independent inflammasome activation induced by simultaneous TLR priming and triggering stimulus was measured by caspase-1 (CASP1) activity and interleukin release. Both classical and intermediate monocytes from healthy donors exhibited robust CASP1 activation, but only classical monocytes produced high mature interleukin-18 (IL18) release. We also recruited a limited number of coronary artery disease (CAD, n=31) and AMI (n=29) patients to evaluate their inflammasome function and expression profiles. Surprisingly, monocyte subpopulations isolated from blood collected during percutaneous coronary intervention (PCI) from AMI patients presented diminished CASP1 activity and abrogated IL18 release despite increased NLRP3 gene expression. This unexpected attenuated rapid inflammasome activation was accompanied by a significant increase of TNFAIP3 and IRAKM expression. Moreover, TNFAIP3 protein levels of circulating monocytes showed positive correlation with high sensitive troponin T (hsTnT), implying an association between TNFAIP3 upregulation and the severity of tissue injury. We suggest this monocyte attenuation to be a protective phenotype aftermath following a very early inflammatory wave in the ischemic area. Damage-associated molecular patterns (DAMPs) or other signals trigger a transitory negative feedback loop within newly recruited circulating monocytes as a mechanism to reduce post-injury tissue damage.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Cell-type specialization is encoded by specific chromatin topologies

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

Lymphocyte access to lymphoma is impaired by high endothelial venule regression

Blood endothelial cells display remarkable plasticity depending on the demands of a malignant microenvironment. While studies in solid tumors focus on their role in metabolic adaptations, formation of high endothelial venules (HEVs) in lymph nodes extends their role to the organization of immune cell interactions. As a response to lymphoma growth, blood vessel density increases; however, the fate of HEVs remains elusive. Here, we report that lymphoma causes severe HEV regression in mouse models that phenocopies aggressive human B cell lymphomas. HEV dedifferentiation occurrs as a consequence of a disrupted lymph-carrying conduit system. Mechanosensitive fibroblastic reticular cells then deregulate CCL21 migration paths, followed by deterioration of dendritic cell proximity to HEVs. Loss of this crosstalk deprives HEVs of lymphotoxin-β-receptor (LTβR) signaling, which is indispensable for their differentiation and lymphocyte transmigration. Collectively, this study reveals a remodeling cascade of the lymph node microenvironment that is detrimental for immune cell trafficking in lymphoma.

PHF3 regulates neuronal gene expression through the Pol II CTD reader domain SPOC

The C-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) is a regulatory hub for transcription and RNA processing. Here, we identify PHD-finger protein 3 (PHF3) as a regulator of transcription and mRNA stability that docks onto Pol II CTD through its SPOC domain. We characterize SPOC as a CTD reader domain that preferentially binds two phosphorylated Serine-2 marks in adjacent CTD repeats. PHF3 drives liquid-liquid phase separation of phosphorylated Pol II, colocalizes with Pol II clusters and tracks with Pol II across the length of genes. PHF3 knock-out or SPOC deletion in human cells results in increased Pol II stalling, reduced elongation rate and an increase in mRNA stability, with marked derepression of neuronal genes. Key neuronal genes are aberrantly expressed in Phf3 knock-out mouse embryonic stem cells, resulting in impaired neuronal differentiation. Our data suggest that PHF3 acts as a prominent effector of neuronal gene regulation by bridging transcription with mRNA decay.

Abstract 196: Redefining cancer subtypes using multi-omics and deep learning

Cancer is a heterogeneous collection of diseases traditionally classified by the tissue of origin. The diversity of the molecular profiles of cancers has a big impact on the way patients are diagnosed and treated, how they respond to their prescribed treatments, the duration of survival after diagnosis, and factors such as remission, recurrence, or spread (metastasis) of the disease. While such diagnostic and prognostic outcomes are potentially predictable by taking a closer look into the changes of the genome, epigenome, transcriptome, proteome, and various other omics platforms, the contemporary cancer treatments still predominantly don't make the best use of such multi-omics profiling of patient samples. Therefore, multi-omics profiling of cancers holds great potential to define a molecularly coherent subtype definition of cancers in order to achieve the eventual goal of matching the best possible treatment to the subgroup of patients.

However, the current subtypes from consortiums such as TCGA have been defined by heterogeneous methods and molecular markers by different teams. A subset of these studies have not attempted to characterize molecular subtypes, but rather taken histopathologically defined subtypes as the gold standard and tried to characterize molecular features of these subtypes. Here we evaluate TCGA cancer subtypes based on the molecular profile coherence score. This novel metric combines survival statistics, pathways information, tumor purity estimates, and mutational signatures. We expect that subtypes that are patient subgroups should display molecular signature homogeneity. We evaluate TCGA subtypes from 21 cancers using these criteria and compare the subtypes with our own definition using multi-omics data in a deep learning framework. We have refined the several subtypes from multiple cancers towards more molecularly coherent patient subgroups.

Citation Format: Altuna Akalin, Bora Uyar, Jonathan Ronen, Vedran Franke. Redefining cancer subtypes using multi-omics and deep learning [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 196.

Transcriptomic profiling of SARS-CoV-2 infected human cell lines identifies HSP90 as target for COVID-19 therapy

Detailed knowledge of the molecular biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is crucial for understanding of viral replication, host responses, and disease progression. Here, we report gene expression profiles of three SARS-CoV- and SARS-CoV-2-infected human cell lines. SARS-CoV-2 elicited an approximately two-fold higher stimulation of the innate immune response compared to SARS-CoV in the human epithelial cell line Calu-3, including induction of miRNA-155. Single-cell RNA sequencing of infected cells showed that genes induced by virus infections were broadly upregulated, whereas interferon beta/lambda genes, a pro-inflammatory cytokines such as IL-6, were expressed only in small subsets of infected cells. Temporal analysis suggested that transcriptional activities of interferon regulatory factors precede those of nuclear factor κB. Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 activity resulted in a reduction of viral replication and pro-inflammatory cytokine expression in primary human airway epithelial cells.

Single-nucleus transcriptomics reveals functional compartmentalization in syncytial skeletal muscle cells

Syncytial skeletal muscle cells contain hundreds of nuclei in a shared cytoplasm. We investigated nuclear heterogeneity and transcriptional dynamics in the uninjured and regenerating muscle using single-nucleus RNA-sequencing (snRNAseq) of isolated nuclei from muscle fibers. This revealed distinct nuclear subtypes unrelated to fiber type diversity, previously unknown subtypes as well as the expected ones at the neuromuscular and myotendinous junctions. In fibers of the Mdx dystrophy mouse model, distinct subtypes emerged, among them nuclei expressing a repair signature that were also abundant in the muscle of dystrophy patients, and a nuclear population associated with necrotic fibers. Finally, modifications of our approach revealed the compartmentalization in the rare and specialized muscle spindle. Our data identifies nuclear compartments of the myofiber and defines a molecular roadmap for their functional analyses; the data can be freely explored on the MyoExplorer server ( https://shiny.mdc-berlin.de/MyoExplorer/ ).

Non-coding RNA transcriptome analysis of circulating human monocytes reveals subpopulation-enriched transcripts and selective isoforms differentially expressed during acute myocardial infarction

Background

Acute inflammation leads to increased myelopoiesis in the bone marrow and recruitment of immune cells to respective sites of injury. Yet, in certain cases of high demand such as acute tissue damage as in acute myocardial infarction (AMI), innate immune cells can be additionally recruited from lymphoid organs such as the spleen. Monocytes are important in both onset and resolution of inflammation comprising a huge repertoire of coding and non-coding RNA molecules regulating a fast-acting immune response. Long non-coding RNAs (lncRNA) emerge as novel regulatory cell-specific molecules, which are able to influence cell physiology in numerous ways.

Purpose

Delineate monocyte subpopulations from patients with AMI compared to healthy controls based on their non-coding transcriptome, different isoform usage and associated source of origin.

Methods and results

Next generation bulk RNA sequencing has been performed from human monocytes after pre-sorting for CD14 expression and subsequent separation into three major subpopulations based on their surface markers CD14 and CD16 using a FACS Aria II cell sorter. Non-coding RNA expression shows clear differences between cells from healthy controls versus patients with AMI and among subpopulations. Intensive bioinformatics analysis revealed numerous annotated lncRNA, antisense, pseudogene and circular RNAs displaying significant differences in their expression profiles. Some lncRNA show positive or negative correlations with their respective coding antisense gene or their convergent or divergent sense neighbouring gene, such as linc02207 vs MCTP2, SAP30-AS vs SAP30, BID-AS vs BID, HLX-AS vs HLX or TNFAIP3-AS vs TNFAIP3. We also identified the existence of different isoforms of certain lncRNA and a distinct preference thereof, such as linc02207- 217 isoform being most abundant.

Real-time PCR was applied for validation of lncRNA expression profile and their nucleocytoplasmic distribution observing preferential nuclear expression for most lncRNA. In vitro assays for silencing and overexpression approaches as well as treatment with inflammatory stimuli and in silico bioinformatics analysis will help to unravel their functionality. Additionally mapping our bulk RNA data to deconvoluted single cell RNAseq data from publicly available human cell atlas (HCA), is applied to additionally unravel the origin of particular subpopulations within the AMI monocytes, associating their origin to either spleen or bone marrow.

Conclusions

High read depth bulk RNA sequencing revealed a human monocyte-specific long non-coding RNA transcriptome differentially expressed in three predefined subpopulations and disease stages. Distinct expression of lncRNA and associated coding genes as well as distinct isoform preference implies specific functionalities. Studying these molecules and their interaction partners will unravel functional participants in the development of therapeutic approaches for cardiovascular diseases.

Funding Acknowledgement

Type of funding source: Public grant(s) – National budget only. Main funding source(s): DZHK, SNF

HOT or not: examining the basis of high-occupancy target regions

High-occupancy target (HOT) regions are segments of the genome with unusually high number of transcription factor binding sites. These regions are observed in multiple species and thought to have biological importance due to high transcription factor occupancy. Furthermore, they coincide with house-keeping gene promoters and consequently associated genes are stably expressed across multiple cell types. Despite these features, HOT regions are solely defined using ChIP-seq experiments and shown to lack canonical motifs for transcription factors that are thought to be bound there. Although, ChIP-seq experiments are the golden standard for finding genome-wide binding sites of a protein, they are not noise free. Here, we show that HOT regions are likely to be ChIP-seq artifacts and they are similar to previously proposed ‘hyper-ChIPable’ regions. Using ChIP-seq data sets for knocked-out transcription factors, we demonstrate presence of false positive signals on HOT regions. We observe sequence characteristics and genomic features that are discriminatory of HOT regions, such as GC/CpG-rich k-mers, enrichment of RNA–DNA hybrids (R-loops) and DNA tertiary structures (G-quadruplex DNA). The artificial ChIP-seq enrichment on HOT regions could be associated to these discriminatory features. Furthermore, we propose strategies to deal with such artifacts for the future ChIP-seq studies.

Single-cell RNA-sequencing of herpes simplex virus 1-infected cells connects NRF2 activation to an antiviral program

Herpesvirus infection initiates a range of perturbations in the host cell, which remain poorly understood at the level of individual cells. Here, we quantify the transcriptome of single human primary fibroblasts during the first hours of lytic infection with HSV-1. By applying a generalizable analysis scheme, we define a precise temporal order of early viral gene expression and propose a set-wise emergence of viral genes. We identify host cell genes and pathways relevant for infection by combining three different computational approaches: gene and pathway overdispersion analysis, prediction of cell-state transition probabilities, as well as future cell states. One transcriptional program, which correlates with increased resistance to infection, implicates the transcription factor NRF2. Consequently, Bardoxolone methyl and Sulforaphane, two known NRF2 agonists, impair virus production, suggesting that NRF2 activation restricts viral infection. Our study provides insights into early stages of HSV-1 infection and serves as a general blueprint for the investigation of heterogeneous cell states in virus infection.

Alternative 3' UTRs direct localization of functionally diverse protein isoforms in neuronal compartments

The proper subcellular localization of RNAs and local translational regulation is crucial in highly compartmentalized cells, such as neurons. RNA localization is mediated by specific cis-regulatory elements usually found in mRNA 3′UTRs. Therefore, processes that generate alternative 3′UTRs—alternative splicing and polyadenylation—have the potential to diversify mRNA localization patterns in neurons. Here, we performed mapping of alternative 3′UTRs in neurites and soma isolated from mESC-derived neurons. Our analysis identified 593 genes with differentially localized 3′UTR isoforms. In particular, we have shown that two isoforms of Cdc42 gene with distinct functions in neuronal polarity are differentially localized between neurites and soma of mESC-derived and mouse primary cortical neurons, at both mRNA and protein level. Using reporter assays and 3′UTR swapping experiments, we have identified the role of alternative 3′UTRs and mRNA transport in differential localization of alternative CDC42 protein isoforms. Moreover, we used SILAC to identify isoform-specific Cdc42 3′UTR-bound proteome with potential role in Cdc42 localization and translation. Our analysis points to usage of alternative 3′UTR isoforms as a novel mechanism to provide for differential localization of functionally diverse alternative protein isoforms.

P5394Whole transcriptome sequencing to define monocyte subpopulation-specific lncRNA profiles in acute myocardial infarction patients and healthy controls

Background

Monocytes are important immune cells in both onset and resolution of inflammation during pathologies such as acute myocardial infarction (AMI) and atherosclerosis. Long non-coding RNAs (lncRNAs) have emerged as novel regulatory and highly cell-specific molecules that can modulate cell physiology in numerous ways such as mRNA (de-)stabilization, micro RNA sponging or scaffolding of RNA binding proteins.

Purpose

Define a complete human monocyte subpopulation-specific transcriptome of long non-coding RNAs and characterize the difference in the profile of these RNA molecules in AMI patients.

Methods and results

Human monocyte subpopulations, defined as classical, intermediate and non-classical based on the expression of the surface markers CD14 and CD16, were collected on a FACS Aria II. Ribosomal-depleted cDNA libraries generated from total RNA were processed for Next Generation Sequencing on a HiSeq Illumina 2000. Computationally intensive bioinformatics revealedannotated lncRNAs, antisense, pseudogene and circular RNAs with significant difference in their expression profiles within subpopulations of healthy donors such as MEG3 or TERC, potential role players in cardiovascular disease. Our data also unraveled novel non-annotated ncRNAs not yet reported to reference databases, which are expected to be monocyte-specific. We applied certain criteria to identify potential candidate molecules such as annotation with existing Ensembl ID and a pre-determined expression level. Thereupon we selected differentially regulated long non-coding RNAs differentially expressed in cardiovascular disease and discovered 18 annotated potential lncRNAs dysregulated in classical monocytes of AMI patients such as HLX antisense, which might be involved in monocyte differentiation. Additional 5 targets appeared specific only for differences in intermediate and 3 with additional specific differences only in non-classical monocytes.Real-time PCR was applied for validation of long non-coding linear and circular RNAs differential expression and also to determine their nucleocytoplasmic distribution. We observed preferential nuclear expression for most lncRNAs in contrast to cytoplasmic circRNAs. In vitro assays for silencing and overexpressing certain target molecules as well astreatment withinflammatory stimuli and in silicoanalysis with different bioinformatics tools such as FANTOM and UCSC browser will help to unravel their functionality.

Conclusions

Next generation sequencing allowed us to define a human monocyte subpopulation-specific long non-coding transcriptome that presented significant differences in both lncRNA and circRNAs within monocyte subpopulations of healthy subjects and AMI patients. Studying functional mechanisms of identified lncRNAs and their interaction with the coding genome will help to unravel novel regulatory means of monocytes in acute myocardial infarction providing new opportunities for therapeutic approaches.

Acknowledgement/Funding

DZHK (German Centre for Cardiovascular Research), partner site Berlin, Germany;Berlin Institute of Health (BIH);Swiss National Science Foundation (SNF)

Surgical removal of the index node marked using magnetic seed localization to assess response to neoadjuvant immunotherapy in patients with stage III melanoma

This pilot study explored the value of localized index node removal after neoadjuvant immunotherapy in patients with stage III melanoma, for use as a response indicator to guide the extent of completion lymph node dissection. Promising technology.

PiGx: reproducible genomics analysis pipelines with GNU Guix

In bioinformatics, as well as other computationally intensive research fields, there is a need for workflows that can reliably produce consistent output, from known sources, independent of the software environment or configuration settings of the machine on which they are executed. Indeed, this is essential for controlled comparison between different observations and for the wider dissemination of workflows. However, providing this type of reproducibility and traceability is often complicated by the need to accommodate the myriad dependencies included in a larger body of software, each of which generally comes in various versions. Moreover, in many fields (bioinformatics being a prime example), these versions are subject to continual change due to rapidly evolving technologies, further complicating problems related to reproducibility. Here, we propose a principled approach for building analysis pipelines and managing their dependencies with GNU Guix. As a case study to demonstrate the utility of our approach, we present a set of highly reproducible pipelines called PiGx for the analysis of RNA sequencing, chromatin immunoprecipitation sequencing, bisulfite-treated DNA sequencing, and single-cell resolution RNA sequencing. All pipelines process raw experimental data and generate reports containing publication-ready plots and figures, with interactive report elements and standard observables. Users may install these highly reproducible packages and apply them to their own datasets without any special computational expertise beyond the use of the command line. We hope such a toolkit will provide immediate benefit to laboratory workers wishing to process their own datasets or bioinformaticians seeking to automate all, or parts of, their analyses. In the long term, we hope our approach to reproducibility will serve as a blueprint for reproducible workflows in other areas. Our pipelines, along with their corresponding documentation and sample reports, are available at http://bioinformatics.mdc-berlin.de/pigx.

Development and validation of a nomogram to predict recurrence and melanoma-specific mortality in patients with negative sentinel lymph nodes

Background

Patients with melanoma and negative sentinel nodes (SNs) have varying outcomes, dependent on several prognostic factors. Considering all these factors in a prediction model might aid in identifying patients who could benefit from a personalized treatment strategy. The objective was to construct and validate a nomogram for recurrence and melanoma‐specific mortality (MSM) in patients with melanoma and negative SNs.

Methods

A total of 3220 patients with negative SNs were identified from a cohort of 4124 patients from four EORTC Melanoma Group centres who underwent sentinel lymph node biopsy. Prognostic factors for recurrence and MSM were studied with Cox regression analysis. Significant factors were incorporated in the models. Performance was assessed by discrimination (c‐index) and calibration in cross‐validation across the four centres. A nomogram was developed for graphical presentation.

Results

There were 3180 eligible patients. The final prediction model for recurrence and the calibrated model for MSM included three independent prognostic factors: ulceration, anatomical location and Breslow thickness. The c‐index was 0·74 for recurrence and 0·76 for the calibrated MSM model. Cross‐validation across the four centres showed reasonable model performance. A nomogram was developed based on these models. One‐third of the patients had a 5‐year recurrence probability of 8·2 per cent or less, and one‐third had a recurrence probability of 23·0 per cent or more.

Conclusion

A nomogram for predicting recurrence and MSM in patients with melanoma and negative SNs was constructed and validated. It could provide personalized estimates useful for tailoring surveillance strategies (reduce or increase intensity), and selection of patients for adjuvant therapy or clinical trials.

Could personalize care

Long non-coding RNA exchange during the oocyte-to-embryo transition in mice

The oocyte-to-embryo transition (OET) transforms a differentiated gamete into pluripotent blastomeres. The accompanying maternal-zygotic RNA exchange involves remodeling of the long non-coding RNA (lncRNA) pool. Here, we used next generation sequencing and de novo transcript assembly to define the core population of 1,600 lncRNAs expressed during the OET (lncRNAs). Relative to mRNAs, OET lncRNAs were less expressed and had shorter transcripts, mainly due to fewer exons and shorter 5′ terminal exons. Approximately half of OET lncRNA promoters originated in retrotransposons suggesting their recent emergence. Except for a small group of ubiquitous lncRNAs, maternal and zygotic lncRNAs formed two distinct populations. The bulk of maternal lncRNAs was degraded before the zygotic genome activation. Interestingly, maternal lncRNAs seemed to undergo cytoplasmic polyadenylation observed for dormant mRNAs. We also identified lncRNAs giving rise to trans-acting short interfering RNAs, which represent a novel lncRNA category. Altogether, we defined the core OET lncRNA transcriptome and characterized its remodeling during early development. Our results are consistent with the notion that rapidly evolving lncRNAs constitute signatures of cells-of-origin while a minority plays an active role in control of gene expression across OET. Our data presented here provide an excellent source for further OET lncRNA studies.

Widespread activation of antisense transcription of the host genome during herpes simplex virus 1 infection

Background

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

Results

Here, we show that HSV-1 induces the expression of about 1000 antisense transcripts from the human host cell genome. A subset of these is also activated by the closely related varicella zoster virus. Antisense transcripts originate either at gene promoters or within the gene body, and they show different susceptibility to the inhibition of early and immediate early viral gene expression. Overexpression of the major viral transcription factor ICP4 is sufficient to turn on a subset of antisense transcripts. Histone marks around transcription start sites of HSV-1-induced and constitutively transcribed antisense transcripts are highly similar, indicating that the genetic loci are already poised to transcribe these novel RNAs. Furthermore, an antisense transcript overlapping with the BBC3 gene (also known as PUMA) transcriptionally silences this potent inducer of apoptosis in cis.

Conclusions

We show for the first time that a virus induces widespread antisense transcription of the host cell genome. We provide evidence that HSV-1 uses this to downregulate a strong inducer of apoptosis. Our findings open new perspectives on global and specific alterations of host cell transcription by viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1329-5) contains supplementary material, which is available to authorized users.

RNA localization is a key determinant of neurite-enriched proteome

Protein subcellular localization is fundamental to the establishment of the body axis, cell migration, synaptic plasticity, and a vast range of other biological processes. Protein localization occurs through three mechanisms: protein transport, mRNA localization, and local translation. However, the relative contribution of each process to neuronal polarity remains unknown. Using neurons differentiated from mouse embryonic stem cells, we analyze protein and RNA expression and translation rates in isolated cell bodies and neurites genome-wide. We quantify 7323 proteins and the entire transcriptome, and identify hundreds of neurite-localized proteins and locally translated mRNAs. Our results demonstrate that mRNA localization is the primary mechanism for protein localization in neurites that may account for half of the neurite-localized proteome. Moreover, we identify multiple neurite-targeted non-coding RNAs and RNA-binding proteins with potential regulatory roles. These results provide further insight into the mechanisms underlying the establishment of neuronal polarity.

Subcellular localization of RNAs and proteins is important for polarized cells such as neurons. Here the authors differentiate mouse embryonic stem cells into neurons, and analyze the local transcriptome, proteome, and translated transcriptome in their cell bodies and neurites, providing a unique resource for future studies on neuronal polarity.

Long terminal repeats power evolution of genes and gene expression programs in mammalian oocytes and zygotes

Retrotransposons are “copy-and-paste” insertional mutagens that substantially contribute to mammalian genome content. Retrotransposons often carry long terminal repeats (LTRs) for retrovirus-like reverse transcription and integration into the genome. We report an extraordinary impact of a group of LTRs from the mammalian endogenous retrovirus-related ERVL retrotransposon class on gene expression in the germline and beyond. In mouse, we identified more than 800 LTRs from ORR1, MT, MT2, and MLT families, which resemble mobile gene-remodeling platforms that supply promoters and first exons. The LTR-mediated gene remodeling also extends to hamster, human, and bovine oocytes. The LTRs function in a stage-specific manner during the oocyte-to-embryo transition by activating transcription, altering protein-coding sequences, producing noncoding RNAs, and even supporting evolution of new protein-coding genes. These functions result, for example, in recycling processed pseudogenes into mRNAs or lncRNAs with regulatory roles. The functional potential of the studied LTRs is even higher, because we show that dormant LTR promoter activity can rescue loss of an essential upstream promoter. We also report a novel protein-coding gene evolution—D6Ertd527e—in which an MT LTR provided a promoter and the 5′ exon with a functional start codon while the bulk of the protein-coding sequence evolved through a CAG repeat expansion. Altogether, ERVL LTRs provide molecular mechanisms for stochastically scanning, rewiring, and recycling genetic information on an extraordinary scale. ERVL LTRs thus offer means for a comprehensive survey of the genome's expression potential, tightly intertwining with gene expression and evolution in the germline.

Maternal SIN3A regulates reprogramming of gene expression during mouse preimplantation development

The oocyte-to-embryo transition entails genome activation and a dramatic reprogramming of gene expression that is required for continued development. Superimposed on genome activation and reprogramming is development of a transcriptionally repressive state at the level of chromatin structure. Inducing global histone hyperacetylation relieves this repression and histone deacetylases 1 and 2 (HDAC1 and HDAC2) are involved in establishing the repressive state. Because SIN3A is an HDAC1/2-containing complex, we investigated whether it is involved in reprogramming gene expression during the course of genome activation. We find that Sin3a mRNA is recruited during maturation and that inhibiting its recruitment not only inhibits development beyond the 2-cell stage but also compromises the fidelity of reprogramming gene expression. The SIN3A that is synthesized during oocyte maturation reaches a maximum level in the mid-1-cell embryo and is essentially absent by the mid-2-cell stage. Overexpressing SIN3A in 1-cell embryos has no obvious effect on pre- and postimplantation development. These results provide a mechanism by which reprogramming can occur using a maternally inherited transcription machinery, namely, recruitment of mRNAs encoding transcription factors and chromatin remodelers, such as SIN3A.

Sculpting the Transcriptome During the Oocyte-to-Embryo Transition in Mouse

In mouse, the oocyte-to-embryo transition entails converting a highly differentiated oocyte to totipotent blastomeres. This transition is driven by degradation of maternal mRNAs, which results in loss of oocyte identity, and reprogramming of gene expression during the course of zygotic gene activation, which occurs primarily during the two-cell stage and confers blastomere totipotency. Full-grown oocytes are transcriptionally quiescent and mRNAs are remarkably stable in oocytes due to the RNA-binding protein MSY2, which stabilizes mRNAs, and low activity of the 5' and 3' RNA degradation machinery. Oocyte maturation initiates a transition from mRNA stability to instability due to phosphorylation of MSY2, which makes mRNAs more susceptible to the RNA degradation machinery, and recruitment of dormant maternal mRNAs that encode for critical components of the 5' and 3' RNA degradation machinery. Small RNAs (miRNA, siRNA, and piRNA) play little, if any, role in mRNA degradation that occurs during maturation. Many mRNAs are totally degraded but a substantial fraction is only partially degraded, their degradation completed by the end of the two-cell stage. Genome activation initiates during the one-cell stage, is promiscuous, low level, and genome wide (and includes both inter- and intragenic regions) and produces transcripts that are inefficiently spliced and polyadenylated. The major wave of genome activation in two-cell embryos involves expression of thousands of new genes. This unique pattern of gene expression is the product of maternal mRNAs recruited during maturation that encode for transcription factors and chromatin remodelers, as well as dramatic changes in chromatin structure due to incorporation of histone variants and modified histones.

The first murine zygotic transcription is promiscuous and uncoupled from splicing and 3' processing

Initiation of zygotic transcription in mammals is poorly understood. In mice, zygotic transcription is first detected shortly after pronucleus formation in 1-cell embryos, but the identity of the transcribed loci and mechanisms regulating their expression are not known. Using total RNA-Seq, we have found that transcription in 1-cell embryos is highly promiscuous, such that intergenic regions are extensively expressed and thousands of genes are transcribed at comparably low levels. Striking is that transcription can occur in the absence of defined core-promoter elements. Furthermore, accumulation of translatable zygotic mRNAs is minimal in 1-cell embryos because of inefficient splicing and 3' processing of nascent transcripts. These findings provide novel insights into regulation of gene expression in 1-cell mouse embryos that may confer a protective mechanism against precocious gene expression that is the product of a relaxed chromatin structure present in 1-cell embryos. The results also suggest that the first zygotic transcription itself is an active component of chromatin remodeling in 1-cell embryos.

Genomation: a toolkit to summarize, annotate and visualize genomic intervals

Biological insights can be obtained through computational integration of genomics data sets consisting of diverse types of information. The integration is often hampered by a large variety of existing file formats, often containing similar information, and the necessity to use complicated tools to achieve the desired results. We have built an R package, genomation, to expedite the extraction of biological information from high throughput data. The package works with a variety of genomic interval file types and enables easy summarization and annotation of high throughput data sets with given genomic annotations.

AVAILABILITY AND IMPLEMENTATION

The software is currently distributed under MIT artistic license and freely available at http://bioinformatics.mdc-berlin.de/genomation, and through the Bioconductor framework.

Archaeal aminoacyl-tRNA synthetases interact with the ribosome to recycle tRNAs

Aminoacyl-tRNA synthetases (aaRS) are essential enzymes catalyzing the formation of aminoacyl-tRNAs, the immediate precursors for encoded peptides in ribosomal protein synthesis. Previous studies have suggested a link between tRNA aminoacylation and high-molecular-weight cellular complexes such as the cytoskeleton or ribosomes. However, the structural basis of these interactions and potential mechanistic implications are not well understood. To biochemically characterize these interactions we have used a system of two interacting archaeal aaRSs: an atypical methanogenic-type seryl-tRNA synthetase and an archaeal ArgRS. More specifically, we have shown by thermophoresis and surface plasmon resonance that these two aaRSs bind to the large ribosomal subunit with micromolar affinities. We have identified the L7/L12 stalk and the proteins located near the stalk base as the main sites for aaRS binding. Finally, we have performed a bioinformatics analysis of synonymous codons in the Methanothermobacter thermautotrophicus genome that supports a mechanism in which the deacylated tRNAs may be recharged by aaRSs bound to the ribosome and reused at the next occurrence of a codon encoding the same amino acid. These results suggest a mechanism of tRNA recycling in which aaRSs associate with the L7/L12 stalk region to recapture the tRNAs released from the preceding ribosome in polysomes.

A cohesin-independent role for NIPBL at promoters provides insights in CdLS

The cohesin complex is crucial for chromosome segregation during mitosis and has recently also been implicated in transcriptional regulation and chromatin architecture. The NIPBL protein is required for the loading of cohesin onto chromatin, but how and where cohesin is loaded in vertebrate cells is unclear. Heterozygous mutations of NIPBL were found in 50% of the cases of Cornelia de Lange Syndrome (CdLS), a human developmental syndrome with a complex phenotype. However, no defects in the mitotic function of cohesin have been observed so far and the links between NIPBL mutations and the observed developmental defects are unclear. We show that NIPBL binds to chromatin in somatic cells with a different timing than cohesin. Further, we observe that high-affinity NIPBL binding sites localize to different regions than cohesin and almost exclusively to the promoters of active genes. NIPBL or cohesin knockdown reduce transcription of these genes differently, suggesting a cohesin-independent role of NIPBL for transcription. Motif analysis and comparison to published data show that NIPBL co-localizes with a specific set of other transcription factors. In cells derived from CdLS patients NIPBL binding levels are reduced and several of the NIPBL-bound genes have previously been observed to be mis-expressed in CdLS. In summary, our observations indicate that NIPBL mutations might cause developmental defects in different ways. First, defects of NIPBL might lead to cohesin-loading defects and thereby alter gene expression and second, NIPBL deficiency might affect genes directly via its role at the respective promoters.

The male germ cell gene regulator CTCFL is functionally different from CTCF and binds CTCF-like consensus sites in a nucleosome composition-dependent manner

Background

CTCF is a highly conserved and essential zinc finger protein expressed in virtually all cell types. In conjunction with cohesin, it organizes chromatin into loops, thereby regulating gene expression and epigenetic events. The function of CTCFL or BORIS, the testis-specific paralog of CTCF, is less clear.

Results

Using immunohistochemistry on testis sections and fluorescence-based microscopy on intact live seminiferous tubules, we show that CTCFL is only transiently present during spermatogenesis, prior to the onset of meiosis, when the protein co-localizes in nuclei with ubiquitously expressed CTCF. CTCFL distribution overlaps completely with that of Stra8, a retinoic acid-inducible protein essential for the propagation of meiosis. We find that absence of CTCFL in mice causes sub-fertility because of a partially penetrant testicular atrophy. CTCFL deficiency affects the expression of a number of testis-specific genes, including Gal3st1 and Prss50. Combined, these data indicate that CTCFL has a unique role in spermatogenesis. Genome-wide RNA expression studies in ES cells expressing a V5- and GFP-tagged form of CTCFL show that genes that are downregulated in CTCFL-deficient testis are upregulated in ES cells. These data indicate that CTCFL is a male germ cell gene regulator. Furthermore, genome-wide DNA-binding analysis shows that CTCFL binds a consensus sequence that is very similar to that of CTCF. However, only ~3,700 out of the ~5,700 CTCFL- and ~31,000 CTCF-binding sites overlap. CTCFL binds promoters with loosely assembled nucleosomes, whereas CTCF favors consensus sites surrounded by phased nucleosomes. Finally, an ES cell-based rescue assay shows that CTCFL is functionally different from CTCF.

Conclusions

Our data suggest that nucleosome composition specifies the genome-wide binding of CTCFL and CTCF. We propose that the transient expression of CTCFL in spermatogonia and preleptotene spermatocytes serves to occupy a subset of promoters and maintain the expression of male germ cell genes.

Using cloud computing infrastructure with CloudBioLinux, CloudMan, and Galaxy

Cloud computing has revolutionized availability and access to computing and storage resources, making it possible to provision a large computational infrastructure with only a few clicks in a Web browser. However, those resources are typically provided in the form of low-level infrastructure components that need to be procured and configured before use. In this unit, we demonstrate how to utilize cloud computing resources to perform open-ended bioinformatic analyses, with fully automated management of the underlying cloud infrastructure. By combining three projects, CloudBioLinux, CloudMan, and Galaxy, into a cohesive unit, we have enabled researchers to gain access to more than 100 preconfigured bioinformatics tools and gigabytes of reference genomes on top of the flexible cloud computing infrastructure. The protocol demonstrates how to set up the available infrastructure and how to use the tools via a graphical desktop interface, a parallel command-line interface, and the Web-based Galaxy interface.

Prediction of interacting protein residues using sequence and structure data

Identifying hotspots responsible for protein interactions with other macromolecules or drugs provides insight into functional aspects of the protein network, and is a pivotal task in systems biology and drug discovery. Here, we present the protocol for the application of a machine-learning method - Random Forest - to prediction of interacting residues in proteins, based on either the structural parameters or the primary sequence alone.

Characterisation of deep arterio-venous anastomoses within monochorionic placentae by vascular casting

OBJECTIVE

To characterise arterio-venous anastomoses (AVA) in monochorionic (MC) placentae and determine (i) whether shared cotyledons lie beneath the co-termination of an artery from one twin and a vein to the contralateral twin and (ii) whether all AVA can be detected by visual inspection of the chorionic plate.

METHODS

Vascular casts were made of 15 MC placentae. The number of typical AVAs suspected visually before digestion was compared with the number of AVAs identified after acid digestion.

RESULTS

Thirty-three of 67 (49%) suspected typical AVAs were confirmed as typical after casting. There were five false positives and no false negatives. The remainder were classified as atypical AVAs, found in > or =90% of MC placentae. Type I (small vascular connections between two apparently normal cotyledons not seen before casting) and Type II (shared cotyledons arising within larger apparently normal cotyledons) atypical AVAs were found in 53% and 73% of placentae, respectively.

CONCLUSIONS

Only half the shared cotyledons in MC placentae are characterised by co-termination of an artery and vein on the chorionic plate. We report the existence of deep anastomoses beneath the chorionic plate that cannot be visualised by chorionic plate inspection. These findings have implications for laser treatment of twin-twin transfusion syndrome.