Iron supplementation is sufficient to rescue skeletal muscle mass and function in cancer cachexia.

Cachexia is a wasting syndrome characterized by devastating skeletal muscle atrophy that dramatically increases mortality in various diseases, most notably in cancer patients with a penetrance of up to 80%. Knowledge regarding the mechanism of cancer-induced cachexia remains very scarce, making cachexia an unmet medical need. In this study, we discovered strong alterations of iron metabolism in the skeletal muscle of both cancer patients and tumor-bearing mice, characterized by decreased iron availability in mitochondria. We found that modulation of iron levels directly influences myotube size in vitro and muscle mass in otherwise healthy mice. Furthermore, iron supplementation was sufficient to preserve both muscle function and mass, prolong survival in tumor-bearing mice, and even rescues strength in human subjects within an unexpectedly short time frame. Importantly, iron supplementation refuels mitochondrial oxidative metabolism and energy production. Overall, our findings provide new mechanistic insights in cancer-induced skeletal muscle wasting, and support targeting iron metabolism as a potential therapeutic option for muscle wasting diseases.

The acetyltransferase p300 is recruited in trans to multiple enhancer sites by lncSmad7.

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

Hair follicle stem cell progeny heal blisters while pausing skin development.

Injury in adult tissue generally reactivates developmental programs to foster regeneration, but it is not known whether this paradigm applies to growing tissue. Here, by employing blisters, we show that epidermal wounds heal at the expense of skin development. The regenerated epidermis suppresses the expression of tissue morphogenesis genes accompanied by delayed hair follicle (HF) growth. Lineage tracing experiments, cell proliferation dynamics, and mathematical modeling reveal that the progeny of HF junctional zone stem cells, which undergo a morphological transformation, repair the blisters while not promoting HF development. In contrast, the contribution of interfollicular stem cell progeny to blister healing is small. These findings demonstrate that HF development can be sacrificed for the sake of epidermal wound regeneration. Our study elucidates the key cellular mechanism of wound healing in skin blistering diseases.

Identification of altered biological processes in heterogeneous RNA-sequencing data by discretization of expression profiles.

Heterogeneity is a fundamental feature of complex phenotypes. So far, genomic screenings have profiled thousands of samples providing insights into the transcriptome of the cell. However, disentangling the heterogeneity of these transcriptomic Big Data to identify defective biological processes remains challenging. Here we present GSECA, a method exploiting the bimodal behavior of RNA-sequencing gene expression profiles to identify altered gene sets in heterogeneous patient cohorts. Using simulated and experimental RNA-sequencing data sets, we show that GSECA provides higher performances than other available algorithms in detecting truly altered biological processes in large cohorts. Applied to 5941 samples from 14 different cancer types, GSECA correctly identified the alteration of the PI3K/AKT signaling pathway driven by the somatic loss of PTEN and verified the emerging role of PTEN in modulating immune-related processes. In particular, we showed that, in prostate cancer, PTEN loss appears to establish an immunosuppressive tumor microenvironment through the activation of STAT3, and low PTEN expression levels have a detrimental impact on patient disease-free survival. GSECA is available at https://github.com/matteocereda/GSECA.

Genome-Wide Analysis of Smad7-Mediated Transcription in Mouse Embryonic Stem Cells.

Smad7 has been identified as a negative regulator of the transforming growth factor TGF-ß pathway by direct interaction with the TGF-ß type I receptor (TßR-I). Although Smad7 has also been shown to play TGF-ß unrelated functions in the cytoplasm and in the nucleus, a comprehensive analysis of its nuclear function has not yet been performed. Here, we show that in ESCs Smad7 is mainly nuclear and acts as a general transcription factor regulating several genes unrelated to the TGF-ß pathway. Loss of Smad7 results in the downregulation of several key stemness master regulators, including Pou5f1 and Zfp42, and in the upregulation of developmental genes, with consequent loss of the stem phenotype. Integrative analysis of genome-wide mapping data for Smad7 and ESC self-renewal and pluripotency transcriptional regulators revealed that Smad7 co-occupies promoters of highly expressed key stemness regulators genes, by binding to a specific consensus response element NCGGAAMM. Altogether, our data establishes Smad7 as a new, integral component of the regulatory circuitry that controls ESC identity.

Expression and Function of C1orf132 Long-Noncoding RNA in Breast Cancer Cell Lines and Tissues.

