Recurrent chromosomal translocations in sarcomas create a megacomplex that mislocalizes NuA4/TIP60 to Polycomb target loci.

Chromosomal translocations frequently promote carcinogenesis by producing gain-of-function fusion proteins. Recent studies have identified highly recurrent chromosomal translocations in patients with endometrial stromal sarcomas (ESSs) and ossifying fibromyxoid tumors (OFMTs), leading to an in-frame fusion of PHF1 (PCL1) to six different subunits of the NuA4/TIP60 complex. While NuA4/TIP60 is a coactivator that acetylates chromatin and loads the H2A.Z histone variant, PHF1 is part of the Polycomb repressive complex 2 (PRC2) linked to transcriptional repression of key developmental genes through methylation of histone H3 on lysine 27. In this study, we characterize the fusion protein produced by the EPC1-PHF1 translocation. The chimeric protein assembles a megacomplex harboring both NuA4/TIP60 and PRC2 activities and leads to mislocalization of chromatin marks in the genome, in particular over an entire topologically associating domain including part of the HOXD cluster. This is linked to aberrant gene expression-most notably increased expression of PRC2 target genes. Furthermore, we show that JAZF1-implicated with a PRC2 component in the most frequent translocation in ESSs, JAZF1-SUZ12-is a potent transcription activator that physically associates with NuA4/TIP60, its fusion creating outcomes similar to those of EPC1-PHF1 Importantly, the specific increased expression of PRC2 targets/HOX genes was also confirmed with ESS patient samples. Altogether, these results indicate that most chromosomal translocations linked to these sarcomas use the same molecular oncogenic mechanism through a physical merge of NuA4/TIP60 and PRC2 complexes, leading to mislocalization of histone marks and aberrant Polycomb target gene expression.

Non-canonical transcriptional regulation of the poor prognostic factor UGT2B17 in chronic lymphocytic leukemic and normal B cells.

BACKGROUND: High expression of the glycosyltransferase UGT2B17 represents an independent adverse prognostic marker in chronic lymphocytic leukemia (CLL). It also constitutes a predictive marker for therapeutic response and a drug resistance mechanism. The key determinants driving expression of the UGT2B17 gene in normal and leukemic B-cells remain undefined. The UGT2B17 transcriptome is complex and is comprised of at least 10 alternative transcripts, identified by previous RNA-sequencing of liver and intestine. We hypothesized that the transcriptional program regulating UGT2B17 in B-lymphocytes is distinct from the canonical expression previously characterized in the liver. RESULTS: RNA-sequencing and genomics data revealed a specific genomic landscape at the UGT2B17 locus in normal and leukemic B-cells. RNA-sequencing and quantitative PCR data indicated that the UGT2B17 enzyme is solely encoded by alternative transcripts expressed in CLL patient cells and not by the canonical transcript widely expressed in the liver and intestine. Chromatin accessible regions (ATAC-Seq) in CLL cells mapped with alternative promoters and non-coding exons, which may be derived from endogenous retrotransposon elements. By luciferase reporter assays, we identified key cis-regulatory STAT3, RELA and interferon regulatory factor (IRF) binding sequences driving the expression of UGT2B17 in lymphoblastoid and leukemic B-cells. Electrophoretic mobility shift assays and pharmacological inhibition demonstrated key roles for the CLL prosurvival transcription factors STAT3 and NF-κB in the leukemic expression of UGT2B17. CONCLUSIONS: UGT2B17 expression in B-CLL is driven by key regulators of CLL progression. Our data suggest that a NF-κB/STAT3/IRF/UGT2B17 axis may represent a novel B-cell pathway promoting disease progression and drug resistance.

MRG Proteins Are Shared by Multiple Protein Complexes With Distinct Functions.

