The Impact of miR-155-5p on Myotube Differentiation: Elucidating Molecular Targets in Skeletal Muscle Disorders.

MicroRNAs are small regulatory molecules that control gene expression. An emerging property of muscle miRNAs is the cooperative regulation of transcriptional and epitranscriptional events controlling muscle phenotype. miR-155 has been related to muscular dystrophy and muscle cell atrophy. However, the function of miR-155 and its molecular targets in muscular dystrophies remain poorly understood. Through in silico and in vitro approaches, we identify distinct transcriptional profiles induced by miR-155-5p in muscle cells. The treated myotubes changed the expression of 359 genes (166 upregulated and 193 downregulated). We reanalyzed muscle transcriptomic data from dystrophin-deficient patients and detected overlap with gene expression patterns in miR-155-treated myotubes. Our analysis indicated that miR-155 regulates a set of transcripts, including Aldh1l, Nek2, Bub1b, Ramp3, Slc16a4, Plce1, Dync1i1, and Nr1h3. Enrichment analysis demonstrates 20 targets involved in metabolism, cell cycle regulation, muscle cell maintenance, and the immune system. Moreover, digital cytometry confirmed a significant increase in M2 macrophages, indicating miR-155's effects on immune response in dystrophic muscles. We highlight a critical miR-155 associated with disease-related pathways in skeletal muscle disorders.

The Transcriptomic Landscape of Age-Induced Changes in Human Visceral Fat and the Predicted Omentum-Liver Connectome in Males.

Aging causes alterations in body composition. Specifically, visceral fat mass increases with age and is associated with age-related diseases. The pathogenic potential of visceral fat accumulation has been associated with its anatomical location and metabolic activity. Visceral fat may control systemic metabolism by secreting molecules that act in distal tissues, mainly the liver, through the portal vein. Currently, little is known about age-related changes in visceral fat in humans. Aiming to identify molecular and cellular changes occurring with aging in the visceral fat of humans, we analyzed publicly available transcriptomic data of 355 omentum samples from the Genotype-Tissue Expression portal (GTEx) of 20-79-year-old males and females. We identified the functional enrichment of genes associated with aging, inferred age-related changes in visceral fat cellularity by deconvolution analysis, profiled the senescence-associated secretory phenotype of visceral adipose tissue, and predicted the connectivity of the age-induced visceral fat secretome with the liver. We demonstrate that age induces alterations in visceral fat cellularity, synchronous to changes in metabolic pathways and a shift toward a pro-inflammatory secretory phenotype. Furthermore, our approach identified candidates such as ADIPOQ-ADIPOR1/ADIPOR2, FCN2-LPR1, and TF-TFR2 to mediate visceral fat-liver crosstalk in the context of aging. These findings cast light on how alterations in visceral fat with aging contribute to liver dysfunction and age-related disease etiology.

Low muscle mass in lung cancer is associated with an inflammatory and immunosuppressive tumor microenvironment.

BACKGROUND: Computed tomographies (CT) are useful for identifying muscle loss in non-small lung cancer (NSCLC) cachectic patients. However, we lack consensus on the best cutoff point for pectoralis muscle loss. We aimed to characterize NSCLC patients based on muscularity, clinical data, and the transcriptional profile from the tumor microenvironment to build a cachexia classification model. METHODS: We used machine learning to generate a muscle loss prediction model, and the tumor's cellular and transcriptional profile was characterized in patients with low muscularity. First, we measured the pectoralis muscle area (PMA) of 211 treatment-naive NSCLC patients using CT available in The Cancer Imaging Archive. The cutoffs were established using machine learning algorithms (CART and Cutoff Finder) on PMA, clinical, and survival data. We evaluated the prediction model in a validation set (36 NSCLC). Tumor RNA-Seq (GSE103584) was used to profile the transcriptome and cellular composition based on digital cytometry. RESULTS: CART demonstrated that a lower PMA was associated with a high risk of death (HR = 1.99). Cutoff Finder selected PMA cutoffs separating low-muscularity (LM) patients based on the risk of death (P-value = 0.003; discovery set). The cutoff presented 84% of success in classifying low muscle mass. The high risk of LM patients was also found in the validation set. Tumor RNA-Seq revealed 90 upregulated secretory genes in LM that potentially interact with muscle cell receptors. The LM upregulated genes enriched inflammatory biological processes. Digital cytometry revealed that LM patients presented high proportions of cytotoxic and exhausted CD8+ T cells. CONCLUSIONS: Our prediction model identified cutoffs that distinguished patients with lower PMA and survival with an inflammatory and immunosuppressive TME enriched with inflammatory factors and CD8+ T cells.

