Search

Non-coding variation in complex human disease has been well established by genome-wide association studies, and it is thought to involve regulatory elements, such as enhancers, whose variation affects the expression of the gene responsible for the disease. The regulatory elements often lie far from the gene they regulate, or within introns of genes differing from the regulated gene, making it difficult to identify the gene whose function is affected by a given enhancer variation.Enhancers are connected to their target gene promoters via long-range physical interactions (loops). In our study, we re-mapped, onto the human genome, more than 10,000 enhancers connected to promoters via long-range interactions, that we had previously identified in mouse brain-derived neural stem cells by RNApolII-ChIA-PET analysis, coupled to ChIP-seq mapping of DNA/chromatin regions carrying epigenetic enhancer marks. These interactions are thought to be functionally relevant. We discovered, in the human genome, thousands of DNA regions syntenic with the interacting mouse DNA regions (enhancers+connected promoters). We further annotated these human regions regarding their overlap with sequence variants (single nucleotide polymorphisms, SNPs; copy number variants, CNVs), that were previously associated with neurodevelopmental disease in human. We document various cases in which the genetic variant, associated in human to neurodevelopmental disease, affects an enhancer involved in long-range interactions: SNPs, previously identified by genome-wide association studies to be associated with schizophrenia, bipolar disorder, and intelligence, are located within our human syntenic enhancers, and alter transcription factors recognition sites. Similarly, CNVs associated to autism spectrum disease and other neurodevelopmental disorders overlap with our human syntenic enhancers. Some of these enhancers are connected (in mouse) to homologs of genes already associated to the human disease, strengthening the hypothesis that the gene is indeed involved in the disease. Other enhancers are connected to genes not previously associated with the disease, pointing to their possible pathogenetic involvement.Our observations provide a resource for further exploration of neural disease, in parallel with the now widespread genome-wide identification of DNA variants in patients with neural disease.

Circulating miRNAs (c-miRNAs) are promising biomarkers for HF diagnosis and prognosis. There are no studies on HF pediatric patients undergoing VAD-implantation. Aims of this study were: to examine the c-miRNAs profile in HF children; to evaluate the effects of VAD on c-miRNAs levels; to in vitro validate putative c-miRNA targets. c-miRNA profile was determined in serum of HF children by NGS before and one month after VAD-implant. The c-miRNA differentially expressed were analyzed by real time-PCR, before and at 4 hrs,1,3,7,14,30 days after VAD-implant. A miRNA mimic transfection study in HepG2 cells was performed to validate putative miRNA targets selected through miRWalk database. Thirteen c-miRNAs were modified at 30 days after VAD-implant compared to pre-VAD at NSG, and, among them, six c-miRNAs were confirmed by Real-TimePCR. Putative targets of the validated c-miRNAs are involved in the hemostatic process. The in vitro study confirmed a down-regulatory effect of hsa-miR-409-3p towards coagulation factor 7 (F7) and F2. Of note, all patients had thrombotic events requiring pump change. In conclusion, in HF children, the level of six c-miRNAs involved in the regulation of hemostatic events changed after 30 days of VAD-treatment. In particular, the lowering of c-miR-409-3p regulating both F7 and F2 could reflect a pro-thrombotic state after VAD-implant.

DNA sequence variants (single nucleotide polymorphisms or variants, SNPs/SNVs; copy number variants, CNVs) associated to neurodevelopmental disorders (NDD) and traits often map on putative transcriptional regulatory elements, including, in particular, enhancers. However, the genes controlled by these enhancers remain poorly defined. Traditionally, the activity of a given enhancer, and the effect of its possible alteration associated to the sequence variants, has been thought to influence the nearest gene promoter. However, the obtainment of genome-wide long-range interaction maps in neural cells chromatin challenged this view, showing that a given enhancer is very frequently not connected to the nearest promoter, but to a more distant one, skipping genes in between. In this Perspective, we review some recent papers, who generated long-range interaction maps (by HiC, RNApolII ChIA-PET, Capture-HiC, or PLACseq), and overlapped the identified long-range interacting DNA segments with DNA sequence variants associated to NDD (such as schizophrenia, bipolar disorder and autism) and traits (intelligence). This strategy allowed to attribute the function of enhancers, hosting the NDD-related sequence variants, to a connected gene promoter lying far away on the linear chromosome map. Some of these enhancer-connected genes had indeed been already identified as contributive to the diseases, by the identification of mutations within the gene’s protein-coding regions (exons), validating the approach. Significantly, however, the connected genes also include many genes that were not previously found mutated in their exons, pointing to novel candidate contributors to NDD and traits. Thus, long-range interaction maps, in combination with DNA variants detected in association with NDD, can be used as “pointers” to identify novel candidate disease-relevant genes. Functional manipulation of the long-range interaction network involving enhancers and promoters by CRISPR-Cas9-based approaches is beginning to probe for the functional significance of the identified interactions, and the enhancers and the genes involved, improving our understanding of neural development and its pathology.