MRG15/MORF4L1 is a highly conserved protein in eukaryotes that contains a chromodomain (CHD) recognizing methylation of lysine 36 on histone H3 (H3K36me3) in chromatin. Intriguingly, it has been reported in the literature to interact with several different factors involved in chromatin modifications, gene regulation, alternative mRNA splicing, and DNA repair by homologous recombination. To get a complete and reliable picture of associations in physiological conditions, we used genome editing and tandem affinity purification to analyze the stable native interactome of human MRG15, its paralog MRGX/MORF4L2 that lacks the CHD, and MRGBP (MRG-binding protein) in isogenic K562 cells. We found stable interchangeable association of MRG15 and MRGX with the NuA4/TIP60 histone acetyltransferase/chromatin remodeler, Sin3B histone deacetylase/demethylase, ASH1L histone methyltransferase, and PALB2-BRCA2 DNA repair protein complexes. These associations were further confirmed and analyzed by CRISPR tagging of endogenous proteins and comparison of expressed isoforms. Importantly, based on structural information, point mutations could be introduced that specifically disrupt MRG15 association with some complexes but not others. Most interestingly, we also identified a new abundant native complex formed by MRG15/X-MRGBP-BRD8-EP400NL (EP400 N-terminal like) that is functionally similar to the yeast TINTIN (Trimer Independent of NuA4 for Transcription Interactions with Nucleosomes) complex. Our results show that EP400NL, being homologous to the N-terminal region of NuA4/TIP60 subunit EP400, creates TINTIN by competing for BRD8 association. Functional genomics indicate that human TINTIN plays a role in transcription of specific genes. This is most likely linked to the H4ac-binding bromodomain of BRD8 along the H3K36me3-binding CHD of MRG15 on the coding region of transcribed genes. Taken together, our data provide a complete detailed picture of human MRG proteins-associated protein complexes, which are essential to understand and correlate their diverse biological functions in chromatin-based nuclear processes.

Multi-Omics Data Integration Reveals Key Variables Contributing to Subgingival Microbiome Dysbiosis-Induced Inflammatory Response in a Hyperglycemic Microenvironment.

Subgingival microbiome dysbiosis promotes the development of periodontitis, an irreversible chronic inflammatory disease associated with metabolic diseases. However, studies regarding the effects of a hyperglycemic microenvironment on host-microbiome interactions and host inflammatory response during periodontitis are still scarce. Here, we investigated the impacts of a hyperglycemic microenvironment on the inflammatory response and transcriptome of a gingival coculture model stimulated with dysbiotic subgingival microbiomes. HGF-1 cells overlaid with U937 macrophage-like cells were stimulated with subgingival microbiomes collected from four healthy donors and four patients with periodontitis. Pro-inflammatory cytokines and matrix metalloproteinases were measured while the coculture RNA was submitted to a microarray analysis. Subgingival microbiomes were submitted to 16s rRNA gene sequencing. Data were analyzed using an advanced multi-omics bioinformatic data integration model. Our results show that the genes krt76, krt27, pnma5, mansc4, rab41, thoc6, tm6sf2, and znf506 as well as the pro-inflammatory cytokines IL-1ß, GM-CSF, FGF2, IL-10, the metalloproteinases MMP3 and MMP8, and bacteria from the ASV 105, ASV 211, ASV 299, Prevotella, Campylobacter and Fretibacterium genera are key intercorrelated variables contributing to periodontitis-induced inflammatory response in a hyperglycemic microenvironment. In conclusion, our multi-omics integration analysis unveiled the complex interrelationships involved in the regulation of periodontal inflammation in response to a hyperglycemic microenvironment.

Differential Inhibition of HIV Replication by the 12 Interferon Alpha Subtypes.

