Genomic and Transcriptomic Analyses of Breast Cancer Primaries and Matched Metastases in AURORA, the Breast International Group (BIG) Molecular Screening Initiative.

AURORA aims to study the processes of relapse in metastatic breast cancer (MBC) by performing multi-omics profiling on paired primary tumors and early-course metastases. Among 381 patients (primary tumor and metastasis pairs: 252 targeted gene sequencing, 152 RNA sequencing, 67 single nucleotide polymorphism arrays), we found a driver role for GATA1 and MEN1 somatic mutations. Metastases were enriched in ESR1, PTEN, CDH1, PIK3CA, and RB1 mutations; MDM4 and MYC amplifications; and ARID1A deletions. An increase in clonality was observed in driver genes such as ERBB2 and RB1. Intrinsic subtype switching occurred in 36% of cases. Luminal A/B to HER2-enriched switching was associated with TP53 and/or PIK3CA mutations. Metastases had lower immune score and increased immune-permissive cells. High tumor mutational burden correlated to shorter time to relapse in HR+/HER2- cancers. ESCAT tier I/II alterations were detected in 51% of patients and matched therapy was used in 7%. Integration of multi-omics analyses in clinical practice could affect treatment strategies in MBC. SIGNIFICANCE: The AURORA program, through the genomic and transcriptomic analyses of matched primary and metastatic samples from 381 patients with breast cancer, coupled with prospectively collected clinical data, identified genomic alterations enriched in metastases and prognostic biomarkers. ESCAT tier I/II alterations were detected in more than half of the patients.This article is highlighted in the In This Issue feature, p. 2659.

The AURORA pilot study for molecular screening of patients with advanced breast cancer-a study of the breast international group.

Several studies have demonstrated the feasibility of molecular screening of tumour samples for matching patients with cancer to targeted therapies. However, most of them have been carried out at institutional or national level. Herein, we report on the pilot phase of AURORA (NCT02102165), a European multinational collaborative molecular screening initiative for advanced breast cancer patients. Forty-one patients were prospectively enroled at four participating centres across Europe. Metastatic tumours were biopsied and profiled using an Ion Torrent sequencing platform at a central facility. Sequencing results were obtained for 63% of the patients in real-time with variable turnaround time stemming from delays between patient consent and biopsy. At least one clinically actionable mutation was identified in 73% of patients. We used the Illumina sequencing technology for orthogonal validation and achieved an average of 66% concordance of substitution calls per patient. Additionally, copy number aberrations inferred from the Ion Torrent sequencing were compared to single nucleotide polymorphism arrays and found to be 59% concordant on average. Although this study demonstrates that powerful next generation genomic techniques are logistically ready for international molecular screening programs in routine clinical settings, technical challenges remain to be addressed in order to ensure the accuracy and clinical utility of the genomic data.

Clinical management of breast cancer heterogeneity.

Traditionally, intertumour heterogeneity in breast cancer has been documented in terms of different histological subtypes, treatment sensitivity profiles, and clinical outcomes among different patients. Results of high-throughput molecular profiling studies have subsequently revealed the true extent of this heterogeneity. Further complicating this scenario, the heterogeneous expression of the oestrogen receptor (ER), progesterone receptor (PR), and HER2 has been reported in different areas of the same tumour. Furthermore, discordance, in terms of ER, PR and HER2 expression, has also been reported between primary tumours and their matched metastatic lesions. High-throughput molecular profiling studies have confirmed that spatial and temporal intratumour heterogeneity of breast cancers exist at a level beyond common expectations. We describe the different levels of tumour heterogeneity, and discuss the strategies that can be adopted by clinicians to tackle treatment response and resistance issues associated with such heterogeneity, including a rationally selected combination of agents that target driver mutations, the targeting of deleterious passenger mutations, identifying and eradicating the 'lethal' clone, targeting the tumour microenvironment, or using adaptive treatments and immunotherapy. The identification of the most-appropriate strategies and their implementation in the clinic will prove highly challenging and necessitate the adoption of radically new practices for the optimal clinical management of breast malignancies.

Identification of a microRNA landscape targeting the PI3K/Akt signaling pathway in inflammation-induced colorectal carcinogenesis.

