Motivation: Statistically assessing the relation between a set of genomic regions and other genomic features is a common challenging task in genomic and epigenomic analyses. Randomization based approaches implicitly take into account the complexity of the genome without the need of assuming an underlying statistical model. Summary: regioneR is an R package that implements a permutation test framework specifically designed to work with genomic regions. In addition to the predefined randomization and evaluation strategies, regioneR is fully customizable allowing the use of custom strategies to adapt it to specific questions. Finally, it also implements a novel function to evaluate the local specificity of the detected association. Availability and implementation: regioneR is an R package released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor ( http://www.bioconductor.org/packages/regioneR ). Contact:rmalinverni@carrerasresearch.org

Porcine reproductive and respiratory syndrome is a major cause of economic loss for the swine industry worldwide. Porcine reproductive and respiratory syndrome virus (PRRSV) triggers weak and atypical innate immune responses, but key genes and mechanisms by which the virus interferes with the host innate immunity have not yet been elucidated. In this study, genes that control the response of the main target of PRRSV, porcine alveolar macrophages (PAMs), were profiledin vitrowith a time-course experiment spanning the first round of virus replication. PAMs were obtained from six piglets and challenged with the Lelystad PRRSV strain, and gene expression was investigated using Affymetrix microarrays and real-time PCR. Of the 1409 differentially expressed transcripts identified by analysis of variance, two, five, 25, 16 and 100 differed from controls by a minimum of 1.5-fold at 1, 3, 6, 9 and 12 h post-infection (p.i.), respectively. A PRRSV infection effect was detectable between 3 and 6 h p.i., and was characterized by a consistent downregulation of gene expression, followed by the start of the host innate immune response at 9 h p.i. The expression of beta interferon 1 (IFN-β), but not ofIFN-α, was strongly upregulated, whilst few genes commonly expressed in response to viral infections and/or induced by interferons were found to be differentially expressed. A predominance of anti-apoptotic transcripts (e.g. interleukin-10), a shift towards a T-helper cell type 2 response and a weak upregulation of tumour necrosis factor-αexpression were observed within 12 h p.i., reinforcing the hypotheses that PRRSV has developed sophisticated mechanisms to escape the host defence.

Genetic loss-of-function studies on development, cancer and somatic cell reprogramming have suggested that the group of macroH2A histone variants might function through stabilizing the differentiated state by a yet unknown mechanism. Here, we present results demonstrating that macroH2A variants have a major function in maintaining nuclear organization and heterochromatin architecture. Specifically, we find that a substantial amount of macroH2A is associated with heterochromatic repeat sequences. We further identify macroH2A on sites of interstitial heterochromatin decorated by histone H3 trimethylated on K9 (H3K9me3). Loss of macroH2A leads to major defects in nuclear organization, including reduced nuclear circularity, disruption of nucleoli and a global loss of dense heterochromatin. Domains formed by DNA repeat sequences are disorganized, expanded and fragmented, and mildly re-expressed when depleted of macroH2A. At the molecular level, we find that macroH2A is required for the interaction of repeat sequences with the nucleostructural protein lamin B1. Taken together, our results argue that a major function of macroH2A histone variants is to link nucleosome composition to higher-order chromatin architecture.

Polycomb proteins play an essential role in maintaining the repression of developmental genes in self-renewing embryonic stem cells. The exact mechanism allowing the derepression of polycomb target genes during cell differentiation remains unclear. Our project aimed to identify Cbx8 binding sites in differentiating mouse embryonic stem cells. Therefore, we used a genome-wide chromatin immunoprecipitation of endogenous Cbx8 coupled to direct massive parallel sequencing (ChIP-Seq). Our analysis identified 171 high confidence peaks. By crossing our data with previously published microarray analysis, we show that several differentiation genes transiently recruit Cbx8 during their early activation. Depletion of Cbx8 partially impairs the transcriptional activation of these genes. Both interaction analysis, as well as chromatin immunoprecipitation experiments support the idea that activating Cbx8 acts in the context of an intact PRC1 complex. Prolonged gene activation results in eviction of PRC1 despite persisting H3K27me3 and H2A ubiquitination. The composition of PRC1 is highly modular and changes when embryonic stem cells commit to differentiation. We further demonstrate that the exchange of Cbx7 for Cbx8 is required for the effective activation of differentiation genes. Taken together, our results establish a function for a Cbx8-containing complex in facilitating the transition from a Polycomb-repressed chromatin state to an active state. As this affects several key regulatory differentiation genes this mechanism is likely to contribute to the robust execution of differentiation programs. Cell fate transitions have long been known to be accompanied by alterations in chromatin structure. But only during the last few years has it become clear that chromatin modifications form the molecular basis of an epigenetic memory that defines cell identity. The Polycomb Group Proteins (PcGs) form two major protein complexes known as polycomb repressive complexes 1 and 2 (PRC1 and PRC2). Their function is essential for the maintenance of transcriptional repression during embryogenesis through the methylation of the lysine 27 on histone H3 and the subsequent ubiquitination of histone H2A. The chromobox homolog 8, Cbx8, which is part of the PRC1 complex, is therefore generally defined as a repressor of gene transcription. The genome wide profiling of Cbx8 during the early steps of mouse embryonic stem (mES) cells differentiation provided us with surprising results involving Cbx8 in gene activation. Our results point out that Cbx8 is part of a PRC1 complex involved in the transition from a Polycomb repressed state to an active state.

Histone variants are structural components of eukaryotic chromatin that can replace replication-coupled histones in the nucleosome. The histone variant macroH2A1.1 contains a macrodomain capable of binding NAD⁺-derived metabolites. Here we report that macroH2A1.1 is rapidly induced during myogenic differentiation through a switch in alternative splicing, and that myotubes that lack macroH2A1.1 have a defect in mitochondrial respiratory capacity. We found that the metabolite-binding macrodomain was essential for sustained optimal mitochondrial function but dispensable for gene regulation. Through direct binding, macroH2A1.1 inhibits basal poly-ADP ribose polymerase 1 (PARP-1) activity and thus reduces nuclear NAD⁺ consumption. The resultant accumulation of the NAD⁺ precursor NMN allows for maintenance of mitochondrial NAD⁺ pools that are critical for respiration. Our data indicate that macroH2A1.1-containing chromatin regulates mitochondrial respiration by limiting nuclear NAD⁺ consumption and establishing a buffer of NAD⁺ precursors in differentiated cells.

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.