Type I interferons (IFNs) are a family of cytokines that represent a first line of defense against virus infections. The 12 different IFN-α subtypes share a receptor on target cells and trigger similar signaling cascades. Several studies have collectively shown that this apparent redundancy conceals qualitatively different responses induced by individual subtypes, which display different efficacies of inhibition of HIV replication. Some studies, however, provided evidence that the disparities are quantitative rather than qualitative. Since RNA expression analyses show a large but incomplete overlap of the genes induced, they may support both models. To explore if the IFN-α subtypes induce functionally relevant different anti-HIV activities, we have compared the efficacies of inhibition of all 12 subtypes on HIV spread and on specific steps of the viral replication cycle, including viral entry, reverse transcription, protein synthesis, and virus release. Finding different hierarchies of inhibition would validate the induction of qualitatively different responses. We found that while most subtypes similarly inhibit virus entry, they display distinctive potencies on other early steps of HIV replication. In addition, only some subtypes were able to target effectively the late steps. The extent of induction of known anti-HIV factors helps to explain some, but not all differences observed, confirming the participation of additional IFN-induced anti-HIV effectors. Our findings support the notion that different IFN-α subtypes can induce the expression of qualitatively different antiviral activities. IMPORTANCE The initial response against viruses relies in large part on type I interferons, which include 12 subtypes of IFN-α. These cytokines bind to a common receptor on the cell surface and trigger the expression of incompletely overlapping sets of genes. Whether the anti-HIV responses induced by IFN-α subtypes differ in the extent of expression or in the nature of the genes involved remains debated. Also, RNA expression profiles led to opposite conclusions, depending on the importance attributed to the induction of common or distinctive genes. To explore if relevant anti-HIV activities can be differently induced by the IFN-α subtypes, we compared their relative efficacies on specific steps of the replication cycle. We show that the hierarchy of IFN potencies depends on the step analyzed, supporting qualitatively different responses. This work will also prompt the search for novel IFN-induced anti-HIV factors acting on specific steps of the replication cycle.

Extensive and Accurate Benchmarking of DIA Acquisition Methods and Software Tools Using a Complex Proteomic Standard.

Over the past decade, the data-independent acquisition mode has gained popularity for broad coverage of complex proteomes by LC-MS/MS and quantification of low-abundance proteins. However, there is no consensus in the literature on the best data acquisition parameters and processing tools to use for this specific application. Here, we present the most comprehensive comparison of DIA workflows on Orbitrap instruments published so far in the field of proteomics. Using a standard human 48 proteins mixture (UPS1-Sigma) at 8 different concentrations in an E. coli proteome background, we tested 36 workflows including 4 different DIA window acquisition schemes and 6 different software tools (DIA-NN, DIA-Umpire, OpenSWATH, ScaffoldDIA, Skyline, and Spectronaut) with or without the use of a DDA spectral library. On the basis of the number of proteins identified, quantification linearity and reproducibility, as well as sensitivity and specificity in 28 pairwise comparisons of different UPS1 concentrations, we summarize the major considerations and propose guidelines for choosing the DIA workflow best suited for LC-MS/MS proteomic analyses. Our 96 DIA raw files and software outputs have been deposited on ProteomeXchange for testing or developing new DIA processing tools.

Hedgehog signaling pathway regulates gene expression profile of epididymal principal cells through the primary cilium.

Primary cilia (PC) are organelles that sense and respond to dynamic changes of the extracellular milieu through the regulation of target genes. By using the epididymis as a model system, we determined the contribution of primary cilia in the regulation of epithelial cell functions through the transduction of the Hedgehog (Hh) signaling pathway. Both Sonic (SHH) and Indian Hedgehog (IHH) ligands were detected in epididymal epithelial cells by confocal microscopy and found secreted in the extracellular space. Gene expression profiling preformed on ciliated epithelial cells indicated that 153 and 1052 genes were differentially expressed following treatment with the Hh agonist SAG or the Hh antagonist cyclopamine (Cyclo), respectively. Strikingly, gene ontology analysis indicated that genes associated with immune response were the most affected following Hh modulation. The contribution of epididymal PC to canonical Hh pathway transduction was validated by ciliobrevin D treatment, which induced a significant decrease in PC length and a reduction in the expression Hh signaling targets. Such findings bring us closer to a molecular understanding of the subtle immune balance observed in some epithelia, including the epididymis and the intestine, which are organs featuring both tolerance toward autoimmune spermatozoa (or commensal bacteria) and defense against pathogens.

Identification of Abundant and Functional dodecaRNAs (doRNAs) Derived from Ribosomal RNA.

