TULIPs decorate the three-dimensional genome of PFA ependymoma.

Posterior fossa group A (PFA) ependymoma is a lethal brain cancer diagnosed in infants and young children. The lack of driver events in the PFA linear genome led us to search its 3D genome for characteristic features. Here, we reconstructed 3D genomes from diverse childhood tumor types and uncovered a global topology in PFA that is highly reminiscent of stem and progenitor cells in a variety of human tissues. A remarkable feature exclusively present in PFA are type B ultra long-range interactions in PFAs (TULIPs), regions separated by great distances along the linear genome that interact with each other in the 3D nuclear space with surprising strength. TULIPs occur in all PFA samples and recur at predictable genomic coordinates, and their formation is induced by expression of EZHIP. The universality of TULIPs across PFA samples suggests a conservation of molecular principles that could be exploited therapeutically.

Deep Intronic FGF14 GAA Repeat Expansion in Late-Onset Cerebellar Ataxia.

BACKGROUND: The late-onset cerebellar ataxias (LOCAs) have largely resisted molecular diagnosis. METHODS: We sequenced the genomes of six persons with autosomal dominant LOCA who were members of three French Canadian families and identified a candidate pathogenic repeat expansion. We then tested for association between the repeat expansion and disease in two independent case-control series - one French Canadian (66 patients and 209 controls) and the other German (228 patients and 199 controls). We also genotyped the repeat in 20 Australian and 31 Indian index patients. We assayed gene and protein expression in two postmortem cerebellum specimens and two induced pluripotent stem-cell (iPSC)-derived motor-neuron cell lines. RESULTS: In the six French Canadian patients, we identified a GAA repeat expansion deep in the first intron of FGF14, which encodes fibroblast growth factor 14. Cosegregation of the repeat expansion with disease in the families supported a pathogenic threshold of at least 250 GAA repeats ([GAA]≥250). There was significant association between FGF14 (GAA)≥250 expansions and LOCA in the French Canadian series (odds ratio, 105.60; 95% confidence interval [CI], 31.09 to 334.20; P<0.001) and in the German series (odds ratio, 8.76; 95% CI, 3.45 to 20.84; P<0.001). The repeat expansion was present in 61%, 18%, 15%, and 10% of French Canadian, German, Australian, and Indian index patients, respectively. In total, we identified 128 patients with LOCA who carried an FGF14 (GAA)≥250 expansion. Postmortem cerebellum specimens and iPSC-derived motor neurons from patients showed reduced expression of FGF14 RNA and protein. CONCLUSIONS: A dominantly inherited deep intronic GAA repeat expansion in FGF14 was found to be associated with LOCA. (Funded by Fondation Groupe Monaco and others.).

MAFG-driven astrocytes promote CNS inflammation.

Multiple sclerosis is a chronic inflammatory disease of the CNS1. Astrocytes contribute to the pathogenesis of multiple sclerosis2, but little is known about the heterogeneity of astrocytes and its regulation. Here we report the analysis of astrocytes in multiple sclerosis and its preclinical model experimental autoimmune encephalomyelitis (EAE) by single-cell RNA sequencing in combination with cell-specific Ribotag RNA profiling, assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), genome-wide analysis of DNA methylation and in vivo CRISPR-Cas9-based genetic perturbations. We identified astrocytes in EAE and multiple sclerosis that were characterized by decreased expression of NRF2 and increased expression of MAFG, which cooperates with MAT2α to promote DNA methylation and represses antioxidant and anti-inflammatory transcriptional programs. Granulocyte-macrophage colony-stimulating factor (GM-CSF) signalling in astrocytes drives the expression of MAFG and MAT2α and pro-inflammatory transcriptional modules, contributing to CNS pathology in EAE and, potentially, multiple sclerosis. Our results identify candidate therapeutic targets in multiple sclerosis.

ATR-X syndrome protein targets tandem repeats and influences allele-specific expression in a size-dependent manner.

