Cell state dependent effects of Bmal1 on melanoma immunity and tumorigenicity

The circadian clock regulator Bmal1 modulates tumorigenesis, but its reported effects are inconsistent. Here, we show that Bmal1 has a context-dependent role in mouse melanoma tumor growth. Loss of Bmal1 in YUMM2.1 or B16-F10 melanoma cells eliminates clock function and diminishes hypoxic gene expression and tumorigenesis, which could be rescued by ectopic expression of HIF1α in YUMM2.1 cells. By contrast, over-expressed wild-type or a transcriptionally inactive mutant Bmal1 non-canonically sequester myosin heavy chain 9 (Myh9) to increase MRTF-SRF activity and AP-1 transcriptional signature, and shift YUMM2.1 cells from a Sox10high to a Sox9high immune resistant, mesenchymal cell state that is found in human melanomas. Our work describes a link between Bmal1, Myh9, mouse melanoma cell plasticity, and tumor immunity. This connection may underlie cancer therapeutic resistance and underpin the link between the circadian clock, MRTF-SRF and the cytoskeleton.

Common activities and predictive gene signature identified for genetic hypomorphs of TP53

Missense mutations that inactivate p53 occur commonly in cancer, and germline mutations in TP53 cause Li Fraumeni syndrome, which is associated with early-onset cancer. In addition, there are over two hundred germline missense variants of p53 that remain uncharacterized. In some cases, these germline variants have been shown to encode lesser-functioning, or hypomorphic, p53 protein, and these alleles are associated with increased cancer risk in humans and mouse models. However, most hypomorphic p53 variants remain un- or mis-classified in clinical genetics databases. There thus exists a significant need to better understand the behavior of p53 hypomorphs and to develop a functional assay that can distinguish hypomorphs from wild-type p53 or benign variants. We report the surprising finding that two different African-centric genetic hypomorphs of p53 that occur in distinct functional domains of the protein share common activities. Specifically, the Pro47Ser variant, located in the transactivation domain, and the Tyr107His variant, located in the DNA binding domain, both share increased propensity to misfold into a conformation specific for mutant, misfolded p53. Additionally, cells and tissues containing these hypomorphic variants show increased NF-κB activity. We identify a common gene expression signature from unstressed lymphocyte cell lines that is shared between multiple germline hypomorphic variants of TP53, and which successfully distinguishes wild-type p53 and a benign variant from lesser-functioning hypomorphic p53 variants. Our findings will allow us to better understand the contribution of p53 hypomorphs to disease risk and should help better inform cancer risk in the carriers of p53 variants.

Integrated microRNA–mRNA Expression Profiling Identifies Novel Targets and Networks Associated with Autism

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder, with mutations in hundreds of genes contributing to its risk. Herein, we studied lymphoblastoid cell lines (LCLs) from children diagnosed with autistic disorder (n = 10) and controls (n = 7) using RNA and miRNA sequencing profiles. The sequencing analysis identified 1700 genes and 102 miRNAs differentially expressed between the ASD and control LCLs (p ≤ 0.05). The top upregulated genes were GABRA4, AUTS2, and IL27, and the top upregulated miRNAs were hsa-miR-6813-3p, hsa-miR-221-5p, and hsa-miR-21-5p. The RT-qPCR analysis confirmed the sequencing results for randomly selected candidates: AUTS2, FMR1, PTEN, hsa-miR-15a-5p, hsa-miR-92a-3p, and hsa-miR-125b-5p. The functional enrichment analysis showed pathways involved in ASD control proliferation of neuronal cells, cell death of immune cells, epilepsy or neurodevelopmental disorders, WNT and PTEN signaling, apoptosis, and cancer. The integration of mRNA and miRNA sequencing profiles by miRWalk2.0 identified correlated changes in miRNAs and their targets’ expression. The integration analysis found significantly dysregulated miRNA–gene pairs in ASD. Overall, these findings suggest that mRNA and miRNA expression profiles in ASD are greatly altered in LCLs and reveal numerous miRNA–gene interactions that regulate critical pathways involved in the proliferation of neuronal cells, cell death of immune cells, and neuronal development.

