Explore the vast range of titles available.

Background Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous (“multimappers”). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. Results In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. Conclusion Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines –currently restricted to specific contexts or communities– to ensure the reliability of genomic and transcriptomic studies.

In most metazoan nuclei, heterochromatin is located at the nuclear periphery in contact with the nuclear lamina, which provides mechanical stability to the nucleus. We show that in cultured cells, chromatin decompaction by the nucleosome binding protein HMGN5 decreases the sturdiness, elasticity and rigidity of the nucleus. Mice overexpressing HMGN5, either globally or only in the heart, are normal at birth but develop hypertrophic heart with large cardiomyoctyes, deformed nuclei and disrupted lamina and die of cardiac malfunction. Chromatin decompaction is seen in cardiomyocytes of newborn mice but misshaped nuclei with disrupted lamina are seen only in adult cardiomyocytes, suggesting that loss of heterochromatin diminishes the ability of the nucleus to withstand the mechanical forces of the contracting heart. Thus, heterochromatin enhances the ability of the nuclear lamina to maintain the sturdiness and shape of the eukaryotic nucleus; a structural role for chromatin that is distinct from its genetic functions. Whether heterochromatin affects the physical properties of the nucleus is not known. Here, Furusawa et al . show that chromatin decompaction decreases the sturdiness of the nucleus in cultured cells and leads to lamina disruption and cardiac abnormalities in adult mice, suggesting a structural, non-genetic function for heterochromatin.

Nadav Ahituv and colleagues use a massively parallel reporter assay to test 4,970 synthetic regulatory element sequences, containing patterns of 12 known liver transcription factor binding sites, in mice and in HepG2 cells. They systematically test the impact of binding site copy number, spacing, combination and order on gene expression. Despite continual progress in the cataloging of vertebrate regulatory elements, little is known about their organization and regulatory architecture. Here we describe a massively parallel experiment to systematically test the impact of copy number, spacing, combination and order of transcription factor binding sites on gene expression. A complex library of ∼5,000 synthetic regulatory elements containing patterns from 12 liver-specific transcription factor binding sites was assayed in mice and in HepG2 cells. We find that certain transcription factors act as direct drivers of gene expression in homotypic clusters of binding sites, independent of spacing between sites, whereas others function only synergistically. Heterotypic enhancers are stronger than their homotypic analogs and favor specific transcription factor binding site combinations, mimicking putative native enhancers. Exhaustive testing of binding site permutations suggests that there is flexibility in binding site order. Our findings provide quantitative support for a flexible model of regulatory element activity and suggest a framework for the design of synthetic tissue-specific enhancers.

MicroRNAs (miRNAs) have been reported to contribute to the pathophysiology of multiple sclerosis (MS), an inflammatory disorder of the central nervous system. Here, we propose a new consensus-based strategy to analyse and integrate miRNA and gene expression data in MS as well as other publically available data to gain a deeper understanding of the role of miRNAs in MS and to overcome the challenges posed by studies with limited patient sample sizes. We processed and analysed microarray datasets and compared the expression of genes and miRNAs in the blood of MS patients and controls. We then used our consensus and integration approach to construct two molecular networks dysregulated in MS: a miRNA- and a gene-based network. We identified 18 differentially expressed (DE) miRNAs and 128 DE genes that may contribute to the regulatory alterations behind MS. The miRNAs were linked to immunological and neurological pathways and we exposed let-7b-5p and miR-345-5p as promising blood-derived disease biomarkers in MS. The results suggest that DE miRNAs are more informative than DE genes in uncovering pathways potentially involved in MS. Our findings provide novel insights into the regulatory mechanisms and networks underlying MS.