ATRX is an X-linked gene of the SWI/SNF family, mutations in which cause syndromal mental retardation and downregulation of α-globin expression. Here we show that ATRX binds to tandem repeat (TR) sequences in both telomeres and euchromatin. Genes associated with these TRs can be dysregulated when ATRX is mutated, and the change in expression is determined by the size of the TR, producing skewed allelic expression. This reveals the characteristics of the affected genes, explains the variable phenotypes seen with identical ATRX mutations, and illustrates a new mechanism underlying variable penetrance. Many of the TRs are G rich and predicted to form non-B DNA structures (including G-quadruplex) in vivo. We show that ATRX binds G-quadruplex structures in vitro, suggesting a mechanism by which ATRX may play a role in various nuclear processes and how this is perturbed when ATRX is mutated.

Apposition of Fibroblasts With Metaplastic Gastric Cells Promotes Dysplastic Transition.

BACKGROUND & AIMS: Elements of field cancerization, including atrophic gastritis, metaplasia, and dysplasia, promote gastric cancer development in association with chronic inflammation. However, it remains unclear how stroma changes during carcinogenesis and how the stroma contributes to progression of gastric preneoplasia. Here we investigated heterogeneity of fibroblasts, one of the most important elements in the stroma, and their roles in neoplastic transformation of metaplasia. METHODS: We used single-cell transcriptomics to evaluate the cellular heterogeneity of mucosal cells from patients with gastric cancer. Tissue sections from the same cohort and tissue microarrays were used to identify the geographical distribution of distinct fibroblast subsets. We further evaluated the role of fibroblasts from pathologic mucosa in dysplastic progression of metaplastic cells using patient-derived metaplastic gastroids and fibroblasts. RESULTS: We identified 4 subsets of fibroblasts within stromal cells defined by the differential expression of PDGFRA, FBLN2, ACTA2, or PDGFRB. Each subset was distributed distinctively throughout stomach tissues with different proportions at each pathologic stage. The PDGFRα+ subset expanded in metaplasia and cancer compared with normal, maintaining a close proximity with the epithelial compartment. Co-culture of metaplasia- or cancer-derived fibroblasts with gastroids showing the characteristics of spasmolytic polypeptide-expressing metaplasia-induced disordered growth, loss of metaplastic markers, and increases in markers of dysplasia. Culture of metaplastic gastroids with conditioned media from metaplasia- or cancer-derived fibroblasts also promoted dysplastic transition. CONCLUSIONS: These findings indicate that fibroblast associations with metaplastic epithelial cells can facilitate direct transition of metaplastic spasmolytic polypeptide-expressing metaplasia cell lineages into dysplastic lineages.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Principles of dimer-specific gene regulation revealed by a comprehensive characterization of NF-κB family DNA binding.

The unique DNA-binding properties of distinct NF-κB dimers influence the selective regulation of NF-κB target genes. To more thoroughly investigate these dimer-specific differences, we combined protein-binding microarrays and surface plasmon resonance to evaluate DNA sites recognized by eight different NF-κB dimers. We observed three distinct binding-specificity classes and clarified mechanisms by which dimers might regulate distinct sets of genes. We identified many new nontraditional NF-κB binding site (κB site) sequences and highlight the plasticity of NF-κB dimers in recognizing κB sites with a single consensus half-site. This study provides a database that can be used in efforts to identify NF-κB target sites and uncover gene regulatory circuitry.

Lessons learned from understanding chemotherapy resistance in epithelial tubo-ovarian carcinoma from BRCA1and BRCA2mutation carriers.

BRCA1 and BRCA2 are multi-functional proteins and key factors for maintaining genomic stability through their roles in DNA double strand break repair by homologous recombination, rescuing stalled or damaged DNA replication forks, and regulation of cell cycle DNA damage checkpoints. Impairment of any of these critical roles results in genomic instability, a phenotypic hallmark of many cancers including breast and epithelial ovarian carcinomas (EOC). Damaging, usually loss of function germline and somatic variants in BRCA1 and BRCA2, are important drivers of the development, progression, and management of high-grade serous tubo-ovarian carcinoma (HGSOC). However, mutations in these genes render patients particularly sensitive to platinum-based chemotherapy, and to the more innovative targeted therapies with poly-(ADP-ribose) polymerase inhibitors (PARPis) that are targeted to BRCA1/BRCA2 mutation carriers. Here, we reviewed the literature on the responsiveness of BRCA1/2-associated HGSOC to platinum-based chemotherapy and PARPis, and propose mechanisms underlying the frequent development of resistance to these therapeutic agents.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly.

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Detection and genomic analysis of BRAF fusions in Juvenile Pilocytic Astrocytoma through the combination and integration of multi-omic data.