Using a modified RNA-sequencing (RNA-seq) approach, we discovered a new family of unusually short RNAs mapping to ribosomal RNA 5.8S, which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. Using a new quantitative detection method that we developed, we confirmed our RNA-seq data and determined that the minimal core doRNA sequence and its 13-nt variant C-doRNA (doRNA with a 5' Cytosine) are the two most abundant doRNAs, which, together, may outnumber microRNAs. The C-doRNA/doRNA ratio is stable within species but differed between species. doRNA and C-doRNA are mainly cytoplasmic and interact with heterogeneous nuclear ribonucleoproteins (hnRNP) A0, A1 and A2B1, but not Argonaute 2. Reporter gene activity assays suggest that C-doRNA may function as a regulator of Annexin II receptor (AXIIR) expression. doRNAs are differentially expressed in prostate cancer cells/tissues and may control cell migration. These findings suggest that unusually short RNAs may be more abundant and important than previously thought.

UGT2B17 modifies drug response in chronic lymphocytic leukaemia.

BACKGROUND: High UGT2B17 is associated with poor prognosis in untreated chronic lymphocytic leukaemia (CLL) patients and its expression is induced in non-responders to fludarabine-containing regimens. We examined whether UGT2B17, the predominant lymphoid glucuronosyltransferase, affects leukaemic drug response and is involved in the metabolic inactivation of anti-leukaemic agents. METHODS: Functional enzymatic assays and patients' plasma samples were analysed by mass-spectrometry to evaluate drug inactivation by UGT2B17. Cytotoxicity assays and RNA sequencing were used to assess drug response and transcriptome changes associated with high UGT2B17 levels. RESULTS: High UGT2B17 in B-cell models led to reduced sensitivity to fludarabine, ibrutinib and idelalisib. UGT2B17 expression in leukaemic cells involved a non-canonical promoter and was induced by short-term treatment with these anti-leukaemics. Glucuronides of both fludarabine and ibrutinib were detected in CLL patients on respective treatment, however UGT2B17 conjugated fludarabine but not ibrutinib. AMP-activated protein kinase emerges as a pathway associated with high UGT2B17 in fludarabine-treated patients and drug-treated cell models. The expression changes linked to UGT2B17 exposed nuclear factor kappa B as a key regulatory hub. CONCLUSIONS: Data imply that UGT2B17 represents a mechanism altering drug response in CLL through direct inactivation but would also involve additional mechanisms for drugs not inactivated by UGT2B17.

Adipose Tissue Transcriptome Is Related to Pollutant Exposure in Polar Bear Mother-Cub Pairs from Svalbard, Norway.

Being at the food chain apex, polar bears (Ursus maritimus) are highly contaminated with persistent organic pollutants (POPs). Females transfer POPs to their offspring through gestation and lactation; therefore, young cubs present higher POPs concentrations than their mothers. Recent studies suggest that POPs affect the lipid metabolism in female polar bears; however, the mechanisms and impact on their offspring remain unknown. Here, we hypothesized that exposure to POPs differentially alters genome-wide gene transcription in the adipose tissue from mother polar bears and their cubs, highlighting molecular differences in response between adults and young. Adipose tissue biopsies were collected from 13 adult female polar bears and their twin cubs in Svalbard, Norway, in April 2011, 2012, and 2013. Total RNA extracted from biopsies was subjected to next-generation RNA sequencing. Plasma concentrations of summed polychlorinated biphenyls, organochlorine pesticides, and polybrominated diphenyl ethers in mothers ranged from 897 to 13620 ng/g wet weight and were associated with altered adipose tissue gene expression in both mothers and cubs. In mothers, 2502 and 2586 genes in total were positively and negatively, respectively, correlated to POP exposure, whereas in cubs, 2585 positively and 1690 negatively genes. Between mothers and cubs, 743 positively and negatively genes overlapped between mothers and cubs suggesting partially shared molecular responses to ΣPOPs. ΣPOP-associated genes were involved in numerous metabolic pathways in mothers and cubs, indicating that POP exposure alters the energy metabolism, which, in turn, may be linked to metabolic dysfunction.

VEXOR: an integrative environment for prioritization of functional variants in fine-mapping analysis.

Motivation: The identification of the functional variants responsible for observed genome-wide association studies (GWAS) signals is one of the most challenging tasks of the post-GWAS research era. Several tools have been developed to annotate genetic variants by their genomic location and potential functional implications. Each of these tools has its own requirements and internal logic, which forces the user to become acquainted with each interface. Results: From an awareness of the amount of work needed to analyze a single locus, we have built a flexible, versatile and easy-to-use web interface designed to help in prioritizing variants and predicting their potential functional implications. This interface acts as a single-point of entry linking association results with reference tools and relevant experiments. Availability and Implementation: VEXOR is an integrative web application implemented through the Shiny framework and available at: http://romix.genome.ulaval.ca/vexor. Contact: arnaud.droit@crchuq.ulaval.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns.

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a "gradient effect" where the regulatory factor occupancy levels follow transcription and ii) a "threshold effect" where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.

