Transmembrane Protein TMEM230, Regulator of Glial Cell Vascular Mimicry and Endothelial Cell Angiogenesis in High-Grade Heterogeneous Infiltrating Gliomas and Glioblastoma.

High-grade gliomas (HGGs) and glioblastoma multiforme (GBM) are characterized by a heterogeneous and aggressive population of tissue-infiltrating cells that promote both destructive tissue remodeling and aberrant vascularization of the brain. The formation of defective and permeable blood vessels and microchannels and destructive tissue remodeling prevent efficient vascular delivery of pharmacological agents to tumor cells and are the significant reason why therapeutic chemotherapy and immunotherapy intervention are primarily ineffective. Vessel-forming endothelial cells and microchannel-forming glial cells that recapitulate vascular mimicry have both infiltration and destructive remodeling tissue capacities. The transmembrane protein TMEM230 (C20orf30) is a master regulator of infiltration, sprouting of endothelial cells, and microchannel formation of glial and phagocytic cells. A high level of TMEM230 expression was identified in patients with HGG, GBM, and U87-MG cells. In this study, we identified candidate genes and molecular pathways that support that aberrantly elevated levels of TMEM230 play an important role in regulating genes associated with the initial stages of cell infiltration and blood vessel and microchannel (also referred to as tumor microtubule) formation in the progression from low-grade to high-grade gliomas. As TMEM230 regulates infiltration, vascularization, and tissue destruction capacities of diverse cell types in the brain, TMEM230 is a promising cancer target for heterogeneous HGG tumors.

Skeletal Muscles of Sedentary and Physically Active Aged People Have Distinctive Genic Extrachromosomal Circular DNA Profiles.

To bring new extrachromosomal circular DNA (eccDNA) enrichment technologies closer to the clinic, specifically for screening, early diagnosis, and monitoring of diseases or lifestyle conditions, it is paramount to identify the differential pattern of the genic eccDNA signal between two states. Current studies using short-read sequenced purified eccDNA data are based on absolute numbers of unique eccDNAs per sample or per gene, length distributions, or standard methods for RNA-seq differential analysis. Previous analyses of RNA-seq data found significant transcriptomics difference between sedentary and active life style skeletal muscle (SkM) in young people but very few in old. The first attempt using circulomics data from SkM and blood of aged lifelong sedentary and physically active males found no difference at eccDNA level. To improve the capability of finding differences between circulomics data groups, we designed a computational method to identify Differentially Produced per Gene Circles (DPpGCs) from short-read sequenced purified eccDNA data based on the circular junction, split-read signal, of the eccDNA, and implemented it into a software tool DifCir in Matlab. We employed DifCir to find to the distinctive features of the influence of the physical activity or inactivity in the aged SkM that would have remained undetected by transcriptomics methods. We mapped the data from tissue from SkM and blood from two groups of aged lifelong sedentary and physically active males using Circle_finder and subsequent merging and filtering, to find the number and length distribution of the unique eccDNA. Next, we used DifCir to find up-DPpGCs in the SkM of the sedentary and active groups. We assessed the functional enrichment of the DPpGCs using Disease Gene Network and Gene Set Enrichment Analysis. To find genes that produce eccDNA in a group without comparison with another group, we introduced a method to find Common PpGCs (CPpGCs) and used it to find CPpGCs in the SkM of the sedentary and active group. Finally, we found the eccDNA that carries whole genes. We discovered that the eccDNA in the SkM of the sedentary group is not statistically different from that of physically active aged men in terms of number and length distribution of eccDNA. In contrast, with DifCir we found distinctive gene-associated eccDNA fingerprints. We identified statistically significant up-DPpGCs in the two groups, with the top up-DPpGCs shed by the genes AGBL4, RNF213, DNAH7, MED13, and WWTR1 in the sedentary group, and ZBTB7C, TBCD, ITPR2, and DDX11-AS1 in the active group. The up-DPpGCs in both groups carry mostly gene fragments rather than whole genes. Though the subtle transcriptomics difference, we found RYR1 to be both transcriptionally up-regulated and up-DPpGCs gene in sedentary SkM. DifCir emphasizes the high sensitivity of the circulome compared to the transcriptome to detect the molecular fingerprints of exercise in aged SkM. It allows efficient identification of gene hotspots that excise more eccDNA in a health state or disease compared to a control condition.

Consistent DNA Hypomethylations in Prostate Cancer.