To dissect the TBX5 regulatory circuit, we focused on microRNAs (miRNAs) that collectively contribute to make TBX5 a pivotal cardiac regulator. We profiled miRNAs in hearts isolated from wild-type, CRE, Tbx5lox/+and Tbx5del/+ mice using a Next Generation Sequencing (NGS) approach. TBX5 deficiency in cardiomyocytes increased the expression of the miR-183 cluster family that is controlled by Kruppel-like factor 4, a transcription factor repressed by TBX5. MiR-182-5p, the most highly expressed miRNA of this family, was functionally analyzed in zebrafish. Transient overexpression of miR-182-5p affected heart morphology, calcium handling and the onset of arrhythmias as detected by ECG tracings. Accordingly, several calcium channel proteins identified as putative miR-182-5p targets were downregulated in miR-182-5p overexpressing hearts. In stable zebrafish transgenic lines, we demonstrated that selective miRNA-182-5p upregulation contributes to arrhythmias. Moreover, cardiac-specific down-regulation of miR-182-5p rescued cardiac defects in a zebrafish model of Holt–Oram syndrome. In conclusion, miR-182-5p exerts an evolutionarily conserved role as a TBX5 effector in the onset of cardiac propensity for arrhythmia, and constitutes a relevant target for mediating the relationship between TBX5, arrhythmia and heart development. To dissect the TBX5 regulatory circuit, we focused on microRNAs (miRNAs) that collectively contribute to make TBX5 a pivotal cardiac regulator. We profiled miRNAs in hearts isolated from wild-type, CRE, Tbx5lox/+and Tbx5del/+ mice using a Next Generation Sequencing (NGS) approach. TBX5 deficiency in cardiomyocytes increased the expression of the miR-183 cluster family that is controlled by Kruppel-like factor 4, a transcription factor repressed by TBX5. MiR-182-5p, the most highly expressed miRNA of this family, was functionally analyzed in zebrafish. Transient overexpression of miR-182-5p affected heart morphology, calcium handling and the onset of arrhythmias as detected by ECG tracings. Accordingly, several calcium channel proteins identified as putative miR-182-5p targets were downregulated in miR-182-5p overexpressing hearts. In stable zebrafish transgenic lines, we demonstrated that selective miRNA-182-5p upregulation contributes to arrhythmias. Moreover, cardiac-specific down-regulation of miR-182-5p rescued cardiac defects in a zebrafish model of Holt–Oram syndrome. In conclusion, miR-182-5p exerts an evolutionarily conserved role as a TBX5 effector in the onset of cardiac propensity for arrhythmia, and constitutes a relevant target for mediating the relationship between TBX5, arrhythmia and heart development.

Five-year overall survival of stage III colorectal cancer (CRC) patients treated with standard adjuvant chemotherapy (ACHT) is highly variable. Genomic biomarkers and/or transcriptomic profiles identified lack of adequate validation. Aim of our study was to identify and validate molecular biomarkers predictive of ACHT response in stage III CRC patients by a transcriptomic approach. From a series of CRC patients who received ACHT, two stage III extreme cohorts (unfavorable vs. favorable prognosis) were selected. RNA-sequencing was performed from fresh frozen explants. Tumors were characterized for somatic mutations. Validation was performed in stage III CRC patients extracted from two GEO datasets. According to disease-free survival (DFS), 108 differentially expressed genes (104/4 up/downregulated in the unfavorable prognosis group) were identified. Among 104 upregulated genes, 42 belonged to olfactory signaling pathways, 62 were classified as pseudogenes (n = 17), uncharacterized noncoding RNA (n = 10), immune response genes (n = 4), microRNA (n = 1), cancer-related genes (n = 14) and cancer-unrelated genes (n = 16). Three out of four down-regulated genes were cancer-related. Mutational status (i.e., RAS, BRAF, PIK3CA) did not differ among the cohorts. In the validation cohort, multivariate analysis showed high PNN and KCNQ1OT1 expression predictive of shorter DFS in ACHT treated patients (p = 0.018 and p = 0.014, respectively); no difference was observed in untreated patients. This is the first study that identifies by a transcriptomic approach and validates PNN and KCNQ1OT1 as molecular biomarkers predictive of chemotherapy response in stage III CRC patients. After a further validation in an independent cohort, PNN and KCNQ1OT1 evaluation could be proposed to prospectively identify stage III CRC patients benefiting from ACHT.

Here we present miR-CATCHv2.0, an implemented experimental method that allows the identification of the microRNA species directly bound to an RNA of interest. After cross-linking of microRNA::RNA::Ago2 complexes using formaldehyde, the RNA is fragmented using sonication and then subjected to affinity purification using two sets of biotinylated tiling probes (ODD and EVEN). Finally, enriched microRNA species are retrieved by means of small RNA sequencing coupled with an ad hoc analytical workflow. In BRAFV600E mutant A375 melanoma cells, miR-CATCHv2.0 allowed us to identify 20 microRNAs that target X1, the most abundant isoform of BRAF mRNA. These microRNAs fall into different functional classes, according to the effect that they exert (decrease/increase in BRAFV600E mRNA and protein levels) and to the mechanism they use to achieve it (destabilization/stabilization of X1 mRNA or decrease/increase in its translation). microRNA-induced variations in BRAFV600E protein levels are most of the times coupled to consistent variations in pMEK levels, in melanoma cell proliferation in vitro and in sensitivity to the BRAF inhibitor vemurafenib in a xenograft model in zebrafish. However, microRNAs exist that uncouple the degree of activation of the ERK pathway from the levels of BRAFV600E protein. Our study proposes miR-CATCHv2.0 as an effective tool for the identification of direct microRNA-target interactions and, by using such a tool, unveils the complexity of the post-transcriptional regulation to which BRAFV600E and the ERK pathway are subjected in melanoma cells.

Tailored treatments for metastatic colorectal cancer (mCRC) have not yet completely evolved due to the variety in response to drugs. Therefore, artificial intelligence has been recently used to develop prognostic and predictive models of treatment response (either activity/efficacy or toxicity) to aid in clinical decision making. In this systematic review, we have examined the ability of learning methods to predict response to chemotherapy alone or combined with targeted therapy in mCRC patients by targeting specific narrative publications in Medline up to April 2022 to identify appropriate original scientific articles. After the literature search, 26 original articles met inclusion and exclusion criteria and were included in the study. Our results show that all investigations conducted on this field have provided generally promising results in predicting the response to therapy or toxic side-effects. By a meta-analytic approach we found that the overall weighted means of the area under the receiver operating characteristic (ROC) curve (AUC) were 0.90, 95% C.I. 0.80–0.95 and 0.83, 95% C.I. 0.74–0.89 in training and validation sets, respectively, indicating a good classification performance in discriminating response vs. non-response. The calculation of overall HR indicates that learning models have strong ability to predict improved survival. Lastly, the delta-radiomics and the 74 gene signatures were able to discriminate response vs. non-response by correctly identifying up to 99% of mCRC patients who were responders and up to 100% of patients who were non-responders. Specifically, when we evaluated the predictive models with tests reaching 80% sensitivity (SE) and 90% specificity (SP), the delta radiomics showed an SE of 99% and an SP of 94% in the training set and an SE of 85% and SP of 92 in the test set, whereas for the 74 gene signatures the SE was 97.6% and the SP 100% in the training set.

