H3K27me3 spreading organizes canonical PRC1 chromatin architecture to regulate developmental programs.

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization. Highlights: The confinement of H3K27me3 at PRC2 nucleation sites without its spreading correlates with increased 3D chromatin interactions.The H3K27M oncohistone concentrates canonical PRC1 that anchors chromatin loop interactions in gliomas, silencing developmental programs.Stem and progenitor cells require factors promoting H3K27me3 confinement, including H3K36me2, to maintain cPRC1 loop architecture.The cPRC1-H3K27me3 interaction is a targetable driver of aberrant self-renewal in tumor cells.

An Arf/Rab cascade controls the growth and invasiveness of glioblastoma.

Glioblastoma is the most common and deadly malignant brain cancer. We now demonstrate that loss of function of the endosomal GTPase Rab35 in human brain tumor initiating cells (BTICs) increases glioblastoma growth and decreases animal survival following BTIC implantation in mouse brains. Mechanistically, we identify that the GTPase Arf5 interacts with the guanine nucleotide exchange factor (GEF) for Rab35, DENND1/connecdenn, and allosterically enhances its GEF activity toward Rab35. Knockdown of either Rab35 or Arf5 increases cell migration, invasiveness, and self-renewal in culture and enhances the growth and invasiveness of BTIC-initiated brain tumors in mice. RNAseq of the tumors reveals up-regulation of the tumor-promoting transcription factor SPOCD1, and disruption of the Arf5/Rab35 axis in glioblastoma cells leads to strong activation of the epidermal growth factor receptor, with resulting enhancement of SPOCD1 levels. These discoveries reveal an unexpected cascade between an Arf and a Rab and indicate a role for the cascade, and thus endosomal trafficking, in brain tumors.

Histone H3.3G34-Mutant Interneuron Progenitors Co-opt PDGFRA for Gliomagenesis.

Histone H3.3 glycine 34 to arginine/valine (G34R/V) mutations drive deadly gliomas and show exquisite regional and temporal specificity, suggesting a developmental context permissive to their effects. Here we show that 50% of G34R/V tumors (n = 95) bear activating PDGFRA mutations that display strong selection pressure at recurrence. Although considered gliomas, G34R/V tumors actually arise in GSX2/DLX-expressing interneuron progenitors, where G34R/V mutations impair neuronal differentiation. The lineage of origin may facilitate PDGFRA co-option through a chromatin loop connecting PDGFRA to GSX2 regulatory elements, promoting PDGFRA overexpression and mutation. At the single-cell level, G34R/V tumors harbor dual neuronal/astroglial identity and lack oligodendroglial programs, actively repressed by GSX2/DLX-mediated cell fate specification. G34R/V may become dispensable for tumor maintenance, whereas mutant-PDGFRA is potently oncogenic. Collectively, our results open novel research avenues in deadly tumors. G34R/V gliomas are neuronal malignancies where interneuron progenitors are stalled in differentiation by G34R/V mutations and malignant gliogenesis is promoted by co-option of a potentially targetable pathway, PDGFRA signaling.

Integrative analysis of vascular endothelial cell genomic features identifies AIDA as a coronary artery disease candidate gene.

BACKGROUND: Genome-wide association studies (GWAS) have identified hundreds of loci associated with coronary artery disease (CAD) and blood pressure (BP) or hypertension. Many of these loci are not linked to traditional risk factors, nor do they include obvious candidate genes, complicating their functional characterization. We hypothesize that many GWAS loci associated with vascular diseases modulate endothelial functions. Endothelial cells play critical roles in regulating vascular homeostasis, such as roles in forming a selective barrier, inflammation, hemostasis, and vascular tone, and endothelial dysfunction is a hallmark of atherosclerosis and hypertension. To test this hypothesis, we generate an integrated map of gene expression, open chromatin region, and 3D interactions in resting and TNFα-treated human endothelial cells. RESULTS: We show that genetic variants associated with CAD and BP are enriched in open chromatin regions identified in endothelial cells. We identify physical loops by Hi-C and link open chromatin peaks that include CAD or BP SNPs with the promoters of genes expressed in endothelial cells. This analysis highlights 991 combinations of open chromatin regions and gene promoters that map to 38 CAD and 92 BP GWAS loci. We validate one CAD locus, by engineering a deletion of the TNFα-sensitive regulatory element using CRISPR/Cas9 and measure the effect on the expression of the novel CAD candidate gene AIDA. CONCLUSIONS: Our data support an important role played by genetic variants acting in the vascular endothelium to modulate inter-individual risk in CAD and hypertension.

