Chromosome-level organization of the regulatory genome in the Drosophila nervous system.

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

A role for mitotic bookmarking of SOX2 in pluripotency and differentiation.

Mitotic bookmarking transcription factors remain bound to chromosomes during mitosis and were proposed to regulate phenotypic maintenance of stem and progenitor cells at the mitosis-to-G1 (M-G1) transition. However, mitotic bookmarking remains largely unexplored in most stem cell types, and its functional relevance for cell fate decisions remains unclear. Here we screened for mitotic chromosome binding within the pluripotency network of embryonic stem (ES) cells and show that SOX2 and OCT4 remain bound to mitotic chromatin through their respective DNA-binding domains. Dynamic characterization using photobleaching-based methods and single-molecule imaging revealed quantitatively similar specific DNA interactions, but different nonspecific DNA interactions, of SOX2 and OCT4 with mitotic chromatin. Using ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) to assess the genome-wide distribution of SOX2 on mitotic chromatin, we demonstrate the bookmarking activity of SOX2 on a small set of genes. Finally, we investigated the function of SOX2 mitotic bookmarking in cell fate decisions and show that its absence at the M-G1 transition impairs pluripotency maintenance and abrogates its ability to induce neuroectodermal differentiation but does not affect reprogramming efficiency toward induced pluripotent stem cells. Our study demonstrates the mitotic bookmarking property of SOX2 and reveals its functional importance in pluripotency maintenance and ES cell differentiation.

Cohesins functionally associate with CTCF on mammalian chromosome arms.

Cohesins mediate sister chromatid cohesion, which is essential for chromosome segregation and postreplicative DNA repair. In addition, cohesins appear to regulate gene expression and enhancer-promoter interactions. These noncanonical functions remained unexplained because knowledge of cohesin-binding sites and functional interactors in metazoans was lacking. We show that the distribution of cohesins on mammalian chromosome arms is not driven by transcriptional activity, in contrast to S. cerevisiae. Instead, mammalian cohesins occupy a subset of DNase I hypersensitive sites, many of which contain sequence motifs resembling the consensus for CTCF, a DNA-binding protein with enhancer blocking function and boundary-element activity. We find cohesins at most CTCF sites and show that CTCF is required for cohesin localization to these sites. Recruitment by CTCF suggests a rationale for noncanonical cohesin functions and, because CTCF binding is sensitive to DNA methylation, allows cohesin positioning to integrate DNA sequence and epigenetic state.

The process of Lewy body formation, rather than simply α-synuclein fibrillization, is one of the major drivers of neurodegeneration.

Parkinson's disease (PD) is characterized by the accumulation of misfolded and aggregated α-synuclein (α-syn) into intraneuronal inclusions named Lewy bodies (LBs). Although it is widely believed that α-syn plays a central role in the pathogenesis of PD, the processes that govern α-syn fibrillization and LB formation remain poorly understood. In this work, we sought to dissect the spatiotemporal events involved in the biogenesis of the LBs at the genetic, molecular, biochemical, structural, and cellular levels. Toward this goal, we further developed a seeding-based model of α-syn fibrillization to generate a neuronal model that reproduces the key events leading to LB formation, including seeding, fibrillization, and the formation of inclusions that recapitulate many of the biochemical, structural, and organizational features of bona fide LBs. Using an integrative omics, biochemical and imaging approach, we dissected the molecular events associated with the different stages of LB formation and their contribution to neuronal dysfunction and degeneration. In addition, we demonstrate that LB formation involves a complex interplay between α-syn fibrillization, posttranslational modifications, and interactions between α-syn aggregates and membranous organelles, including mitochondria, the autophagosome, and endolysosome. Finally, we show that the process of LB formation, rather than simply fibril formation, is one of the major drivers of neurodegeneration through disruption of cellular functions and inducing mitochondria damage and deficits, and synaptic dysfunctions. We believe that this model represents a powerful platform to further investigate the mechanisms of LB formation and clearance and to screen and evaluate therapeutics targeting α-syn aggregation and LB formation.

AMPK promotes survival of c-Myc-positive melanoma cells by suppressing oxidative stress.