MicroRNAs (miRNAs) can regulate the expression of potentially every transcript in the cell, and the definition of miRNA-target interactions is crucial to understand their role in all biological processes. However, the identification of the miRNAs that target a specific mRNA remains a challenge. Here, we describe an innovative method called miR-CATCHv2.0 for the high-throughput identification of the miRNA species bound to an RNA of interest. We also describe how this method can overcome the limitations of the current computational and experimental methods available in this field.

MicroRNAs (miRNAs) can regulate the expression of potentially every transcript in the cell, and the definition of miRNA-target interactions is crucial to understand their role in all biological processes. However, the identification of the miRNAs that target a specific mRNA remains a challenge. Here, we describe an innovative method called miR-CATCHv2.0 for the high-throughput identification of the miRNA species bound to an RNA of interest. We also describe how this method can overcome the limitations of the current computational and experimental methods available in this field.

Simple Summary DDX3X has been repeatedly proven to exert pro-tumorigenic effects in breast cancer. An inhibitor of DDX3X helicase function, FHP01, exhibited very effective antiproliferative and tumor killing activity in vitro and in vivo against several breast cancer cells, yet showed no systemic toxicity but good bioavailability in vivo. FHP01 may therefore represent a novel therapeutic approach with high potential as a personalized strategy for breast cancer, confirming DDX3X as a new actionable molecular target in these tumors. Abstract Inhibition of DDX3X expression or activity reduces proliferation in cells from various tumor tissues, in particular in breast cancer, and its expression often correlates to tumor aggressiveness. This makes DDX3X a prominent candidate for the design of drugs for novel personalized therapeutic strategies. Starting from an in silico drug discovery approach, a group of molecules has been selected by molecular docking at the RNA binding site of DDX3X. Here, the most promising among them, FHP01, was evaluated in breast cancer preclinical models. Specifically, FHP01 exhibited very effective antiproliferative and killing activity against different breast cancer cell types, among which those from triple-negative breast cancer (TNBC). Interestingly, FHP01 also inhibited WNT signaling, a key tumorigenic pathway already correlated to DDX3X functions in breast cancer model cell lines. Ultimately, FHP01 also caused a significant reduction, in vivo, in the growth of MDA MB 231-derived TNBC xenograft models. Importantly, FHP01 showed good bioavailability and no toxicity on normal peripheral blood mononuclear cells in vitro and on several mouse tissues in vivo. Overall, our data suggest that the use of FHP01 and its related compounds may represent a novel therapeutic approach with high potential against breast cancer, including the triple-negative subtype usually correlated to the most unfavorable outcomes because of the lack of available targeted therapies.

The most popular approach to explain cancer is based on the discovery of oncogenes and tumor suppressor genes as a preliminary step in estimating their impact on altered pathways. The present paper proposes a pipeline which aims at detecting "weak" or "indirect" functions impacted by Copy Number Variations (CNVs) of cancer-related genes, integrating such signals over all known oncogenes/tumor suppressor genes of a cancer type. We applied the pipeline to the task of detecting the aberrant functional effects of these alterations across ovarian cancer patients from The Cancer Genome Atlas (TCGA) data.

Among the Herpesviridae, human cytomegalovirus (HCMV) owns the largest genome and displays a huge coding potential. Here, we characterized the UL5 gene product (pUL5) of the clinical isolate TR strain. The protein was predicted as a 166-amino-acid membrane protein with a theoretical mass of 19 kDa. Recombinant virus expressing pUL5 with a tag allowed the identification of two pUL5 non-glycosylated species of approximately 19 and 9 kDa, expressed with early and late kinetic respectively. Experiments in infection confirmed that the lower molecular weight species was translated from an internal ATG in the UL5 open reading frame. Confocal microscopy analysis showed that pUL5 localized within the assembly compartment, but is not incorporated in the virion, as shown by Western blot on purified viral particles. Finally, pull-down experiments coupled with mass spectrometry analysis identified IQGAP1 as a pUL5 interactor, giving new hints on possible roles of pUL5 during HCMV infection.

Background: TBX5 mutations cause Holt-Oram syndrome (HOS) characterized by upper limb and cardiac malformations, but can also contribute to early-onset of atrial fibrillation. Focusing on miRNAs involved in TBX5 regulatory circuits with a cardiac relevant role, we identified miR-182-5p, belonging to miR-183 cluster, found upregulated in Tbx5-depleted hearts of mouse and zebrafish embryos. Methods: To functionally analyse the miR-182-5p role in developing heart, miR-182-5p was dysregulated in zebrafish zygotes of Tg(Myl7:EGFP) and Tg(myl7:gCaMP) transgenic lines. To stably deregulate miR-182-5p in zebrafish heart we exploited the Gal4/UAS system to restrict the miR-182 expression into cardiac context. For physiological analyses we performed the mechanogram of cardiac contraction and electrocardiogram recording. To understand miR-182-5p downstream regulation, in silico analyses, followed by ddPCR/real-time quantifications on dissected zebrafish hearts and rescue experiments both in transient and stable miR-182-5p overexpressing zebrafish embryos were performed. Results: Depletion of Tbx5 from cardiomyocytes increased the expression of miR-182 cluster family that is controlled by Kruppel-like factor 4 (KLF4), a transcription factor repressed by Tbx5. Both transient and stable upregulation of miR-182 in zebrafish affect heart morphology, calcium handling and the onset of arrhythmia while its cardiac-specific downregulation decreases cardiac defects in zebrafish HOS hearts. Expression analyses on selected miR-182-5p putative targets revealed that several calcium channel proteins resulted downregulated in miR-182-5p overexpressing hearts. Transgenic zebrafish line stably overexpressing miR-182-5p in the heart manifested arrhythmia overtime with or without cardiac structural defects. Conclusion: We identified miR-182-5p as a potential suitable target to interfere in the circuit between upstream genetic abnormalities and downstream effectors leading to arrhythmia occurrence.