Inflammation can contribute to tumor formation; however, markers that predict progression are still lacking. In the present study, the well-established azoxymethane (AOM)/dextran sulfate sodium (DSS)-induced mouse model of colitis-associated cancer was used to analyze microRNA (miRNA) modulation accompanying inflammation-induced tumor development and to determine whether inflammation-triggered miRNA alterations affect the expression of genes or pathways involved in cancer. A miRNA microarray experiment was performed to establish miRNA expression profiles in mouse colon at early and late time points during inflammation and/or tumor growth. Chronic inflammation and carcinogenesis were associated with distinct changes in miRNA expression. Nevertheless, prediction algorithms of miRNA-mRNA interactions and computational analyses based on ranked miRNA lists consistently identified putative target genes that play essential roles in tumor growth or that belong to key carcinogenesis-related signaling pathways. We identified PI3K/Akt and the insulin growth factor-1 (IGF-1) as major pathways being affected in the AOM/DSS model. DSS-induced chronic inflammation downregulates miR-133a and miR-143/145, which is reportedly associated with human colorectal cancer and PI3K/Akt activation. Accordingly, conditioned medium from inflammatory cells decreases the expression of these miRNA in colorectal adenocarcinoma Caco-2 cells. Overexpression of miR-223, one of the main miRNA showing strong upregulation during AOM/DSS tumor growth, inhibited Akt phosphorylation and IGF-1R expression in these cells. Cell sorting from mouse colons delineated distinct miRNA expression patterns in epithelial and myeloid cells during the periods preceding and spanning tumor growth. Hence, cell-type-specific miRNA dysregulation and subsequent PI3K/Akt activation may be involved in the transition from intestinal inflammation to cancer.

Myelin-derived lipids modulate macrophage activity by liver X receptor activation.

Multiple sclerosis is a chronic, inflammatory, demyelinating disease of the central nervous system in which macrophages and microglia play a central role. Foamy macrophages and microglia, containing degenerated myelin, are abundantly found in active multiple sclerosis lesions. Recent studies have described an altered macrophage phenotype after myelin internalization. However, it is unclear by which mechanisms myelin affects the phenotype of macrophages and how this phenotype can influence lesion progression. Here we demonstrate, by using genome wide gene expression analysis, that myelin-phagocytosing macrophages have an enhanced expression of genes involved in migration, phagocytosis and inflammation. Interestingly, myelin internalization also induced the expression of genes involved in liver-X-receptor signaling and cholesterol efflux. In vitro validation shows that myelin-phagocytosing macrophages indeed have an increased capacity to dispose intracellular cholesterol. In addition, myelin suppresses the secretion of the pro-inflammatory mediator IL-6 by macrophages, which was mediated by activation of liver-X-receptor ß. Our data show that myelin modulates the phenotype of macrophages by nuclear receptor activation, which may subsequently affect lesion progression in demyelinating diseases such as multiple sclerosis.

MicroRNAs profiling in murine models of acute and chronic asthma: a relationship with mRNAs targets.

BACKGROUND: miRNAs are now recognized as key regulator elements in gene expression. Although they have been associated with a number of human diseases, their implication in acute and chronic asthma and their association with lung remodelling have never been thoroughly investigated. METHODOLOGY/PRINCIPAL FINDINGS: In order to establish a miRNAs expression profile in lung tissue, mice were sensitized and challenged with ovalbumin mimicking acute, intermediate and chronic human asthma. Levels of lung miRNAs were profiled by microarray and in silico analyses were performed to identify potential mRNA targets and to point out signalling pathways and biological processes regulated by miRNA-dependent mechanisms. Fifty-eight, 66 and 75 miRNAs were found to be significantly modulated at short-, intermediate- and long-term challenge, respectively. Inverse correlation with the expression of potential mRNA targets identified mmu-miR-146b, -223, -29b, -29c, -483, -574-5p, -672 and -690 as the best candidates for an active implication in asthma pathogenesis. A functional validation assay was performed by cotransfecting in human lung fibroblasts (WI26) synthetic miRNAs and engineered expression constructs containing the coding sequence of luciferase upstream of the 3'UTR of various potential mRNA targets. The bioinformatics analysis identified miRNA-linked regulation of several signalling pathways, as matrix metalloproteinases, inflammatory response and TGF-ß signalling, and biological processes, including apoptosis and inflammation. CONCLUSIONS/SIGNIFICANCE: This study highlights that specific miRNAs are likely to be involved in asthma disease and could represent a valuable resource both for biological makers identification and for unveiling mechanisms underlying the pathogenesis of asthma.