Although c-Myc is essential for melanocyte development, its role in cutaneous melanoma, the most aggressive skin cancer, is only partly understood. Here we used the NrasQ61KINK4a-/- mouse melanoma model to show that c-Myc is essential for tumor initiation, maintenance, and metastasis. c-Myc-expressing melanoma cells were preferentially found at metastatic sites, correlated with increased tumor aggressiveness and high tumor initiation potential. Abrogation of c-Myc caused apoptosis in primary murine and human melanoma cells. Mechanistically, c-Myc-positive melanoma cells activated and became dependent on the metabolic energy sensor AMP-activated protein kinase (AMPK), a metabolic checkpoint kinase that plays an important role in energy and redox homeostasis under stress conditions. AMPK pathway inhibition caused apoptosis of c-Myc-expressing melanoma cells, while AMPK activation protected against cell death of c-Myc-depleted melanoma cells through suppression of oxidative stress. Furthermore, TCGA database analysis of early-stage human melanoma samples revealed an inverse correlation between C-MYC and patient survival, suggesting that C-MYC expression levels could serve as a prognostic marker for early-stage disease.

Single-cell transcriptional logic of cell-fate specification and axon guidance in early-born retinal neurons.

Retinal ganglion cells (RGCs), cone photoreceptors (cones), horizontal cells and amacrine cells are the first classes of neurons produced in the retina. However, an important question is how this diversity of cell states is transcriptionally produced. Here, we profiled 6067 single retinal cells to provide a comprehensive transcriptomic atlas showing the diversity of the early developing mouse retina. RNA velocities unveiled the dynamics of cell cycle coordination of early retinogenesis and define the transcriptional sequences at work during the hierarchical production of early cell-fate specification. We show that RGC maturation follows six waves of gene expression, with older-generated RGCs transcribing increasing amounts of guidance cues for young peripheral RGC axons that express the matching receptors. Spatial transcriptionally deduced features in subpopulations of RGCs allowed us to define novel molecular markers that are spatially restricted. Finally, the isolation of such a spatially restricted population, ipsilateral RGCs, allowed us to identify their molecular identity at the time they execute axon guidance decisions. Together, these data represent a valuable resource shedding light on transcription factor sequences and guidance cue dynamics during mouse retinal development.

Clustering of mammalian Hox genes with other H3K27me3 targets within an active nuclear domain.

Embryogenesis requires the precise activation and repression of many transcriptional regulators. The Polycomb group proteins and the associated H3K27me3 histone mark are essential to maintain the inactive state of many of these genes. Mammalian Hox genes are targets of Polycomb proteins and form local 3D clusters centered on the H3K27me3 mark. More distal contacts have also been described, yet their selectivity, dynamics, and relation to other layers of chromatin organization remained elusive. We report that repressed Hox genes form mutual intra- and interchromosomal interactions with other genes located in strong domains labeled by H3K27me3. These interactions occur in a central and active nuclear environment that consists of the HiC compartment A, away from peripheral lamina-associated domains. Interactions are independent of nearby H3K27me3-marked loci and determined by chromosomal distance and cell-type-specific scaling factors, thus inducing a moderate reorganization during embryogenesis. These results provide a simplified view of nuclear organization whereby Polycomb proteins may have evolved to repress genes located in gene-dense regions whose position is restricted to central, active, nuclear environments.

Conservation and divergence of regulatory strategies at Hox Loci and the origin of tetrapod digits.

The evolution of tetrapod limbs from fish fins enabled the conquest of land by vertebrates and thus represents a key step in evolution. Despite the use of comparative gene expression analyses, critical aspects of this transformation remain controversial, in particular the origin of digits. Hoxa and Hoxd genes are essential for the specification of the different limb segments and their functional abrogation leads to large truncations of the appendages. Here we show that the selective transcription of mouse Hoxa genes in proximal and distal limbs is related to a bimodal higher order chromatin structure, similar to that reported for Hoxd genes, thus revealing a generic regulatory strategy implemented by both gene clusters during limb development. We found the same bimodal chromatin architecture in fish embryos, indicating that the regulatory mechanism used to pattern tetrapod limbs may predate the divergence between fish and tetrapods. However, when assessed in mice, both fish regulatory landscapes triggered transcription in proximal rather than distal limb territories, supporting an evolutionary scenario whereby digits arose as tetrapod novelties through genetic retrofitting of preexisting regulatory landscapes. We discuss the possibility to consider regulatory circuitries, rather than expression patterns, as essential parameters to define evolutionary synapomorphies.