With approximately 1.4 million men annually diagnosed with prostate cancer (PCa) worldwide, PCa remains a dreaded threat to life and source of devastating morbidity. In recent decades, a significant decrease in age-specific PCa mortality has been achieved by increasing prostate-specific antigen (PSA) screening and improving treatments. Nevertheless, upcoming, augmented recommendations against PSA screening underline an escalating disproportion between the benefit and harm of current diagnosis/prognosis and application of radical treatment standards. Undoubtedly, new potent diagnostic and prognostic tools are urgently needed to alleviate this tensed situation. They should allow a more reliable early assessment of the upcoming threat, in order to enable applying timely adjusted and personalized therapy and monitoring. Here, we present a basic study on an epigenetic screening approach by Methylated DNA Immunoprecipitation (MeDIP). We identified genes associated with hypomethylated CpG islands in three PCa sample cohorts. By adjusting our computational biology analyses to focus on single CpG-enriched 60-nucleotide-long DNA probes, we revealed numerous consistently differential methylated DNA segments in PCa. They were associated among other genes with NOTCH3, CDK2AP1, KLK4, and ADAM15. These can be used for early discrimination, and might contribute to a new epigenetic tumor classification system of PCa. Our analysis shows that we can dissect short, differential methylated CpG-rich DNA fragments and combinations of them that are consistently present in all tumors. We name them tumor cell-specific differential methylated CpG dinucleotide signatures (TUMS).

Epigenetics in the Diagnosis and Therapy of Malignant Melanoma.

Epigenetic mechanisms are fundamentally important for cancer initiation and development. However, a survey of the literature reveals that, to date, they appear less comprehensively investigated in melanoma than in many other cancers, e.g., prostate, breast, and colon carcinoma. The aim of this review is to provide a short summary of epigenetic aspects of functional relevance for melanoma pathogenesis. In addition, some new perspectives from epigenetic research in other cancers with potential for melanoma diagnosis and therapy are introduced. For example, the PrimeEpiHit hypothesis in urothelial carcinoma, which, similarly to malignant melanoma, can also be triggered by a single exogenous noxa, states that one of the first steps for cancer initiation could be epigenetic changes in key genes of one-carbon metabolism. The application of such insights may contribute to further progress in the diagnosis and therapy of melanoma, a deadly type of cancer.

Therapeutic HNF4A mRNA attenuates liver fibrosis in a preclinical model.

BACKGROUND & AIMS: Therapeutic targeting of injuries that require transient restoration of proteins by mRNA delivery is an attractive approach that, until recently, has remained poorly explored. In this study, we examined the therapeutic utility of mRNA delivery for liver fibrosis and cirrhosis. Specifically, we aimed to demonstrate the therapeutic efficacy of human hepatocyte nuclear factor alpha (HNF4A) mRNA in mouse models of fibrosis and cirrhosis. METHODS: We investigated restoration of hepatocyte functions by HNF4A mRNA transfection in vitro, and analyzed the attenuation of liver fibrosis and cirrhosis in multiple mouse models, by delivering hepatocyte-targeted biodegradable lipid nanoparticles (LNPs) encapsulating HNF4A mRNA. To identify potential mechanisms of action, we performed microarray-based gene expression profiling, single-cell RNA sequencing, and chromatin immunoprecipitation. We used primary liver cells and human liver buds for additional functional validation. RESULTS: Expression of HNF4A mRNA led to restoration of the metabolic activity of fibrotic primary murine and human hepatocytes in vitro. Repeated in vivo delivery of LNP-encapsulated HNF4A mRNA induced a robust inhibition of fibrogenesis in 4 independent mouse models of hepatotoxin- and cholestasis-induced liver fibrosis. Mechanistically, we discovered that paraoxonase 1 is a direct target of HNF4A and it contributes to HNF4A-mediated attenuation of liver fibrosis via modulation of liver macrophages and hepatic stellate cells. CONCLUSION: Collectively, our findings provide the first direct preclinical evidence of the applicability of HNF4A mRNA therapeutics for the treatment of fibrosis in the liver. LAY SUMMARY: Liver fibrosis and cirrhosis remain unmet medical needs and contribute to high mortality worldwide. Herein, we take advantage of a promising therapeutic approach to treat liver fibrosis and cirrhosis. We demonstrate that restoration of a key gene, HNF4A, via mRNA encapsulated in lipid nanoparticles decreased injury in multiple mouse models of fibrosis and cirrhosis. Our study provides proof-of-concept that mRNA therapy is a promising strategy for reversing liver fibrosis and cirrhosis.

GeromiRs Are Downregulated in the Tumor Microenvironment during Colon Cancer Colonization of the Liver in a Murine Metastasis Model.

