Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS.

BACKGROUND: Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions. METHODS: We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing. RESULTS: In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects. CONCLUSIONS: Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Molecular Characterization Reveals Subclasses of 1q Gain in Intermediate Risk Wilms Tumors.

Chromosomal alterations have recurrently been identified in Wilms tumors (WTs) and some are associated with poor prognosis. Gain of 1q (1q+) is of special interest given its high prevalence and is currently actively studied for its prognostic value. However, the underlying mutational mechanisms and functional effects remain unknown. In a national unbiased cohort of 30 primary WTs, we integrated somatic SNVs, CNs and SVs with expression data and distinguished four clusters characterized by affected biological processes: muscle differentiation, immune system, kidney development and proliferation. Combined genome-wide CN and SV profiles showed that tumors profoundly differ in both their types of 1q+ and genomic stability and can be grouped into WTs with co-occurring 1p-/1q+, multiple chromosomal gains or CN neutral tumors. We identified 1q+ in eight tumors that differ in mutational mechanisms, subsequent rearrangements and genomic contexts. Moreover, 1q+ tumors were present in all four expression clusters reflecting activation of various biological processes, and individual tumors overexpress different genes on 1q. In conclusion, by integrating CNs, SVs and gene expression, we identified subgroups of 1q+ tumors reflecting differences in the functional effect of 1q gain, indicating that expression data is likely needed for further risk stratification of 1q+ WTs.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing.

PURPOSE: Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events. METHODS: We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing. RESULTS: We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents. CONCLUSION: We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

Amplicon-Based Targeted Next-Generation Sequencing of Formalin-Fixed, Paraffin-Embedded Tissue.

Next-generation sequencing (NGS) is rapidly becoming the method of choice for mutation analysis in both research and diagnostics. The benefit of targeted NGS compared to whole-genome and whole-exome sequencing is that smaller amounts of input material can be used as well as qualitatively suboptimal tissue samples, like formalin-fixed, paraffin-embedded archival tissue.Here, we describe the protocol for targeted next-generation sequencing using the Ion Torrent PGM platform in combination with Ion Ampliseq NGS gene panels for formalin-fixed, paraffin-embedded tissues. Both the manual and the automated workflow are described as well as the bioinformatics for data analysis.

E-cadherin loss induces targetable autocrine activation of growth factor signalling in lobular breast cancer.

Despite the fact that loss of E-cadherin is causal to the development and progression of invasive lobular carcinoma (ILC), options to treat this major breast cancer subtype are limited if tumours develop resistance to anti-oestrogen treatment regimens. This study aimed to identify clinically targetable pathways that are aberrantly active downstream of E-cadherin loss in ILC. Using a combination of reverse-phase protein array (RPPA) analyses, mRNA sequencing, conditioned medium growth assays and CRISPR/Cas9-based knock-out experiments, we demonstrate that E-cadherin loss causes increased responsiveness to autocrine growth factor receptor (GFR)-dependent activation of phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/Akt signalling. Autocrine activation of GFR signalling and its downstream PI3K/Akt hub was independent of oncogenic mutations in PIK3CA, AKT1 or PTEN. Analyses of human ILC samples confirmed growth factor production and pathway activity. Pharmacological inhibition of Akt using AZD5363 or MK2206 resulted in robust inhibition of cell growth and survival of ILC cells, and impeded tumour growth in a mouse ILC model. Because E-cadherin loss evokes hypersensitisation of PI3K/Akt activation independent of oncogenic mutations in this pathway, we propose clinical intervention of PI3K/Akt in ILC based on functional E-cadherin inactivation, irrespective of activating pathway mutations.

Mutational analysis using Sanger and next generation sequencing in sporadic spindle cell hemangiomas: A study of 19 cases.