rTANDEM, an R/Bioconductor package for MS/MS protein identification.

SUMMARY: rTANDEM is an R/Bioconductor package that interfaces the X!Tandem protein identification algorithm. The package can run the multi-threaded algorithm on proteomic data files directly from R. It also provides functions to convert search parameters and results to/from R as well as functions to manipulate parameters and automate searches. An associated R package, shinyTANDEM, provides a web-based graphical interface to visualize and interpret the results. Together, those two packages form an entry point for a general MS/MS-based proteomic pipeline in R/Bioconductor. AVAILABILITY AND IMPLEMENTATION: rTANDEM and shinyTANDEM are distributed in R/Bioconductor, http://bioconductor.org/packages/release/bioc/. The packages are under open licenses (GPL-3 and Artistice-1.0). CONTACT: frederic.fournier@crchuq.ulaval.ca or arnaud.droit@crchuq.ulaval.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

NGS++: a library for rapid prototyping of epigenomics software tools.

MOTIVATION: The development of computational tools to enable testing and analysis of high-throughput-sequencing data is essential to modern genomics research. However, although multiple frameworks have been developed to facilitate access to these tools, comparatively little effort has been made at implementing low-level programming libraries to increase the speed and ease of their development. RESULTS: We propose NGS++, a programming library in C++11 specialized in manipulating both next-generation sequencing (NGS) datasets and genomic information files. This library allows easy integration of new formats and rapid prototyping of new functionalities with a focus on the analysis of genomic regions and features. It offers a powerful, yet versatile and easily extensible interface to read, write and manipulate multiple genomic file formats. By standardizing the internal data structures and presenting a common interface to the data parser, NGS++ offers an effective framework for epigenomics tool development. AVAILABILITY: NGS++ was written in C++ using the C++11 standard. It requires minimal efforts to build and is well-documented via a complete docXygen guide, online documentation and tutorials. Source code, tests, code examples and documentation are available via the website at http://www.ngsplusplus.ca and the github repository at https://github.com/NGS-lib/NGSplusplus. CONTACT: nicolas.gevry@usherbrooke.ca or arnaud.droit@crchuq.ulaval.ca.

Myelin-reactive B cells exacerbate CD4+ T cell-driven CNS autoimmunity in an IL-23-dependent manner.

B cells and T cells collaborate in multiple sclerosis (MS) pathogenesis. IgH[MOG] mice possess a B cell repertoire skewed to recognize myelin oligodendrocyte glycoprotein (MOG). Here, we show that upon immunization with the T cell-obligate autoantigen, MOG[35-55], IgH[MOG] mice develop rapid and exacerbated experimental autoimmune encephalomyelitis (EAE) relative to wildtype (WT) counterparts, characterized by aggregation of T and B cells in the IgH[MOG] meninges and by CD4+ T helper 17 (Th17) cells in the CNS. Production of the Th17 maintenance factor IL-23 is observed from IgH[MOG] CNS-infiltrating and meningeal B cells, and in vivo blockade of IL-23p19 attenuates disease severity in IgH[MOG] mice. In the CNS parenchyma and dura mater of IgH[MOG] mice, we observe an increased frequency of CD4+PD-1+CXCR5- T cells that share numerous characteristics with the recently described T peripheral helper (Tph) cell subset. Further, CNS-infiltrating B and Tph cells from IgH[MOG] mice show increased reactive oxygen species (ROS) production. Meningeal inflammation, Tph-like cell accumulation in the CNS and B/Tph cell production of ROS were all reduced upon p19 blockade. Altogether, MOG-specific B cells promote autoimmune inflammation of the CNS parenchyma and meninges in an IL-23-dependent manner.