GenPipes: an open-source framework for distributed and scalable genomic analyses.

BACKGROUND: With the decreasing cost of sequencing and the rapid developments in genomics technologies and protocols, the need for validated bioinformatics software that enables efficient large-scale data processing is growing. FINDINGS: Here we present GenPipes, a flexible Python-based framework that facilitates the development and deployment of multi-step workflows optimized for high-performance computing clusters and the cloud. GenPipes already implements 12 validated and scalable pipelines for various genomics applications, including RNA sequencing, chromatin immunoprecipitation sequencing, DNA sequencing, methylation sequencing, Hi-C, capture Hi-C, metagenomics, and Pacific Biosciences long-read assembly. The software is available under a GPLv3 open source license and is continuously updated to follow recent advances in genomics and bioinformatics. The framework has already been configured on several servers, and a Docker image is also available to facilitate additional installations. CONCLUSIONS: GenPipes offers genomics researchers a simple method to analyze different types of data, customizable to their needs and resources, as well as the flexibility to create their own workflows.

Author Correction: Germline HAVCR2 mutations altering TIM-3 characterize subcutaneous panniculitis-like T cell lymphomas with hemophagocytic lymphohistiocytic syndrome.

In the version of this article originally published, the main-text sentence "In three patients of European ancestry, we identified the germline variant encoding p.Ile97Met in TIM-3, which was homozygous in two (P12 and P13) and heterozygous in one (P15) in the germline but with no TIM-3 plasma membrane expression in the tumor" misstated the identifiers of the two homozygous individuals, which should have been P13 and P14. The error has been corrected in the HTML, PDF and print versions of the paper.

Germline HAVCR2 mutations altering TIM-3 characterize subcutaneous panniculitis-like T cell lymphomas with hemophagocytic lymphohistiocytic syndrome.

Subcutaneous panniculitis-like T cell lymphoma (SPTCL), a non-Hodgkin lymphoma, can be associated with hemophagocytic lymphohistiocytosis (HLH), a life-threatening immune activation that adversely affects survival1,2. T cell immunoglobulin mucin 3 (TIM-3) is a modulator of immune responses expressed on subgroups of T and innate immune cells. We identify in ~60% of SPTCL cases germline, loss-of-function, missense variants altering highly conserved residues of TIM-3, c.245A>G (p.Tyr82Cys) and c.291A>G (p.Ile97Met), each with specific geographic distribution. The variant encoding p.Tyr82Cys TIM-3 occurs on a potential founder chromosome in patients with East Asian and Polynesian ancestry, while p.Ile97Met TIM-3 occurs in patients with European ancestry. Both variants induce protein misfolding and abrogate TIM-3's plasma membrane expression, leading to persistent immune activation and increased production of inflammatory cytokines, including tumor necrosis factor-α and interleukin-1ß, promoting HLH and SPTCL. Our findings highlight HLH-SPTCL as a new genetic entity and identify mutations causing TIM-3 alterations as a causative genetic defect in SPTCL. While HLH-SPTCL patients with mutant TIM-3 benefit from immunomodulation, therapeutic repression of the TIM-3 checkpoint may have adverse consequences.

Characterization of a FOXG1:TLE1 transcriptional network in glioblastoma-initiating cells.

