Multiplex chromatin interactions with single-molecule precision.

The genomes of multicellular organisms are extensively folded into 3D chromosome territories within the nucleus1. Advanced 3D genome-mapping methods that combine proximity ligation and high-throughput sequencing (such as chromosome conformation capture, Hi-C)2, and chromatin immunoprecipitation techniques (such as chromatin interaction analysis by paired-end tag sequencing, ChIA-PET)3, have revealed topologically associating domains4 with frequent chromatin contacts, and have identified chromatin loops mediated by specific protein factors for insulation and regulation of transcription5-7. However, these methods rely on pairwise proximity ligation and reflect population-level views, and thus cannot reveal the detailed nature of chromatin interactions. Although single-cell Hi-C8 potentially overcomes this issue, this method may be limited by the sparsity of data that is inherent to current single-cell assays. Recent advances in microfluidics have opened opportunities for droplet-based genomic analysis9 but this approach has not yet been adapted for chromatin interaction analysis. Here we describe a strategy for multiplex chromatin-interaction analysis via droplet-based and barcode-linked sequencing, which we name ChIA-Drop. We demonstrate the robustness of ChIA-Drop in capturing complex chromatin interactions with single-molecule precision, which has not been possible using methods based on population-level pairwise contacts. By applying ChIA-Drop to Drosophila cells, we show that chromatin topological structures predominantly consist of multiplex chromatin interactions with high heterogeneity; ChIA-Drop also reveals promoter-centred multivalent interactions, which provide topological insights into transcription.

Picky comprehensively detects high-resolution structural variants in nanopore long reads.