RNA-seq reveals that anti-obesity irisin and triiodothyronine (T3) hormones differentially affect the purinergic signaling transcriptomics in differentiated human adipocytes.

The anti-obesity thyroid hormone, triiodothyronine (T3), and irisin, an exercise- and/or cold-induced myokine, stimulate thermogenesis and energy consumption while decreasing lipid accumulation. The involvement of ATP signaling in adipocyte cell function and obesity has attracted increasing attention, but the crosstalk between the purinergic signaling cascade and anti-obesity hormones lacks experimental evidence. In this study, we investigated the effects of T3 and irisin in the transcriptomics of membrane-bound purinoceptors, ectonucleotidase enzymes and nucleoside transporters participating in the purinergic signaling in cultured human adipocytes. The RNA-seq analysis revealed that differentiated adipocytes express high amounts of ADORA1, P2RY11, P2RY12, and P2RX6 gene transcripts, along with abundant levels of transcriptional products encoding to purine metabolizing enzymes (ENPP2, ENPP1, NT5E, ADA and ADK) and transporters (SLC29A1, SCL29A2). The transcriptomics of purinergic signaling markers changed in parallel to the upsurge of "browning" adipocyte markers, like UCP1 and P2RX5, after treatment with T3 and irisin. Upregulation of ADORA1, ADORA2A and P2RX4 gene transcription was obtained with irisin, whereas T3 preferentially upregulated NT5E, SLC29A2 and P2RY11 genes. Irisin was more powerful than T3 towards inhibition of the leptin gene transcription, the SCL29A1 gene encoding for the ENT1 transporter, the E-NPP2 (autotaxin) gene, and genes that encode for two ADP-sensitive P2Y receptors, P2RY1 and P2RY12. These findings indicate that anti-obesity irisin and T3 hormones differentially affect the purinergic signaling transcriptomics, which might point towards new directions for the treatment of obesity and related metabolic disorders that are worth to be pursued in future functional studies.

Identification of Novel Genes Associated with Fish Skeletal Muscle Adaptation during Fasting and Refeeding Based on a Meta-Analysis.

The regulation of the fish phenotype and muscle growth is influenced by fasting and refeeding periods, which occur in nature and are commonly applied in fish farming. However, the regulators associated with the muscle responses to these manipulations of food availability have not been fully characterized. We aimed to identify novel genes associated with fish skeletal muscle adaptation during fasting and refeeding based on a meta-analysis. Genes related to translational and proliferative machinery were investigated in pacus (Piaractus mesopotamicus) subjected to fasting (four and fifteen days) and refeeding (six hours, three and fifteen days). Our results showed that different fasting and refeeding periods modulate the expression of the genes mtor, rps27a, eef1a2, and cdkn1a. These alterations can indicate the possible protection of the muscle phenotype, in addition to adaptive responses that prioritize energy and substrate savings over cell division, a process regulated by ccnd1. Our study reveals the potential of meta-analysis for the identification of muscle growth regulators and provides new information on muscle responses to fasting and refeeding in fish that are of economic importance to aquaculture.

Interplay Between Immune and Cancer-Associated Fibroblasts: A Path to Target Metalloproteinases in Penile Cancer.

