Biological Random Walks: multi-omics integration for disease gene prioritization

Motivation: Over the past decade, network-based approaches have proven useful in identifying disease modules within the human interactome, often providing insights into key mechanisms and guiding the quest for therapeutic targets. This is all the more important, since experimental investigation of potential gene candidates is an expensive task, thus not always a feasible option. On the other hand, many sources of biological information exist beyond the interactome and an important research direction is the design of effective techniques for their integration. Results: In this work, we introduce the Biological Random Walks (BRW) approach for disease gene prioritization in the human interactome. The proposed framework leverages multiple biological sources within an integrated framework. We perform an extensive, comparative study of BRW's performance against well-established baselines. Availability and implementation: All code is publicly available and can be downloaded at \url{https://github.com/LeoM93/BiologicalRandomWalks}. We used publicly available datasets, details on their retrieval and preprocessing are provided in the supplementary material

Clinical exome performance for reporting secondary genetic findings.

BACKGROUND : Reporting clinically actionable incidental genetic findings in the course of clinical exome testing is recommended by the American College of Medical Genet- ics and Genomics (ACMG). However, the performance of clinical exome methods for reporting small subsets of genes has not been previously reported. METHODS : In this study, 57 exome data sets performed as clinical (n ! 12) or research (n ! 45) tests were retrospec- tively analyzed. Exome sequencing data was examined for adequacy in the detection of potentially pathogenic variant locations in the 56 genes described in the ACMG incidental findings recommendation. All exons of the 56 genes were examined for adequacy of sequencing coverage. In addition, nucleotide positions annotated in HGMD (Human Gene Mutation Database) were examined. RESULTS : The 56 ACMG genes have 18336 nucleotide variants annotated in HGMD. None of the 57 exome data sets possessed a HGMD variant. The clinical exome test had inadequate coverage for " 50% of HGMD vari- ant locations in 7 genes. Six exons from 6 different genes had consistent failure across all 3 test methods; these exons had high GC content (76%–84%). CONCLUSIONS : The use of clinical exome sequencing for the interpretation and reporting of subsets of genes requires recognition of the substantial possibility of inadequate depth and breadth of sequencing coverage at clinically relevant locations. Inadequate depth of coverage may contribute to false-negative clinical ex- ome results

Update on Fundamental Mechanisms of Thyroid Cancer

The incidence of thyroid cancer (TC) has increased worldwide over the past four decades. TC is divided into three main histological types: differentiated (papillary and follicular TC), undifferentiated (poorly differentiated and anaplastic TC), and medullary TC, arising from TC cells. This review discusses the molecular mechanisms associated to the pathogenesis of different types of TC and their clinical relevance. In the last years, progresses in the genetic characterization of TC have provided molecular markers for diagnosis, risk stratification, and treatment targets. Recently, papillary TC, the most frequent form of TC, has been reclassified into two molecular subtypes, named BRAF-like and RAS-like, associated to a different range of cancer risks. Similarly, the genetic characterization of follicular TC has been proposed to complement the new histopathological classification in order to estimate the prognosis. New analyses characterized a comprehensive molecular profile of medullary TC, raising the role of RET mutations. More recent evidences suggested that immune microenvironment associated to TC may play a critical role in tumor invasion, with potential immunotherapeutic implications in advanced and metastatic TC. Several types of ancillary approaches have been developed to improve the diagnostic value of fine needle aspiration biopsies in indeterminate thyroid nodules. Finally, liquid biopsy, as a non-invasive diagnostic tool for body fluid genotyping, brings a new prospective of disease and therapy monitoring. Despite all these novelties, much work remains to be done to fully understand the pathogenesis and biological behaviors of the different types of TC and to transfer this knowledge in clinical practice

Risk stratification of neck lesions detected sonographically during the follow-up of differentiated thyroid cancer

Context: The European Thyroid Association (ETA) has classified post-treatment cervical ultrasound findings in thyroid cancer patients based on their association with disease persistence/recurrence. Objective: To assess this classification's ability to predict the growth and persistence of such lesions during active post-treatment surveillance of patients with differentiated thyroid cancer (DTC). Design: Retrospective, observational study Setting: Thyroid cancer center, large Italian teaching hospital. Patients: Center referrals (2005–2014) were reviewed and patients selected with pathologically confirmed DTC; total thyroidectomy, with or without neck dissection and/or radioiodine remnant ablation; abnormal findings on ≥2 consecutive post-treatment neck sonograms; subsequent follow-up consisting of active surveillance. Baseline ultrasound abnormalities (thyroid bed masses, lymph nodes) were classified according to the ETA system. Patients were divided into group S (those with ≥1 lesion classified as ‘suspicious’) and group I (‘indeterminate’ lesions only). We recorded baseline and follow-up clinical data through 30 June 2015. Main Outcomes: Patients with growth (> 3 mm, largest diameter) of ≥1 lesion during follow-up, patients with ≥1 persistent lesion at the final visit. Results: The cohort included 58 (9%) of the 637 DTC cases screened. A total of 113 lesions were followed (18 thyroid bed masses, 95 lymph nodes). During surveillance (median 3.7 years), group I had significantly lower rates than group S of lesion growth (8% vs. 36%, p=0.01) and persistence (64% vs. 97%, p=0.014). Median time to scan normalization: 2.9 years. Conclusions: The ETA's evidence-based classification of sonographically detected neck abnormalities can help identify PTC patients eligible for more relaxed follow-up