Here we present miR-CATCHv2.0, an implemented experimental method that allows the identification of the microRNA species directly bound to an RNA of interest. After cross-linking of microRNA::RNA::Ago2 complexes using formaldehyde, the RNA is fragmented using sonication and then subjected to affinity purification using two sets of biotinylated tiling probes (ODD and EVEN). Finally, enriched microRNA species are retrieved by means of small RNA sequencing coupled with an ad hoc analytical workflow. In BRAFV600E mutant A375 melanoma cells, miR-CATCHv2.0 allowed us to identify 20 microRNAs that target X1, the most abundant isoform of BRAF mRNA. These microRNAs fall into different functional classes, according to the effect that they exert (decrease/increase in BRAFV600E mRNA and protein levels) and to the mechanism they use to achieve it (destabilization/stabilization of X1 mRNA or decrease/increase in its translation). microRNA-induced variations in BRAFV600E protein levels are most of the times coupled to consistent variations in pMEK levels, in melanoma cell proliferation in vitro and in sensitivity to the BRAF inhibitor vemurafenib in a xenograft model in zebrafish. However, microRNAs exist that uncouple the degree of activation of the ERK pathway from the levels of BRAFV600E protein. Our study proposes miR-CATCHv2.0 as an effective tool for the identification of direct microRNA-target interactions and, by using such a tool, unveils the complexity of the post-transcriptional regulation to which BRAFV600E and the ERK pathway are subjected in melanoma cells.

Metabolomics/lipidomics are important tools to identify novel biomarkers associated with liver damage. Patients with chronic liver disease (CLD) and hepatitis C virus (HCV) infection often have alterations in glucose, lipid and protein metabolism. The aim of this study was to evaluate if dysfunctional lipid and amino acid metabolism was associated with fibrosis severity and insulin resistance in CLD/HCV patients. We analyzed the baseline sera of 75 subjects with CLD/HCV infection HCV genotype-1, with proven liver biopsy prior to antiviral treatment. We measured amino acid (AA) and lipid concentration by gas and liquid chromatography-mass spectrometry respectively. Alterations in peripheral glucose metabolism due to insulin resistance (IR) were assesed by HOMA-IR (Glucose x Insulin/22.5), while adipose tissue IR was estimated as (Adipo-IR = Free Fatty Acids x Insulin). Baseline HOMA-IR and Adipo-IR were related to the degree of liver fibrosis. Reduction in ceramides 18:1/22:0, 18:1/24:0, diacylglycerol 42:6 and increased phosphocholine 40:6 were associated with higher fibrosis. Adipo-IR was related to lower levels of lysophosphatidylcholine 14:0 and 18:2 and with higher levels of sphingomyelin 18:2/24:0 and 18:2/24:1. Almost all AA were positively associated with Adipo-IR but not with HOMA-IR. We further confirmed the potential use of metabolomics and lipidomics in CLD/HCV subjects finding novel biomarkers of hepatic fibrosis and show that the adipose tissue IR is associated with more severe liver disease and is an important marker not only of altered lipid but also AA metabolism.

Copy Number Variants (CNVs) are structural rearrangements contributing to phenotypic variation but also associated with many disease states. In recent years, the identification of CNVs from high-throughput sequencing experiments has become a common practice for both research and clinical purposes. Several computational methods have been developed so far. In this unit, we describe and give instructions on how to run two read count-based tools, XCAVATOR and EXCAVATOR2, which are tailored for the detection of both germline and somatic CNVs from different sequencing experiments (whole-genome, whole-exome, and targeted) in various disease contexts and population genetic studies. © 2018 by John Wiley & Sons, Inc.

In this short report, we pinpoint some technical and conceptual flaws that we found in the article entitled "miR-204-5p and miR-211-5p contribute to BRAF inhibitor resistance in melanoma" (Diaz-Martinez et al., Cancer Research 2018). We also discuss how, in our opinion, these flaws led Diaz-Martinez and colleagues to incorrect conclusions about the biological role that miR-204 and miR-211 play in melanoma and about the terms of their involvement in the phenomenon of resistance to BRAF inhibitors.

Thyroid hormone (T3) dyshomeostasis in the cardiac ischemia-reperfusion (IR) setting negatively impacts on mitochondria function and extracellular matrix remodeling. The modulation of cardiac miRNAs may represent the underlying molecular mechanisms, but a systems biology perspective investigating this critical issue in depth is still lacking. A rat model of myocardial IR, with or without an early short-term T3-replacement, was used to predict putative T3-dependent miRNA-gene interactions targeted to mitochondria quality control and wound healing repair. As evidenced by mRNA and miRNA expression profiling, the T3 supplementation reverted the expression of 87 genes and 11 miRNAs that were dysregulated in the untreated group. In silico crossing and functional analysis of the T3-associated differentially expressed transcripts, identified a signature of interconnected miRNA-gene regulatory circuits that confer resistance to noxious cascades of acute stress. In this network the T3-down-regulated Tp53, Jun and Sp1 transcription factors emerge as critical nodes linking intrinsic cell death and oxidative stress pathways to adverse remodeling cascades. The data presented here provide a novel insight into the molecular basis of T3 cardioprotection in the early post-IR phase and highlight the contribution of a previously unappreciated complex T3-regulatory network that may be helpful in translating T3 replacement into clinical practice.

