Search

Abstract Background As part of a National Institute of Allergy and Infectious Diseases funded collaborative project, we have performed over 150 microarray experiments measuring the response of C57/BL6 mouse bone marrow macrophages to toll-like receptor stimuli. These microarray expression profiles are available freely from our project web site http://www.innateImmunity-systemsbiology.org. Here, we report the development of a database of computationally predicted transcription factor binding sites and related genomic features for a set of over 2000 murine immune genes of interest. Our database, which includes microarray co-expression clusters and a host of web-based query, analysis and visualization facilities, is available freely via the internet. It provides a broad resource to the research community, and a stepping stone towards the delineation of the network of transcriptional regulatory interactions underlying the integrated response of macrophages to pathogens. Description We constructed a database indexed on genes and annotations of the immediate surrounding genomic regions. To facilitate both gene-specific and systems biology oriented research, our database provides the means to analyze individual genes or an entire genomic locus. Although our focus to-date has been on mammalian toll-like receptor signaling pathways, our database structure is not limited to this subject, and is intended to be broadly applicable to immunology. By focusing on selected immune-active genes, we were able to perform computationally intensive expression and sequence analyses that would currently be prohibitive if applied to the entire genome. Using six complementary computational algorithms and methodologies, we identified transcription factor binding sites based on the Position Weight Matrices available in TRANSFAC. For one example transcription factor (ATF3) for which experimental data is available, over 50% of our predicted binding sites coincide with genome-wide chromatin immnuopreciptation (ChIP-chip) results. Our database can be interrogated via a web interface. Genomic annotations and binding site predictions can be automatically viewed with a customized version of the Argo genome browser. Conclusion We present the Innate Immune Database (IIDB) as a community resource for immunologists interested in gene regulatory systems underlying innate responses to pathogens. The database website can be freely accessed at http://db.systemsbiology.net/IIDB.

Background: Microsatellite instability (MSI) due to mismatch repair deficiency (dMMR) is common in colorectal cancer (CRC). These cancers are associated with somatic coding events, but the noncoding pathophysiological impact of this genomic instability is yet poorly understood. Here, we perform an analysis of coding and noncoding MSI events at the different steps of colorectal tumorigenesis using whole exome sequencing and search for associated splicing events via RNA sequencing at the bulk-tumor and single-cell levels., Results: Our results demonstrate that MSI leads to hundreds of noncoding DNA mutations, notably at polypyrimidine U2AF RNA-binding sites which are endowed with cis-activity in splicing, while higher frequency of exon skipping events are observed in the mRNAs of MSI compared to non-MSI CRC. At the DNA level, these noncoding MSI mutations occur very early prior to cell transformation in the dMMR colonic crypt, accounting for only a fraction of the exon skipping in MSI CRC. At the RNA level, the aberrant exon skipping signature is likely to impair colonic cell differentiation in MSI CRC affecting the expression of alternative exons encoding protein isoforms governing cell fate, while also targeting constitutive exons, making dMMR cells immunogenic in early stage before the onset of coding mutations. This signature is characterized by its similarity to the oncogenic U2AF1-S34F splicing mutation observed in several other non-MSI cancer., Conclusions: Overall, these findings provide evidence that a very early RNA splicing signature partly driven by MSI impairs cell differentiation and promotes MSI CRC initiation, far before coding mutations which accumulate later during MSI tumorigenesis., (© 2024. The Author(s).)

Background: BMP9 and BMP10 are two major regulators of vascular homeostasis. These two ligands bind with high affinity to the endothelial type I kinase receptor ALK1, together with a type II receptor, leading to the direct phosphorylation of the SMAD transcription factors. Apart from this canonical pathway, little is known. Interestingly, mutations in this signaling pathway have been identified in two rare cardiovascular diseases, hereditary hemorrhagic telangiectasia and pulmonary arterial hypertension., Methods: To get an overview of the signaling pathways modulated by BMP9 and BMP10 stimulation in endothelial cells, we employed an unbiased phosphoproteomic-based strategy. Identified phosphosites were validated by western blot analysis and regulated targets by RT-qPCR. Cell cycle analysis was analyzed by flow cytometry., Results: Large-scale phosphoproteomics revealed that BMP9 and BMP10 treatment induced a very similar phosphoproteomic profile. These BMPs activated a non-canonical transcriptional SMAD-dependent MAPK pathway (MEKK4/P38). We were able to validate this signaling pathway and demonstrated that this activation required the expression of the protein GADD45β. In turn, activated P38 phosphorylated the heat shock protein HSP27 and the endocytosis protein Eps15 (EGF receptor pathway substrate), and regulated the expression of specific genes (E-selectin, hyaluronan synthase 2 and cyclooxygenase 2). This study also highlighted the modulation in phosphorylation of proteins involved in transcriptional regulation (phosphorylation of the endothelial transcription factor ERG) and cell cycle inhibition (CDK4/6 pathway). Accordingly, we found that BMP10 induced a G1 cell cycle arrest and inhibited the mRNA expression of E2F2, cyclinD1 and cyclinA1., Conclusions: Overall, our phosphoproteomic screen identified numerous proteins whose phosphorylation state is impacted by BMP9 and BMP10 treatment, paving the way for a better understanding of the molecular mechanisms regulated by BMP signaling in vascular diseases., (© 2024. The Author(s).)