A bovine miRNA, bta-miR-154c, withstands in vitro human digestion but does not affect cell viability of colorectal human cell lines after transfection

Colorectal cancer (CRC) is the third most frequent human cancer with over 1.3 million new cases globally. CRC is a complex disease caused by interactions between genetic and environmental factors; in particular, high consumption of red meat, including beef, is considered a risk factor for CRC initiation and progression. Recent data demonstrate that exogenous microRNAs (miRNAs) entering the body via ingestion could pose an effect on the consumer. In this study, we focused on bovine miRNAs that do not share a seed sequence with humans and mice. We identified bta-miR-154c, a bovine miRNA found in edible parts of beef and predicted via cross-species bioinformatic analysis to affect cancer-related pathways in human cells. When bovine tissue was subjected to cooking and a simulation of human digestion, bta-miR-154c was still detected after all procedures, albeit at reduced concentrations. However, lipofection of bta-miR-154c in three different colorectal human cell lines did not affect their viability as evaluated at various time points and concentrations. These data indicate that bta-miR-154c (a) may affect cancer-related pathways in human cells, (b) can withstand digestion and be detected after all stages of an in vitro digestion protocol, but (c) it does not appear to alter epithelial cell viability after entering human enterocytes, even at supraphysiological amounts. Further experiments will elucidate whether bta-miR-154c exerts a different functional effect on the human gut epithelium, which may cause it to contribute to CRC progression through its consumption.

Functional Annotation of MicroRNAs Using Existing Resources

MicroRNAs (miRNAs) are endogenous small noncoding RNAs that are involved in most biological signaling pathways, including the cell cycle, apoptosis, proliferation, immune response, metabolism as well as in biological processes including organ development and in human diseases like cancers. During the past two decades, high-throughput transcriptomic profiling using next generation sequencing and microarrays have been extensively utilized to identify differentially expressed miRNAs across different conditions and diseases. A natural extension of miRNA identification is to the process of functionally annotating known or predicted gene targets of those miRNAs and, by inference, revealing their potential influences on diverse biological pathways and functions. In this chapter, we provide a stepwise guideline on how to perform functional enrichment analyses on miRNAs of interest using publicly available resources such as miRWalk2.0.

Molecular Dysregulation in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) comprises a heterogeneous group of neurodevelopmental disorders with a strong heritable genetic component. At present, ASD is diagnosed solely by behavioral criteria. Advances in genomic analysis have contributed to numerous candidate genes for the risk of ASD, where rare mutations and s common variants contribute to its susceptibility. Moreover, studies show rare de novo variants, copy number variation and single nucleotide polymorphisms (SNPs) also impact neurodevelopment signaling. Exploration of rare and common variants involved in common dysregulated pathways can provide new diagnostic and therapeutic strategies for ASD. Contributions of current innovative molecular strategies to understand etiology of ASD will be explored which are focused on whole exome sequencing (WES), whole genome sequencing (WGS), microRNA, long non-coding RNAs and CRISPR/Cas9 models. Some promising areas of pharmacogenomic and endophenotype directed therapies as novel personalized treatment and prevention will be discussed.

Plasma Markers of Disrupted Gut Permeability in Severe COVID-19 Patients

A disruption of the crosstalk between the gut and the lung has been implicated as a driver of severity during respiratory-related diseases. Lung injury causes systemic inflammation, which disrupts gut barrier integrity, increasing the permeability to gut microbes and their products. This exacerbates inflammation, resulting in positive feedback. We aimed to test whether severe Coronavirus disease 2019 (COVID-19) is associated with markers of disrupted gut permeability. We applied a multi-omic systems biology approach to analyze plasma samples from COVID-19 patients with varying disease severity and SARS-CoV-2 negative controls. We investigated the potential links between plasma markers of gut barrier integrity, microbial translocation, systemic inflammation, metabolome, lipidome, and glycome, and COVID-19 severity. We found that severe COVID-19 is associated with high levels of markers of tight junction permeability and translocation of bacterial and fungal products into the blood. These markers of disrupted intestinal barrier integrity and microbial translocation correlate strongly with higher levels of markers of systemic inflammation and immune activation, lower levels of markers of intestinal function, disrupted plasma metabolome and glycome, and higher mortality rate. Our study highlights an underappreciated factor with significant clinical implications, disruption in gut functions, as a potential force that may contribute to COVID-19 severity.