Spindle cell hemangioma (SCH) is a distinct vascular soft-tissue lesion characterized by cavernous blood vessels and a spindle cell component mainly occurring in the distal extremities of young adults. The majority of cases harbor heterozygous mutations in IDH1/2 sporadically or rarely in association with Maffucci syndrome. However, based on mosaicism and accordingly a low percentage of lesional cells harboring a mutant allele, detection can be challenging. We tested 19 sporadic SCHs by Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), conventional next generation sequencing (NGS), and NGS using a single molecule molecular inversion probes (smMIP)-based library preparation to compare their diagnostic value. Out of 10 cases tested by Sanger sequencing and 2 analyzed using MLPA, 4 and 1, respectively, revealed a mutation in IDH1 (p.R132C). The 7 remaining negative cases and additional 6 cases were investigated using smMIP/NGS, showing hot spot mutations in IDH1 (p.R132C) (8 cases) and IDH2 (3 cases; twice p.R172S and once p.R172G, respectively). One case was negative. Owing to insufficient DNA quality and insufficient coverage, 2 cases were excluded. In total, in 16 out of 17 cases successfully tested, an IDH1/2 mutation was found. Given that IDH1/2 mutations were absent in 161 other vascular lesions tested by smMIP/NGS, the mutation can be considered as highly specific for SCH.

D-amino acid aberrations in cerebrospinal fluid and plasma of smokers.