miR-29b2 and miR-29c play a suppressive role in breast cancer progression. C1orf132 (also named MIR29B2CHG) is the host gene for generating both microRNAs. However, the region also expresses longer transcripts with unknown functions. We employed bioinformatics and experimental approaches to decipher C1orf132 expression and function in breast cancer tissues. We also used the CRISPR/Cas9 technique to excise a predicted C1orf132 distal promoter and followed the behavior of the edited cells by real-time PCR, flow cytometry, migration assay, and RNA-seq techniques. We observed that C1orf132 long transcript is significantly downregulated in triple-negative breast cancer. We also identified a promoter for the longer transcripts of C1orf132 whose functionality was demonstrated by transfecting MCF7 cells with a C1orf132 promoter-GFP construct. Knocking-out the promoter by means of CRISPR/Cas9 revealed no alterations in the expression of the neighboring genes CD46 and CD34, while the expression of miR-29c was reduced by half. Furthermore, the promoter knockout elevated the migration ability of the edited cells. RNA sequencing revealed many up- and downregulated genes involved in various cellular pathways, including epithelial to mesenchymal transition and mammary gland development pathways. Altogether, we are reporting here the existence of an additional/distal promoter with an enhancer effect on miR-29 generation and an inhibitory effect on cell migration.

Involvement of N4BP2L1, PLEKHA4, and BEGAIN genes in breast cancer and muscle cell development.

Patients with breast cancer show altered expression of genes within the pectoralis major skeletal muscle cells of the breast. Through analyses of The Cancer Genome Atlas (TCGA)-breast cancer (BRCA), we identified three previously uncharacterized putative novel tumor suppressor genes expressed in normal muscle cells, whose expression was downregulated in breast tumors. We found that NEDD4 binding protein 2-like 1 (N4BP2L1), pleckstrin homology domain-containing family A member 4 (PLEKHA4), and brain-enriched guanylate kinase-associated protein (BEGAIN) that are normally highly expressed in breast myoepithelial cells and smooth muscle cells were significantly downregulated in breast tumor tissues of a cohort of 50 patients with this cancer. Our data revealed that the low expression of PLEKHA4 in patients with menopause below 50 years correlated with a higher risk of breast cancer. Moreover, we identified N4BP2L1 and BEGAIN as potential biomarkers of HER2-positive breast cancer. Furthermore, low BEGAIN expression in breast cancer patients with blood fat, heart problems, and diabetes correlated with a higher risk of this cancer. In addition, protein and RNA expression analysis of TCGA-BRCA revealed N4BP2L1 as a promising diagnostic protein biomarker in breast cancer. In addition, the in silico data of scRNA-seq showed high expression of these genes in several cell types of normal breast tissue, including breast myoepithelial cells and smooth muscle cells. Thus, our results suggest their possible tumor-suppressive function in breast cancer and muscle development.

3plex enables deep computational investigation of triplex forming lncRNAs.

Long non-coding RNAs (lncRNAs) regulate gene expression through different molecular mechanisms, including DNA binding via the formation of RNA:DNA:DNA triple helices (TPXs). Despite the increasing amount of experimental evidence, TPXs investigation remains challenging. Here we present 3plex, a software able to predict TPX interactions in silico. Given an RNA sequence and a set of DNA sequences, 3plex integrates 1) Hoogsteen pairing rules that describe the biochemical interactions between RNA and DNA nucleotides, 2) RNA secondary structure prediction and 3) determination of the TPX thermal stability derived from a collection of TPX experimental evidences. We systematically collected and uniformly re-analysed published experimental lncRNA binding sites on human and mouse genomes. We used these data to evaluate 3plex performance and showed that its specific features allow a reliable identification of TPX interactions. We compared 3plex with the other available software and obtained comparable or even better accuracy at a fraction of the computation time. Interestingly, by inspecting collected data with 3plex we found that TPXs tend to be shorter and more degenerated than previously expected and that the majority of analysed lncRNAs can directly bind to the genome by TPX formation. Those results suggest that an important fraction of lncRNAs can exert its biological function through this mechanism. The software is available at https://github.com/molinerisLab/3plex.

DNMT3B supports meso-endoderm differentiation from mouse embryonic stem cells.

The correct establishment of DNA methylation patterns during mouse early development is essential for cell fate specification. However, the molecular targets as well as the mechanisms that determine the specificity of the de novo methylation machinery during differentiation are not completely elucidated. Here we show that the DNMT3B-dependent DNA methylation of key developmental regulatory regions at epiblast-like cells (EpiLCs) provides an epigenetic priming that ensures flawless commitment at later stages. Using in vitro stem cell differentiation and loss of function experiments combined with high-throughput genome-wide bisulfite-, bulk-, and single cell RNA-sequencing we dissected the specific role of DNMT3B in cell fate. We identify DNMT3B-dependent regulatory elements on the genome which, in Dnmt3b knockout (3BKO), impair the differentiation into meso-endodermal (ME) progenitors and redirect EpiLCs towards the neuro-ectodermal lineages. Moreover, ectopic expression of DNMT3B in 3BKO re-establishes the DNA methylation of the master regulator Sox2 super-enhancer, downmodulates its expression, and restores the expression of ME markers. Taken together, our data reveal that DNMT3B-dependent methylation at the epiblast stage is essential for the priming of the meso-endodermal lineages and provide functional characterization of the de novo DNMTs during EpiLCs lineage determination.