Analysis of classical neutrophils and polymorphonuclear myeloid-derived suppressor cells in cancer patients and tumor-bearing mice

In this study, using single-cell RNA-seq, cell mass spectrometry, flow cytometry, and functional analysis, we characterized the heterogeneity of polymorphonuclear neutrophils (PMNs) in cancer. We describe three populations of PMNs in tumor-bearing mice: classical PMNs, polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs), and activated PMN-MDSCs with potent immune suppressive activity. In spleens of mice, PMN-MDSCs gradually replaced PMNs during tumor progression. Activated PMN-MDSCs were found only in tumors, where they were present at the very early stages of the disease. These populations of PMNs in mice could be separated based on the expression of CD14. In peripheral blood of cancer patients, we identified two distinct populations of PMNs with characteristics of classical PMNs and PMN-MDSCs. The gene signature of tumor PMN-MDSCs was similar to that in mouse activated PMN-MDSCs and was closely associated with negative clinical outcome in cancer patients. Thus, we provide evidence that PMN-MDSCs are a distinct population of PMNs with unique features and potential for selective targeting opportunities.

Tumor-infiltrating mast cells are associated with resistance to anti-PD-1 therapy

Anti-PD-1 therapy is used as a front-line treatment for many cancers, but mechanistic insight into this therapy resistance is still lacking. Here we generate a humanized (Hu)-mouse melanoma model by injecting fetal liver-derived CD34⁺ cells and implanting autologous thymus in immune-deficient NOD-scid IL2Rγ^null (NSG) mice. Reconstituted Hu-mice are challenged with HLA-matched melanomas and treated with anti-PD-1, which results in restricted tumor growth but not complete regression. Tumor RNA-seq, multiplexed imaging and immunohistology staining show high expression of chemokines, as well as recruitment of FOXP3⁺ Treg and mast cells, in selective tumor regions. Reduced HLA-class I expression and CD8⁺/Granz B⁺ T cells homeostasis are observed in tumor regions where FOXP3⁺ Treg and mast cells co-localize, with such features associated with resistance to anti-PD-1 treatment. Combining anti-PD-1 with sunitinib or imatinib results in the depletion of mast cells and complete regression of tumors. Our results thus implicate mast cell depletion for improving the efficacy of anti-PD-1 therapy.

Immune checkpoint therapies (ICT) are promising for treating various cancers, but response rates vary. Here the authors show, in mouse models, that tumor-infiltrating mast cells colocalize with regulatory T cells, coincide with local reduction of MHC-I and CD8 T cells, and is associated with resistance to ICT, which can be reversed by c-kit inhibitor treatment.

Sialyl-LewisX Glycoantigen Is Enriched on Cells with Persistent HIV Transcription during Therapy

A comprehensive understanding of the phenotype of persistent HIV-infected cells, transcriptionally active and/or transcriptionally inactive, is imperative for developing a cure. The relevance of cell-surface glycosylation to HIV persistence has never been explored. We characterize the relationship between cell-surface glycomic signatures and persistent HIV transcription in vivo. We find that the cell surface of CD4+ T cells actively transcribing HIV, despite suppressive therapy, harbors high levels of fucosylated carbohydrate ligands, including the cell extravasation mediator Sialyl-LewisX (SLeX), compared with HIV-infected transcriptionally inactive cells. These high levels of SLeX are induced by HIV transcription in vitro and are maintained after therapy in vivo. Cells with high-SLeX are enriched with markers associated with HIV susceptibility, signaling pathways that drive HIV transcription, and pathways involved in leukocyte extravasation. We describe a glycomic feature of HIV-infected transcriptionally active cells that not only differentiates them from their transcriptionally inactive counterparts but also may affect their trafficking abilities.

Granzyme B and miR-378a Interaction in Acetaminophen Toxicity in Children

BACKGROUND AND AIM

Hepatic phase I drug-metabolizing enzymes CYP2E1, CYP1A2 and CYP3A4 catalyze the biotransformation of Acetaminophen (APAP) and are important in the mediation of toxicity. The potential role of other hepatic and non-hepatic Phase I enzymes in APAP toxicity has not been established.