Extracellular matrix (ECM) remodeling and inflammation have been reported in penile carcinomas (PeCa). However, the cell types and cellular crosstalk involved in PeCa are unexplored. We aimed to characterize the complexity of cells and pathways involved in the tumor microenvironment (TME) in PeCa and propose target molecules associated with the TME. We first investigated the prognostic impact of cell types with a secretory profile to identify drug targets that modulate TME-enriched cells. The secretome analysis using the PeCa transcriptome revealed the enrichment of inflammation and extracellular matrix pathways. Twenty-three secreted factors were upregulated, mainly collagens and matrix metalloproteinases (MMPs). The deregulation of collagens and MMPs was confirmed by Quantitative reverse transcription - polymerase chain reaction (RT-qPCR). Further, the deconvolution method (digital cytometry) of the bulk samples revealed a high proportion of macrophages and dendritic cells (DCs) and B cells. Increased DCs and B cells were associated with better survival. A high proportion of cancer-associated fibroblasts (CAFs) was observed in low-survival patients. Patients with increased CAFs had decreased immune cell proportions. The treatment with the MMP inhibitor GM6001 in CAF cells derived from PeCa resulted in altered cell viability. We reported a crosstalk between immune cells and CAFs, and the proportion of these cell populations was associated with prognosis. We demonstrate that a drug targeting MMPs modulates CAFs, expanding the therapeutic options of PeCa.

Fibronectin Modulates the Expression of miRNAs in Prostate Cancer Cell Lines.

Prostate cancer (PCa) is a significant cause of cancer-related deaths among men and companion animals, such as dogs. However, despite its high mortality and incidence rates, the molecular mechanisms underlying this disease remain to be fully elucidated. Among the many factors involved in prostate carcinogenesis, the extracellular matrix (ECM) plays a crucial role. This ECM in the prostate is composed mainly of collagen fibers, reticular fibers, elastic fibers, proteoglycans and glycoproteins, such as fibronectin. Fibronectin is a glycoprotein whose dysregulation has been implicated in the development of multiple types of cancer, and it has been associated with cell migration, invasion, and metastasis. Furthermore, our research group has previously shown that fibronectin induces transcriptional changes by modulating the expression of protein coding genes in LNCaP cells. However, potential changes at the post-transcriptional level are still not well understood. This study investigated the impact of exposure to fibronectin on the expression of a key class of regulatory RNAs, the microRNAs (miRNAs), in prostate cancer cell lines LNCaP and PC-3. Five mammalian miRNAs (miR-21, miR-29b, miR-125b, miR-221, and miR-222) were differentially expressed after fibronectin exposure in prostate cell lines. The expression profile of hundreds of mRNAs predicted to be targeted by these miRNAs was analyzed using publicly available RNA-Sequencing data (GSE64025, GSE68645, GSE29155). Also, protein-protein interaction networks and enrichment analysis were performed to gain insights into miRNA biological functions. Altogether, these functional analyzes revealed that fibronectin exposure impacts the expression of miRNAs potentially involved in PCa causing changes in critical signaling pathways such as PI3K-AKT, and response to cell division, death, proliferation, and migration. The relationship here demonstrated between fibronectin exposure and altered miRNA expression improves the comprehension of PCa in both men and other animals, such as dogs, which naturally develop prostate cancer.

Effects of Triiodothyronine on Human Osteoblast-Like Cells: Novel Insights From a Global Transcriptome Analysis.