Congenital myotonic dystrophy type 1 (CDM1) is characterized by severe symptoms that affect patients from birth, with 40% mortality in the neonatal period and impaired skeletal muscle development. In this paper, we examined the relationship between autophagy and abnormal myogenic differentiation of CDM1 myoblasts. We investigated these pathological features at both ultrastructural and molecular levels, utilizing two CDM1 foetal myoblasts, CDM13 and CDM15, with 1800 and 3200 repeats, respectively. The congenital nature of these CDM1 myoblasts was confirmed by the high methylation level at the DMPK locus. Our results indicated that abnormal autophagy was independent of myogenic differentiation, as CDM13 myoblasts differentiated as well as control myoblasts but underwent autophagy like CDM15, displaying impaired differentiation. miRNA expression profiles revealed that CDM15 myoblasts failed to upregulate the complex network of myo-miRNAs under MYOD and MEF2A control, while this network was upregulated in CDM13 myoblasts. Interestingly, the abnormal differentiation of CDM15 myoblasts was associated with cellular stress accompanied by the induction of the interferon type 1 pathway (innate immune response). Indeed, inhibition of the interferon (IFN) type I pathway restores myogenic differentiation of CDM15 myoblasts, suggesting that the inappropriate activation of the innate immune response might contribute to impaired myogenic differentiation and severe muscle symptoms observed in some CDM1 patients. These findings open up the possibility of new therapeutic approaches to treat CDM1.

MicroRNAs, small molecules of 22-25 nt, inhibit translation of target mRNAs and with transcription factors comprise two major layers of gene regulatory networks with strictly intercon- nected activities. Tbx5, a dosage sensitive gene, is a pivotal player involved in heart/limbs development and its mutations are responsi- ble of the Holt-Oram syndrome (HOS) in human, characterized by upper limb malformations and congenital heart defects (CHD)s both in morphology and electrophysiology.

On the grounds that miRNAs present in the blood of prostate cancer (PCa) patients are released in the growth medium by PCa cells, it is conceivable that PCa cells resistant to docetaxel (DCT) (DCTR) will release miRNAs that may be found in PCa patients under DCT therapy if resistant PCa cells appear. We isolated DCTR clones respectively from 22Rv1 and DU-145 PCa cell lines and performed through next-generation sequencing (NGS) the miRNAs profiles of the released miRNAs. The analysis of the NGS data identified 105 and 1 miRNAs which were differentially released in the growth medium of the 22Rv1/DCTR and DU-145/DCTR clones, respectively. Using additional filters, we selected 12 and 1 miRNA more released by all 22Rv1/DCTR and DU-145/DCTR clones, respectively. Moreover, we showed that 6 of them were more represented in the growth medium of the DCTR cells than the ones of DCT-treated cells. We speculated that they have the pre-requisite to be tested as predictive biomarkers of the DCT resistance in PCa patients under DCT therapy. We propose the utilization of clones resistant to a given drug as in vitro model to identify the differentially released miRNAs, which in perspective could be tested as predictive biomarkers of drug resistance in tumor patients under therapy.

Gene drives have enormous potential for the control of insect populations of medical and agricultural relevance. By preferentially biasing their own inheritance gene drives can rapidly introduce genetic traits even if these confer a negative fitness on the population.We have recently developed gene drives based on a CRISPR nuclease construct that is designed to disrupt key genes essential for female fertility in the malaria mosquito. The construct copies itself and the associated genetic disruption from one homologous chromosome to another during gamete formation, in a process called homing that ensures the majority of offspring inherit the drive. Such drives have the potential to cause long-lasting, sustainable population suppression though they are also expected to impose a large selection pressure for resistance in the mosquito. One of these population suppression gene drives showed rapid invasion of a caged population over 4 generations, establishing proof of principle for this technology. In order to assess the potential for the emergence of resistance to the gene drive in this population we allowed it to run for 25 generations and monitored the frequency of the gene drive over time. Following the initial increase of the gene drive we noticed a gradual decrease in its frequency that was accompanied emergence of small, nuclease-induced mutations at the target gene that are resistant to further cleavage and restore its functionality. Such mutations show rates of increase consistent with positive selection in the face of the gene drive. Our findings represent the first documented example of selection for resistance to a synthetic gene drive and lead to important design recommendations and considerations in order to mitigate for resistance for future gene drive applications.

Objective: Ventricular assist device (VAD) utilization is increasing in children with heart failure (HF) unresponsive to medical therapy, allowing for bridge to transplantation. The use of circulating biomarkers is crucial for the management of adult patients with HF, but their role in pediatric setting is unknown. In this study we assessed the effects of VAD implant on plasma levels of inflammatory and cardiac biomarkers as well as of circulating micro RNAs (miRNAs) in pediatric patients with HF. Methods: Blood samples were collected from 12 HF pediatric patients [73.91±27.94 (mean±SD) months, 6 males, 19±3.8 LVEF%, Interagency Registry for Mechanically Assisted Circulatory Support (INTERMACS) profiles 1/2] before and at 4 hrs, 1, 3, 7, 14 and 30 days after VAD implant. Indications for support was idiopathic dilated (9 patients), non compaction (2 patients) and hypertrophic cardiomyopathy (1 patient). Inflammatory and cardiac biomarkers (Figure 1A) were measured in frozen plasma through the use of immunoassay tests. The circulating miRNAs were extracted from serum samples and miRNA profile was determined by NGS. and validated by real-time PCR. Results: Inflammatory and cardiac biomarker levels during time-course were reported in Figure 1A. Cytokines increased significantly up to 4 hrs after VAD implant, while cardiac biomarkers up to 1 day. Only NT-proBNP decreased significantly up to 1 month of VAD implant compared to pre-VAD values. After sequencing, a total of 169 miRNAs were detected in serum samples. Among them, 13 miRNAs were simultaneously modulated at 1 month after VAD implant compared to pre-VAD levels, i.e. 9 miRNAs decrease and 4 increase (Figure 1B). Real-time PCR validation confirmed a reduction for miR-409-3p and miR-432-5p. Finally, significant correlations with bio-humoral markers were observed for these miRNAs (Figure 1 C). Conclusions: VAD implant induced modifications in circulating levels of cytokine, cardiac biomarkers and miRNAs. NT-proBNP levels and miRNAs significantly decreased after 1 month compared to the values pre-VAD, suggesting the involvement of these biomarkers in molecular mechanisms triggered by VAD implant in children.