METHODS

PCR array containing 84 genes involved in phase I drug metabolism was examined in subgroups of hospitalized children for APAP overdose, categorized as no toxicity (ALT ≤ 45 IU/L, n=5) and moderate toxicity (ALT ≥ 500 IU/L, n=5).

RESULTS

Significant downregulation was observed for ALDH6A1, CYP4F12 and GZMB in the no toxicity subgroup and ALDH1A1, CYP27A1 and GZMB in the moderate toxicity subgroup. qRTPCR confirmed significant downregulation for ALDH1A1, CYP4F12, and GZMB. In-silico analysis identified GZMB 3'UTR to be a target of miR-378a-5p. Overexpression of miR-378a-5p reduced the luciferase activity of GZMB 3'UTR reporter plasmid reportedly by 50%. NK-92 cells transfected with the miR-378a-5p mimic extended the effect of APAP on GZMB protein expression compared to mimic controls. In addition, miR-378a-5p was significantly upregulated in blood samples of children with APAP overdose undergoing NAC treatment.

CONCLUSION

Overall, our study suggests the presence of a novel signaling pathway, whereby miR- 378a-5p inhibits GZMB expression in children with APAP overdose.

TLR9 agonist MGN1703 enhances B cell differentiation and function in lymph nodes

Background

TLR9 agonists are being developed as immunotherapy against malignancies and infections. TLR9 is primarily expressed in B cells and plasmacytoid dendritic cells (pDCs). TLR9 signalling may be critically important for B cell activity in lymph nodes but little is known about the in vivo impact of TLR9 agonism on human lymph node B cells. As a pre-defined sub-study within our clinical trial investigating TLR9 agonist MGN1703 (lefitolimod) treatment in the context of developing HIV cure strategies (NCT02443935), we assessed TLR9 agonist-mediated effects in lymph nodes.

Methods

Participants received MGN1703 for 24 weeks concurrent with antiretroviral therapy. Seven participants completed the sub-study including lymph node resection at baseline and after 24 weeks of treatment. A variety of tissue-based immunologic and virologic parameters were assessed.

Findings

MGN1703 dosing increased B cell differentiation; activated pDCs, NK cells, and T cells; and induced a robust interferon response in lymph nodes. Expression of Activation-Induced cytidine Deaminase, an essential regulator of B cell diversification and somatic hypermutation, was highly elevated. During MGN1703 treatment IgG production increased and antibody glycosylation patterns were changed.

Interpretation

Our data present novel evidence that the TLR9 agonist MGN1703 modulates human lymph node B cells in vivo. These findings warrant further considerations in the development of TLR9 agonists as immunotherapy against cancers and infectious diseases.

Fund

This work was supported by Aarhus University Research Foundation, the Danish Council for Independent Research and the NovoNordisk Foundation. Mologen AG provided study drug free of charge.

Correction to: Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

In the original publication of this article [1], the Figure 1 and Figure 2 were wrong. The Figure 1 "Heat map showing the quantity of DNA repair genes, from red to blue in ascending order, per species' genome (numbers at the top of the figure represent the species code that is found in Table 1). Each DNA repair gene pathway was analyzed separately in rows. Radiated species' genomes are richer in DNA repair genes. Analytical data can be found in Additional file 2: Table S2. M mammals, B&R birds and reptiles, BF bony fishes" should be the picture of Figure 2. The figure 2 "Linear regression analysis. The number of DNA repair genes is linearly related to genome size and protein number. As a negative control, we show that genome size is not linearly related with protein number" should be the picture of figure 1.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

miR-16 and miR-103 impact 5-HT4 receptor signalling and correlate with symptom profile in irritable bowel syndrome