Background: Thyroid hormones play a significant role in bone development and maintenance, with triiodothyronine (T3) particularly being an important modulator of osteoblast differentiation, proliferation, and maintenance. However, details of the biological processes (BPs) and molecular pathways affected by T3 in osteoblasts remain unclear. Methods: To address this issue, primary cultures of human adipose-derived mesenchymal stem cells were subjected to our previously established osteoinduction protocol, and the resultant osteoblast-like cells were treated with 1 nm or 10 nm T3 for 72 h. RNA sequencing (RNA-Seq) was performed using the Illumina platform, and differentially expressed genes (DEGs) were identified from the raw data using Kallisto and DESeq2. Enrichment analysis of DEGs was performed against the Gene Ontology Consortium database for BP terms using the R package clusterProfiler and protein network analysis by STRING. Results: Approximately 16,300 genes were analyzed by RNA-Seq, with 343 DEGs regulated in the 1 nm T3 group and 467 upregulated in the 10 nm T3 group. Several independent BP terms related to bone metabolism were significantly enriched, with a number of genes shared among them (FGFR2, WNT5A, WNT3, ROR2, VEGFA, FBLN1, S1PR1, PRKCZ, TGFB3, and OSR1 for 1nM T3; and FZD1, SMAD6, NOG, NEO1, and ENG for 10 nm T3). An osteoblast-related search in the literature regarding this set of genes suggests that both T3 doses are unfavorable for osteoblast development, mainly hindering BMP and canonical and non-canonical WNT signaling. Conclusions: Therefore, this study provides new directions toward the elucidation of the mechanisms of T3 action on osteoblast metabolism, with potential future implications for the treatment of endocrine-related bone pathologies.

Prediction of Non-canonical Routes for SARS-CoV-2 Infection in Human Placenta Cells.

The SARS-CoV-2 is the causative agent of the COVID-19 pandemic. The data available about COVID-19 during pregnancy have demonstrated placental infection; however, the mechanisms associated with intrauterine transmission of SARS-CoV-2 is still debated. Intriguingly, while canonical SARS-CoV-2 cell entry mediators are expressed at low levels in placental cells, the receptors for viruses that cause congenital infections such as the cytomegalovirus and Zika virus are highly expressed in these cells. Here we analyzed the transcriptional profile (microarray and single-cell RNA-Seq) of proteins potentially interacting with coronaviruses to identify non- canonical mediators of SARS-CoV-2 infection and replication in the placenta. Despite low levels of the canonical cell entry mediators ACE2 and TMPRSS2, we show that cells of the syncytiotrophoblast, villous cytotrophoblast, and extravillous trophoblast co-express high levels of the potential non-canonical cell-entry mediators DPP4 and CTSL. We also found changes in the expression of DAAM1 and PAICS genes during pregnancy, which are translated into proteins also predicted to interact with coronaviruses proteins. These results provide new insight into the interaction between SARS-CoV-2 and host proteins that may act as non-canonical routes for SARS-CoV-2 infection and replication in the placenta cells.

Multidimensional Single-Nuclei RNA-Seq Reconstruction of Adipose Tissue Reveals Adipocyte Plasticity Underlying Thermogenic Response.

Adipose tissue has been classified based on its morphology and function as white, brown, or beige/brite. It plays an essential role as a regulator of systemic metabolism through paracrine and endocrine signals. Recently, multiple adipocyte subtypes have been revealed using RNA sequencing technology, going beyond simply defined morphology but also by their cellular origin, adaptation to metabolic stress, and plasticity. Here, we performed an in-depth analysis of publicly available single-nuclei RNAseq from adipose tissue and utilized a workflow template to characterize adipocyte plasticity, heterogeneity, and secretome profiles. The reanalyzed dataset led to the identification of different subtypes of adipocytes including three subpopulations of thermogenic adipocytes, and provided a characterization of distinct transcriptional profiles along the adipocyte trajectory under thermogenic challenges. This study provides a useful resource for further investigations regarding mechanisms related to adipocyte plasticity and trans-differentiation.

Drug Repositioning Based on the Reversal of Gene Expression Signatures Identifies TOP2A as a Therapeutic Target for Rectal Cancer.