Acquired genomic structural variants (SVs) are major hallmarks of cancer genomes, but they are challenging to reconstruct from short-read sequencing data. Here we exploited the long reads of the nanopore platform using our customized pipeline, Picky ( https://github.com/TheJacksonLaboratory/Picky ), to reveal SVs of diverse architecture in a breast cancer model. We identified the full spectrum of SVs with superior specificity and sensitivity relative to short-read analyses, and uncovered repetitive DNA as the major source of variation. Examination of genome-wide breakpoints at nucleotide resolution uncovered micro-insertions as the common structural features associated with SVs. Breakpoint density across the genome is associated with the propensity for interchromosomal connectivity and was found to be enriched in promoters and transcribed regions of the genome. Furthermore, we observed an over-representation of reciprocal translocations from chromosomal double-crossovers through phased SVs. We demonstrate that Picky analysis is an effective tool for comprehensive detection of SVs in cancer genomes from long-read data.

Chromatin connectivity maps reveal dynamic promoter-enhancer long-range associations.

In multicellular organisms, transcription regulation is one of the central mechanisms modelling lineage differentiation and cell-fate determination. Transcription requires dynamic chromatin configurations between promoters and their corresponding distal regulatory elements. It is believed that their communication occurs within large discrete foci of aggregated RNA polymerases termed transcription factories in three-dimensional nuclear space. However, the dynamic nature of chromatin connectivity has not been characterized at the genome-wide level. Here, through a chromatin interaction analysis with paired-end tagging approach using an antibody that primarily recognizes the pre-initiation complexes of RNA polymerase II, we explore the transcriptional interactomes of three mouse cells of progressive lineage commitment, including pluripotent embryonic stem cells, neural stem cells and neurosphere stem/progenitor cells. Our global chromatin connectivity maps reveal approximately 40,000 long-range interactions, suggest precise enhancer-promoter associations and delineate cell-type-specific chromatin structures. Analysis of the complex regulatory repertoire shows that there are extensive colocalizations among promoters and distal-acting enhancers. Most of the enhancers associate with promoters located beyond their nearest active genes, indicating that the linear juxtaposition is not the only guiding principle driving enhancer target selection. Although promoter-enhancer interactions exhibit high cell-type specificity, promoters involved in interactions are found to be generally common and mostly active among different cells. Chromatin connectivity networks reveal that the pivotal genes of reprogramming functions are transcribed within physical proximity to each other in embryonic stem cells, linking chromatin architecture to coordinated gene expression. Our study sets the stage for the full-scale dissection of spatial and temporal genome structures and their roles in orchestrating development.

Transcriptional responses of the marine green alga Micromonas pusilla and an infecting prasinovirus under different phosphate conditions.

Prasinophytes are widespread marine algae for which responses to nutrient limitation and viral infection are not well understood. We studied the picoprasinophyte, Micromonas pusilla, grown under phosphate-replete (0.65 ± 0.07 d-1 ) and 10-fold lower (low)-phosphate (0.11 ± 0.04 d-1 ) conditions, and infected by the phycodnavirus MpV-SP1. Expression of 17% of Micromonas genes in uninfected cells differed by >1.5-fold (q < 0.01) between nutrient conditions, with genes for P-metabolism and the uniquely-enriched Sel1-like repeat (SLR) family having higher relative transcript abundances, while phospholipid-synthesis genes were lower in low-P than P-replete. Approximately 70% (P-replete) and 30% (low-P) of cells were lysed 24 h post-infection, and expression of ≤5.8% of host genes changed relative to uninfected treatments. Host genes for CAZymes and glycolysis were activated by infection, supporting importance in viral production, which was significantly lower in slower growing (low-P) hosts. All MpV-SP1 genes were expressed, and our analyses suggest responses to differing host-phosphate backgrounds involve few viral genes, while the temporal program of infection involves many more, and is largely independent of host-phosphate background. Our study (i) identifies genes previously unassociated with nutrient acclimation or viral infection, (ii) provides insights into the temporal program of prasinovirus gene expression by hosts and (iii) establishes cell biological aspects of an ecologically important host-prasinovirus system that differ from other marine algal-virus systems.

A statistical method for detecting differentially expressed SNVs based on next-generation RNA-seq data.

In this article, we propose a new statistical method-MutRSeq-for detecting differentially expressed single nucleotide variants (SNVs) based on RNA-seq data. Specifically, we focus on nonsynonymous mutations and employ a hierarchical likelihood approach to jointly model observed mutation events as well as read count measurements from RNA-seq experiments. We then introduce a likelihood ratio-based test statistic, which detects changes not only in overall expression levels, but also in allele-specific expression patterns. In addition, this method can jointly test multiple mutations in one gene/pathway. The simulation studies suggest that the proposed method achieves better power than a few competitors under a range of different settings. In the end, we apply this method to a breast cancer data set and identify genes with nonsynonymous mutations differentially expressed between the triple negative breast cancer tumors and other subtypes of breast cancer tumors.

Marine algae and land plants share conserved phytochrome signaling systems.

Phytochrome photosensors control a vast gene network in streptophyte plants, acting as master regulators of diverse growth and developmental processes throughout the life cycle. In contrast with their absence in known chlorophyte algal genomes and most sequenced prasinophyte algal genomes, a phytochrome is found in Micromonas pusilla, a widely distributed marine picoprasinophyte (<2 µm cell diameter). Together with phytochromes identified from other prasinophyte lineages, we establish that prasinophyte and streptophyte phytochromes share core light-input and signaling-output domain architectures except for the loss of C-terminal response regulator receiver domains in the streptophyte phytochrome lineage. Phylogenetic reconstructions robustly support the presence of phytochrome in the common progenitor of green algae and land plants. These analyses reveal a monophyletic clade containing streptophyte, prasinophyte, cryptophyte, and glaucophyte phytochromes implying an origin in the eukaryotic ancestor of the Archaeplastida. Transcriptomic measurements reveal diurnal regulation of phytochrome and bilin chromophore biosynthetic genes in Micromonas. Expression of these genes precedes both light-mediated phytochrome redistribution from the cytoplasm to the nucleus and increased expression of photosynthesis-associated genes. Prasinophyte phytochromes perceive wavelengths of light transmitted farther through seawater than the red/far-red light sensed by land plant phytochromes. Prasinophyte phytochromes also retain light-regulated histidine kinase activity lost in the streptophyte phytochrome lineage. Our studies demonstrate that light-mediated nuclear translocation of phytochrome predates the emergence of land plants and likely represents a widespread signaling mechanism in unicellular algae.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants.

BACKGROUND: Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification. RESULTS: We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages. CONCLUSIONS: Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

ST2 as a marker for risk of therapy-resistant graft-versus-host disease and death.

BACKGROUND: No plasma biomarkers are associated with the response of acute graft-versus-host disease (GVHD) to therapy after allogeneic hematopoietic stem-cell transplantation. METHODS: We compared 12 biomarkers in plasma obtained a median of 16 days after therapy initiation from 10 patients with a complete response by day 28 after therapy initiation and in plasma obtained from 10 patients with progressive GVHD during therapy. The lead biomarker, suppression of tumorigenicity 2 (ST2), was measured at the beginning of treatment for GVHD in plasma from 381 patients and during the first month after transplantation in three independent sets totaling 673 patients to determine the association of this biomarker with treatment-resistant GVHD and 6-month mortality after treatment or transplantation. RESULTS: Of the 12 markers, ST2 had the most significant association with resistance to GVHD therapy and subsequent death without relapse. As compared with patients with low ST2 values at therapy initiation, patients with high ST2 values were 2.3 times as likely to have treatment-resistant GVHD (95% confidence interval [CI], 1.5 to 3.6) and 3.7 times as likely to die within 6 months after therapy (95% CI, 2.3 to 5.9). Patients with low ST2 values had lower mortality without relapse than patients with high ST2 values, regardless of the GVHD grade (11% vs. 31% among patients with grade I or II GVHD and 14% vs. 67% among patients with grade III or IV GVHD, P<0.001 for both comparisons). Plasma ST2 values at day 14 after transplantation were associated with 6-month mortality without relapse, regardless of the intensity of the conditioning regimen. CONCLUSIONS: ST2 levels measured at the initiation of therapy for GVHD and during the first month after transplantation improved risk stratification for treatment-resistant GVHD and death without relapse after transplantation. (Funded by the National Institutes of Health.)

Extrachromosomal DNA Associates with Nuclear Condensates and Reorganizes Chromatin Structures to Enhance Oncogenic Transcription.

Extrachromosomal, circular DNA (ecDNA) is a prevalent oncogenic alteration in cancer genomes, often associated with aggressive tumor behavior and poor patient outcome. While previous studies proposed a chromatin-based mobile enhancer model for ecDNA-driven oncogenesis, its precise mechanism and impact remains unclear across diverse cancer types. Our study, utilizing advanced multi-omics profiling, epigenetic editing, and imaging approaches in three cancer models, reveals that ecDNA hubs are an integrated part of nuclear condensates and exhibit cancer-type specific chromatin connectivity. Epigenetic silencing of the ecDNA-specific regulatory modules or chemically disrupting liquid-liquid phase separation breaks down ecDNA hubs, displaces MED1 co-activator binding, inhibits oncogenic transcription, and promotes cell death. These findings substantiate the trans -activator function of ecDNA and underscore a structural mechanism driving oncogenesis. This refined understanding expands our views of oncogene regulation and opens potential avenues for novel therapeutic strategies in cancer treatment.

A novel soluble form of Tim-3 associated with severe graft-versus-host disease.

The T cell Ig and mucin domain 3 (Tim-3) receptor has been implicated as a negative regulator of adaptive immune responses. We have utilized a proteomic strategy to identify novel proteins associated with graft versus host disease (GVHD) after allogeneic hematopoietic cell transplantation (HCT). Mass spectrometry analysis of plasma from subjects with mid-gut and upper-gut GVHD compared with those without GVHD identified increased levels of a protein identified with high confidence as Tim-3. A follow-up validation study using an immunoassay to measure Tim-3 levels in individual plasma samples from 127 patients demonstrated significantly higher plasma Tim-3 concentrations in patients with the more severe mid-gut GVHD, compared with those with upper-gut GVHD (P = .005), patients without GVHD (P = .002), and normal controls (P < .0001). Surface expression of Tim-3 was increased on CD8(+) T cells from patients with grade 2 to 4 acute GVHD (P = .01). Mass spectrometry-based profiling of plasma from multiple subjects diagnosed with common diseases provided evidence for restricted release of soluble Tim-3 in the context of GVHD. These findings have mechanistic implications for the development of novel strategies for targeting the Tim-3 immune regulatory pathway as an approach to improving control of GVHD.

ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs.

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection.

Background: It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms. Methods: To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq. Results: Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity. Conclusions: Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription.

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

Targeted quantitative analysis of Streptococcus pyogenes virulence factors by multiple reaction monitoring.

In many studies, particularly in the field of systems biology, it is essential that identical protein sets are precisely quantified in multiple samples such as those representing differentially perturbed cell states. The high degree of reproducibility required for such experiments has not been achieved by classical mass spectrometry-based proteomics methods. In this study we describe the implementation of a targeted quantitative approach by which predetermined protein sets are first identified and subsequently quantified at high sensitivity reliably in multiple samples. This approach consists of three steps. First, the proteome is extensively mapped out by multidimensional fractionation and tandem mass spectrometry, and the data generated are assembled in the PeptideAtlas database. Second, based on this proteome map, peptides uniquely identifying the proteins of interest, proteotypic peptides, are selected, and multiple reaction monitoring (MRM) transitions are established and validated by MS2 spectrum acquisition. This process of peptide selection, transition selection, and validation is supported by a suite of software tools, TIQAM (Targeted Identification for Quantitative Analysis by MRM), described in this study. Third, the selected target protein set is quantified in multiple samples by MRM. Applying this approach we were able to reliably quantify low abundance virulence factors from cultures of the human pathogen Streptococcus pyogenes exposed to increasing amounts of plasma. The resulting quantitative protein patterns enabled us to clearly define the subset of virulence proteins that is regulated upon plasma exposure.

Chromatin interaction analyses elucidate the roles of PRC2-bound silencers in mouse development.

Lineage-specific gene expression is modulated by a balance between transcriptional activation and repression during animal development. Knowledge about enhancer-centered transcriptional activation has advanced considerably, but silencers and their roles in normal development remain poorly understood. Here, we performed chromatin interaction analyses of Polycomb repressive complex 2 (PRC2), a key inducer of transcriptional gene silencing, to uncover silencers, their molecular identity and associated chromatin connectivity. Systematic analysis of cis-regulatory silencer elements reveals their chromatin features and gene-targeting specificity. Deletion of certain PRC2-bound silencers in mice results in transcriptional derepression of their interacting genes and pleiotropic developmental phenotypes, including embryonic lethality. While some PRC2-bound elements function as silencers in pluripotent cells, they can transition into active tissue-specific enhancers during development, highlighting their regulatory versatility. Our study characterizes the molecular profile of silencers and their associated chromatin architectures, and suggests the possibility of targeted reactivation of epigenetically silenced genes.

Quantitative serum proteomics using dual stable isotope coding and nano LC-MS/MSMS.

Stable isotope coding technique in combination with mass spectrometry has emerged as a powerful tool to accurately identify and differentially quantify proteins within complex protein mixtures. We present a novel methodology to increase the yield of quantified proteins while maintaining a high stable-isotopic labeling efficacy. With this approach, intact proteins in complex biological sample such as sera are labeled with the designated dual stable isotope coding (DSIC) systems. In brief, intact proteins are coded sequentially with acrylamide to label Cysteine residues (Cys) and with succinic anhydride to label Lysine residues (Lys). Protein samples coded with this dual stable isotope are subjected to an online 2D-HPLC fractionation. The resolved protein fractions are individually digested with trypsin and analyzed with nano LC-MS/MSMS. Our results show that the DSIC labeling efficiency is 100% for Cysteine (Cys) labeled with acrylamide and 98% for Lysine (Lys) labeled with succinic anhydride. A comparative analysis of DSIC labeling and single labeling of Cysteine residues was made. Analysis of an entire anion-exchange chromatography subfraction of sera yielded 165 identified proteins (criteria: error rate <5% and unique peptides >or=2), 104 of which were quantified using the single labeling method (i.e., Cysteine acrylamide labeling only). In contrast, using same criteria for identification, a total 185 proteins were identified and 174 proteins were quantified using the DSIC labeling technique.

Ultra-long Read Sequencing for Whole Genomic DNA Analysis.

Third generation single-molecule DNA sequencing technologies offer significantly longer read length that can facilitate the assembly of complex genomes and analysis of complex structural variants. Nanopore platforms perform single-molecule sequencing by directly measuring the current changes mediated by DNA passage through the pores and can generate hundreds of kilobase (kb) reads with minimal capital cost. This platform has been adopted by many researchers for a variety of applications. Achieving longer sequencing read lengths is the most critical factor to leverage the value of nanopore sequencing platforms. To generate ultra-long reads, special consideration is required to avoid DNA breakages and gain efficiency to generate productive sequencing templates. Here, we provide the detailed protocol of ultra-long DNA sequencing including high molecular weight (HMW) DNA extraction from fresh or frozen cells, library construction by mechanical shearing or transposase fragmentation, and sequencing on a nanopore device. From 20-25 µg of HMW DNA, the method can achieve N50 read length of 50-70 kb with mechanical shearing and N50 of 90-100 kb read length with transposase mediated fragmentation. The protocol can be applied to DNA extracted from mammalian cells to perform whole genome sequencing for the detection of structural variants and genome assembly. Additional improvements on the DNA extraction and enzymatic reactions will further increase the read length and expand its utility.

Complex Microbiome in Brain Abscess Revealed by Whole-Genome Culture-Independent and Culture-Based Sequencing.

Brain abscess is a severe infectious disease with high mortality and mobility. Although culture-based techniques have been widely used for the investigation of microbial composition of brain abscess, these approaches are inherent biased. Recent studies using 16S ribosomal sequencing approaches revealed high complexity of the bacterial community involved in brain abscess but fail to detect fungal and viral composition. In the study, both culture-independent nanopore metagenomic sequencing and culture-based whole-genome sequencing using both the Illumina and the Nanopore platforms were conducted to investigate the microbial composition and genomic characterization in brain abscess. Culture-independent metagenomic sequencing revealed not only a larger taxonomic diversity of bacteria but also the presence of fungi and virus communities. The culture-based whole-genome sequencing identified a novel species in Prevotella and reconstructs a Streptococcus constellatus with a high GC-skew genome. Antibiotic-resistance genes CfxA and ErmF associated with resistance to penicillin and clindamycin were also identified in culture-based and culture-free sequencing. This study implies current understanding of brain abscess need to consider the broader diversity of microorganisms.

Linked-read Sequencing Analysis Reveals Tumor-specific Genome Variation Landscapes in Neurofibromatosis Type 2 (NF2) Patients.

HYPOTHESIS: We hypothesize that genomic variants including deletions, insertions, inversions, and tandem duplications beyond the changes in tumor suppressor NF2 gene affect gene expression of tumor-specific pathways in vestibular schwannomas (VS) patients with Neurofibromatosis type 2 (NF2), thus contributing to their clinical behavior. BACKGROUND: Genomic variation could reconfigure transcription in NF2 transformation process. Therefore, genome-wide high-resolution characterization of structural variants (SV) landscapes in NF2 tumors can expand our understanding of the genes regulating the clinical phenotypes in NF2-associated VS. METHODS: We performed whole-genome haplotype-specific structural variation analysis using synthetic linked reads generated through microfluidics-based barcoding of high molecular weight DNA followed by high-coverage Illumina paired-end whole-genome sequencing from 10 patients' tumors of different growth rates and their matching blood samples. RESULTS: NF2 tumor-specific deletions and large SVs were detected and can be classified based on their association with tumor growth rates. Through detailed annotation of these mutations, we uncover common alleles affected by these deletions and large SVs that can be associated with signaling pathways implicated in cell proliferation and tumorigenesis. CONCLUSION: The genomic variation landscape of NF2-related VS was investigated through whole-genome linked-read sequencing. Large SVs, in addition to deletions, were identified and may serve as modulators of clinical behavior.

Mapping the Global Chromatin Connectivity Network for Sox2 Function in Neural Stem Cell Maintenance.

The SOX2 transcription factor is critical for neural stem cell (NSC) maintenance and brain development. Through chromatin immunoprecipitation (ChIP) and chromatin interaction analysis (ChIA-PET), we determined genome-wide SOX2-bound regions and Pol II-mediated long-range chromatin interactions in brain-derived NSCs. SOX2-bound DNA was highly enriched in distal chromatin regions interacting with promoters and carrying epigenetic enhancer marks. Sox2 deletion caused widespread reduction of Pol II-mediated long-range interactions and decreased gene expression. Genes showing reduced expression in Sox2-deleted cells were significantly enriched in interactions between promoters and SOX2-bound distal enhancers. Expression of one such gene, Suppressor of Cytokine Signaling 3 (Socs3), rescued the self-renewal defect of Sox2-ablated NSCs. Our work identifies SOX2 as a major regulator of gene expression through connections to the enhancer network in NSCs. Through the definition of such a connectivity network, our study shows the way to the identification of genes and enhancers involved in NSC maintenance and neurodevelopmental disorders.