Irritable bowel syndrome (IBS) is a gut-brain disorder involving alterations in intestinal sensitivity and motility. Serotonin 5-HT4 receptors are promising candidates in IBS pathophysiology since they regulate gut motor function and stool consistency, and targeted 5-HT4R selective drug intervention has been proven beneficial in subgroups of patients. We identified a single nucleotide polymorphism (SNP) (rs201253747) c.*61 T > C within the 5-HT4 receptor gene HTR4 to be predominantly present in diarrhoea-IBS patients (IBS-D). It affects a binding site for the miR-16 family and miR-103/miR-107 within the isoforms HTR4b/i and putatively impairs HTR4 expression. Subsequent miRNA-profiling revealed downregulation of miR-16 and miR-103 in the jejunum of IBS-D patients correlating with symptoms. In vitro assays confirmed expression regulation via three 3'UTR binding sites. The novel isoform HTR4b_2 lacking two of the three miRNA binding sites escapes miR-16/103/107 regulation in SNP carriers. We provide the first evidence that HTR4 expression is fine-tuned by miRNAs, and that this regulation is impaired either by the SNP c.*61 T > C or by diminished levels of miR-16 and miR-103 suggesting that HTR4 might be involved in the development of IBS-D.

MicroRNA regulation of CYP 1A2, CYP3A4 and CYP2E1 expression in acetaminophen toxicity

MicroRNAs (miRNAs) that regulate the cytochrome P-450 isoforms involved in acetaminophen (APAP) toxicity were examined in HepaRG cells treated with APAP (20 mM). In-vitro studies found that APAP protein adducts were increased at 1 h, followed by ALT increases at 12 and 24 h. CYP1A2, CYP3A4 and CYP2E1 mRNA levels were decreased, while miRNAs were increased for miR-122-5p, miR-378a-5p, miR-27b-3p at 6 h and miR-125b-5p at 12 h and miR-27b-3p at 24 h. Putative miRNA binding sites on the 3′UTRs of the CYPs were identified in-silico. Overexpression of miR-122-5p and miR-378a-5p in cells suppressed protein expression of CYP1A2, CYP3A4 and CYP2E1. Luciferase reporter assays confirmed the interaction between miR-122 and the 3′UTR of the CYP1A2 and CYP3A4. Thus, the in-vitro experiments showed that miR-122-5p and miR-378a-5p upregulation were associated with translational repression of CYPs. Serum samples of children with APAP overdose had significant elevation of miR-122-5p, miR-378a-5p, miR-125b-5p and miR-27b-3p, compared to healthy controls and receiver operator curves of the miRNAs had AUCs of 91 to 100%. Collectively, the data suggest that miRNA elevations in APAP toxicity represent a regulatory response to modify CYP1A2, CYP3A4 and CYP2E1 translation due to cellular stress and injury.

Toxicity study of ochratoxin A using HEK293 and HepG2 cell lines based on microRNA profiling

Ochratoxin A (OTA) induced DNA damage, cytotoxicity, and apoptosis in mammalian cell lines. Micro RNAs (miRNAs) are involved in physiological and developmental processes and contribute to cancer development and progression. In our study, high-throughput miRNA profiling and Kyoto Encyclopedia of Genes and Genomes analysis were applied to comparatively study the toxicity of OTA in HEK293 cells and HepG2 cells treated with 25 μM OTA for 24 h. In these two cells, the same changing miRNAs were mostly related to signal transduction pathways, whereas the different changing miRNAs were mostly related to human cancer pathways. DGCR8, Dicer1, and Drosha were significantly suppressed in HEK293 cells, indicating an impairment of miRNA biogenesis. The damage seemed more extensive in HEK293 cells. Cell models and in vivo models were also compared. Many miRNAs in vitro were markedly different from those in vivo; however, OTA toxicity was observed both in vitro and in vivo. The classification of deregulated pathways is similar. The biogenesis of miRNA was impaired in both lines. In conclusion, deregulated miRNAs in vitro are mostly related to human cancer and signal transduction pathways. The deregulated pathways in vivo are similar to those in vitro.

A Comprehensive Review on the Genetic Regulation of Cisplatin-induced Nephrotoxicity