BACKGROUND: Juvenile Pilocytic Astrocytomas (JPAs) are one of the most common pediatric brain tumors, and they are driven by aberrant activation of the mitogen-activated protein kinase (MAPK) signaling pathway. RAF-fusions are the most common genetic alterations identified in JPAs, with the prototypical KIAA1549-BRAF fusion leading to loss of BRAF's auto-inhibitory domain and subsequent constitutive kinase activation. JPAs are highly vascular and show pervasive immune infiltration, which can lead to low tumor cell purity in clinical samples. This can result in gene fusions that are difficult to detect with conventional omics approaches including RNA-Seq. METHODS: To this effect, we applied RNA-Seq as well as linked-read whole-genome sequencing and in situ Hi-C as new approaches to detect and characterize low-frequency gene fusions at the genomic, transcriptomic and spatial level. RESULTS: Integration of these datasets allowed the identification and detailed characterization of two novel BRAF fusion partners, PTPRZ1 and TOP2B, in addition to the canonical fusion with partner KIAA1549. Additionally, our Hi-C datasets enabled investigations of 3D genome architecture in JPAs which showed a high level of correlation in 3D compartment annotations between JPAs compared to other pediatric tumors, and high similarity to normal adult astrocytes. We detected interactions between BRAF and its fusion partners exclusively in tumor samples containing BRAF fusions. CONCLUSIONS: We demonstrate the power of integrating multi-omic datasets to identify low frequency fusions and characterize the JPA genome at high resolution. We suggest that linked-reads and Hi-C could be used in clinic for the detection and characterization of JPAs.

Single-Cell RNA-Seq Analysis of Cells from Degenerating and Non-Degenerating Intervertebral Discs from the Same Individual Reveals New Biomarkers for Intervertebral Disc Degeneration.

In this study, we used single-cell transcriptomic analysis to identify new specific biomarkers for nucleus pulposus (NP) and inner annulus fibrosis (iAF) cells, and to define cell populations within non-degenerating (nD) and degenerating (D) human intervertebral discs (IVD) of the same individual. Cluster analysis based on differential gene expression delineated 14 cell clusters. Gene expression profiles at single-cell resolution revealed the potential functional differences linked to degeneration, and among NP and iAF subpopulations. GO and KEGG analyses discovered molecular functions, biological processes, and transcription factors linked to cell type and degeneration state. We propose two lists of biomarkers, one as specific cell type, including C2orf40, MGP, MSMP, CD44, EIF1, LGALS1, RGCC, EPYC, HILPDA, ACAN, MT1F, CHI3L1, ID1, ID3 and TMED2. The second list proposes predictive IVD degeneration genes, including MT1G, SPP1, HMGA1, FN1, FBXO2, SPARC, VIM, CTGF, MGST1, TAF1D, CAPS, SPTSSB, S100A1, CHI3L2, PLA2G2A, TNRSF11B, FGFBP2, MGP, SLPI, DCN, MT-ND2, MTCYB, ADIRF, FRZB, CLEC3A, UPP1, S100A2, PRG4, COL2A1, SOD2 and MT2A. Protein and mRNA expression of MGST1, vimentin, SOD2 and SYF2 (p29) genes validated our scRNA-seq findings. Our data provide new insights into disc cells phenotypes and biomarkers of IVD degeneration that could improve diagnostic and therapeutic options.

Whole-exome sequencing of non-BRCA1/BRCA2 mutation carrier cases at high-risk for hereditary breast/ovarian cancer.