Spatial separation of Xist RNA and polycomb proteins revealed by superresolution microscopy.

In female mammals, one of the two X chromosomes is transcriptionally silenced to equalize X-linked gene dosage relative to XY males, a process termed X chromosome inactivation. Mechanistically, this is thought to occur via directed recruitment of chromatin modifying factors by the master regulator, X-inactive specific transcript (Xist) RNA, which localizes in cis along the entire length of the chromosome. A well-studied example is the recruitment of polycomb repressive complex 2 (PRC2), for which there is evidence of a direct interaction involving the PRC2 proteins Enhancer of zeste 2 (Ezh2) and Supressor of zeste 12 (Suz12) and the A-repeat region located at the 5' end of Xist RNA. In this study, we have analyzed Xist-mediated recruitment of PRC2 using two approaches, microarray-based epigenomic mapping and superresolution 3D structured illumination microscopy. Making use of an ES cell line carrying an inducible Xist transgene located on mouse chromosome 17, we show that 24 h after synchronous induction of Xist expression, acquired PRC2 binding sites map predominantly to gene-rich regions, notably within gene bodies. Paradoxically, these new sites of PRC2 deposition do not correlate with Xist-mediated gene silencing. The 3D structured illumination microscopy was performed to assess the relative localization of PRC2 proteins and Xist RNA. Unexpectedly, we observed significant spatial separation and absence of colocalization both in the inducible Xist transgene ES cell line and in normal XX somatic cells. Our observations argue against direct interaction between Xist RNA and PRC2 proteins and, as such, prompt a reappraisal of the mechanism for PRC2 recruitment in X chromosome inactivation.

ROS-induced ribosome impairment underlies ZAKα-mediated metabolic decline in obesity and aging.

The ribotoxic stress response (RSR) is a signaling pathway in which the p38- and c-Jun N-terminal kinase (JNK)-activating mitogen-activated protein kinase kinase kinase (MAP3K) ZAKα senses stalling and/or collision of ribosomes. Here, we show that reactive oxygen species (ROS)-generating agents trigger ribosomal impairment and ZAKα activation. Conversely, zebrafish larvae deficient for ZAKα are protected from ROS-induced pathology. Livers of mice fed a ROS-generating diet exhibit ZAKα-activating changes in ribosomal elongation dynamics. Highlighting a role for the RSR in metabolic regulation, ZAK-knockout mice are protected from developing high-fat high-sugar (HFHS) diet-induced blood glucose intolerance and liver steatosis. Finally, ZAK ablation slows animals from developing the hallmarks of metabolic aging. Our work highlights ROS-induced ribosomal impairment as a physiological activation signal for ZAKα that underlies metabolic adaptation in obesity and aging.

The Transcriptional Regulator Prdm1 Is Essential for the Early Development of the Sensory Whisker Follicle and Is Linked to the Beta-Catenin First Dermal Signal.

Prdm1 mutant mice are one of the rare mutant strains that do not develop whisker hair follicles while still displaying a pelage. Here, we show that Prdm1 is expressed at the earliest stage of whisker development in clusters of mesenchymal cells before placode formation. Its conditional knockout in the murine soma leads to the loss of expression of Bmp2, Shh, Bmp4, Krt17, Edar, and Gli1, though leaving the ß-catenin-driven first dermal signal intact. Furthermore, we show that Prdm1 expressing cells not only act as a signaling center but also as a multipotent progenitor population contributing to the several lineages of the adult whisker. We confirm by genetic ablation experiments that the absence of macro vibrissae reverberates on the organization of nerve wiring in the mystacial pads and leads to the reorganization of the barrel cortex. We demonstrate that Lef1 acts upstream of Prdm1 and identify a primate-specific deletion of a Lef1 enhancer named Leaf. This loss may have been significant in the evolutionary process, leading to the progressive defunctionalization and disappearance of vibrissae in primates.

T cell receptor signaling controls Foxp3 expression via PI3K, Akt, and mTOR.