Cisplatin (CDDP) is a well-known antineoplastic drug which has been extensively utilized over the last decades in the treatment of numerous kinds of tumors. However, CDDP induces a wide range of toxicities in a dose-dependent manner, among which nephrotoxicity is of particular importance. Still, the mechanism of CDDP-induced renal damage is not completely understood; moreover, the knowledge about the role of microRNAs (miRNAs) in the nephrotoxic response is still unknown. miRNAs are known to interact with the representative members of a diverse range of regulatory pathways (including postnatal development, proliferation, inflammation and fibrosis) and pathological conditions, including kidney diseases: polycystic kidney diseases (PKDs), diabetic nephropathy (DN), kidney cancer, and drug-induced kidney injury. In this review, we shed light on the following important aspects: (i) information on genes/proteins and their interactions with previously known pathways engaged with CDDP-induced nephrotoxicity, (ii) information on newly discovered biomarkers, especially, miRNAs for detecting CDDP-induced nephrotoxicity and (iii) information to improve our understanding on CDDP. This information will not only help the researchers belonging to nephrotoxicity field, but also supply an indisputable help for oncologists to better understand and manage the side effects induced by CDDP during cancer treatment. Moreover, we provide up-to-date information about different in vivo and in vitro models that have been utilized over the last decades to study CDDP-induced renal injury. Taken together, this review offers a comprehensive network on genes, miRNAs, pathways and animal models which will serve as a useful resource to understand the molecular mechanism of CDDP-induced nephrotoxicity.

Alteration of mRNA and microRNA expression profiles in rat muscular type vasculature in early postnatal development

The vascular system is characterized by a high degree of plasticity. In particular, functional and structural remodeling of the arterial system takes place during early postnatal development. However, the mechanisms providing such alterations in the rapidly growing organisms are poorly understood, especially for the peripheral vasculature. To explore this, we performed mRNA- and miRNA microarray analysis on muscular type saphenous arteries of young (10-12 days) and adult (2-3 months) rats. Thirty-eight significant pathways (such as oxidative phosphorylation, MAPK signaling, metabolism, cell cycle, DNA replication and focal adhesion) were obtained on differentially regulated genes during postnatal development. Many differentially regulated genes were determined as target- and miRNA-hubs. We also found 92 miRNAs differentially expressed in arteries of young and adult rats. Several significantly regulated pathways were found on these regulated miRNAs. Interestingly, these biological cascades also contain those significantly enriched pathways that were previously identified based on the differently expressed genes. Our data indicate that the expression of many genes involved in the regulation of pathways that are relevant for different functions in arteries may be under the control of miRNAs and these miRNAs regulate the functional, and structural remodeling occurring in the vascular system during early postnatal development.

MicroRNA profiling of rats with ochratoxin A nephrotoxicity

Background

Nephrotoxicity is the most prominent one among the various toxicities of ochratoxin A (OTA). MicroRNAs (miRNAs) are small non-coding RNAs that have an impact on a wide range of biological processes by regulating gene expression at post-transcriptional level or protein systhesis level. The objective of this study is to analyze miRNA profiling in the kidneys of rats gavaged with OTA.

Results

To profile miRNAs in the kidneys of rats with OTA nephrotoxicity, high-throughput sequencing and bioinformatics approaches were applied to analyze the miRNAs in the kidney of rats following OTA treatment. A total of 409 known miRNAs and 8 novel miRNAs were identified in the kidney and the levels of the novel miRNAs were varied in response to different doses of OTA. Expression of miR-129, miR-130a, miR-130b, miR-141, miR-218b and miR-3588 were uniquely suppressed in mid dose but then elevated in high dose, with opposite expression to their target genes. The expression pattern was closely related with the “MAPK signaling pathway”. Dicer1 and Drosha were significantly suppressed, indicating an impairment of miRNA biogenesis in response to OTA.

Conclusions

The abrogation of miRNA maturation process suggests a new target of OTA toxicity. Moreover, the identification of the differentially expressed miRNAs provides us a molecular insight into the nephrtoxicity of OTA.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-333) contains supplementary material, which is available to authorized users.

miRWalk database for miRNA-target interactions