Inferring regulatory networks from expression data using tree-based methods.

One of the pressing open problems of computational systems biology is the elucidation of the topology of genetic regulatory networks (GRNs) using high throughput genomic data, in particular microarray gene expression data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) challenge aims to evaluate the success of GRN inference algorithms on benchmarks of simulated data. In this article, we present GENIE3, a new algorithm for the inference of GRNs that was best performer in the DREAM4 In Silico Multifactorial challenge. GENIE3 decomposes the prediction of a regulatory network between p genes into p different regression problems. In each of the regression problems, the expression pattern of one of the genes (target gene) is predicted from the expression patterns of all the other genes (input genes), using tree-based ensemble methods Random Forests or Extra-Trees. The importance of an input gene in the prediction of the target gene expression pattern is taken as an indication of a putative regulatory link. Putative regulatory links are then aggregated over all genes to provide a ranking of interactions from which the whole network is reconstructed. In addition to performing well on the DREAM4 In Silico Multifactorial challenge simulated data, we show that GENIE3 compares favorably with existing algorithms to decipher the genetic regulatory network of Escherichia coli. It doesn't make any assumption about the nature of gene regulation, can deal with combinatorial and non-linear interactions, produces directed GRNs, and is fast and scalable. In conclusion, we propose a new algorithm for GRN inference that performs well on both synthetic and real gene expression data. The algorithm, based on feature selection with tree-based ensemble methods, is simple and generic, making it adaptable to other types of genomic data and interactions.

Hereditary cutaneomucosal venous malformations are caused by TIE2 mutations with widely variable hyper-phosphorylating effects.

Mutations in the angiopoietin receptor TIE2/TEK have been identified as the cause for autosomal dominantly inherited cutaneomucosal venous malformation (VMCM). Thus far, two specific germline substitutions (R849W and Y897S), located in the kinase domain of TIE2, have been reported in five families. The mutations result in a fourfold increase in ligand-independent phosphorylation of the receptor. Here, we report 12 new families with TEK mutations. Although the phenotype is primarily characterized by small multifocal cutaneous vascular malformations, many affected members also have mucosal lesions. In addition, cardiac malformations are observed in some families. Six of the identified mutations are new, with three located in the tyrosine kinase domain, two in the kinase insert domain, and another in the carboxy terminal tail. The remaining six are R849W substitutions. Overexpression of the new mutants resulted in ligand-independent hyperphosphorylation of the receptor, suggesting this is a general feature of VMCM-causative TIE2 mutations. Moreover, variation in the level of activation demonstrates, to the best of our knowledge for the first time, that widely differing levels of chronic TIE2 hyperphosphorylation are tolerated in the heterozygous state, and are compatible with normal endothelial cell function except in the context of highly localized areas of lesion pathogenesis.

Supervised learning with decision tree-based methods in computational and systems biology.

At the intersection between artificial intelligence and statistics, supervised learning allows algorithms to automatically build predictive models from just observations of a system. During the last twenty years, supervised learning has been a tool of choice to analyze the always increasing and complexifying data generated in the context of molecular biology, with successful applications in genome annotation, function prediction, or biomarker discovery. Among supervised learning methods, decision tree-based methods stand out as non parametric methods that have the unique feature of combining interpretability, efficiency, and, when used in ensembles of trees, excellent accuracy. The goal of this paper is to provide an accessible and comprehensive introduction to this class of methods. The first part of the review is devoted to an intuitive but complete description of decision tree-based methods and a discussion of their strengths and limitations with respect to other supervised learning methods. The second part of the review provides a survey of their applications in the context of computational and systems biology.

Genotype-phenotype associations in Sotos syndrome: an analysis of 266 individuals with NSD1 aberrations.