// Marianna Vitiello 1, 2 , Andrea Tuccoli 1 , Romina D’Aurizio 3 , Samanta Sarti 1, 4 , Laura Giannecchini 1 , Simone Lubrano 1, 4 , Andrea Marranci 1, 4 , Monica Evangelista 2 , Silvia Peppicelli 5 , Chiara Ippolito 6 , Ivana Barravecchia 7 , Elena Guzzolino 7 , Valentina Montagnani 8 , Michael Gowen 9 , Elisa Mercoledi 1 , Alberto Mercatanti 2 , Laura Comelli 2 , Salvatore Gurrieri 1 , Lawrence W. Wu 10 , Omotayo Ope 10 , Keith Flaherty 11 , Genevieve M. Boland 11 , Marc R. Hammond 11 , Lawrence Kwong 12 , Mario Chiariello 2, 13 , Barbara Stecca 8 , Gao Zhang 10 , Alessandra Salvetti 14 , Debora Angeloni 7 , Letizia Pitto 2 , Lido Calorini 5 , Giovanna Chiorino 15 , Marco Pellegrini 3 , Meenhard Herlyn 10 , Iman Osman 9 , Laura Poliseno 1, 2 1 Oncogenomics Unit, Core Research Laboratory, Istituto Toscano Tumori (ITT), AOUP, Pisa, Italy 2 Institute of Clinical Physiology (IFC), CNR, Pisa, Italy 3 Laboratory of Integrative Systems Medicine (LISM), Institute of Informatics and Telematics (IIT), CNR, Pisa, Italy 4 University of Siena, Italy 5 Section of Experimental Pathology and Oncology, Department of Experimental and Clinical Biomedical Sciences, University of Firenze, Italy 6 Unit of Histology, Department of Clinical and Experimental Medicine, University of Pisa, Italy 7 Scuola Superiore Sant’Anna, Pisa, Italy 8 Tumor Cell Biology Unit, Core Research Laboratory, Istituto Toscano Tumori (ITT), AOUC, Firenze, Italy 9 New York University, New York, NY, USA 10 The Wistar Institute, Philadelphia, PA, USA 11 Massachusetts General Hospital, Boston, MA, USA 12 MD Anderson Cancer Center, Houston, TX, USA 13 Signal Transduction Unit, Core Research Laboratory, Istituto Toscano Tumori (ITT), AOUS, Siena, Italy 14 Unit of Experimental Biology and Genetics, Department of Clinical and Experimental Medicine, University of Pisa, Italy 15 Fondazione Edo and Elvo Tempia, Biella, Italy Correspondence to: Laura Poliseno, email: laura.poliseno@gmail.com Keywords: melanoma, BRAFV600E, ERK pathway, miR-204 family, context-dependency Received: March 11, 2016 Accepted: February 06, 2017 Published: March 06, 2017 ABSTRACT Despite increasing amounts of experimental evidence depicting the involvement of non-coding RNAs in cancer, the study of BRAFV600E-regulated genes has thus far focused mainly on protein-coding ones. Here, we identify and study the microRNAs that BRAFV600E regulates through the ERK pathway. By performing small RNA sequencing on A375 melanoma cells and a vemurafenib-resistant clone that was taken as negative control, we discover miR-204 and miR-211 as the miRNAs most induced by vemurafenib. We also demonstrate that, although belonging to the same family, these two miRNAs have distinctive features. miR-204 is under the control of STAT3 and its expression is induced in amelanotic melanoma cells, where it acts as an effector of vemurafenib’s anti-motility activity by targeting AP1S2. Conversely, miR-211, a known transcriptional target of MITF, is induced in melanotic melanoma cells, where it targets EDEM1 and consequently impairs the degradation of TYROSINASE (TYR) through the ER-associated degradation (ERAD) pathway. In doing so, miR-211 serves as an effector of vemurafenib’s pro-pigmentation activity. We also show that such an increase in pigmentation in turn represents an adaptive response that needs to be overcome using appropriate inhibitors in order to increase the efficacy of vemurafenib. In summary, we unveil the distinct and context-dependent activities exerted by miR-204 family members in melanoma cells. Our work challenges the widely accepted “same miRNA family = same function” rule and provides a rationale for a novel treatment strategy for melanotic melanomas that is based on the combination of ERK pathway inhibitors with pigmentation inhibitors.

Safe recombinant vaccines, based on a small number of antigenic proteins, are emerging as the most attractive, cost-effective solution against infectious diseases. In the present work, we confirmed previous data from our laboratory showing that whole viable bacterial cell treatment with proteases followed by the identification of released peptides by mass spectrometry is the method of choice for the rapid and reliable identification of vaccine candidates in Gram-positive bacteria. When applied to the Group B Streptococcus COH1 strain, 43 surface-associated proteins were identified, including all the protective antigens described in the literature as well as a new protective antigen, the cell wall-anchored protein SAN_1485 belonging to the serine-rich repeat protein family. This strategy overcomes the difficulties so far encountered in the identification of novel vaccine candidates and speeds up the entire vaccine discovery process by reducing the number of recombinant proteins to be tested in the animal model.