Regulatory T (Treg) cells safeguard against autoimmunity and immune pathology. Because determinants of the Treg cell fate are not completely understood, we have delineated signaling events that control the de novo expression of Foxp3 in naive peripheral CD4 T cells and in thymocytes. We report that premature termination of TCR signaling and inibition of phosphatidyl inositol 3-kinase (PI3K) p110alpha, p110delta, protein kinase B (Akt), or mammalian target of rapamycin (mTOR) conferred Foxp3 expression and Treg-like gene expression profiles. Conversely, continued TCR signaling and constitutive PI3K/Akt/mTOR activity antagonised Foxp3 induction. At the chromatin level, di- and trimethylation of lysine 4 of histone H3 (H3K4me2 and -3) near the Foxp3 transcription start site (TSS) and within the 5' untranslated region (UTR) preceded active Foxp3 expression and, like Foxp3 inducibility, was lost upon continued TCR stimulation. These data demonstrate that the PI3K/Akt/mTOR signaling network regulates Foxp3 expression.

Three-dimensional chromatin interactions remain stable upon CAG/CTG repeat expansion.

Expanded CAG/CTG repeats underlie 13 neurological disorders, including myotonic dystrophy type 1 (DM1) and Huntington's disease (HD). Upon expansion, disease loci acquire heterochromatic characteristics, which may provoke changes to chromatin conformation and thereby affect both gene expression and repeat instability. Here, we tested this hypothesis by performing 4C sequencing at the DMPK and HTT loci from DM1 and HD-derived cells. We find that allele sizes ranging from 15 to 1700 repeats displayed similar chromatin interaction profiles. This was true for both loci and for alleles with different DNA methylation levels and CTCF binding. Moreover, the ectopic insertion of an expanded CAG repeat tract did not change the conformation of the surrounding chromatin. We conclude that CAG/CTG repeat expansions are not enough to alter chromatin conformation in cis. Therefore, it is unlikely that changes in chromatin interactions drive repeat instability or changes in gene expression in these disorders.

Multicenter analysis of sputum microbiota in tuberculosis patients.

The impact of tuberculosis and of anti-tuberculosis therapy on composition and modification of human lung microbiota has been the object of several investigations. However, no clear outcome has been presented so far and the relationship between M. tuberculosis pulmonary infection and the resident lung microbiota remains vague. In this work we describe the results obtained from a multicenter study of the microbiota of sputum samples from patients with tuberculosis or unrelated lung diseases and healthy donors recruited in Switzerland, Italy and Bangladesh, with the ultimate goal of discovering a microbiota-based biomarker associated with tuberculosis. Bacterial 16S rDNA amplification, high-throughput sequencing and extensive bioinformatic analyses revealed patient-specific flora and high variability in taxon abundance. No common signature could be identified among the individuals enrolled except for minor differences which were not consistent among the different geographical settings. Moreover, anti-tuberculosis therapy did not cause any important variation in microbiota diversity, thus precluding its exploitation as a biomarker for the follow up of tuberculosis patients undergoing treatment.

Integrated proteogenomic deep sequencing and analytics accurately identify non-canonical peptides in tumor immunopeptidomes.

Efforts to precisely identify tumor human leukocyte antigen (HLA) bound peptides capable of mediating T cell-based tumor rejection still face important challenges. Recent studies suggest that non-canonical tumor-specific HLA peptides derived from annotated non-coding regions could elicit anti-tumor immune responses. However, sensitive and accurate mass spectrometry (MS)-based proteogenomics approaches are required to robustly identify these non-canonical peptides. We present an MS-based analytical approach that characterizes the non-canonical tumor HLA peptide repertoire, by incorporating whole exome sequencing, bulk and single-cell transcriptomics, ribosome profiling, and two MS/MS search tools in combination. This approach results in the accurate identification of hundreds of shared and tumor-specific non-canonical HLA peptides, including an immunogenic peptide derived from an open reading frame downstream of the melanoma stem cell marker gene ABCB5. These findings hold great promise for the discovery of previously unknown tumor antigens for cancer immunotherapy.

Large scale genomic reorganization of topological domains at the HoxD locus.