Background Despite the advent of Chromatin Immunoprecipitation Sequencing (ChIP-seq) having revolutionised our understanding of the mammalian genome’s regulatory landscape, many challenges remain. In particular, because of their repetitive nature, the sequencing reads derived from transposable elements (TEs) pose a real bioinformatics challenge, to the point that standard analysis pipelines typically ignore reads whose genomic origin cannot be unambiguously ascertained. Results We show that discarding ambiguously mapping reads may lead to a systematic underestimation of the number of reads associated with young TE families/subfamilies. We also provide evidence suggesting that the strategy of randomly permuting the location of the read mappings (or the TEs) that is often used to compute the background for enrichment calculations at TE families/subfamilies can result in both false positive and negative enrichments. To address these problems, we present the Transposable Element Enrichment Estimator (T3E), a tool that makes use of ChIP-seq data to characterise the epigenetic profile of associated TE families/subfamilies. T3E weights the number of read mappings assigned to the individual TE copies of a family/subfamily by the overall number of genomic loci to which the corresponding reads map, and this is done at the single nucleotide level. In addition, T3E computes ChIP-seq enrichment relative to a background estimated based on the distribution of the read mappings in the input control DNA. We demonstrated the capabilities of T3E on 23 different ChIP-seq libraries. T3E identified enrichments that were consistent with previous studies. Furthermore, T3E detected context-specific enrichments that are likely to pinpoint unexplored TE families/subfamilies with individual TE copies that have been frequently exapted as cis -regulatory elements during the evolution of mammalian regulatory networks. Conclusions T3E is a novel open-source computational tool (available for use at: https://github.com/michelleapaz/T3E ) that overcomes some of the pitfalls associated with the analysis of ChIP-seq data arising from the repetitive mammalian genome and provides a framework to shed light on the epigenetics of entire TE families/subfamilies.

A key problem in development is to understand how genes turn on or off at the right place and right time during embryogenesis. Such decisions are made by non-coding sequences called ‘enhancers.’ Much of our models of how enhancers work rely on the assumption that genes are activated de novo as stable domains across embryonic tissues. Such a view has been strengthened by the intensive landmark studies of the early patterning of the anterior-posterior (AP) axis of the Drosophila embryo, where indeed gene expression domains seem to arise more or less stably. However, careful analysis of gene expression patterns in other model systems (including the AP patterning in vertebrates and short-germ insects like the beetle Tribolium castaneum ) painted a different, very dynamic view of gene regulation, where genes are oftentimes expressed in a wavelike fashion. How such gene expression waves are mediated at the enhancer level is so far unclear. Here, we establish the AP patterning of the short-germ beetle Tribolium as a model system to study dynamic and temporal pattern formation at the enhancer level. To that end, we established an enhancer prediction system in Tribolium based on time- and tissue-specific ATAC-seq and an enhancer live reporter system based on MS2 tagging. Using this experimental framework, we discovered several Tribolium enhancers, and assessed the spatiotemporal activities of some of them in live embryos. We found our data consistent with a model in which the timing of gene expression during embryonic pattern formation is mediated by a balancing act between enhancers that induce rapid changes in gene expression patterns (that we call ‘dynamic enhancers’) and enhancers that stabilize gene expression patterns (that we call ‘static enhancers’). However, more data is needed for a strong support for this or any other alternative models.

Superovulation is the epitome for generating oocytes for molecular embryology in mice, and it is used to model medically assisted reproduction in humans. However, whether a superovulated oocyte is normal, is an open question. This study establishes for the first time that superovulation is associated with proteome changes that affect phenotypic traits in mice, whereas the transcriptome is far less predictive. The proteins that were differentially expressed in superovulated mouse oocytes and embryos compared to their naturally ovulated counterparts were enriched in ontology terms describing abnormal mammalian phenotypes: a thinner zona pellucida, a smaller oocyte diameter, increased frequency of cleavage arrest, and defective blastocyst formation, which could all be verified functionally. Moreover, our findings indicate that embryos with such abnormalities are negatively selected during preimplantation, and ascribe these abnormalities to incomplete ovarian maturation during the time of the conventional superovulation, since they could be corrected upon postponement of the ovulatory stimulus by 24 h. Our data place constraints on the common view that superovulated oocytes are suitable for drawing general conclusions about developmental processes, and underscore the importance of including the proteins in a modern molecular definition of oocyte quality.

Innate lymphoid cells (ILC) not only are responsible for shaping the innate immune response but also actively modulate T cell responses. However, the molecular processes regulating ILC-T cell interaction are not yet completely understood. The protein butyrophilin 2a2 (Btn2a2), a co-stimulatory molecule first identified on antigen-presenting cells, has a pivotal role in the maintenance of T cell homeostasis, but the main effector cell and the respective ligands remain elusive. We analyzed the role of Btn2a2 in the ILC-T cell cross talk. We found that the expression of Btn2a2 is upregulated in ILC2 following stimulation with IL-33/IL-25/TSLP. In vitro and in vivo experiments indicated that lack of Btn2a2 expression on ILC2 resulted in elevated T cell responses. We observed an enhanced proliferation of T cells as well as increased secretion of the type 2 cytokines IL-4/IL-5/IL-13 following cocultures with Btn2a2-deficient ILC2. In vivo transfer experiments confirmed the regulatory role of Btn2a2 on ILC2 as Btn2a2-deficient ILC2 induced stronger T cell responses and prevented chronic helminth infections. Taken together, we identified Btn2a2 as a significant player in the regulation of ILC2–T cell interactions.

The molecular heterogeneity of feline mammary carcinomas (FMCs) represents a prognostic and therapeutic challenge. RNA-Seq-based comparative transcriptomic profiling serves to identify recurrent and exclusive differentially expressed genes (DEGs) across sample types and molecular subtypes. Using mass-parallel RNA-Seq, we identified DEGs and performed comparative function-based analysis across 15 tumours (four basal-like triple-negative [TN], eight normal-like TN, and three luminal B f HER2 negative [LB f HER2−]), two cell lines (CL, TiHo-0906, and TiHo-1403) isolated from the primary tumours (LB f HER2−) of two cats included in this study, and 13 healthy mammary tissue controls. DEGs in tumours were predominantly upregulated; dysregulation of CLs transcriptome was more extensive, including mostly downregulated genes. Cell-cycle and metabolic-related DEGs were upregulated in both tumours and CLs, including therapeutically-targetable cell cycle regulators (e.g. CCNB1 , CCNB2 , CDK1 , CDK4 , GTSE1 , MCM4 , and MCM5) , metabolic-related genes (e.g. FADS2 and SLC16A3 ), heat-shock proteins (e.g. HSPH1, HSP90B1 , and HSPA5 ), genes controlling centrosome disjunction (e.g. RACGAP1 and NEK2 ), and collagen molecules (e.g. COL2A1 ). DEGs specifically upregulated in basal-like TN tumours were involved in antigen processing and presentation, in normal-like TN tumours encoded G protein-coupled receptors (GPCRs), and in LB f HER2− tumours were associated with lysosomes, phagosomes, and endosomes formation. Downregulated DEGs in CLs were associated with structural and signalling cell surface components. Hence, our results suggest that upregulation of genes enhancing proliferation and metabolism is a common feature among FMCs and derived CLs. In contrast, the dissimilarities observed in dysregulation of membrane components highlight CLs’ disconnection with the tumour microenvironment. Furthermore, recurrent and exclusive DEGs associated with dysregulated pathways might be useful for the development of prognostically and therapeutically-relevant targeted panels.

Detailed information about stage-specific changes in gene expression is crucial for understanding the gene regulatory networks underlying development and the various signal transduction pathways contributing to morphogenesis. Here we describe the global gene expression dynamics during early murine limb development, when cartilage, tendons, muscle, joints, vasculature and nerves are specified and the musculoskeletal system of limbs is established. We used whole-genome microarrays to identify genes with differential expression at 5 stages of limb development (E9.5 to 13.5), during fore- and hind-limb patterning. We found that the onset of limb formation is characterized by an up-regulation of transcription factors, which is followed by a massive activation of genes during E10.5 and E11.5 which levels off at later time points. Among the 3520 genes identified as significantly up-regulated in the limb, we find ~30% to be novel, dramatically expanding the repertoire of candidate genes likely to function in the limb. Hierarchical and stage-specific clustering identified expression profiles that are likely to correlate with functional programs during limb development and further characterization of these transcripts will provide new insights into specific tissue patterning processes. Here, we provide for the first time a comprehensive analysis of developmentally regulated genes during murine limb development, and provide some novel insights into the expression dynamics governing limb morphogenesis.