The development of vascular networks in microfluidic chips is crucial for the long-term culture of three-dimensional cell aggregates such as spheroids, organoids, tumoroids, or tissue explants. Despite rapid advancement in microvascular network systems and organoid technologies, vascularizing organoids-on-chips remains a challenge in tissue engineering. Most existing microfluidic devices poorly reflect the complexity of in vivo flows and require complex technical set-ups. Considering these constraints, we develop a platform to establish and monitor the formation of endothelial networks around mesenchymal and pancreatic islet spheroids, as well as blood vessel organoids generated from pluripotent stem cells, cultured for up to 30 days on-chip. We show that these networks establish functional connections with the endothelium-rich spheroids and vascular organoids, as they successfully provide intravascular perfusion to these structures. We find that organoid growth, maturation, and function are enhanced when cultured on-chip using our vascularization method. This microphysiological system represents a viable organ-on-chip model to vascularize diverse biological 3D tissues and sets the stage to establish organoid perfusions using advanced microfluidics., (© 2024. The Author(s).)

Although the involvement of protein kinase CK2 in cancer is well-documented, there is a need for selective CK2 inhibitors suitable for investigating CK2 specific roles in cancer-related biological pathways and further exploring its therapeutic potential. Here, we report the discovery of AB668, an outstanding selective inhibitor that binds CK2 through a bivalent mode, interacting both at the ATP site and an allosteric αD pocket unique to CK2. Using caspase activation assay, live-cell imaging, and transcriptomic analysis, we have compared the effects of this bivalent inhibitor to representative ATP-competitive inhibitors, CX-4945, and SGC-CK2-1. Our results show that in contrast to CX-4945 or SGC-CK2-1, AB668, by targeting the CK2 αD pocket, has a distinct mechanism of action regarding its anti-cancer activity, inducing apoptotic cell death in several cancer cell lines and stimulating distinct biological pathways in renal cell carcinoma., Competing Interests: The authors A.B. and I.K. declare the following competing financial interest(s): A.B. and I.K. acknowledge greater than 5% ownership of KAIROS DISCOVERY, a company that is developing small molecules as targeted cancer therapies. The other authors declare no competing interests. The authors A.B., T.L., C.C., O.F., and I.K. declare that they are authors of a related patent WO2022008475. The other authors and their immediate family members, have no related patents to declare., (© 2024 The Authors.)

Clear cell Renal Cell Carcinoma (ccRCC) is one of the most prevalent kidney cancers, which is often asymptomatic and thus discovered at a metastatic state (mRCC). mRCC are highly heterogeneous tumors composed of subclonal populations that lead to poor treatment response rate. Several recent works explored the potential of ccRCC tumoroids culture derived from patients. However, these models were produced following a scaffold-based method using collagen I or Matrigel that exhibit lot variability and whose complexity could induce treatment response modifications and phenotypic alterations. Following the observation that ccRCC tumoroids can create their own niche by secreting extracellular matrix components, we developed the first scaffold-free tumoroid model of ccRCC tumors. Tumoroids from mice as well as from human tumors were generated with high success rate (≥90%) using a magnetic suspension method and standard culture media. Immunofluorescence analysis revealed their self-organization capacities to maintain multiple tumor-resident cell types, including endothelial progenitor cells. Transcriptomic analysis showed the reproducibility of the method highlighting that the majority of gene expression patterns was conserved in tumoroids compared to their matching tumor tissue. Moreover, this model enables to evaluate drug effects and invasiveness of renal cancer cells in a 3D context, providing a robust preclinical tool for drug screening and biomarker assessment in line with alternative ex vivo methods like tumor tissue slice culture or in vivo xenograft models., (© 2023. The Author(s).)