Holt-Oram Syndrome (HOS) is an autosomal dominant heart-hand syndrome caused by mutations in the TBX5 gene, a transcription factor capable of regulating hundreds of cardiac-specific genes through complex transcriptional networks. Here we show that, in zebrafish, modulation of a single miRNA is sufficient to rescue the morphogenetic defects generated by HOS. The analysis of miRNA-seq profiling revealed a decreased expression of miR-19a in Tbx5-depleted zebrafish embryos compared to the wild type. We revealed that the transcription of the miR-17-92 cluster, which harbors miR-19a, is induced by Tbx5 and that a defined dosage of miR-19a is essential for the correct development of the heart. Importantly, we highlighted that miR-19a replacement is able to rescue cardiac and pectoral fin defects and to increase the viability of HOS zebrafish embryos. We further observed that miR-19a replacement shifts the global gene expression profile of HOS-like zebrafish embryos towards the wild type condition, confirming the ability of miR-19a to rescue the Tbx5 phenotype. In conclusion our data demonstrate the importance of Tbx5/miR-19a regulatory circuit in heart development and provide a proof of principle that morphogenetic defects associated with HOS can be rescued by transient miRNA modulation.

Structural variants are genomic rearrangements larger than 50 bp accounting for around 1% of the variation among human genomes. They impact on phenotypic diversity and play a role in various diseases including neurological/neurocognitive disorders and cancer development and progression. Dissecting structural variants from next-generation sequencing data presents several challenges and a number of approaches have been proposed in the literature. In this mini review, we describe and summarize the latest tools - and their underlying algorithms - designed for the analysis of whole-genome sequencing, whole-exome sequencing, custom captures, and amplicon sequencing data, pointing out the major advantages/drawbacks. We also report a summary of the most recent applications of third-generation sequencing platforms. This assessment provides a guided indication - with particular emphasis on human genetics and copy number variants - for researchers involved in the investigation of these genomic events.

Background By examining the genotype calls generated by the 1000 Genomes Project we discovered that the human reference genome GRCh37 contains almost 20,000 loci in which the reference allele has never been observed in healthy individuals and around 70,000 loci in which it has been observed only in the heterozygous state. Results We show that a large fraction of this rare reference allele (RRA) loci belongs to coding, functional and regulatory elements of the genome and could be linked to rare Mendelian disorders as well as cancer. We also demonstrate that classical germline and somatic variant calling tools are not capable to recognize the rare allele when present in these loci. To overcome such limitations, we developed a novel tool, named RAREVATOR, that is able to identify and call the rare allele in these genomic positions. By using a small cancer dataset we compared our tool with two state-of-the-art callers and we found that RAREVATOR identified more than 1,500 germline and 22 somatic RRA variants missed by the two methods and which belong to significantly mutated pathways. Conclusions These results show that, to date, the investigation of around 100,000 loci of the human genome has been missed by re-sequencing experiments based on the GRCh37 assembly and that our tool can fill the gap left by other methods. Moreover, the investigation of the latest version of the human reference genome, GRCh38, showed that although the GRC corrected almost all insertions and a small part of SNVs and deletions, a large number of functionally relevant RRAs still remain unchanged. For this reason, also future resequencing experiments, based on GRCh38, will benefit from RAREVATOR analysis results. RAREVATOR is freely available at http://sourceforge.net/projects/rarevator. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1481-9) contains supplementary material, which is available to authorized users.

Mechanical unloading by left ventricular assist devices (LVADs) in advanced heart failure (HF), in addition to improving symptoms and end-organ perfusion, is supposed to stimulate cellular and molecular responses which can reverse maladaptive cardiac remodeling. As microRNAs (miRNAs) are key regulators in remodeling processes, a comparative miRNA profiling in transplanted hearts of HF patients with/without LVAD assistance could aid to comprehend underlying molecular mechanisms. Next generation sequencing (NGS) was used to analyze miRNA differential expression in left ventricles of HF patients who underwent heart transplantation directly (n=9) or following a period of LVAD support (n=8). After data validation by quantitative real-time PCR, association with functional clinical parameters was investigated. Bioinformatics’ tools were then used for prediction of putative targets of modulated miRNAs and relative pathway enrichment. The analysis revealed 13 upregulated and 10 downregulated miRNAs in failing hearts subjected to LVAD assistance. In particular, the expression level of some of them (miR-338-3p, miR-142-5p and -3p, miR-216a-5p, miR-223-3p, miR-27a-5p, and miR-378g) showed correlation with off-pump cardiac index values. Predicted targets of these miRNAs were involved in focal adhesion/integrin pathway and in actin cytoskeleton regulation. The identified miRNAs might contribute to molecular regulation of reverse remodeling and heart recovery mechanisms.

Safe recombinant vaccines, based on a small number of antigenic proteins, are emerging as the most attractive, cost-effective solution against infectious diseases. In the present work, we confirmed previous data from our laboratory showing that whole viable bacterial cell treatment with proteases followed by the identification of released peptides by mass spectrometry is the method of choice for the rapid and reliable identification of vaccine candidates in Gram-positive bacteria. When applied to the Group B Streptococcus COH1 strain, 43 surface-associated proteins were identified, including all the protective antigens described in the literature as well as a new protective antigen, the cell wall-anchored protein SAN_1485 belonging to the serine-rich repeat protein family. This strategy overcomes the difficulties so far encountered in the identification of novel vaccine candidates and speeds up the entire vaccine discovery process by reducing the number of recombinant proteins to be tested in the animal model.