Rapid and Technically Simple Detection of SARS-CoV-2 Variants Using CRISPR Cas12 and Cas13.

The worldwide proliferation of the SARS-CoV-2 virus in the past 3 years has allowed the virus to accumulate numerous mutations. Dangerous lineages have emerged one after another, each leading to a new wave of the pandemic. In this study, we have developed the THRASOS pipeline to rapidly discover lineage-specific mutation signatures and thus advise the development of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-based diagnostic tests. We also optimized a strategy to modify loop-mediated isothermal amplification amplicons for downstream use with Cas12 and Cas13 for future multiplexing. The close ancestry of the BA.1 and BA.2 variants of SARS-CoV-2 (Omicron) made these excellent candidates for the development of a first test using this workflow. With a quick turnaround time and low requirements for laboratory equipment, the test we have created is ideally suited for settings such as mobile clinics lacking equipment such as Next-Generation Sequencers or Sanger Sequencers and the personnel to run these devices.

Early ART reduces viral seeding and innate immunity in liver and lungs of SIV-infected macaques.

Identifying immune cells and anatomical tissues that contribute to the establishment of viral reservoirs is of central importance in HIV-1 cure research. Herein, we used rhesus macaques (RMs) infected with SIVmac251 to analyze viral seeding in the liver and lungs of either untreated or early antiretroviral therapy-treated (ART-treated) RMs. Consistent with viral replication and sensing, transcriptomic analyses showed higher levels of inflammation, pyroptosis, and chemokine genes as well as of interferon-stimulating gene (ISG) transcripts, in the absence of ART. Our results highlighted the infiltration of monocyte-derived macrophages (HLA-DR+CD11b+CD14+CD16+) in inflamed liver and lung tissues associated with the expression of CD183 and CX3CR1 but also with markers of tissue-resident macrophages (CD206+ and LYVE+). Sorting of myeloid cell subsets demonstrated that CD14+CD206-, CD14+CD206+, and CD14-CD206+ cell populations were infected, in the liver and lungs, in SIVmac251-infected RMs. Of importance, early ART drastically reduced viral seeding consistent with the absence of ISG detection but also of genes related to inflammation and tissue damage. Viral DNA was only detected in CD206+HLA-DR+CD11b+ cells in ART-treated RMs. The observation of pulmonary and hepatic viral rebound after ART interruption reinforces the importance of early ART implementation to limit viral seeding and inflammatory reactions.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Cis-regulatory hubs: a new 3D model of complex disease genetics with an application to schizophrenia.

The 3D conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases such as schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs, and schizophrenia heritability compared with equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared with promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

DIA proteomics data from a UPS1-spiked E.coli protein mixture processed with six software tools.

In this article, we provide a proteomic reference dataset that has been initially generated for a benchmarking of software tools for Data-Independent Acquisition (DIA) analysis. This large dataset includes 96 DIA .raw files acquired from a complex proteomic standard composed of an E.coli protein background spiked-in with 8 different concentrations of 48 human proteins (UPS1 Sigma). These 8 samples were analyzed in triplicates on an Orbitrap mass spectrometer with 4 different DIA window schemes. We also provide the spectral libraries and FASTA file used for their analysis and the software outputs of the six tools used in this study: DIA-NN, Spectronaut, ScaffoldDIA, DIA-Umpire, Skyline and OpenSWATH. This dataset also contains post-processed quantification tables where the peptides and proteins have been validated, their intensities normalized and the missing values imputed with a noise value. All the files are available on ProteomeXchange. Altogether, these files represent the most comprehensive DIA reference dataset acquired on an Orbitrap instrument ever published. It will be a very useful resource to the proteomic scientists in order to assess the performance of DIA software tools or to test their processing pipelines, to the software developers to improve their tools or develop new ones and to the students for their training on proteomics data analysis.