Rectal cancer is a common disease with high mortality rates and limited therapeutic options. Here we combined the gene expression signatures of rectal cancer patients with the reverse drug-induced gene-expression profiles to identify drug repositioning candidates for cancer therapy. Among the predicted repurposable drugs, topoisomerase II inhibitors (doxorubicin, teniposide, idarubicin, mitoxantrone, and epirubicin) presented a high potential to reverse rectal cancer gene expression signatures. We showed that these drugs effectively reduced the growth of colorectal cancer cell lines closely representing rectal cancer signatures. We also found a clear correlation between topoisomerase 2A (TOP2A) gene copy number or expression levels with the sensitivity to topoisomerase II inhibitors. Furthermore, CRISPR-Cas9 and shRNA screenings confirmed that loss-of-function of the TOP2A has the highest efficacy in reducing cellular proliferation. Finally, we observed significant TOP2A copy number gains and increased expression in independent cohorts of rectal cancer patients. These findings can be translated into clinical practice to evaluate TOP2A status for targeted and personalized therapies based on topoisomerase II inhibitors in rectal cancer patients.

Germline variants in DNA repair genes are associated with young-onset head and neck cancer.

The genetic predisposition to head and neck carcinomas (HNSCC) and how the known risk factors (papillomavirus infection, alcohol, and tobacco consumption) contribute to the early-onset disease are barely explored. Although HNSCC at early onset is rare, its frequency is increasing in recent years. Germline and somatic variants were assessed to build a comprehensive genetic influence pattern in HNSCC predisposition and patient outcome. Whole-exome sequencing was performed in 45 oral and oropharynx carcinomas paired with normal samples of young adults (≤49 years). We found FANCG, CDKN2A, and TPP germline variants previously associated with HNSCC risk. At least one germline variant in DNA repair pathway genes was detected in 67% of cases. Germline and somatic variants (including copy number variations) in FAT1 gene were identified in 9 patients (20%) and 12 tumors (30%), respectively. Somatic variants were found in HNSCC associated genes, such as TP53, CDKN2A, and PIK3CA. To date, 55 of 521 cases from the large cohort of TCGA presented < 49 years old. A comparison between the somatic alterations of TCGA-HNSCC at early onset and our dataset revealed strong similarities. Protein-protein interaction analysis between somatic and germline altered genes revealed a central role of TP53. Altogether, germline alterations in DNA repair genes potentially contribute to an increased risk of developing HNSCC at early-onset, while FAT1 could impact the prognosis.

Preventive training does not interfere with mRNA-encoding myosin and collagen expression during pulmonary arterial hypertension.

To gain insight on the impact of preventive exercise during pulmonary arterial hypertension (PAH), we evaluated the gene expression of myosins and gene-encoding proteins associated with the extracellular matrix remodeling of right hypertrophied ventricles. We used 32 male Wistar rats, separated in four groups: Sedentary Control (S, n = 8); Control with Training (T, n = 8); Sedentary with Pulmonary Arterial Hypertension (SPAH, n = 8); and Pulmonary Arterial Hypertension with Training (TPAH, n = 8). All rats underwent a two-week adaptation period; T and TPAH group rats then proceeded to an eight-week training period on a treadmill. At the beginning of the 11th week, S and T groups received an intraperitoneal injection of saline, and SPAH and TPAH groups received an injection of monocrotaline (60 mg/kg). Rats in the T and TPAH groups then continued with the training protocol until the 13th week. We assessed exercise capacity, echocardiography analysis, Fulton's index, cross-sectional areas of cardiomyocytes, collagen content and types, and fractal dimension (FD). Transcript abundance of myosins and extracellular matrix genes were estimated through reverse transcription-quantitative PCR (RT-qPCR). When compared to the SPAH group, the TPAH group showed increases in functional capacity and pulmonary artery acceleration time/pulmonary ejection time ratio and decreases in Fulton's index and cross-sectional areas of myocyte cells. However, preventive exercise did not induce alterations in col1a1 and myh7 gene expression. Our findings demonstrate that preventive exercise improved functional capacity, reduced cardiac hypertrophy, and attenuated PH development without interfering in mRNA-encoding myosin and collagen expression during PAH.

Prediction of SARS-CoV Interaction with Host Proteins during Lung Aging Reveals a Potential Role for TRIB3 in COVID-19.