The current study aimed to identify new breast and/or ovarian cancer predisposition genes. For that, whole-exome sequencing (WES) was performed in the germline DNA of 52 non-BRCA1/BRCA2/TP53 mutation carrier women at high-risk for hereditary breast and ovarian cancer (HBOC). All variants were classified using information from population and disease specific databases, in silico prediction tools and the American College of Medical Genetics and Genomics (ACMG) criteria. Loss of heterozygosity (LOH) of tumor samples and segregation analyses were performed whenever possible. The variants identified were investigated in a second, independent cohort of 17 BC cases. Pathogenic/Likely Pathogenic variants were identified in known cancer genes such as CHEK2, MUTYH, PMS2, and RAD51C. Rare and potentially pathogenic variants were identified in DNA repair genes (FAN1, POLQ, and RAD54L) and other cancer-related genes such as DROSHA and SLC34A2. Interestingly, the variant c.149T>G in the FAN1 gene was identified in two unrelated families, and exhibited LOH in the tumor tissue of one of them. In conclusion, this is the largest Brazilian WES study involving families at high-risk for HBOC which has brought novel insights into the role of potentially new genetic risk factors for hereditary breast and ovarian cancer.

Molecular characterization of DICER1-mutated pituitary blastoma.

Pituitary blastoma (PitB) has recently been identified as a rare and potentially lethal pediatric intracranial tumor. All cases that have been studied molecularly possess at least one DICER1 pathogenic variant. Here, we characterized nine pituitary samples, including three fresh frozen PitBs, three normal fetal pituitary glands and three normal postnatal pituitary glands using small-RNA-Seq, RNA-Seq, methylation profiling, whole genome sequencing and Nanostring® miRNA analyses; an extended series of 21 pituitary samples was used for validation purposes. These analyses demonstrated that DICER1 RNase IIIb hotspot mutations in PitBs induced improper processing of miRNA precursors, resulting in aberrant 5p-derived miRNA products and a skewed distribution of miRNAs favoring mature 3p over 5p miRNAs. This led to dysregulation of hundreds of 5p and 3p miRNAs and concomitant dysregulation of numerous mRNA targets. Gene expression analysis revealed PRAME as the most significantly upregulated gene (500-fold increase). PRAME is a member of the Retinoic Acid Receptor (RAR) signaling pathway and in PitBs, the RAR, WNT and NOTCH pathways are dysregulated. Cancer Hallmarks analysis showed that PI3K pathway is activated in the tumors. Whole genome sequencing demonstrated a quiet genome with very few somatic alterations. The comparison of methylation profiles to publicly available data from ~ 3000 other central nervous system tumors revealed that PitBs have a distinct methylation profile compared to all other tumors, including pituitary adenomas. In conclusion, this comprehensive characterization of DICER1-related PitB revealed key molecular underpinnings of PitB and identified pathways that could potentially be exploited in the treatment of this tumor.

Developmental trajectory of oligodendrocyte progenitor cells in the human brain revealed by single cell RNA sequencing.

Characterizing the developmental trajectory of oligodendrocyte progenitor cells (OPC) is of great interest given the importance of these cells in the remyelination process. However, studies of human OPC development remain limited by the availability of whole cell samples and material that encompasses a wide age range, including time of peak myelination. In this study, we apply single cell RNA sequencing to viable whole cells across the age span and link transcriptomic signatures of oligodendrocyte-lineage cells with stage-specific functional properties. Cells were isolated from surgical tissue samples of second-trimester fetal, 2-year-old pediatric, 13-year-old adolescent, and adult donors by mechanical and enzymatic digestion, followed by percoll gradient centrifugation. Gene expression was analyzed using droplet-based RNA sequencing (10X Chromium). Louvain clustering analysis identified three distinct cellular subpopulations based on 5,613 genes, comprised of an early OPC (e-OPC) group, a late OPC group (l-OPC), and a mature OL (MOL) group. Gene ontology terms enriched for e-OPCs included cell cycle and development, for l-OPCs included extracellular matrix and cell adhesion, and for MOLs included myelination and cytoskeleton. The e-OPCs were mostly confined to the premyelinating fetal group, and the l-OPCs were most highly represented in the pediatric age group, corresponding to the peak age of myelination. Cells expressing a signature characteristic of l-OPCs were identified in the adult brain in situ using RNAScope. These findings highlight the transcriptomic variability in OL-lineage cells before, during, and after peak myelination and contribute to identifying novel pathways required to achieve remyelination.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly.