Glioblastoma (GBM) is the most common and deadly malignant brain cancer of glial cell origin, with a median patient survival of less than 20 months. Transcription factors FOXG1 and TLE1 promote GBM propagation by supporting maintenance of brain tumour-initiating cells (BTICs) with stem-like properties. Here, we characterize FOXG1 and TLE1 target genes in GBM patient-derived BTICs using ChIP-Seq and RNA-Seq approaches. These studies identify 150 direct FOXG1 targets, several of which are also TLE1 targets, involved in cell proliferation, differentiation, survival, chemotaxis and angiogenesis. Negative regulators of NOTCH signalling, including CHAC1, are among the transcriptional repression targets of FOXG1:TLE1 complexes, suggesting a crosstalk between FOXG1:TLE1 and NOTCH-mediated pathways in GBM. These results provide previously unavailable insight into the transcriptional programs underlying the tumour-promoting functions of FOXG1:TLE1 in GBM.

A critical assessment of topologically associating domain prediction tools.

Topologically associating domains (TADs) have been proposed to be the basic unit of chromosome folding and have been shown to play key roles in genome organization and gene regulation. Several different tools are available for TAD prediction, but their properties have never been thoroughly assessed. In this manuscript, we compare the output of seven different TAD prediction tools on two published Hi-C data sets. TAD predictions varied greatly between tools in number, size distribution and other biological properties. Assessed against a manual annotation of TADs, individual TAD boundary predictions were found to be quite reliable, but their assembly into complete TAD structures was much less so. In addition, many tools were sensitive to sequencing depth and resolution of the interaction frequency matrix. This manuscript provides users and designers of TAD prediction tools with information that will help guide the choice of tools and the interpretation of their predictions.

Transcriptional regulation of mouse hypoglossal motor neuron somatotopic map formation.

Somatic motor neurons in the hypoglossal nucleus innervate tongue muscles controlling vital functions such as chewing, swallowing and respiration. Formation of functional hypoglossal nerve circuits depends on the establishment of precise hypoglossal motor neuron maps correlating with specific tongue muscle innervations. Little is known about the molecular mechanisms controlling mammalian hypoglossal motor neuron topographic map formation. Here we show that combinatorial expression of transcription factors Runx1, SCIP and FoxP1 defines separate mouse hypoglossal motor neuron groups with different topological, neurotransmitter and calcium-buffering phenotypes. Runx1 and SCIP are coexpressed in ventromedial hypoglossal motor neurons involved in control of tongue protrusion whereas FoxP1 is expressed in dorsomedial motor neurons associated with tongue retraction. Establishment of separate hypoglossal motor neuron maps depends in part on Runx1-mediated suppression of ventrolateral and dorsomedial motor neuron phenotypes and regulation of FoxP1 expression pattern. These findings suggest that combinatorial actions of Runx1, SCIP and FoxP1 are important for mouse hypoglossal nucleus somatotopic map formation.

Cyclin D2 is a GATA4 cofactor in cardiogenesis.

The G1 cyclins play a pivotal role in regulation of cell differentiation and proliferation. The mechanisms underlying their cell-specific roles are incompletely understood. Here, we show that a G1 cyclin, cyclin D2 (CycD2), enhances the activity of transcription factor GATA4, a key regulator of cardiomyocyte growth and differentiation. GATA4 recruits CycD2 to its target promoters, and their interaction results in synergistic activation of GATA-dependent transcription. This effect is specific to CycD2 because CycD1 is unable to potentiate activity of GATA4 and is CDK-independent. GATA4 physically interacts with CycD2 through a discreet N-terminal activation domain that is essential for the cardiogenic activity of GATA4. Human mutations in this domain that are linked to congenital heart disease interfere with CycD2-GATA4 synergy. Cardiogenesis assays in Xenopus embryos indicate that CycD2 enhances the cardiogenic function of GATA4. Together, our data uncover a role for CycD2 as a cardiogenic coactivator of GATA4 and suggest a paradigm for cell-specific effects of cyclin Ds.

Transcription factors FOXG1 and Groucho/TLE promote glioblastoma growth.