Tumor antigens can emerge through multiple mechanisms, including translation of noncoding genomic regions. This noncanonical category of tumor antigens has recently gained attention; however, our understanding of how they recur within and between cancer types is still in its infancy. Therefore, we developed a proteogenomic pipeline based on deep learning de novo mass spectrometry (MS) to enable the discovery of noncanonical MHC class I-associated peptides (ncMAP) from noncoding regions. Considering that the emergence of tumor antigens can also involve posttranslational modifications (PTM), we included an open search component in our pipeline. Leveraging the wealth of MS-based immunopeptidomics, we analyzed data from 26 MHC class I immunopeptidomic studies across 11 different cancer types. We validated the de novo identified ncMAPs, along with the most abundant PTMs, using spectral matching and controlled their FDR to 1%. The noncanonical presentation appeared to be 5 times enriched for the A03 HLA supertype, with a projected population coverage of 55%. The data reveal an atlas of 8,601 ncMAPs with varying levels of cancer selectivity and suggest 17 cancer-selective ncMAPs as attractive therapeutic targets according to a stringent cutoff. In summary, the combination of the open-source pipeline and the atlas of ncMAPs reported herein could facilitate the identification and screening of ncMAPs as targets for T-cell therapies or vaccine development., (©2023 The Authors; Published by the American Association for Cancer Research.)

T cells have the potential to maintain immunological memory and self-tolerance by recognizing antigens from pathogens or tumors. In pathological situations, failure to generate de novo T cells causes immunodeficiency resulting in acute infections and complications. Hematopoietic stem cells (HSC) transplantation constitutes a valuable option to restore proper immune function. However, delayed T cell reconstitution is observed compared to other lineages. To overcome this difficulty, we developed a new approach to identify populations with efficient lymphoid reconstitution properties. To this end, we use a DNA barcoding strategy based on the insertion into a cell chromosome of a lentivirus (LV) carrying a non-coding DNA fragment named barcode (BC). These will segregate through cell divisions and be present in cells' progeny. The remarkable characteristic of the method is that different cell types can be tracked simultaneously in the same mouse. Thus, we in vivo barcoded LMPP and CLP progenitors to test their ability to reconstitute the lymphoid lineage. Barcoded progenitors were co-grafted in immuno-compromised mice and their fate analyzed by evaluating the BC composition in transplanted mice. The results highlight the predominant role of LMPP progenitors for lymphoid generation and reveal valuable novel insights to be reconsidered in clinical transplantation assays.

CK2 is a hetero-tetrameric serine/threonine protein kinase made up of two CK2α/α' catalytic subunits and two CK2β regulatory subunits. The free CK2α subunit and the tetrameric holoenzyme have distinct substrate specificity profiles, suggesting that the spatiotemporal organization of the individual CK2 subunits observed in living cells is crucial in the control of the many cellular processes that are governed by this pleiotropic kinase. Indeed, previous studies reported that the unbalanced expression of CK2 subunits is sufficient to drive epithelial to mesenchymal transition (EMT), a process involved in cancer invasion and metastasis. Moreover, sub-stoichiometric expression of CK2β compared to CK2α in a subset of breast cancer tumors was correlated with the induction of EMT markers and increased epithelial cell plasticity in breast carcinoma progression. Phenotypic changes of epithelial cells are often associated with the activation of phosphotyrosine signaling. Herein, using phosphotyrosine enrichment coupled with affinity capture and proteomic analysis, we show that decreased expression of CK2β in MCF10A mammary epithelial cells triggers the phosphorylation of a number of proteins on tyrosine residues and promotes the striking activation of the FAK1-Src-PAX1 signaling pathway. Moreover, morphometric analyses also reveal that CK2β loss increases the number and the spatial distribution of focal adhesion signaling complexes that coordinate the adhesive and migratory processes. Together, our findings allow positioning CK2β as a gatekeeper for cell spreading by restraining focal adhesion formation and invasion of mammary epithelial cells., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Filhol, Hesse, Bouin, Albigès-Rizo, Jeanneret, Battail, Pflieger and Cochet.)