The glutamatergic neurotransmission system and the N-methyl-D-aspartate receptor (NMDAR) have been implicated in smoking and alcohol consumption behavior. Preclinical studies have demonstrated that nicotine and ethanol influence NMDAR functionality, which may have a role in tendencies to consume these substances. Nonetheless, little is known about concentrations of NMDAR coagonists in the cerebrospinal fluid (CSF) and plasma of individuals who smoke or consume alcohol. Glycine and L- and D-stereoisomers of alanine, serine, and proline were therefore measured using ultra-high-performance liquid chromatography-tandem mass spectrometry in 403 healthy subjects. Nicotine and alcohol consumption were quantified using questionnaires. Possible differences in NMDAR coagonist concentrations in plasma and CSF were investigated using ANCOVA with age, body mass index, and storage duration as covariates. The significance threshold was Bonferroni corrected (α=0.00625). Compared with non-smokers, smokers displayed lower levels of D-proline in plasma (p=0.0027, Cohen's d=-0.41) and D-proline in CSF (p=0.0026, Cohen's d=-0.43). D-Serine in CSF was higher in smokers than in non-smokers (p=0.0052, Cohen's d=0.41). After subdividing participants based on smoking quantity, dose-dependent decreases were demonstrated in smokers for D-proline in plasma (F=5.65, p=0.0039) and D-proline in CSF (F=5.20, p=0.0060). No differences in NMDAR coagonist levels between alcohol consumption groups were detected. To our knowledge, this is the first report to implicate D-amino acids in smoking behavior of humans. Whether such concentration differences lie at the root of or result from smoking habits may be addressed in prospective studies.

A common variant in ERBB4 regulates GABA concentrations in human cerebrospinal fluid.

The neuregulin 1 (NRG1) receptor ErbB4 is involved in the development of cortical inhibitory GABAergic circuits and NRG1-ErbB4 signaling has been implicated in schizophrenia (SCZ). A magnetic resonance spectroscopy ((1)H-MRS) study has demonstrated that a single-nucleotide polymorphism in ERBB4, rs7598440, influences human cortical GABA concentrations. Other work has highlighted the significant impact of this genetic variant on expression of ERBB4 in the hippocampus and dorsolateral prefrontal cortex in human post mortem tissue. Our aim was to examine the association of rs7598440 with cerebrospinal fluid (CSF) GABA levels in healthy volunteers (n=155). We detected a significant dose-dependent association of the rs7598440 genotype with CSF GABA levels (G-allele standardized ß=-0.23; 95% CIs: -0.39 to -0.07; P=0.0066). GABA concentrations were highest in A homozygous, intermediate in heterozygous, and lowest in G homozygous subjects. When excluding subjects on psychotropic medication (three subjects using antidepressants), the results did not change (G-allele standardized ß=-0.23; 95% CIs: -0.40 to -0.07; P=0.0051). The explained variance in CSF GABA by rs7598440 in our model is 5.2% (P=0.004). The directionality of our findings agrees with the aforementioned (1)H-MRS and gene expression studies. Our observation therefore strengthens the evidence that the A-allele of rs7598440 in ERBB4 is associated with increased GABA concentrations in the human central nervous system (CNS). To our knowledge, our finding constitutes the first confirmation that CSF can be used to study genotype-phenotype correlations of GABA levels in the CNS. Such quantitative genetic analyses may be extrapolated to other CSF constituents relevant to SCZ in future studies.

The association of the alpha-5 subunit of the nicotinic acetylcholine receptor gene and the brain-derived neurotrophic factor gene with different aspects of smoking behavior.

Recent studies show that different aspects of smoking behavior are associated with the α-5 subunit of the nicotinic acetylcholine receptor (CHRNA5) gene and the gene coding for brain-derived neurotrophic factor (BDNF). This raises the question whether the amount of cigarettes smoked per day has a different genetic background than smoking initiation and what other smoking phenotypes may be relevant. The aim of this study was to replicate these associations in a large population-based sample. We investigated the association with smoking initiation and the number of cigarettes used per day and additional smoking phenotypes in a population-based sample of 2166 participants of Dutch origin. Rs6265 in BDNF was not associated with smoking initiation. This single nucleotide polymorphism was associated with smoking cessation. Rs16969968 in CHRNA5 was associated with the amount of nicotine used and in particular smoking 25 cigarettes or more per day. Overall, the results confirm the involvement of the CHRNA5 gene in the amount of nicotine use and further suggest involvement of the BDNF gene in smoking behavior.

Neutrophil recruitment and barrier impairment in celiac disease: a genomic study.

BACKGROUND & AIMS: Celiac disease is an enteropathy featuring villous atrophy, crypt hyperplasia, and lymphocytosis. Tissue remodeling is driven by an inflammatory reaction to gluten in genetically susceptible individuals. The adaptive pathway is considered the major immune response but recent evidence has indicated the involvement of innate immunity as well. To assess the contribution of either immune response we performed global gene expression profiling of the regenerating mucosa. METHODS: Microarray hybridizations were performed with biopsy samples from 13 untreated patients, 31 patients on a gluten-free diet in various stages of remission, and 21 controls. Additional data were generated using low-density array and conventional quantitative reverse-transcription polymerase chain reaction, and immunohistochemistry. RESULTS: A total of 108 differentially expressed immune-related genes were identified (50 innate, 43 adaptive, 9 both innate/adaptive, and 6 immunoregulatory). Expression levels showed a gradual change as opposed to the discrete histological transitions. In addition to details provided on the adaptive and innate immune pathways used, we observed a chronic recruitment of activated neutrophils. Neutrophil involvement was unabated in otherwise completely normalized remission patients. CONCLUSIONS: We observed a contribution of both the innate and adaptive immune response in celiac disease pathogenesis. The discrepancy between the histological classification and the observed incremental change in immune-gene expression may have consequences for current diagnostic inclusion criteria. Enhanced neutrophil infiltration in both active and remission patients points to a genetic impairment of the intestinal barrier that may contribute to the cause rather than the consequence of celiac disease.

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect.

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.

Gene expression profiling and phenotype analyses of S. cerevisiae in response to changing copper reveals six genes with new roles in copper and iron metabolism.

Exhaustive microarray time course analyses of Saccharomyces cerevisiae during copper starvation and copper excess reveal new aspects of metal-induced gene regulation. Aside from identifying targets of established copper- and iron-responsive transcription factors, we find that genes encoding mitochondrial proteins are downregulated and that copper-independent iron transport genes are preferentially upregulated, both during prolonged copper deprivation. The experiments also suggest the presence of a small regulatory iron pool that links copper and iron responses. One hundred twenty-eight genes with putative roles in metal metabolism were further investigated by several systematic phenotype screens. Of the novel phenotypes uncovered, hsp12-Delta and arn1-Delta display increased sensitivity to copper, cyc1-Delta and crr1-Delta show resistance to high copper, vma13-Delta exhibits increased sensitivity to iron deprivation, and pep12-Delta results in reduced growth in high copper and low iron. Besides revealing new components of eukaryotic metal trafficking pathways, the results underscore the previously determined intimate links between iron and copper metabolism and mitochondrial and vacuolar function in metal trafficking. The analyses further suggest that copper starvation can specifically lead to downregulation of respiratory function to preserve iron and copper for other cellular processes.