Reduced expression of THRβ in papillary thyroid carcinomas: relationship with BRAF mutation, aggressiveness and miR expression

Purpose Down-regulation of thyroid hormone receptor beta (THRβ) gene has been described in several human malignancies, including thyroid cancer. In this study, we analyzed THRβ mRNA expression in surgical specimens from a series of human papillary thyroid carcinomas (PTCs), characterized by their genotypic and clinical–biological features. Methods Thirty-six PTCs were divided into two groups according to the 2009 American Thyroid Association risk classification (17 low, 19 intermediate), and each group was divided into subgroups based on the presence or absence of the BRAFV600E mutation (21 BRAF mutated, 15 BRAF wild type). Gene expression was analyzed using fluidic cards containing probes and primers specific for the THRβ gene, as well as for genes of thyroperoxidase (TPO), sodium/iodide symporter (NIS), thyroglobulin (Tg) and thyroid stimulating hormone receptor (TSH-R) and for some miRNAs involved in thyroid neoplasia and targeting THRβ. The mRNA levels of each tumor tissue were compared with their correspondent normal counterpart. Results THRβ transcript was down-regulated in all PTCs examined. No significant differences were found between intermediate- vs low-risk PTCs patients, and BRAF-mutated vs BRAF wild-type groups. THRβ expression was directly correlated with NIS, TPO, Tg and TSH-R, and inversely correlated to miR-21, -146a, -181a and -221 expression. Conclusions Our results demonstrate that down-regulation of THRβ is a common feature of PTCs. While it is not associated with a more aggressive phenotype of PTC, it correlates with the reduction of all the markers of differentiation and is associated with overexpression of some miRNAs supposed to play a role in thyroid tumorigenesis

Fibronectin-1 expression is increased in aggressive thyroid cancer and favors the migration and invasion of cancer cells

n this study we analyzed the expression levels of markers of epithelial-to-mesenchymal transition (EMT) in several papillary thyroid carcinomas (PTCs) and the relation with tumor genotypes and clinicopathological characteristics. The role of fibronectin-1 (FN1) was investigated by analyzing the effects of FN1 silencing in two human thyroid cancer cell lines. Most of EMT markers were significantly over-expressed in a group of 36 PTCs. In particular, FN1 mRNA levels were higher in tumor vs non-tumor tissue (117.3, p < 0.001) and also in aggressive and BRAF(V600E) samples. Similar results were observed (and confirmed at the protein level) when FN1 expression was analyzed in a validation group of 50 PTCs and six lymph node (LN) metastases. Silencing of FN1 in TPC-1 and BCPAP thyroid cancer cells significantly reduced proliferation, adhesion, migration, and invasion in both cell lines. Collectively, our data indicate that FN1 overexpression is an important determinant of thyroid cancer aggressiveness

RET mutation and increased angiogenesis in medullary thyroid carcinomas

Advanced medullary thyroid cancers (MTCs) are now being treated with drugs that inhibit receptor tyrosine kinases, many of which involved in angiogenesis. Response rates vary widely, and toxic effects are common, so treatment should be reserved for MTCs likely to be responsive to these drugs. RET mutations are common in MTCs, but it is unclear how they influence the microvascularization of these tumors. We examined 45 MTCs with germ-line or somatic RET mutations (RETmut group) and 34 with wild-type RET (RETwt). Taqman Low-Density Arrays were used to assess proangiogenic gene expression. Immunohistochemistry was used to assess intratumoral, peritumoral and nontumoral expression levels of VEGFR1, R2, R3, PDGFRa, PDGFB and NOTCH3. We also assessed microvessel density (MVD) and lymphatic vessel density (LVD) based on CD31-positive and podoplanin-positive vessel counts, respectively, and vascular pericyte density based on staining for a-smooth muscle actin (a-SMA), a pericyte marker. Compared with RETwt tumors, RETmut tumors exhibited upregulated expression of proangiogenic genes (mRNA and protein), especially VEGFR1, PDGFB and NOTCH3. MVDs and LVDs were similar in the two groups. However, microvessels in RETmut tumors were more likely to be a-SMA positive, indicating enhanced coverage by pericytes, which play key roles in vessel sprouting, maturation and stabilization. These data suggest that angiogenesis in RETmut MTCs may be more intense and complete than that found in RETwt tumors, a feature that might increase their susceptibility to antiangiogenic therapy. Given their increased vascular pericyte density, RETmut MTCs might also benefit from combined or preliminary treatment with PDGF inhibitors