Aims: BMP9 and BMP10 mutations were recently identified in patients with pulmonary arterial hypertension, but their specific roles in the pathogenesis of the disease are still unclear. We aimed to study the roles of BMP9 and BMP10 in cardiovascular homeostasis and pulmonary hypertension using transgenic mouse models deficient in Bmp9 and/or Bmp10., Methods and Results: Single- and double-knockout mice for Bmp9 (constitutive) and/or Bmp10 (tamoxifen inducible) were generated. Single-knock-out (KO) mice developed no obvious age-dependent phenotype when compared with their wild-type littermates. However, combined deficiency in Bmp9 and Bmp10 led to vascular defects resulting in a decrease in peripheral vascular resistance and blood pressure and the progressive development of high-output heart failure and pulmonary hemosiderosis. RNAseq analysis of the lungs of the double-KO mice revealed differential expression of genes involved in inflammation and vascular homeostasis. We next challenged these mice to chronic hypoxia. After 3 weeks of hypoxic exposure, Bmp10-cKO mice showed an enlarged heart. However, although genetic deletion of Bmp9 in the single- and double-KO mice attenuated the muscularization of pulmonary arterioles induced by chronic hypoxia, we observed no differences in Bmp10-cKO mice. Consistent with these results, endothelin-1 levels were significantly reduced in Bmp9 deficient mice but not Bmp10-cKO mice. Furthermore, the effects of BMP9 on vasoconstriction were inhibited by bosentan, an endothelin receptor antagonist, in a chick chorioallantoic membrane assay., Conclusions: Our data show redundant roles for BMP9 and BMP10 in cardiovascular homeostasis under normoxic conditions (only combined deletion of both Bmp9 and Bmp10 was associated with severe defects) but highlight specific roles under chronic hypoxic conditions. We obtained evidence that BMP9 contributes to chronic hypoxia-induced pulmonary vascular remodelling, whereas BMP10 plays a role in hypoxia-induced cardiac remodelling in mice., (© The Author(s) 2021. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Cutaneous involvement of chronic graft-versus-host disease (cGVHD) has a wide range of manifestations including a lichenoid form with a currently assumed mixed Th1/Th17 signature and a sclerotic form with Th1 signature. Despite substantial heterogeneity of innate and adaptive immune cells recruited to the skin and of the different clinical manifestations, treatment depends mainly on the severity of the skin involvement and relies on systemic, high-dose glucocorticoids alone or in combination with a calcineurin inhibitor. We performed the first study using RNA sequencing to profile and compare the transcriptome of lichen planus cGVHD (n = 8), morphea cGVHD (n = 5), and healthy controls (n = 6). Our findings revealed shared and unique inflammatory pathways to each cGVHD subtype that are both pathogenic and targetable. In particular, the deregulation of IFN signaling pathway was strongly associated with cutaneous cGVHD, whereas the triggering receptor expressed on myeloid cells 1 pathway was found to be specific of lichen planus and likely contributes to its pathogenesis. The results were confirmed at a protein level by performing immunohistochemistry staining and at a transcriptomic level using real-time quantitative polymerase chain reaction., (© 2022 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

Understanding the emergence of lymphoid committed cells from multipotent progenitors (MPP) is a great challenge in hematopoiesis. To gain deeper insight into the dynamic expression changes associated with these transitions, we report the quantitative transcriptome of two MPP subsets and the common lymphoid progenitor (CLP). While the transcriptome is rather stable between MPP2 and MPP3, expression changes increase with differentiation. Among those, we found that pioneer lymphoid genes such as Rag1 , Mpeg1 , and Dntt are expressed continuously from MPP2. Others, such as CD93, are CLP specific, suggesting their potential use as new markers to improve purification of lymphoid populations. Notably, a six-transcription factor network orchestrates the lymphoid differentiation program. Additionally, we pinpointed 24 long intergenic-non-coding RNA (lincRNA) differentially expressed through commitment and further identified seven novel forms. Collectively, our approach provides a comprehensive landscape of coding and non-coding transcriptomes expressed during lymphoid commitment.

The human placenta shares properties with solid tumors, such as rapid growth, tissue invasion, cell migration, angiogenesis, and immune evasion. However, the mechanisms that drive the evolution from premalignant proliferative placental diseases-called hydatidiform moles-to their malignant counterparts, gestational choriocarcinoma, as well as the factors underlying the increased aggressiveness of choriocarcinoma arising after term delivery compared to those developing from hydatidiform moles, are unknown. Using a 730-gene panel covering 13 cancer-associated canonical pathways, we compared the transcriptomic profiles of complete moles to those of postmolar choriocarcinoma samples and those of postmolar to post-term delivery choriocarcinoma. We identified 33 genes differentially expressed between complete moles and postmolar choriocarcinoma, which revealed TGF-β pathway dysregulation. We found the strong expression of SALL4, an upstream regulator of TGF-β, in postmolar choriocarcinoma, compared to moles, in which its expression was almost null. Finally, there were no differentially expressed genes between postmolar and post-term delivery choriocarcinoma samples. To conclude, the TGF-β pathway appears to be a crucial step in the progression of placental malignancies. Further studies should investigate the value of TGF- β family members as biomarkers and new therapeutic targets.

The inflammatory gene NLRP7 is the major gene responsible for recurrent complete hydatidiform moles (CHM), an abnormal pregnancy that can develop into gestational choriocarcinoma (CC). However, the role of NLRP7 in the development and immune tolerance of CC has not been investigated. Three approaches were employed to define the role of NLRP7 in CC development: (i) a clinical study that analyzed human placenta and sera collected from women with normal pregnancies, CHM or CC; (ii) an in vitro study that investigated the impact of NLRP7 knockdown on tumor growth and organization; and (iii) an in vivo study that used two CC mouse models, including an orthotopic model. NLRP7 and circulating inflammatory cytokines were upregulated in tumor cells and in CHM and CC. In tumor cells, NLRP7 functions in an inflammasome-independent manner and promoted their proliferation and 3D organization. Gravid mice placentas injected with CC cells invalidated for NLRP7 , exhibited higher maternal immune response, developed smaller tumors, and displayed less metastases. Our data characterized the critical role of NLRP7 in CC and provided evidence of its contribution to the development of an immunosuppressive maternal microenvironment that not only downregulates the maternal immune response but also fosters the growth and progression of CC.

Background: The duplication of genes is one of the main genetic mechanisms that led to the gain in complexity of biological tissue. Although the implication of duplicated gene expression in brain evolution was extensively studied through comparisons between organs, their role in the regional specialization of the adult human central nervous system has not yet been well described., Results: Our work explored intra-organ expression properties of paralogs through multiple territories of the human central nervous system (CNS) using transcriptome data generated by the Genotype-Tissue Expression (GTEx) consortium. Interestingly, we found that paralogs were associated with region-specific expression in CNS, suggesting their involvement in the differentiation of these territories. Beside the influence of gene expression level on region-specificity, we observed the contribution of both duplication age and duplication type to the CNS region-specificity of paralogs. Indeed, we found that small scale duplicated genes (SSDs) and in particular ySSDs (SSDs younger than the 2 rounds of whole genome duplications) were more CNS region-specific than other paralogs. Next, by studying the two paralogs of ySSD pairs, we observed that when they were region-specific, they tend to be specific to the same region more often than for other paralogs, showing the high co-expression of ySSD pairs. The extension of this analysis to families of paralogs showed that the families with co-expressed gene members (i.e. homogeneous families) were enriched in ySSDs. Furthermore, these homogeneous families tended to be region-specific families, where the majority of their gene members were specifically expressed in the same region., Conclusions: Overall, our study suggests the involvement of ySSDs in the differentiation of human central nervous system territories. Therefore, we show the relevance of exploring region-specific expression of paralogs at the intra-organ level.

Kinase-targeted agents demonstrate antitumor activity in advanced metastatic clear cell renal cell carcinoma (ccRCC), which remains largely incurable. Integration of genomic approaches through small-molecules and genetically based high-throughput screening holds the promise of improved discovery of candidate targets for cancer therapy. The 786-O cell line represents a model for most ccRCC that have a loss of functional pVHL (von Hippel-Lindau). A multiplexed assay was used to study the cellular fitness of a panel of engineered ccRCC isogenic 786-O VHL - cell lines in response to a collection of targeted cancer therapeutics including kinase inhibitors, allowing the interrogation of over 2880 drug-gene pairs. Among diverse patterns of drug sensitivities, investigation of the mechanistic effect of one selected drug combination on tumor spheroids and ex vivo renal tumor slice cultures showed that VHL-defective ccRCC cells were more vulnerable to the combined inhibition of the CK2 and ATM kinases than wild-type VHL cells. Importantly, we found that HIF-2α acts as a key mediator that potentiates the response to combined CK2/ATM inhibition by triggering ROS-dependent apoptosis. Importantly, our findings reveal a selective killing of VHL-deficient renal carcinoma cells and provide a rationale for a mechanism-based use of combined CK2/ATM inhibitors for improved patient care in metastatic VHL-ccRCC.

Transmission distorters (TDs) are genetic elements that favor their own transmission to the detriments of others. Slx/Slxl1 (Sycp3-like-X-linked and Slx-like1) and Sly (Sycp3-like-Y-linked) are TDs, which have been coamplified on the X and Y chromosomes of Mus species. They are involved in an intragenomic conflict in which each favors its own transmission, resulting in sex ratio distortion of the progeny when Slx/Slxl1 versus Sly copy number is unbalanced. They are specifically expressed in male postmeiotic gametes (spermatids) and have opposite effects on gene expression: Sly knockdown leads to the upregulation of hundreds of spermatid-expressed genes, whereas Slx/Slxl1-deficiency downregulates them. When both Slx/Slxl1 and Sly are knocked down, sex ratio distortion and gene deregulation are corrected. Slx/Slxl1 and Sly are, therefore, in competition but the molecular mechanism remains unknown. By comparing their chromatin-binding profiles and protein partners, we show that SLX/SLXL1 and SLY proteins compete for interaction with H3K4me3-reader SSTY1 (Spermiogenesis-specific-transcript-on-the-Y1) at the promoter of thousands of genes to drive their expression, and that the opposite effect they have on gene expression is mediated by different abilities to recruit SMRT/N-Cor transcriptional complex. Their target genes are predominantly spermatid-specific multicopy genes encoded by the sex chromosomes and the autosomal Speer/Takusan. Many of them have coamplified with not only Slx/Slxl1/Sly but also Ssty during muroid rodent evolution. Overall, we identify Ssty as a key element of the X versus Y intragenomic conflict, which may have influenced gene content and hybrid sterility beyond Mus lineage since Ssty amplification on the Y predated that of Slx/Slxl1/Sly., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

A fundamental shortcoming in the current treatment of schizophrenia is the lack of valid criteria to predict who will respond to antipsychotic treatment. The identification of blood-based biological markers of the therapeutic response would enable clinicians to identify the subgroup of patients in whom conventional antipsychotic treatment is ineffective and offer alternative treatments. As part of the Optimisation of Treatment and Management of Schizophrenia in Europe (OPTiMiSE) programme, we conducted an RNA-Seq analysis on 188 subjects with first episode psychosis, all of whom were subsequently treated with amisulpride for 4 weeks. We compared gene expression on total RNA from patients' blood before and after treatment and identified 32 genes for which the expression changed after treatment in good responders only. These findings were replicated in an independent sample of 24 patients with first episode psychosis. Six genes showed a significant difference in expression level between good and poor responders before starting treatment, allowing to predict treatment outcome with a predictive value of 93.8% when combined with clinical features. Collectively, these findings identified new mechanisms to explain symptom improvement after amisulpride medication and highlight the potential of combining gene expression profiling with clinical data to predict treatment response in first episode psychoses.

Epigenetic regulation of gene expression is tightly controlled by the dynamic modification of histones by chemical groups, the diversity of which has largely expanded over the past decade with the discovery of lysine acylations, catalyzed from acyl-coenzymes A. We investigated the dynamics of lysine acetylation and crotonylation on histones H3 and H4 during mouse spermatogenesis. Lysine crotonylation appeared to be of significant abundance compared to acetylation, particularly on Lys27 of histone H3 (H3K27cr) that accumulates in sperm in a cleaved form of H3. We identified the genomic localization of H3K27cr and studied its effects on transcription compared to the classical active mark H3K27ac at promoters and distal enhancers. The presence of both marks was strongly associated with highest gene expression. Assessment of their co-localization with transcription regulators (SLY, SOX30) and chromatin-binding proteins (BRD4, BRDT, BORIS and CTCF) indicated systematic highest binding when both active marks were present and different selective binding when present alone at chromatin. H3K27cr and H3K27ac finally mark the building of some sperm super-enhancers. This integrated analysis of omics data provides an unprecedented level of understanding of gene expression regulation by H3K27cr in comparison to H3K27ac, and reveals both synergistic and specific actions of each histone modification., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

High-throughput RNA-sequencing has become the gold standard method for whole-transcriptome gene expression analysis, and is widely used in numerous applications to study cell and tissue transcriptomes. It is also being increasingly used in a number of clinical applications, including expression profiling for diagnostics and alternative transcript detection. However, despite its many advantages, RNA sequencing can be challenging in some situations, for instance in cases of low input amounts or degraded RNA samples. Several protocols have been proposed to overcome these challenges, and many are available as commercial kits. In this study, we systematically test three recent commercial technologies for RNA-seq library preparation (TruSeq, SMARTer and SMARTer Ultra-Low) on human biological reference materials, using standard (1 mg), low (100 ng and 10 ng) and ultra-low (<1 ng) input amounts, and for mRNA and total RNA, stranded and unstranded. The results are analyzed using read quality and alignment metrics, gene detection and differential gene expression metrics. Overall, we show that the TruSeq kit performs well with an input amount of 100 ng, while the SMARTer kit shows decreased performance for inputs of 100 and 10 ng, and the SMARTer Ultra-Low kit performs relatively well for input amounts <1 ng. All the results are discussed in detail, and we provide guidelines for biologists for the selection of an RNA-seq library preparation kit.

Elucidating population structure and levels of genetic diversity and recombination is necessary to understand the evolution and adaptation of species. Candida albicans is the second most frequent agent of human fungal infections worldwide, causing high-mortality rates. Here we present the genomic sequences of 182 C. albicans isolates collected worldwide, including commensal isolates, as well as ones responsible for superficial and invasive infections, constituting the largest dataset to date for this major fungal pathogen. Although, C. albicans shows a predominantly clonal population structure, we find evidence of gene flow between previously known and newly identified genetic clusters, supporting the occurrence of (para)sexuality in nature. A highly clonal lineage, which experimentally shows reduced fitness, has undergone pseudogenization in genes required for virulence and morphogenesis, which may explain its niche restriction. Candida albicans thus takes advantage of both clonality and gene flow to diversify.

The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima's D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution., (© The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are autosomal dominant neuromuscular diseases caused by microsatellite expansions and belong to the family of RNA-dominant disorders. Availability of cellular models in which the DM mutation is expressed within its natural context is essential to facilitate efforts to identify new therapeutic compounds. Here, we generated immortalized DM1 and DM2 human muscle cell lines that display nuclear RNA aggregates of expanded repeats, a hallmark of myotonic dystrophy. Selected clones of DM1 and DM2 immortalized myoblasts behave as parental primary myoblasts with a reduced fusion capacity of immortalized DM1 myoblasts when compared with control and DM2 cells. Alternative splicing defects were observed in differentiated DM1 muscle cell lines, but not in DM2 lines. Splicing alterations did not result from differentiation delay because similar changes were found in immortalized DM1 transdifferentiated fibroblasts in which myogenic differentiation has been forced by overexpression of MYOD1. As a proof-of-concept, we show that antisense approaches alleviate disease-associated defects, and an RNA-seq analysis confirmed that the vast majority of mis-spliced events in immortalized DM1 muscle cells were affected by antisense treatment, with half of them significantly rescued in treated DM1 cells. Immortalized DM1 muscle cell lines displaying characteristic disease-associated molecular features such as nuclear RNA aggregates and splicing defects can be used as robust readouts for the screening of therapeutic compounds. Therefore, immortalized DM1 and DM2 muscle cell lines represent new models and tools to investigate molecular pathophysiological mechanisms and evaluate the in vitro effects of compounds on RNA toxicity associated with myotonic dystrophy mutations., (© 2017. Published by The Company of Biologists Ltd.)

Bromodomain and Extra-terminal motif (BET) proteins play a central role in transcription regulation and chromatin signalling pathways. They are present in unicellular eukaryotes and in this study, the role of the BET protein Bdf1 has been explored in Saccharomyces cerevisiae. Mutation of Bdf1 bromodomains revealed defects on both the formation of spores and the meiotic progression, blocking cells at the exit from prophase, before the first meiotic division. This phenotype is associated with a massive deregulation of the transcription of meiotic genes and Bdf1 bromodomains are required for appropriate expression of the key meiotic transcription factor NDT80 and almost all the Ndt80-inducible genes, including APC complex components. Bdf1 notably accumulates on the promoter of Ndt80 and its recruitment is dependent on Bdf1 bromodomains. In addition, the ectopic expression of NDT80 during meiosis partially bypasses this dependency. Finally, purification of Bdf1 partners identified two independent complexes with Bdf2 or the SWR complex, neither of which was required to complete sporulation. Taken together, our results unveil a new role for Bdf1 -working independently from its predominant protein partners Bdf2 and the SWR1 complex-as a regulator of meiosis-specific genes., Competing Interests: The authors have declared that no competing interests exist.

LHP1-INTERACTING FACTOR2 (LIF2), a heterogeneous nuclear ribonucleoprotein involved in Arabidopsis thaliana cell fate and stress responses, interacts with LIKE HETEROCHROMATIN PROTEIN1 (LHP1), a Polycomb Repressive Complex1 subunit. To investigate LIF2-LHP1 functional interplay, we mapped their genome-wide distributions in wild-type, lif2 , and lhp1 backgrounds, under standard and stress conditions. Interestingly, LHP1-targeted regions form local clusters, suggesting an underlying functional organization of the plant genome. Regions targeted by both LIF2 and LHP1 were enriched in stress-responsive genes, the H2A.Z histone variant, and antagonistic histone marks. We identified specific motifs within the targeted regions, including a G-box-like motif, a GAGA motif, and a telo -box. LIF2 and LHP1 can operate both antagonistically and synergistically. In response to methyl jasmonate treatment, LIF2 was rapidly recruited to chromatin, where it mediated transcriptional gene activation. Thus, LIF2 and LHP1 participate in transcriptional switches in stress-response pathways., (© 2016 American Society of Plant Biologists. All rights reserved.)

Major physiological changes are governed by alternative splicing of RNA, and its misregulation may lead to specific diseases. With the use of a genome-wide approach, we show here that this splicing step can be modified by medication and demonstrate the effects of the biguanide metformin, on alternative splicing. The mechanism of action involves AMPK activation and downregulation of the RBM3 RNA-binding protein. The effects of metformin treatment were tested on myotonic dystrophy type I (DM1), a multisystemic disease considered to be a spliceopathy. We show that this drug promotes a corrective effect on several splicing defects associated with DM1 in derivatives of human embryonic stem cells carrying the causal mutation of DM1 as well as in primary myoblasts derived from patients. The biological effects of metformin were shown to be compatible with typical therapeutic dosages in a clinical investigation involving diabetic patients. The drug appears to act as a modifier of alternative splicing of a subset of genes and may therefore have novel therapeutic potential for many more diseases besides those directly linked to defective alternative splicing.

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement., (Copyright © 2014, American Association for the Advancement of Science.)

The duplication of mammalian genomes is under the control of a spatiotemporal program that orchestrates the positioning and the timing of firing of replication origins. The molecular mechanisms coordinating the activation of about [Formula: see text] predicted origins remain poorly understood, partly due to the intrinsic rarity of replication bubbles, making it difficult to purify short nascent strands (SNS). The precise identification of origins based on the high-throughput sequencing of SNS constitutes a new methodological challenge. We propose a new statistical method with a controlled resolution, adapted to the detection of replication origins from SNS data. We detected an average of 80,000 replication origins in different cell lines. To evaluate the consistency between different protocols, we compared SNS detections with bubble trapping detections. This comparison demonstrated a good agreement between genome-wide methods, with 65% of SNS-detected origins validated by bubble trapping, and 44% of bubble trapping origins validated by SNS origins, when compared at the same resolution. We investigated the interplay between the spatial and the temporal programs of replication at fine scales. We show that most of the origins detected in regions replicated in early S phase are shared by all the cell lines investigated whereas cell-type-specific origins tend to be replicated in late S phase. We shed a new light on the key role of CpG islands, by showing that 80% of the origins associated with CGIs are constitutive. Our results further show that at least 76% of CGIs are origins of replication. The analysis of associations with chromatin marks at different timing of cell division revealed new potential epigenetic regulators driving the spatiotemporal activity of replication origins. We highlight the potential role of H4K20me1 and H3K27me3, the coupling of which is correlated with increased efficiency of replication origins, clearly identifying those marks as potential key regulators of replication origins.

Background: As part of a National Institute of Allergy and Infectious Diseases funded collaborative project, we have performed over 150 microarray experiments measuring the response of C57/BL6 mouse bone marrow macrophages to toll-like receptor stimuli. These microarray expression profiles are available freely from our project web site http://www.innateImmunity-systemsbiology.org. Here, we report the development of a database of computationally predicted transcription factor binding sites and related genomic features for a set of over 2000 murine immune genes of interest. Our database, which includes microarray co-expression clusters and a host of web-based query, analysis and visualization facilities, is available freely via the internet. It provides a broad resource to the research community, and a stepping stone towards the delineation of the network of transcriptional regulatory interactions underlying the integrated response of macrophages to pathogens., Description: We constructed a database indexed on genes and annotations of the immediate surrounding genomic regions. To facilitate both gene-specific and systems biology oriented research, our database provides the means to analyze individual genes or an entire genomic locus. Although our focus to-date has been on mammalian toll-like receptor signaling pathways, our database structure is not limited to this subject, and is intended to be broadly applicable to immunology. By focusing on selected immune-active genes, we were able to perform computationally intensive expression and sequence analyses that would currently be prohibitive if applied to the entire genome. Using six complementary computational algorithms and methodologies, we identified transcription factor binding sites based on the Position Weight Matrices available in TRANSFAC. For one example transcription factor (ATF3) for which experimental data is available, over 50% of our predicted binding sites coincide with genome-wide chromatin immnuopreciptation (ChIP-chip) results. Our database can be interrogated via a web interface. Genomic annotations and binding site predictions can be automatically viewed with a customized version of the Argo genome browser., Conclusion: We present the Innate Immune Database (IIDB) as a community resource for immunologists interested in gene regulatory systems underlying innate responses to pathogens. The database website can be freely accessed at http://db.systemsbiology.net/IIDB.