Cancer is a phenomenon broadly related to ageing in various ways such as cell cycle deregulation, metabolic defects or telomerases dysfunction as principal processes. Although the tumor cell is the main actor in cancer progression, it is not the only element of the disease. Cells and the matrix surrounding the tumor, called the tumor microenvironment (TME), play key roles in cancer progression. Phenotypic changes of the TME are indispensable for disease progression and a few of these transformations are produced by epigenetic changes including miRNA dysregulation. In this study, we found that a specific group of miRNAs in the liver TME produced by colon cancer called geromiRs, which are miRNAs related to the ageing process, are significantly downregulated. The three principal cell types involved in the liver TME, namely, liver sinusoidal endothelial cells, hepatic stellate (Ito) cells and Kupffer cells, were isolated from a murine hepatic metastasis model, and the miRNA and gene expression profiles were studied. From the 115 geromiRs and their associated hallmarks of aging, which we compiled from the literature, 75 were represented in the used microarrays, 26 out of them were downregulated in the TME cells during colon cancer colonization of the liver, and none of them were upregulated. The histone modification hallmark of the downregulated geromiRs is significantly enriched with the geromiRs miR-15a, miR-16, miR-26a, miR-29a, miR-29b and miR-29c. We built a network of all of the geromiRs downregulated in the TME cells and their gene targets from the MirTarBase database, and we analyzed the expression of these geromiR gene targets in the TME. We found that Cercam and Spsb4, identified as prognostic markers in a few cancer types, are associated with downregulated geromiRs and are upregulated in the TME cells.

P3BSseq: parallel processing pipeline software for automatic analysis of bisulfite sequencing data.

Motivation: Bisulfite sequencing (BSseq) processing is among the most cumbersome next generation sequencing (NGS) applications. Though some BSseq processing tools are available, they are scattered, require puzzling parameters and are running-time and memory-usage demanding. Results: We developed P3BSseq, a parallel processing pipeline for fast, accurate and automatic analysis of BSseq reads that trims, aligns, annotates, records the intermediate results, performs bisulfite conversion quality assessment, generates BED methylome and report files following the NIH standards. P3BSseq outperforms the known BSseq mappers regarding running time, computer hardware requirements (processing power and memory use) and is optimized to process the upcoming, extended BSseq reads. We optimized the P3BSseq parameters for directional and non-directional libraries, and for single-end and paired-end reads of Whole Genome and Reduced Representation BSseq. P3BSseq is a user-friendly streamlined solution for BSseq upstream analysis, requiring only basic computer and NGS knowledge. Availability and Implementation: P3BSseq binaries and documentation are available at: http://sourceforge.net/p/p3bsseq/wiki/Home/ Contact: mararabra@yahoo.co.uk Supplimentary Information: Supplementary data are available at Bioinformatics online.

DNA methylation regulates discrimination of enhancers from promoters through a H3K4me1-H3K4me3 seesaw mechanism.

BACKGROUND: DNA methylation at promoters is largely correlated with inhibition of gene expression. However, the role of DNA methylation at enhancers is not fully understood, although a crosstalk with chromatin marks is expected. Actually, there exist contradictory reports about positive and negative correlations between DNA methylation and H3K4me1, a chromatin hallmark of enhancers. RESULTS: We investigated the relationship between DNA methylation and active chromatin marks through genome-wide correlations, and found anti-correlation between H3K4me1 and H3K4me3 enrichment at low and intermediate DNA methylation loci. We hypothesized "seesaw" dynamics between H3K4me1 and H3K4me3 in the low and intermediate DNA methylation range, in which DNA methylation discriminates between enhancers and promoters, marked by H3K4me1 and H3K4me3, respectively. Low methylated regions are H3K4me3 enriched, while those with intermediate DNA methylation levels are progressively H3K4me1 enriched. Additionally, the enrichment of H3K27ac, distinguishing active from primed enhancers, follows a plateau in the lower range of the intermediate DNA methylation level, corresponding to active enhancers, and decreases linearly in the higher range of the intermediate DNA methylation. Thus, the decrease of the DNA methylation switches smoothly the state of the enhancers from a primed to an active state. We summarize these observations into a rule of thumb of one-out-of-three methylation marks: "In each genomic region only one out of these three methylation marks {DNA methylation, H3K4me1, H3K4me3} is high. If it is the DNA methylation, the region is inactive. If it is H3K4me1, the region is an enhancer, and if it is H3K4me3, the region is a promoter". To test our model, we used available genome-wide datasets of H3K4 methyltransferases knockouts. Our analysis suggests that CXXC proteins, as readers of non-methylated CpGs would regulate the "seesaw" mechanism that focuses H3K4me3 to unmethylated sites, while being repulsed from H3K4me1 decorated enhancers and CpG island shores. CONCLUSIONS: Our results show that DNA methylation discriminates promoters from enhancers through H3K4me1-H3K4me3 seesaw mechanism, and suggest its possible function in the inheritance of chromatin marks after cell division. Our analyses suggest aberrant formation of promoter-like regions and ectopic transcription of hypomethylated regions of DNA. Such mechanism process can have important implications in biological process in where it has been reported abnormal DNA methylation status such as cancer and aging.

Does mouse embryo primordial germ cell activation start before implantation as suggested by single-cell transcriptomics dynamics?