We identified 266 individuals with intragenic NSD1 mutations or 5q35 microdeletions encompassing NSD1 (referred to as "NSD1-positive individuals"), through analyses of 530 subjects with diverse phenotypes. Truncating NSD1 mutations occurred throughout the gene, but pathogenic missense mutations occurred only in functional domains (P < 2 x 10(-16)). Sotos syndrome was clinically diagnosed in 99% of NSD1-positive individuals, independent of the molecular analyses, indicating that NSD1 aberrations are essentially specific to this condition. Furthermore, our data suggest that 93% of patients who have been clinically diagnosed with Sotos syndrome have identifiable NSD1 abnormalities, of which 83% are intragenic mutations and 10% are 5q35 microdeletions. We reviewed the clinical phenotypes of 239 NSD1-positive individuals. Facial dysmorphism, learning disability, and childhood overgrowth were present in 90% of the individuals. However, both the height and head circumference of 10% of the individuals were within the normal range, indicating that overgrowth is not obligatory for the diagnosis of Sotos syndrome. A broad spectrum of associated clinical features was also present, the occurrence of which was largely independent of genotype, since individuals with identical mutations had different phenotypes. We compared the phenotypes of patients with intragenic NSD1 mutations with those of patients with 5q35 microdeletions. Patients with microdeletions had less-prominent overgrowth (P = .0003) and more-severe learning disability (P = 3 x 10(-9)) than patients with mutations. However, all features present in patients with microdeletions were also observed in patients with mutations, and there was no correlation between deletion size and the clinical phenotype, suggesting that the deletion of additional genes in patients with 5q35 microdeletions has little specific effect on phenotype. We identified only 13 familial cases. The reasons for the low vertical transmission rate are unclear, although familial cases were more likely than nonfamilial cases (P = .005) to carry missense mutations, suggesting that the underlying NSD1 mutational mechanism in Sotos syndrome may influence reproductive fitness.

Confirmation of a hereditary motor and sensory neuropathy IIC locus at chromosome 12q23-q24.

Hereditary motor and sensory neuropathy type IIC (HMSN IIC) is an autosomal dominant axonal neuropathy. The cardinal features include distal muscle wasting and weakness, vocal cord paralysis, and mild sensory impairment. Recently, HMSN IIC locus was mapped to chromosome 12q23-24. Two families affected by HMSN IIC were identified and evaluated for linkage to this region. Segregation analysis in both families was consistent with linkage to chromosome 12q23-24. Combined analysis generated a multipoint LOD score of 2.1 at marker D12S1583 and refined the HMSN IIC gene interval to The clinical and molecular findings are discussed.

Constitutional aneuploidy and cancer predisposition caused by biallelic mutations in BUB1B.

Mosaic variegated aneuploidy is a rare recessive condition characterized by growth retardation, microcephaly, childhood cancer and constitutional mosaicism for chromosomal gains and losses. In five families with mosaic variegated aneuploidy, including two with embryonal rhabdomyosarcoma, we identified truncating and missense mutations of BUB1B, which encodes BUBR1, a key protein in the mitotic spindle checkpoint. These data are the first to relate germline mutations in a spindle checkpoint gene with a human disorder and strongly support a causal link between aneuploidy and cancer development.

Mutations in the transcription factor gene SOX18 underlie recessive and dominant forms of hypotrichosis-lymphedema-telangiectasia.

Hereditary lymphedema is a developmental disorder characterized by chronic swelling of the extremities due to dysfunction of the lymphatic vessels. Two responsible genes have been identified: the vascular endothelial growth factor receptor 3 (VEGFR3) gene, implicated in congenital lymphedema, or Milroy disease, and the forkhead-related transcription factor gene FOXC2, causing lymphedema-distichiasis. We describe three families with an unusual association of hypotrichosis, lymphedema, and telangiectasia. Using microsatellite analysis, we first excluded both VEGFR3 and FOXC2 as causative genes; we then considered the murine ragged phenotype, caused by mutations in the Sox18 transcription factor, as a likely counterpart to the human disease, because it presents a combination of hair and cardiovascular anomalies, including symptoms of lymphatic dysfunction. Two of the families were consanguineous; in affected members of these families, we identified homozygous missense mutations in the SOX18 gene, located in 20q13. The two amino acid substitutions, W95R and A104P, affect conserved residues in the first alpha helix of the DNA-binding domain of the transcription factor. In the third family, the parents were nonconsanguineous, and both the affected child and his brother, who died in utero with hydrops fetalis, showed a heterozygous nonsense mutation that truncates the SOX18 protein in its transactivation domain; this substitution was not found in genomic DNA from either parent and hence constitutes a de novo germline mutation. Thus, we show that SOX18 mutations in humans cause both recessive and dominant hypotrichosis-lymphedema-telangiectasia, suggesting that, in addition to its established role in hair and blood vessel development, the SOX18 transcription factor plays a role in the development and/or maintenance of lymphatic vessels.