BACKGROUND: The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use. RESULTS: Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. CONCLUSIONS: The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

RNA Sequencing Reveals that Kaposi Sarcoma-Associated Herpesvirus Infection Mimics Hypoxia Gene Expression Signature.

Kaposi sarcoma-associated herpesvirus (KSHV) causes several tumors and hyperproliferative disorders. Hypoxia and hypoxia-inducible factors (HIFs) activate latent and lytic KSHV genes, and several KSHV proteins increase the cellular levels of HIF. Here, we used RNA sequencing, qRT-PCR, Taqman assays, and pathway analysis to explore the miRNA and mRNA response of uninfected and KSHV-infected cells to hypoxia, to compare this with the genetic changes seen in chronic latent KSHV infection, and to explore the degree to which hypoxia and KSHV infection interact in modulating mRNA and miRNA expression. We found that the gene expression signatures for KSHV infection and hypoxia have a 34% overlap. Moreover, there were considerable similarities between the genes up-regulated by hypoxia in uninfected (SLK) and in KSHV-infected (SLKK) cells. hsa-miR-210, a HIF-target known to have pro-angiogenic and anti-apoptotic properties, was significantly up-regulated by both KSHV infection and hypoxia using Taqman assays. Interestingly, expression of KSHV-encoded miRNAs was not affected by hypoxia. These results demonstrate that KSHV harnesses a part of the hypoxic cellular response and that a substantial portion of hypoxia-induced changes in cellular gene expression are induced by KSHV infection. Therefore, targeting hypoxic pathways may be a useful way to develop therapeutic strategies for KSHV-related diseases.

Identification and characterization of enhancers controlling the inflammatory gene expression program in macrophages.

Enhancers determine tissue-specific gene expression programs. Enhancers are marked by high histone H3 lysine 4 mono-methylation (H3K4me1) and by the acetyl-transferase p300, which has allowed genome-wide enhancer identification. However, the regulatory principles by which subsets of enhancers become active in specific developmental and/or environmental contexts are unknown. We exploited inducible p300 binding to chromatin to identify, and then mechanistically dissect, enhancers controlling endotoxin-stimulated gene expression in macrophages. In these enhancers, binding sites for the lineage-restricted and constitutive Ets protein PU.1 coexisted with those for ubiquitous stress-inducible transcription factors such as NF-kappaB, IRF, and AP-1. PU.1 was required for maintaining H3K4me1 at macrophage-specific enhancers. Reciprocally, ectopic expression of PU.1 reactivated these enhancers in fibroblasts. Thus, the combinatorial assembly of tissue- and signal-specific transcription factors determines the activity of a distinct group of enhancers. We suggest that this may represent a general paradigm in tissue-restricted and stimulus-responsive gene regulation.

Estrogen receptor-α directly regulates the hypoxia-inducible factor 1 pathway associated with antiestrogen response in breast cancer.

A majority of breast cancers are driven by estrogen via estrogen receptor-α (ERα). Our previous studies indicate that hypoxia-inducible factor 1α (HIF-1α) cooperates with ERα in breast cancer cells. However, whether ERα is implicated in the direct regulation of HIF-1α and the role of HIF-1α in endocrine therapy response are unknown. In this study we found that a subpopulation of HIF-1α targets, many of them bearing both hypoxia response elements and estrogen response elements, are regulated by ERα in normoxia and hypoxia. Interestingly, the HIF-1α gene itself also bears an estrogen response element, and its expression is directly regulated by ERα. Clinical data revealed that expression of the HIF-1α gene or a hypoxia metagene signature is associated with a poor outcome to endocrine treatment in ERα(+) breast cancer. HIF-1α was able to confer endocrine therapy resistance to ERα(+) breast cancer cells. Our findings define, for the first time to our knowledge, a direct regulatory pathway between ERα and HIF-1α, which might modulate hormone response in treatment.