BACKGROUND: The transcriptional activation of HoxD genes during mammalian limb development involves dynamic interactions with two topologically associating domains (TADs) flanking the HoxD cluster. In particular, the activation of the most posterior HoxD genes in developing digits is controlled by regulatory elements located in the centromeric TAD (C-DOM) through long-range contacts. RESULTS: To assess the structure-function relationships underlying such interactions, we measured compaction levels and TAD discreteness using a combination of chromosome conformation capture (4C-seq) and DNA FISH. We assessed the robustness of the TAD architecture by using a series of genomic deletions and inversions that impact the integrity of this chromatin domain and that remodel long-range contacts. We report multi-partite associations between HoxD genes and up to three enhancers. We find that the loss of native chromatin topology leads to the remodeling of TAD structure following distinct parameters. CONCLUSIONS: Our results reveal that the recomposition of TAD architectures after large genomic re-arrangements is dependent on a boundary-selection mechanism in which CTCF mediates the gating of long-range contacts in combination with genomic distance and sequence specificity. Accordingly, the building of a recomposed TAD at this locus depends on distinct functional and constitutive parameters.

Transcriptome-wide co-expression analysis identifies LRRC2 as a novel mediator of mitochondrial and cardiac function.

Mitochondrial dysfunction contributes to myriad monogenic and complex pathologies. To understand the underlying mechanisms, it is essential to define the full complement of proteins that modulate mitochondrial function. To identify such proteins, we performed a meta-analysis of publicly available gene expression data. Gene co-expression analysis of a large and heterogeneous compendium of microarray data nominated a sub-population of transcripts that whilst highly correlated with known mitochondrial protein-encoding transcripts (MPETs), are not themselves recognized as generating proteins either localized to the mitochondrion or pertinent to functions therein. To focus the analysis on a medically-important condition with a strong yet incompletely understood mitochondrial component, candidates were cross-referenced with an MPET-enriched module independently generated via genome-wide co-expression network analysis of a human heart failure gene expression dataset. The strongest uncharacterized candidate in the analysis was Leucine Rich Repeat Containing 2 (LRRC2). LRRC2 was found to be localized to the mitochondria in human cells and transcriptionally-regulated by the mitochondrial master regulator Pgc-1α. We report that Lrrc2 transcript abundance correlates with that of ß-MHC, a canonical marker of cardiac hypertrophy in humans and experimentally demonstrated an elevation in Lrrc2 transcript in in vitro and in vivo rodent models of cardiac hypertrophy as well as in patients with dilated cardiomyopathy. RNAi-mediated Lrrc2 knockdown in a rat-derived cardiomyocyte cell line resulted in enhanced expression of canonical hypertrophic biomarkers as well as increased mitochondrial mass in the context of increased Pgc-1α expression. In conclusion, our meta-analysis represents a simple yet powerful springboard for the nomination of putative mitochondrially-pertinent proteins relevant to cardiac function and enabled the identification of LRRC2 as a novel mitochondrially-relevant protein and regulator of the hypertrophic response.

Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics.

BACKGROUND: Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors. METHODS: We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes. RESULTS: Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1 -/- mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1 -/- mouse model. CONCLUSIONS: These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

Temporal dynamics and developmental memory of 3D chromatin architecture at Hox gene loci.

Hox genes are essential regulators of embryonic development. Their step-wise transcriptional activation follows their genomic topology and the various states of activation are subsequently memorized into domains of progressively overlapping gene products. We have analyzed the 3D chromatin organization of Hox clusters during their early activation in vivo, using high-resolution circular chromosome conformation capture. Initially, Hox clusters are organized as single chromatin compartments containing all genes and bivalent chromatin marks. Transcriptional activation is associated with a dynamic bi-modal 3D organization, whereby the genes switch autonomously from an inactive to an active compartment. These local 3D dynamics occur within a framework of constitutive interactions within the surrounding Topological Associated Domains, indicating that this regulation process is mostly cluster intrinsic. The step-wise progression in time is fixed at various body levels and thus can account for the chromatin architectures previously described at a later stage for different anterior to posterior levels.DOI: http://dx.doi.org/10.7554/eLife.02557.001.

HTSstation: a web application and open-access libraries for high-throughput sequencing data analysis.

The HTSstation analysis portal is a suite of simple web forms coupled to modular analysis pipelines for various applications of High-Throughput Sequencing including ChIP-seq, RNA-seq, 4C-seq and re-sequencing. HTSstation offers biologists the possibility to rapidly investigate their HTS data using an intuitive web application with heuristically pre-defined parameters. A number of open-source software components have been implemented and can be used to build, configure and run HTS analysis pipelines reactively. Besides, our programming framework empowers developers with the possibility to design their own workflows and integrate additional third-party software. The HTSstation web application is accessible at http://htsstation.epfl.ch.