A synonymous RET substitution enhances the oncogenic effect of an in-cis missense mutation by increasing constitutive splicing efficiency

Synonymous mutations continue to be filtered out from most large-scale cancer genome studies, but several lines of evidence suggest they can play driver roles in neoplastic disease. We investigated a case of an aggressive, apparently sporadic medullary thyroid carcinoma (MTC) harboring a somatic RET p.Cys634Arg mutation (a known MTC driver). A germ-line RET substitution (p.Cys630=) had also been found but was considered clinically irrelevant because of its synonymous nature. Next generation sequencing (NGS) of the tumor tissues revealed that the RET mutations were in cis. There was no evidence of gene amplification. Expression analysis found an increase of RET transcript in p.Cys630=;p.Cys634Arg patient compared with that found in 7 MTCs harboring p.Cys634 mutations. Minigene expression assays demonstrated that the presence of the synonymous RET mutation was sufficient to explain the increased RET mRNA level. In silico analyses and RNA immunoprecipitation experiments showed that the p.Cys630 = variant created new exonic splicing enhancer motifs that enhanced SRp55 recruitment to the mutant allele, leading to more efficient maturation of its pre-mRNA and an increased abundance of mature mRNA encoding a constitutively active RET receptor. These findings document a novel mechanism by which synonymous mutations can contribute to cancer progression

Identification of Exosomal microRNAs and Their Targets in Papillary Thyroid Cancer Cells

The release of molecules in exosomal cargoes is involved in tumor development and progression. We compared the profiles of exosomal microRNAs released by two thyroid cancer cell lines (TPC-1 and K1) with that of non-tumorigenic thyroid cells (Nthy-ori-3-1), and we explored the network of miRNA–target interaction. After extraction and characterization of exosomes, expression levels of microRNAs were investigated using custom TaqMan Advanced array cards, and compared with those expressed in the total cell extracts. The functional enrichment and network-based analysis of the miRNAs’ targets was also performed. Five microRNAs (miR-21-5p, miR-31-5p, miR-221-3p, miR-222-3p, and let-7i-3p) were significantly deregulated in the exosomes of tumor cells vs. non-tumorigenic cells, and three of them (miR-31-5p, miR-222-3p, and let-7i-3p) in the more aggressive K1 compared to TPC-1 cells. The network analysis of the five miRNAs identified some genes as targets of more than one miRNAs. These findings permitted the identification of exosomal microRNAs secreted by aggressive PTC cells, and indicated that their main targets are regulators of the tumor microenvironment. A deeper analysis of the functional role of the targets of exosomal miRNAs will provide further information on novel targets of molecular treatments for these neoplasms

Increased expression of pro-angiogenic factors and vascularization in thyroid hyperfunctioning adenomas with and without TSH receptor activating mutations

Autonomously functioning thyroid nodules (AFTN) are known to receive an increased blood influx necessary to sustain their high rate of growth and hormone production. Here, we investigated the expression of hematic and lymphatic vases in a series of 20 AFTN compared with the contralateral non-tumor tissues of the same patients, and the transcript levels of proteins involved in the control of vascular proliferation, including the vascular endothelial growth factor (VEGF) and platelet-derived growth factors (PDGF) and their receptors and the endothelial nitric oxide synthase (eNOS). In parallel, the expression of the differentiation markers sodium/iodide symporter (NIS), thyroperoxidase (TPO), thyroglobulin (Tg), and TSH receptor (TSHR) was also investigated. The data were further analyzed comparing subgroups of tumors with or without mutations in the TSHR gene. Analysis by means of CD31 and D2-40 immunostaining showed in AFTN an increased number of hematic, but not lymphatic, vessels in parallel with an enhanced proliferation rate shown by increased Ki67 staining. Quantitative RT-PCR analysis revealed an increase of VEGF, VEGFR1 and 2, PDGF-A, PDGF-B, and eNOS expression in tumor versus normal tissues. Also, higher transcript levels of NIS, TPO, and Tg were detected. Comparison of the two subgroups of samples revealed only few differences in the expression of the genes examined. In conclusion, these data demonstrate an increased expression of angiogenesis-related factors associated with an enhanced proliferation of hematic, but not lymphatic, vessels in AFTNs. In this context, the presence of TSHR mutations may only slightly influence the expression of pro-angiogenic growth factors