Glioblastoma (GBM) is the most common and deadly malignant brain cancer, with a median survival of <2 years. GBM displays a cellular complexity that includes brain tumour-initiating cells (BTICs), which are considered as potential key targets for GBM therapies. Here we show that the transcription factors FOXG1 and Groucho/TLE are expressed in poorly differentiated astroglial cells in human GBM specimens and in primary cultures of GBM-derived BTICs, where they form a complex. FOXG1 knockdown in BTICs causes downregulation of neural stem/progenitor and proliferation markers, increased replicative senescence, upregulation of astroglial differentiation genes and decreased BTIC-initiated tumour growth after intracranial transplantation into host mice. These effects are phenocopied by Groucho/TLE knockdown or dominant inhibition of the FOXG1:Groucho/TLE complex. These results provide evidence that transcriptional programmes regulated by FOXG1 and Groucho/TLE are important for BTIC-initiated brain tumour growth, implicating FOXG1 and Groucho/TLE in GBM tumourigenesis.

Multivariate phenotypic divergence due to the fixation of beneficial mutations in experimentally evolved lineages of a filamentous fungus.

The potential for evolutionary change is limited by the availability of genetic variation. Mutations are the ultimate source of new alleles, yet there have been few experimental investigations of the role of novel mutations in multivariate phenotypic evolution. Here, we evaluated the degree of multivariate phenotypic divergence observed in a long-term evolution experiment whereby replicate lineages of the filamentous fungus Aspergillus nidulans were derived from a single genotype and allowed to fix novel (beneficial) mutations while maintained at two different population sizes. We asked three fundamental questions regarding phenotypic divergence following approximately 800 generations of adaptation: (1) whether divergence was limited by mutational supply, (2) whether divergence proceeded in relatively many (few) multivariate directions, and (3) to what degree phenotypic divergence scaled with changes in fitness (i.e. adaptation). We found no evidence that mutational supply limited phenotypic divergence. Divergence also occurred in all possible phenotypic directions, implying that pleiotropy was either weak or sufficiently variable among new mutations so as not to constrain the direction of multivariate evolution. The degree of total phenotypic divergence from the common ancestor was positively correlated with the extent of adaptation. These results are discussed in the context of the evolution of complex phenotypes through the input of adaptive mutations.

Antagonistic interactions peak at intermediate genetic distance in clinical and laboratory strains of Pseudomonas aeruginosa.

BACKGROUND: Bacteria excrete costly toxins to defend their ecological niche. The evolution of such antagonistic interactions between individuals is expected to depend on both the social environment and the strength of resource competition. Antagonism is expected to be weak among highly similar genotypes because most individuals are immune to antagonistic agents and among dissimilar genotypes because these are unlikely to be competing for the same resources and antagonism should not yield much benefit. The strength of antagonism is therefore expected to peak at intermediate genetic distance. RESULTS: We studied the ability of laboratory strains of Pseudomonas aeruginosa to prevent growth of 55 different clinical P. aeruginosa isolates derived from cystic fibrosis patients. Genetic distance was determined using genetic fingerprints. We found that the strength of antagonism was maximal among genotypes of intermediate genetic distance and we show that genetic distance and resource use are linked. CONCLUSIONS: Our results suggest that the importance of social interactions like antagonism may be modulated by the strength of resource competition.

Fitness-associated sexual reproduction in a filamentous fungus.

Sex is a long-standing evolutionary enigma. Although the majority of eukaryotes reproduce sexually at least sometimes [1-3], the evolution of sex from an asexual ancestor has been difficult to explain because it requires sexually reproducing lineages to overcome the manifold costs of sex, including the destruction of favorable gene combinations created by selection [4, 5]. Conditions for the evolution of sex are much broader if individuals can reproduce either sexually or asexually (i.e., facultative sex) and allocate disproportionately more resources to sex when their fitness is low (fitness-associated-sex or FAS [6-10]). Although facultatively sexual organisms have been shown to engage in more sex when stressed [11], direct evidence for FAS is lacking. We provide evidence using 53 genotypes of the filamentous fungus Aspergillus nidulans in a reciprocal transplant experiment across three environments. Different genotypes achieved highest fitness in different environments and genotypes invested relatively more in sex in environments in which their fitness was lower, showing that allocation to sexual reproduction is a function of how well-adapted a genotype is to its environment. FAS in A. nidulans is unlikely to have evolved as a strategy to resist or avoid stress because asexual spores are more dispersive and equally resistant [12, 13].