STUDY HYPOTHESIS: Does primordial germ cell (PGC) activation start before mouse embryo implantation, and does the possible regulation of the DNA (cytosine-5-)-methyltransferase 3-like (Dnmt3l) by transcription factor AP-2, gamma (TCFAP2C) have a role in this activation and in the primitive endoderm (PE)-epiblast (EPI) lineage specification? STUDY FINDING: A burst of expression of PGC markers, such as Dppa3/Stella, Ifitm2/Fragilis, Fkbp6 and Prdm4, is observed from embryonic day (E) 3.25, and some of them, together with the late germ cell markers Zp3, Mcf2 and Morc1, become restricted to the EPI subpopulation at E4.5, while the dynamics analysis of the PE-EPI transitions in the single-cell data suggests that TCFAP2C transitorily represses Dnmt3l in EPI cells at E3.5 and such repression is withdrawn with reactivation of Dnmt3l expression in PE and EPI cells at E4.5. WHAT IS KNOWN ALREADY: In the mouse preimplantation embryo, cells with the same phenotype take different fates based on the orchestration between topological clues (cell polarity, positional history and division orientation) and gene regulatory rules (at transcriptomics and epigenomics level), prompting the proposal of positional, stochastic and combined models explaining the specification mechanism. PGC specification starts at E6.0-6.5 post-implantation. In view of the important role of DNA methylation in developmental events, the cross-talk between some transcription factors and DNA methyltransferases is of particular relevance. TCFAP2C has a CpG DNA methylation motif that is not methylated in pluripotent cells and that could potentially bind on DNMT3L, the stimulatory DNA methyltransferase co-factor that assists in the process of de novo DNA methylation. Chromatin-immunoprecipitation analysis has demonstrated that Dnmt3l is indeed a target of TCFAP2C. STUDY DESIGN, SAMPLES/MATERIALS, METHODS: We aimed to assess the timing of early preimplantation events and to understand better the segregation of the inner cell mass (ICM) into PE and EPI. We designed a single-cell transcriptomics dynamics computational study to identify markers of the PE-EPI bifurcation in ICM cells through searching for statistically significant (using the Student's t-test method) differently expressed genes (DEGs) between PE and EPI cells from E3.5 to E4.5. The DEGs common for E3.5 and E4.5 were used as the markers defining the steady states. We collected microarray and next-generation sequencing transcriptomics data from public databases from bulk populations and single cells from mice at E3.25, E3.5 and E4.5. The results are based on three independent single-cell transcriptomics data sets, with a fold change of 3 and P-value <0.01 for the DEG selection. MAIN RESULTS AND THE ROLE OF CHANCE: The dynamics analysis revealed new transitory E3.5 and steady PE and EPI markers. Among the transitory E3.5 PE markers (Dnmt3l, Dusp4, Cpne8, Akap13, Dcaf12l1, Aaed1, B4galt6, BC100530, Rnpc3, Tfpi, Lgalsl, Ckap4 and Fbxl20), several (Dusp4, Akap13, Cpn8, Dcaf12l1 and Tfpi) are related to the extracellular regulated kinase pathway. We also identified new transitory E3.5 EPI markers (Sgk1, Mal, Ubxn2a, Atg16l2, Gm13102, Tcfap2c, Hexb, Slc1a1, Svip, Liph and Mier3), six new stable PE markers (Sdc4, Cpn1, Dkk1, Havcr1, F2r/Par1 and Slc7a6os) as well as three new stable EPI markers (Zp3, Mcf2 and Hexb), which are known to be late stage germ cell markers. We found that mouse PGC marker activation starts at least at E3.25 preimplantation. The transcriptomics dynamics analyses support the regulation of Dnmt3l expression by TCFAP2C. LIMITATIONS, REASONS FOR CAUTION: Since the regulation of Dnmt3l by TCFAP2C is based on computational prediction of DNA methylation motifs, Chip-Seq and transcriptomics data, functional studies are required to validate this result. WIDER IMPLICATIONS OF THE FINDINGS: We identified a collection of previously undescribed E3.5-specific PE and EPI markers, and new steady PE and EPI markers. Identification of these genes, many of which encode cell membrane proteins, will facilitate the isolation and characterization of early PE and EPI populations. Since it is so well established in the literature that mouse PGC specification is a post-implantation event, it was surprising for us to see activation of PGC markers as early as E3.25 preimplantation, and identify the newly found steady EPI markers as late germ cell markers. The discovery of such early activation of PGC markers has important implications in the derivation of germ cells from pluripotent cells (embryonic stem cells or induced pluripotent stem cells), since the initial stages of such derivation resemble early development. The early activation of PGC markers points out the difficulty of separating PGC cells from pluripotent populations. Collectively, our results suggest that the combining of the precision of single-cell omics data with dynamic analysis of time-series data can establish the timing of some developmental stages as earlier than previously thought. LARGE-SCALE DATA: Not applicable. STUDY FUNDING AND COMPETING INTERESTS: This work was supported by grants DFG15/14 and DFG15/020 from Diputación Foral de Gipuzkoa (Spain), and grant II14/00016 from I + D + I National Plan 2013-2016 (Spain) and FEDER funds. The authors declare no conflict of interest.

Galunisertib downregulates mutant type I collagen expression and promotes MSCs osteogenesis in pediatric osteogenesis imperfecta.

Qualitative alterations in type I collagen due to pathogenic variants in the COL1A1 or COL1A2 genes, result in moderate and severe Osteogenesis Imperfecta (OI), a rare disease characterized by bone fragility. The TGF-ß signaling pathway is overactive in OI patients and certain OI mouse models, and inhibition of TGF-ß through anti-TGF-ß monoclonal antibody therapy in phase I clinical trials in OI adults is rendering encouraging results. However, the impact of TGF-ß inhibition on osteogenic differentiation of mesenchymal stem cells from OI patients (OI-MSCs) is unknown. The following study demonstrates that pediatric skeletal OI-MSCs have imbalanced osteogenesis favoring the osteogenic commitment. Galunisertib, a small molecule inhibitor (SMI) that targets the TGF-ß receptor I (TßRI), favored the final osteogenic maturation of OI-MSCs. Mechanistically, galunisertib downregulated type I collagen expression in OI-MSCs, with greater impact on mutant type I collagen, and concomitantly, modulated the expression of unfolded protein response (UPR) and autophagy markers. In vivo, galunisertib improved trabecular bone parameters only in female oim/oim mice. These results further suggest that type I collagen is a tunable target within the bone ECM that deserves investigation and that the SMI, galunisertib, is a promising new candidate for the anti-TGF-ß targeting for the treatment of OI.

rhBMP-2 induces terminal differentiation of human bone marrow mesenchymal stromal cells only by synergizing with other signals.

BACKGROUND: Recombinant human bone morphogenetic protein 2 (rhBMP-2) and human bone marrow mesenchymal stromal cells (hBM-MSCs) have been thoroughly studied for research and translational bone regeneration purposes. rhBMP-2 induces bone formation in vivo, and hBM-MSCs are its target, bone-forming cells. In this article, we studied how rhBMP-2 drives the multilineage differentiation of hBM-MSCs both in vivo and in vitro. METHODS: rhBMP-2 and hBM-MSCs were tested in an in vivo subcutaneous implantation model to assess their ability to form mature bone and undergo multilineage differentiation. Then, the hBM-MSCs were treated in vitro with rhBMP-2 for short-term or long-term cell-culture periods, alone or in combination with osteogenic, adipogenic or chondrogenic media, aiming to determine the role of rhBMP-2 in these differentiation processes. RESULTS: The data indicate that hBM-MSCs respond to rhBMP-2 in the short term but fail to differentiate in long-term culture conditions; these cells overexpress the rhBMP-2 target genes DKK1, HEY-1 and SOST osteogenesis inhibitors. However, in combination with other differentiation signals, rhBMP-2 acts as a potentiator of multilineage differentiation, not only of osteogenesis but also of adipogenesis and chondrogenesis, both in vitro and in vivo. CONCLUSIONS: Altogether, our data indicate that rhBMP-2 alone is unable to induce in vitro osteogenic terminal differentiation of hBM-MSCs, but synergizes with other signals to potentiate multiple differentiation phenotypes. Therefore, rhBMP-2 triggers on hBM-MSCs different specific phenotype differentiation depending on the signalling environment.

Single-cell transcriptomic analysis to identify endomembrane regulation of metalloproteins and motor proteins in autoimmunity.

TMEM230 promotes antigen processing, trafficking, and presentation by regulating the endomembrane system of membrane bound organelles (lysosomes, proteosomes and mitochondria) and phagosomes. Activation of the immune system requires trafficking of various cargos between the endomembrane system and cell plasma membrane. The Golgi apparatus is the hub of the endomembrane system and essential for the generation, maintenance, recycling, and trafficking of the components of the endomembrane system itself and immune system. Intracellular trafficking and secretion of immune system components depend on mitochondrial metalloproteins for ATP synthesis that powers motor protein transport of endomembrane cargo. Glycan modifying enzyme genes and motor proteins are essential for the activation of the immune system and trafficking of antigens between the endomembrane system and the plasma membrane. Recently, TMEM230 was identified as co-regulated with RNASET2 in lysosomes and with metalloproteins in various cell types and organelles, including mitochondria in autoimmune diseases. Aberrant metalloproteinase secretion by motor proteins is a major contributor to tissue remodeling of synovial membrane and joint tissue destruction in rheumatoid arthritis (RA) by promoting infiltration of blood vessels, bone erosion, and loss of cartilage by phagocytes. In this study, we identified that specific glycan processing enzymes are upregulated in certain cell types (fibroblast or endothelial cells) that function in destructive tissue remodeling in rheumatoid arthritis compared to osteoarthritis (OA). TMEM230 was identified as a regulator in the secretion of metaloproteinases and heparanase necessary tissue remodeling in OA and RA. In dendritic (DC), natural killer and T cells, TMEM230 was expressed at low or no levels in RA compared to OA. TMEM230 expression in DC likely is necessary for regulatory or helper T cells to maintain tolerance to self-antigens and prevent susceptibility to autoimmune disease. To identify how TMEM230 and the endomembrane system contribute to autoimmunity we investigated, glycan modifying enzymes, metalloproteinases and motor protein genes co-regulated with or regulated by TMEM230 in synovial tissue by analyzing published single cell transcriptomic datasets from RA patient derived synovial tissue.

Systemic Lupus Erythematosus Patients with DNASE1L3·Deficiency Have a Distinctive and Specific Genic Circular DNA Profile in Plasma.

Cell-free (cf) extrachromosomal circular DNA (eccDNA) has a potential clinical application as a biomarker. Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with a complex immunological pathogenesis, associated with autoantibody synthesis. A previous study found that SLE patients with deoxyribonuclease 1-like 3 (DNASE1L3) deficiency exhibit changes in the frequency of short and long eccDNA in plasma compared to controls. Here, using the DifCir method for differential analysis of short-read sequenced purified eccDNA data based on the split-read signal of the eccDNA on circulomics data, we show that SLE patients with DNASE1L3 deficiency have a distinctive profile of eccDNA excised by gene regions compared to controls. Moreover, this profile is specific; cf-eccDNA from the top 93 genes is detected in all SLE with DNASE1L3 deficiency samples, and none in the control plasma. The top protein coding gene producing eccDNA-carrying gene fragments is the transcription factor BARX2, which is involved in skeletal muscle morphogenesis and connective tissue development. The top gene ontology terms are 'positive regulation of torc1 signaling' and 'chondrocyte development'. The top Harmonizome terms are 'lymphopenia', 'metabolic syndrome x', 'asthma', 'cardiovascular system disease', 'leukemia', and 'immune system disease'. Here, we show that gene associations of cf-eccDNA can serve as a biomarker in the autoimmune rheumatic diseases.

Cell-Free Genic Extrachromosomal Circular DNA Profiles of DNase Knockouts Associated with Systemic Lupus Erythematosus and Relation with Common Fragile Sites.

Cell-free extrachromosomal circular DNA (cf-eccDNA) has been proposed as a promising early biomarker for disease diagnosis, progression and drug response. Its established biomarker features are changes in the number and length distribution of cf-eccDNA. Another novel promising biomarker is a set of eccDNA excised from a panel of genes specific to a condition compared to a control. Deficiencies in two endonucleases that specifically target DNA, Dnase1 and Dnase1l3, are associated with systemic lupus erythematosus (SLE). To study the genic eccDNA profiles in the case of their deficiencies, we mapped sequenced eccDNA data from plasma, liver and buffy coat from Dnase1 and Dnase1l3 knockouts (KOs), and wild type controls in mouse. Next, we performed an eccDNA differential analysis between KO and control groups using our DifCir algorithm. We found a specific genic cf-eccDNA fingerprint of the Dnase1l3 group compared to the wild type controls involving 131 genes; 26% of them were associated with human chromosomal fragile sites (CFSs) and with a statistically significant enrichment of CFS-associated genes. We found six genes in common with the genic cf-eccDNA profile of SLE patients with DNASE1L3 deficiency, namely Rorb, Mvb12b, Osbpl10, Fto, Tnik and Arhgap10; all of them were specific and present in all human plasma samples, and none of them were associated with CFSs. A not so distinctive genic cf-eccDNA difference involving only seven genes was observed in the case of the Dnase1 group compared to the wild type. In tissue-liver and buffy coat-we did not detect the same genic eccDNA difference observed in the plasma samples. These results point to a specific role of a set of genic eccDNA in plasma from DNase KOs, as well as a relation with CFS genes, confirming the promise of the genic cf-eccDNA in studying diseases and the need for further research on the relationship between eccDNA and CFSs.

Circular and Circulating DNA in Inflammatory Bowel Disease: From Pathogenesis to Potential Molecular Therapies.

Inflammatory bowel diseases (IBD), including Crohn's Disease (CD) and Ulcerative Colitis (UC) are chronic multifactorial disorders which affect the gastrointestinal tract with variable extent. Despite extensive research, their etiology and exact pathogenesis are still unknown. Cell-free DNAs (cfDNAs) are defined as any DNA fragments which are free from the origin cell and able to circulate into the bloodstream with or without microvescicles. CfDNAs are now being increasingly studied in different human diseases, like cancer or inflammatory diseases. However, to date it is unclear how IBD etiology is linked to cfDNAs in plasma. Extrachromosomal circular DNA (eccDNA) are non-plasmidic, nuclear, circular and closed DNA molecules found in all eukaryotes tested. CfDNAs appear to play an important role in autoimmune diseases, inflammatory processes, and cancer; recently, interest has also grown in IBD, and their role in the pathogenesis of IBD has been suggested. We now suggest that eccDNAs also play a role in IBD. In this review, we have comprehensively collected available knowledge in literature regarding cfDNA, eccDNA, and structures involving them such as neutrophil extracellular traps and exosomes, and their role in IBD. Finally, we focused on old and novel potential molecular therapies and drug delivery systems, such as nanoparticles, for IBD treatment.

A distinct circular DNA profile intersects with proteome changes in the genotoxic stress-related hSOD1G93A model of ALS.

BACKGROUND: Numerous genes, including SOD1, mutated in familial and sporadic amyotrophic lateral sclerosis (f/sALS) share a role in DNA damage and repair, emphasizing genome disintegration in ALS. One possible outcome of chromosomal instability and repair processes is extrachromosomal circular DNA (eccDNA) formation. Therefore, eccDNA might accumulate in f/sALS with yet unknown function. METHODS: We combined rolling circle amplification with linear DNA digestion to purify eccDNA from the cervical spinal cord of 9 co-isogenic symptomatic hSOD1G93A mutants and 10 controls, followed by deep short-read sequencing. We mapped the eccDNAs and performed differential analysis based on the split read signal of the eccDNAs, referred as DifCir, between the ALS and control specimens, to find differentially produced per gene circles (DPpGC) in the two groups. Compared were eccDNA abundances, length distributions and genic profiles. We further assessed proteome alterations in ALS by mass spectrometry, and matched the DPpGCs with differentially expressed proteins (DEPs) in ALS. Additionally, we aligned the ALS-specific DPpGCs to ALS risk gene databases. RESULTS: We found a six-fold enrichment in the number of unique eccDNAs in the genotoxic ALS-model relative to controls. We uncovered a distinct genic circulome profile characterized by 225 up-DPpGCs, i.e., genes that produced more eccDNAs from distinct gene sequences in ALS than under control conditions. The inter-sample recurrence rate was at least 89% for the top 6 up-DPpGCs. ALS proteome analyses revealed 42 corresponding DEPs, of which 19 underlying genes were itemized for an ALS risk in GWAS databases. The up-DPpGCs and their DEP tandems mainly impart neuron-specific functions, and gene set enrichment analyses indicated an overrepresentation of the adenylate cyclase modulating G protein pathway. CONCLUSIONS: We prove, for the first time, a significant enrichment of eccDNA in the ALS-affected spinal cord. Our triple circulome, proteome and genome approach provide indication for a potential importance of certain eccDNAs in ALS neurodegeneration and a yet unconsidered role as ALS biomarkers. The related functional pathways might open up new targets for therapeutic intervention.

Single-cell RNA sequencing of human non-hematopoietic bone marrow cells reveals a unique set of inter-species conserved biomarkers for native mesenchymal stromal cells.

BACKGROUND: Native bone marrow (BM) mesenchymal stem/stromal cells (BM-MSCs) participate in generating and shaping the skeleton and BM throughout the lifespan. Moreover, BM-MSCs regulate hematopoiesis by contributing to the hematopoietic stem cell niche in providing critical cytokines, chemokines and extracellular matrix components. However, BM-MSCs contain a heterogeneous cell population that remains ill-defined. Although studies on the taxonomy of native BM-MSCs in mice have just started to emerge, the taxonomy of native human BM-MSCs remains unelucidated. METHODS: By using single-cell RNA sequencing (scRNA-seq), we aimed to define a proper taxonomy for native human BM non-hematopoietic subsets including endothelial cells (ECs) and mural cells (MCs) but with a focal point on MSCs. To this end, transcriptomic scRNA-seq data were generated from 5 distinct BM donors and were analyzed together with other transcriptomic data and with computational biology analyses at different levels to identify, characterize and classify distinct native cell subsets with relevant biomarkers. RESULTS: We could ascribe novel specific biomarkers to ECs, MCs and MSCs. Unlike ECs and MCs, MSCs exhibited an adipogenic transcriptomic pattern while co-expressing genes related to hematopoiesis support and multilineage commitment potential. Furthermore, by a comparative analysis of scRNA-seq of BM cells from humans and mice, we identified core genes conserved in both species. Notably, we identified MARCKS, CXCL12, PDGFRA, and LEPR together with adipogenic factors as archetypal biomarkers of native MSCs within BM. In addition, our data suggest some complex gene nodes regulating critical biological functions of native BM-MSCs together with a preferential commitment toward an adipocyte lineage. CONCLUSIONS: Overall, our taxonomy for native BM non-hematopoietic compartment provides an explicit depiction of gene expression in human ECs, MCs and MSCs at single-cell resolution. This analysis helps enhance our understanding of the phenotype and the complexity of biological functions of native human BM-MSCs.

Human skeletal muscle organoids model fetal myogenesis and sustain uncommitted PAX7 myogenic progenitors.

In vitro culture systems that structurally model human myogenesis and promote PAX7+ myogenic progenitor maturation have not been established. Here we report that human skeletal muscle organoids can be differentiated from induced pluripotent stem cell lines to contain paraxial mesoderm and neuromesodermal progenitors and develop into organized structures reassembling neural plate border and dermomyotome. Culture conditions instigate neural lineage arrest and promote fetal hypaxial myogenesis toward limb axial anatomical identity, with generation of sustainable uncommitted PAX7 myogenic progenitors and fibroadipogenic (PDGFRa+) progenitor populations equivalent to those from the second trimester of human gestation. Single-cell comparison to human fetal and adult myogenic progenitor /satellite cells reveals distinct molecular signatures for non-dividing myogenic progenitors in activated (CD44High/CD98+/MYOD1+) and dormant (PAX7High/FBN1High/SPRY1High) states. Our approach provides a robust 3D in vitro developmental system for investigating muscle tissue morphogenesis and homeostasis.

Humans contains around 650 skeletal muscles which allow the body to move around and maintain its posture. Skeletal muscles are made up of individual cells that bundle together into highly organized structures. If this group of muscles fail to develop correctly in the embryo and/or fetus, this can lead to muscular disorders that can make it painful and difficult to move. One way to better understand how skeletal muscles are formed, and how this process can go wrong, is to grow them in the laboratory. This can be achieved using induced pluripotent stem cells (iPSCs), human adult cells that have been 'reprogrammed' to behave like cells in the embryo that can develop in to almost any cell in the body. The iPSCs can then be converted into specific cell types in the laboratory, including the cells that make up skeletal muscle. Here, Mavrommatis et al. created a protocol for developing iPSCs into three-dimensional organoids which resemble how cells of the skeletal muscle look and arrange themselves in the fetus. To form the skeletal muscle organoid, Mavrommatis et al. treated iPSCs that were growing in a three-dimensional environment with various factors that are found early on in development. This caused the iPSCs to organize themselves in to embryonic and fetal structures that will eventually give rise to the parts of the body that contain skeletal muscle, such as the limbs. Within the organoid were cells that produced Pax7, a protein commonly found in myogenic progenitors that specifically mature into skeletal muscle cells in the fetus. Pax 7 is also present in 'satellite cells' that help to regrow damaged skeletal muscle in adults. Indeed, Mavrommatis et al. found that the myogenic progenitors produced by the organoid were able to regenerate muscle when transplanted in to adult mice. These findings suggest that this organoid protocol can generate cells that will give rise to skeletal muscle. In the future, these lab-grown progenitors could potentially be created from cells isolated from patients and used to repair muscle injuries. The organoid model could also provide new insights in to how skeletal muscles develop in the fetus, and how genetic mutations linked with muscular disorders disrupt this process.

The common incidence-age multistep model of neurodegenerative diseases revisited: wider general age range of incidence corresponds to fewer disease steps.

BACKGROUND: Previously, we collected age-stratified incidence data of 404 epidemiological datasets of 10 neurodegenerative diseases (NDs), namely Amyotrophic Lateral Sclerosis (ALS), Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD), Fronto Temporal Dementia (FTD), Dementia with Lewy Bodies (DLB), Parkinsonism (PDM), Parkinson's disease with Dementia (PDD), Creutzfeldt-Jakob disease (CJD), and Multiple Sclerosis (MS). We tested whether each ND follows a multistep model, found the number of steps necessary for the onset of each ND, found the number of common steps with other NDs and the number of specific steps of each ND, and built a parsimony tree of the genealogy of the NDs. The tree disclosed three groups of NDs: the stem NDs with less than 3 steps; the trunk NDs with 5-7 steps; and the crown NDs with more than 7 steps. METHODS: We made a multidimensional reduction of the previously collected age-stratified incidence epidemiological data of the 10 NDs. We studied the general range of incidence of the 10 NDs using the age- and sex-stratified incidence data. First, we calculated the log of the incidence versus the log of the age for each ND. Next, we calculated the age intervals of the spread of the incidence of each ND. We calculated the regression of the steps obtained with the multistep model versus the age of incidence of the NDs. RESULTS: We found that the number of steps of the NDs is inversely correlated with the age of incidence of the NDs, and calculated the number of years required for a single step for each ND. Based on these results, we extended the genealogy tree model of the NDs to account for the time needed for a ND step to occur. CONCLUSION: The extended genealogy tree disclosed three groups of NDs according to the estimated time needed for a step to occur: the stem ND, HD, with 32.5 years/step, the trunk NDs ALS, FTD, PD and CJD, with 6.7-13.7 years/step; and the crown NDs PDM, PDD, AD and DLB, with 2.3-3.8 years/step. Thus, the NDs cluster into three groups according to both the number of steps and the number of years for a step to occur.