COVID-19 is prevalent in the elderly. Old individuals are more likely to develop pneumonia and respiratory failure due to alveolar damage, suggesting that lung senescence may increase the susceptibility to SARS-CoV-2 infection and replication. Considering that human coronavirus (HCoVs; SARS-CoV-2 and SARS-CoV) require host cellular factors for infection and replication, we analyzed Genotype-Tissue Expression (GTEx) data to test whether lung aging is associated with transcriptional changes in human protein-coding genes that potentially interact with these viruses. We found decreased expression of the gene tribbles homolog 3 (TRIB3) during aging in male individuals, and its protein was predicted to interact with HCoVs nucleocapsid protein and RNA-dependent RNA polymerase. Using publicly available lung single-cell data, we found TRIB3 expressed mainly in alveolar epithelial cells that express SARS-CoV-2 receptor ACE2. Functional enrichment analysis of age-related genes, in common with SARS-CoV-induced perturbations, revealed genes associated with the mitotic cell cycle and surfactant metabolism. Given that TRIB3 was previously reported to decrease virus infection and replication, the decreased expression of TRIB3 in aged lungs may help explain why older male patients are related to more severe cases of the COVID-19. Thus, drugs that stimulate TRIB3 expression should be evaluated as a potential therapy for the disease.

Transcriptome of Two Canine Prostate Cancer Cells Treated With Toceranib Phosphate Reveals Distinct Antitumor Profiles Associated With the PDGFR Pathway.

Canine prostate cancer (PC) presents a poor antitumor response, usually late diagnosis and prognosis. Toceranib phosphate (TP) is a nonspecific inhibitor of receptor tyrosine kinases (RTKs), including vascular endothelial growth factor receptor (VEGFR), platelet-derived growth factor receptor (PDGFR), and c-KIT. This study aimed to evaluate VEGFR2, PDGFR-ß, and c-KIT protein expression in two established canine PC cell lines (PC1 and PC2) and the transcriptome profile of the cells after treatment with TP. Immunofluorescence (IF) analysis revealed VEGFR2 and PDGFR-ß protein expression and the absence of c-KIT protein expression in both cell lines. After TP treatment, only the viability of PC1 cells decreased in a dose-dependent manner. Transcriptome and enrichment analyses of treated PC1 cells revealed 181 upregulated genes, which were related to decreased angiogenesis and cell proliferation. In addition, we found upregulated PDGFR-A, PDGFR-ß, and PDGF-D expression in PC1 cells, and the upregulation of PDGFR-ß was also observed in treated PC1 cells by qPCR. PC2 cells had fewer protein-protein interactions (PPIs), with 18 upregulated and 22 downregulated genes; the upregulated genes were involved in the regulation of parallel pathways and mechanisms related to proliferation, which could be associated with the resistance observed after treatment. The canine PC1 cell line but not the PC2 cell line showed decreased viability after treatment with TP, although both cell lines expressed PDGFR and VEGFR receptors. Further studies could explain the mechanism of resistance in PC2 cells and provide a basis for personalized treatment for dogs with PC.

Comparison of microRNA Expression Profile in Chronic Myeloid Leukemia Patients Newly Diagnosed and Treated by Allogeneic Hematopoietic Stem Cell Transplantation.

Chronic myeloid leukemia (CML) results from a translocation between chromosomes 9 and 22, which generates the Philadelphia chromosome. This forms BCR/ABL1, an active tyrosine kinase protein that promotes cell growth and replication. Despite great progress in CML treatment in the form of tyrosine kinase inhibitors, allogeneic-hematopoietic stem cell transplantation (allo-HSCT) is currently used as an important treatment alternative for patients resistant to these inhibitors. Studies have shown that unregulated expression of microRNAs, which act as oncogenes or tumor suppressors, is associated with human cancers. This contributes to tumor formation and development by stimulating proliferation, angiogenesis, and invasion. Research has demonstrated the potential of microRNAs as biomarkers for cancer diagnosis, prognosis, and therapeutic targets. In the present study, we compared the circulating microRNA expression profiles of 14 newly diagnosed patients with chronic phase-CML and 14 Philadelphia chromosome-negative patients after allo-HSCT. For each patient, we tested 758 microRNAs by reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis. The global expression profile of microRNAs revealed 16 upregulated and 30 downregulated microRNAs. Target genes were analyzed, and key pathways were extracted and compared. Bioinformatics tools were used to analyze data. Among the downregulated miRNA target genes, some genes related to cell proliferation pathways were identified. These results reveal the comprehensive microRNA profile of CML patients and the main pathways related to the target genes of these miRNAs in cytogenetic remission after allo-HSCT. These results provide new resources for exploring stem cell transplantation-based CML treatment strategies.