Objective: To identify novel genes implicated in epilepsy with auditory features (EAF) in phenotypically heterogeneous families with unknown molecular basis. Methods: We identified 15 probands with EAF in whom an LGI1 mutation had been excluded. We performed electroclinical phenotyping on all probands and available affected relatives. We used whole-exome sequencing (WES) in 20 individuals with EAF (including all the probands and 5 relatives) to identify single nucleotide variants, small insertions/deletions, and copy number variants. Results: WES revealed likely pathogenic variants in genes that had not been previously associated with EAF: a CNTNAP2 intragenic deletion, 2 truncating mutations of DEPDC5, and a missense SCN1A change. Conclusions: EAF is a clinically and molecularly heterogeneous disease. The association of EAF with CNTNAP2, DEPDC5 ,a ndSCN1A mutations widens the phenotypic spectrum related to these genes. CNTNAP2 encodes CASPR2, a member of the voltage-gated potassium channel complex in which LGI1 plays a role. The finding of a CNTNAP2 deletion emphasizes the importance of this complex in EAF and shows biological convergence. Neurol Genet 2015;1:e5; doi

Purpose: MicroRNAs (miRNAs) are short, single-stranded RNAs that anneal with complementary sequences in mRNAs thereby suppressing protein expression and often decreasing mRNA stability. Due to the high number of genes that are potentially targeted by one miRNA and to the relatively low number of different miRNAs coded in the genome, miRNAs are the best candidates to orchestrate gene expression in embryonic development and in response to inducing signals, by integrating distinct transcription factors pathways. Indeed, the manipulation of miRNA expression or function can have a profound impact on cellular phenotype. miRNA are critically involved in almost all the biological processes in health and disease including several cardiovascular disorders. Tbx5 is a transcription factor crucial for heart development and implicated in cardiac function. In human, TBX5 mutations are associated with Holt-Oram Syndrome (HOS), which is characterized by upper limb and congenital heart defects. In mouse model of HOS, expression profiling revealed that Tbx5 regulates hundreds of genes through complex transcriptional networks, acting both directly on gene expression and indirectly as “regulating of regulators”. Our data indicate that miRNAs can be important effectors of Tbx5. The goal of our project is to better characterize the miRNAs/Tbx5 regulatory network with the aim: a) to understand the molecular roots of the complex HOS phenotype; b) to verify if it might be possible, by miRNA modulation, to restore the aberrant pathways generated by Tbx5 alteration. Methods: Using zebrafish as model systems of HOS, we characterized the miRNA profiles of wild type (WT) and Tbx5-depleted embryos by NGS. The impact of differentially expressed miRNAs on cardiac development of WT embryos and Tbx5 morphants was investigated by gain and loss of function experiments. Results: Among the identified miRNAs responsive to Tbx5 dosage we focused on miR-19a, whose expression decreases in HOS zebrafish embryos. We showed that: a) miR-19a dysregulation affects heart development; b) miR-19a modulates the expression of genes important for epicardial development and retinoic acid metabolism; c) miR-19a replacement was able to rescue cardiac and pectoral fin defects and increase the viability of Tbx5 morphants. Conclusions-Our data demonstrate that in zebrafish the modulation of a single miRNA is sufficient to reduce the defects generated by a complex syndrome like HOS, supporting enthusiasm for the continued exploration of miRNAs as a new class of drug target.

We developed a novel software tool, EXCAVATOR, for the detection of copy number variants (CNVs) from whole-exome sequencing data. EXCAVATOR combines a three-step normalization procedure with a novel heterogeneous hidden Markov model algorithm and a calling method that classifies genomic regions into five copy number states. We validate EXCAVATOR on three datasets and compare the results with three other methods. These analyses show that EXCAVATOR outperforms the other methods and is therefore a valuable tool for the investigation of CNVs in largescale projects, as well as in clinical research and diagnostics. EXCAVATOR is freely available at http://sourceforge.net/projects/excavatortool/. © 2013 Magi et al.; licensee BioMed Central Ltd.

Introduction: Regorafenib is indicated for the treatment of mCRC patients who have failed all other therapies. Nevertheless a substantial percentage of patients experiences rapid disease progression (PD) and serious adverse events may occur. For these reasons, clinical and/or molecular markers able to improve the cost/benefit ratio are urgently needed. Circulating microRNAs (c-miRNAs) have been recognized as possible prognostic and diagnostic markers in mCRC. The aim of this study was to describe the early changes in plasma levels of 10 selected c-miRNAs during the treatment with regorafenib and to investigate their correlation with clinical outcome. Methods: Plasma samples of patients treated with regorafenib at our Institution were collected at baseline (D1) and after 15 days of treatment (D15). Plasma levels of c-miR-17, c-miR-21, c-miR-29, c-miR-34, c-miR-92, c-miR-126, c-miR-141, c-miR-221, c-miR-601, c-miR-760 were analysed by means of real-time PCR. Paired levels at D1 and D15 were compared by means of Wilcoxon test for each c-miRNA. C-miRNAs showing significant changes were further analysed in order to identify possible correlations with outcome. Results: Thirty-four patients were included in the present study. Main characteristics were the following: M/F = 50%/50%; median age = 65 (range 48-78 years); ECOG-PS 0/1-2 = 71%/29%; time from diagnosis of metastases = 18 months 15%/85%. Median PFS and OS were 2.4 and 6.5 months, respectively. One (3%) patient achieved a response and 16 (47%) had disease stabilization (disease control rate: 50%). As compared to D1, the following c-miRNAs increased at D15: c-miR-601 (p = 0.01), c-miR-141 (p = 0.04) and c-miR-21 (p = 0.06). Despite a median increase in the overall population, 12 (35%) out of 34 patients showed reduced level of c-miR-21 at D15. Nine out of 12 (75%) patients with reduced levels of c-miR-21 achieved disease control, as compared to 8 out of 23 (35%) patients with increased levels (Fisher 's Exact Test, p = 0.035). Median PFS of patients with increased and decreased level of levels of c-miR-21 were 2.1 and 3.9 months, respectively (HR = 1.89 95%CI 0.92-4.14 p = 0.08). Data on OS are not yet mature. Early modifications of c-miR-21 levels showed a sensitivity of 82% in predicting benefit from regorafenib. Conclusion: The early modulation of c-miR-21 levels may predict benefit from regorafenib in terms of disease control. These results need validation in independent series.