miRWalk (http://mirwalk.uni-hd.de/) is a publicly available comprehensive resource, hosting the predicted as well as the experimentally validated microRNA (miRNA)-target interaction pairs. This database allows obtaining the possible miRNA-binding site predictions within the complete sequence of all known genes of three genomes (human, mouse, and rat). Moreover, it also integrates many novel features such as a comparative platform of miRNA-binding sites resulting from ten different prediction datasets, a holistic view of genetic networks of miRNA-gene pathway, and miRNA-gene-Online Mendelian Inheritance in Man disorder interactions, and unique experimentally validated information (e.g., cell lines, diseases, miRNA processing proteins). In this chapter, we describe a schematic workflow on how one can access the stored information from miRWalk and subsequently summarize its applications.

A schematic workflow for collecting information about the interaction between copy number variants and microRNAs using existing resources

MicroRNAs (miRNAs) and copy number variations (CNVs) are two extensively studied genomic components in the field of modern biology-as they have been found to be associated with many disorders such as cancer, Alzheimer, pancreatitis, HIV susceptibility, beta-thalassemia, and glomerulonephritis. Several studies suggested that an alteration in CNV-miRNA interaction could result in some human diseases such as cancer. Therefore, the possible miRNA-binding site information within the CNV genes opens new avenues in understanding such disorders. In this chapter, we present a schematic approach for collecting the information on CNV-miRNA interactions using miRWalk and TargetScan databases.

Evidence for activation of the unfolded protein response in collagen IV nephropathies

Thin-basement-membrane nephropathy (TBMN) and Alport syndrome (AS) are progressive collagen IV nephropathies caused by mutations in COL4A3/A4/A5 genes. These nephropathies invariably present with microscopic hematuria and frequently progress to proteinuria and CKD or ESRD during long-term follow-up. Nonetheless, the exact molecular mechanisms by which these mutations exert their deleterious effects on the glomerulus remain elusive. We hypothesized that defective trafficking of the COL4A3 chain causes a strong intracellular effect on the cell responsible for COL4A3 expression, the podocyte. To this end, we overexpressed normal and mutant COL4A3 chains (G1334E mutation) in human undifferentiated podocytes and tested their effects in various intracellular pathways using a microarray approach. COL4A3 overexpression in the podocyte caused chain retention in the endoplasmic reticulum (ER) that was associated with activation of unfolded protein response (UPR)-related markers of ER stress. Notably, the overexpression of normal or mutant COL4A3 chains differentially activated the UPR pathway. Similar results were observed in a novel knockin mouse carrying the Col4a3-G1332E mutation, which produced a phenotype consistent with AS, and in biopsy specimens from patients with TBMN carrying a heterozygous COL4A3-G1334E mutation. These results suggest that ER stress arising from defective localization of collagen IV chains in human podocytes contributes to the pathogenesis of TBMN and AS through activation of the UPR, a finding that may pave the way for novel therapeutic interventions for a variety of collagenopathies.

In-Silico Algorithms for the Screening of Possible microRNA Binding Sites and Their Interactions

MicroRNAs (miRNAs) comprise a recently discovered class of small, non-coding RNA molecules of 21-25 nucleotides in length that regulate the gene expression by base-pairing with the transcripts of their targets i.e. protein-coding genes, leading to down-regulation or repression of the target genes. However, target gene activation has also been described. miRNAs are involved in diverse regulatory pathways, including control of developmental timing, apoptosis, cell proliferation, cell differentiation, modulation of immune response to macrophages, and organ development and are associated with many diseases, such as cancer. Computational prediction of miRNA targets is much more challenging in animals than in plants, because animal miRNAs often perform imperfect base-pairing with their target sites, unlike plant miRNAs which almost always bind their targets with near perfect complementarity. In the past years, a large number of target prediction programs and databases on experimentally validated information have been developed for animal miRNAs to fulfil the need of experimental scientists conducting miRNA research. In this review we first succinctly describe the prediction criteria (rules or principles) adapted by prediction algorithms to generate possible miRNA binding site interactions and introduce most relevant algorithms, and databases. We then summarize their applications with the help of some previously published studies. We further provide experimentally validated functional binding sites outside 3’-UTR region of target mRNAs and the resources which offer such predictions. Finally, the issue of experimental validation of miRNA binding sites will be briefly discussed.

A miR-1207-5p binding site polymorphism abolishes regulation of HBEGF and is associated with disease severity in CFHR5 nephropathy

Heparin binding epidermal growth factor (HBEGF) is expressed in podocytes and was shown to play a role in glomerular physiology. MicroRNA binding sites on the 3′UTR of HBEGF were predicted using miRWalk algorithm and followed by DNA sequencing in 103 patients diagnosed with mild or severe glomerulopathy. A single nucleotide polymorphism, miRSNP C1936T (rs13385), was identified at the 3′UTR of HBEGF that corresponds to the second base of the hsa-miR-1207-5p seed region. When AB8/13 undifferentiated podocytes were transfected with miRNA mimics of hsa-miR-1207-5p, the HBEGF protein levels were reduced by about 50%. A DNA fragment containing the miRSNP allele-1936C was cloned into the pMIR-Report Luciferase vector and co-transfected with miRNA mimics of hsa-miR-1207-5p into AB8/13 podocytes. In agreement with western blot data, this resulted in reduced luciferase expression demonstrating the ability of hsa-miR-1207-5p to directly regulate HBEGF expression. On the contrary, in the presence of the miRSNP 1936T allele, this regulation was abolished. Collectively, these results demonstrate that variant 1936T of this miRSNP prevents hsa-miR-1207-5p from down-regulating HBEGF in podocytes. We hypothesized that this variant has a functional role as a genetic modifier. To this end, we showed that in a cohort of 78 patients diagnosed with CFHR5 nephropathy (also known as C3-glomerulopathy), inheritance of miRSNP 1936T allele was significantly increased in the group demonstrating progression to chronic renal failure on long follow-up. No similar association was detected in a cohort of patients with thin basement membrane nephropathy. This is the first report associating a miRSNP as genetic modifier to a monogenic renal disorder.

miRWalk--database: prediction of possible miRNA binding sites by "walking" the genes of three genomes

MicroRNAs are small, non-coding RNA molecules that can complementarily bind to the mRNA 3'-UTR region to regulate the gene expression by transcriptional repression or induction of mRNA degradation. Increasing evidence suggests a new mechanism by which miRNAs may regulate target gene expression by binding in promoter and amino acid coding regions. Most of the existing databases on miRNAs are restricted to mRNA 3'-UTR region. To address this issue, we present miRWalk, a comprehensive database on miRNAs, which hosts predicted as well as validated miRNA binding sites, information on all known genes of human, mouse and rat. All mRNAs, mitochondrial genes and 10 kb upstream flanking regions of all known genes of human, mouse and rat were analyzed by using a newly developed algorithm named 'miRWalk' as well as with eight already established programs for putative miRNA binding sites. An automated and extensive text-mining search was performed on PubMed database to extract validated information on miRNAs. Combined information was put into a MySQL database. miRWalk presents predicted and validated information on miRNA-target interaction. Such a resource enables researchers to validate new targets of miRNA not only on 3'-UTR, but also on the other regions of all known genes. The 'Validated Target module' is updated every month and the 'Predicted Target module' is updated every 6 months. miRWalk is freely available at http://mirwalk.uni-hd.de/.

Alternative splicing and nonsense-mediated RNA decay contribute to the regulation of SHOX expression

The human SHOX gene is composed of seven exons and encodes a paired-related homeodomain transcription factor. SHOX mutations or deletions have been associated with different short stature syndromes implying a role in growth and bone formation. During development, SHOX is expressed in a highly specific spatiotemporal expression pattern, the underlying regulatory mechanisms of which remain largely unknown. We have analysed SHOX expression in diverse embryonic, fetal and adult human tissues and detected expression in many tissues that were not known to express SHOX before, e.g. distinct brain regions. By using RT-PCR and comparing the results with RNA-Seq data, we have identified four novel exons (exon 2a, 7-1, 7-2 and 7-3) contributing to different SHOX isoforms, and also established an expression profile for the emerging new SHOX isoforms. Interestingly, we found the exon 7 variants to be exclusively expressed in fetal neural tissues, which could argue for a specific role of these variants during brain development. A bioinformatical analysis of the three novel 3′UTR exons yielded insights into the putative role of the different 3′UTRs as targets for miRNA binding. Functional analysis revealed that inclusion of exon 2a leads to nonsense-mediated RNA decay altering SHOX expression in a tissue and time specific manner. In conclusion, SHOX expression is regulated by different mechanisms and alternative splicing coupled with nonsense-mediated RNA decay constitutes a further component that can be used to fine tune the SHOX expression level.