A meta-analysis of microRNA networks regulated by melatonin in cancer: Portrait of potential candidates for breast cancer treatment.

Melatonin is a ubiquitous molecule with a broad spectrum of functions including widespread anti-cancer activities. Identifying how melatonin intervenes in complex molecular signaling at the gene level is essential to guide proper therapies. Using meta-analysis approach, herein we examined the role of melatonin in regulating the expression of 46 microRNAs (miRNAs) and their target genes in breast, oral, gastric, colorectal, and prostate cancers, and glioblastoma. The deregulated miRNA-associated target genes revealed their involvement in the regulation of cellular proliferation, differentiation, apoptosis, senescence, and autophagy. Melatonin changes the expression of miRNA-associated genes in breast, gastric, and oral cancers. These genes are associated with cellular senescence, the hedgehog signaling pathway, cell proliferation, p53 signaling, and the hippo signaling pathway. Conversely, colorectal and prostate cancers as well as glioblastoma and oral carcinoma present a clear pattern of less pronounced changes in the expression of miRNA-associated genes. Most notably, colorectal cancer displayed a unique molecular change in response to melatonin. Considering breast cancer network complexity, we compared the genes found during the meta-analysis with RNA-Seq data from breast cancer-bearing mice treated with melatonin. Mechanistically, melatonin upregulated genes associated with immune responses and apoptotic processes, whereas it downregulated genes involved in cellular aggressiveness/metastasis (eg, mitosis, telomerase activity, and angiogenesis). We further characterized the expression profile of our gene subsets with human breast cancer and found eight upregulated genes and 16 downregulated genes that were appositively correlated with melatonin. Our results pose a multi-dimension network of tumor-associated genes regulated by miRNAs potentially targeted by melatonin.

MicroRNA-mRNA Co-sequencing Identifies Transcriptional and Post-transcriptional Regulatory Networks Underlying Muscle Wasting in Cancer Cachexia.

Cancer cachexia is a metabolic syndrome with alterations in gene regulatory networks that consequently lead to skeletal muscle wasting. Integrating microRNAs-mRNAs omics profiles offers an opportunity to understand transcriptional and post-transcriptional regulatory networks underlying muscle wasting. Here, we used RNA sequencing to simultaneously integrate and explore microRNAs and mRNAs expression profiles in the tibialis anterior (TA) muscles of the Lewis Lung Carcinoma (LLC) model of cancer cachexia. We found 1,008 mRNAs and 18 microRNAs differentially expressed in cachectic mice compared with controls. Although our transcriptomic analysis demonstrated a high heterogeneity in mRNA profiles of cachectic mice, we identified a reduced number of differentially expressed genes that were uniformly regulated within cachectic muscles. This set of uniformly regulated genes is associated with the extracellular matrix (ECM), proteolysis, and inflammatory response. We also used transcriptomic data to perform enrichment analysis of transcriptional factor binding sites in promoter sequences, which revealed activation of the atrophy-related transcription factors NF-κB, Stat3, AP-1, and FoxO. Furthermore, the integration of mRNA and microRNA expression profiles identified post-transcriptional regulation by microRNAs of genes involved in ECM organization, cell migration, transcription factors binding, ion transport, and the FoxO signaling pathway. Our integrative analysis of microRNA-mRNA co-profiles comprehensively characterized regulatory relationships of molecular pathways and revealed microRNAs targeting ECM-associated genes in cancer cachexia.