The genomic profiling and MAMLD1 expression in human and canines with Cushing's disease.

BACKGROUND: Cushing's disease (CD) is defined as hypercortisolemia caused by adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (corticotroph PA) that afflicts humans and dogs. In order to map common aberrant genomic features of CD between humans and dogs, we performed genomic sequencing and immunostaining on corticotroph PA. METHODS: For inclusion, humans and dog were diagnosed with CD. Whole exome sequencing (WES) was conducted on 6 human corticotroph PA. Transcriptome RNA-Seq was performed on 6 human and 7 dog corticotroph PA. Immunohistochemistry (IHC) was complete on 31 human corticotroph PA. Corticotroph PA were compared with normal tissue and between species analysis were also performed. RESULTS: Eight genes (MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, FAM131B) were significantly (P < 0.05) overexpressed across human and canine corticotroph PA. IHC revealed MAMLD1 to be positively (3+) expressed in the nucleus of ACTH-secreting tumor cells of human corticotroph PA (22/31, 70.9%), but absent in healthy human pituitary glands. CONCLUSIONS: In this small exploratory cohort, we provide the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.

Gene integrated set profile analysis: a context-based approach for inferring biological endpoints.

The identification of genes with specific patterns of change (e.g. down-regulated and methylated) as phenotype drivers or samples with similar profiles for a given gene set as drivers of clinical outcome, requires the integration of several genomic data types for which an 'integrate by intersection' (IBI) approach is often applied. In this approach, results from separate analyses of each data type are intersected, which has the limitation of a smaller intersection with more data types. We introduce a new method, GISPA (Gene Integrated Set Profile Analysis) for integrated genomic analysis and its variation, SISPA (Sample Integrated Set Profile Analysis) for defining respective genes and samples with the context of similar, a priori specified molecular profiles. With GISPA, the user defines a molecular profile that is compared among several classes and obtains ranked gene sets that satisfy the profile as drivers of each class. With SISPA, the user defines a gene set that satisfies a profile and obtains sample groups of profile activity. Our results from applying GISPA to human multiple myeloma (MM) cell lines contained genes of known profiles and importance, along with several novel targets, and their further SISPA application to MM coMMpass trial data showed clinical relevance.

Next-generation sequencing and clinical outcomes of patients with lung adenocarcinoma treated with stereotactic body radiotherapy.

BACKGROUND: Genetic aberrations are well characterized in lung adenocarcinomas (LACs) and clinical outcomes have been influenced by targeted therapies in the advanced setting. Stereotactic body radiotherapy (SBRT) is the standard-of-care therapy for patients with nonoperable, early-stage LAC, but to the authors' knowledge, no information is available regarding the impact of genomic changes in these patients. The current study sought to determine the frequency and clinical impact of genetic aberrations in this population. METHODS: Under an Institutional Review Board-approved protocol, the records of 242 consecutive patients with early-stage lung cancers were reviewed; inclusion criteria included LAC histology with an adequate tumor sample for the successful use of next-generation sequencing and fluorescence in situ hybridization testing. Univariate analysis was performed to identify factors associated with clinical outcomes. RESULTS: LAC samples from 98 of the 242 patients were reviewed (40.5%), of whom 45 patients (46.0%) had genetic testing. The following mutations were noted: KRAS in 20.0% of samples, BRAF in 2.2% of samples, SMAD family member 4 (SMAD4) in 4.4% of samples, epidermal growth factor receptor (EGFR) in 15.6% of samples, STK1 in 2.2% of samples, tumor protein 53 (TP53) in 15.6% of samples, and phosphatase and tensin homolog (PTEN) in 2.2% of samples. The following gene rearrangements were observed: anaplastic lymphoma kinase (ALK) in 8.9% of samples, RET in 2.2% of samples, and MET amplification in 17.8% of samples. The median total delivered SBRT dose was 50 grays (range, 48-60 grays) over a median of 5 fractions (range, 3-8 fractions). The KRAS mutation was associated with worse local control (odds ratio [OR], 3.64; P<.05). MET amplification was associated with worse regional (OR, 4.64; P<.05) and distant (OR, 3.73; P<.05) disease control. CONCLUSIONS: To the authors' knowledge, the current series is the first to quantify genetic mutations and their association with clinical outcomes in patients with early-stage LAC treated with SBRT. KRAS mutations were associated with worse local control and MET amplification was associated with worse regional and distant disease control, findings that need to be validated in a prospective setting. Cancer 2017;123:3681-3690. © 2017 American Cancer Society.

Patient-derived xenografts faithfully replicated clinical outcome in a phase II co-clinical trial of arsenic trioxide in relapsed small cell lung cancer.

BACKGROUND: SCLC has limited treatment options and inadequate preclinical models. Promising activity of arsenic trioxide (ASO) recorded in conventional preclinical models of SCLC supported the clinical evaluation of ASO in patients. We assessed the efficacy of ASO in relapsed SCLC patients and in corresponding patient-derived xenografts (PDX). METHODS: Single arm, Simon 2-stage, phase II trial to enroll patients with relapsed SCLC who have failed at least one line of therapy. ASO was administered as an intravenous infusion over 1-2 h daily for 4 days in week 1 and for 2 days in weeks 2-6 of an 8-week cycle. Treatment continued until disease progression. Pretreatment tumor biopsy was employed for PDX generation through direct implantation into subcutaneous pockets of SCID mice without in vitro manipulation and serially propagated for five generations. Ex vivo efficacy of cisplatin (3 mg/kg i.p. weekly) and ASO (3.75 mg/kg i.p. every other day) was tested in PDX representative of platinum sensitive and platinum refractory SCLC. RESULTS: The best response in 17 evaluable patients was stable disease in 2 (12 %), progressive disease in 15 (88 %) patients and median time-to-progression of seven (range 1-7) weeks. PDX was successfully grown in 5 of 9 (56 %) transplanted biopsy samples. Serially-propagated PDXs preserved characteristic small cell histology and genomic stability confirmed by immunohistochemistry, short tandem repeat (STR) profiling and targeted sequencing. ASO showed in vitro cytotoxicity but lacked in vivo efficacy against SCLC PDX tumor growth. CONCLUSIONS: Cisplatin inhibited growth of PDX derived from platinum-sensitive SCLC but was ineffective against PDX from platinum-refractory SCLC. Strong concordance between clinical and ex vivo effects of ASO and cisplatin in SCLC supports the use of PDX models to prescreen promising anticancer agents prior to clinical testing in SCLC patients. Trial Registration The study was registered at http://www.clinicaltrials.gov (NCT01470248).

A new method for temperature-field reconstruction during ultrasound-monitored cryosurgery using potential-field analogy.

The current study aims at developing computational tools in order to gain information about the thermal history in areas invisible to ultrasound imaging during cryosurgery. This invisibility results from the high absorption rate of the ultrasound energy by the frozen region, which leads to an apparent opacity in the cryotreated area and a shadow behind it. A proof-of-concept for freezing-front estimation is demonstrated in the current study, using the new potential-field analogy method (PFAM). This method is further integrated with a recently developed temperature-field reconstruction method (TFRM) to estimate the temperature distribution within the frozen region. This study uses prostate cryosurgery as a developmental model and trans-rectal ultrasound imaging as a choice of practice. Results of this study indicate that the proposed PFAM is a viable and computationally inexpensive solution to estimate the extent of freezing in the acoustic shadow region. Comparison of PFAM estimations and experimental data shows an average mismatch of less than 2 mm in freezing-front location, which is comparable to the uncertainty in ultrasound imaging. Comparison of the integrated PFAM + TFRM scheme with a full-scale finite-elements analysis (FEA) indicates an average mismatch of 0.9 mm for the freezing front location and 0.1 mm for the lethal temperature isotherm of -45 °C. Comparison of the integrated PFAM + TFRM scheme with experimental temperature measurements show a difference in the range of 2 °C and 6 °C for selected points of measurement. Results of this study demonstrate the integrated PFAM + TFRM scheme as a viable and computationally inexpensive means to gain information about the thermal history in the frozen region during ultrasound-monitored cryosurgery.

Gene expression profiling of clear cell papillary renal cell carcinoma: comparison with clear cell renal cell carcinoma and papillary renal cell carcinoma.

Clear cell papillary renal cell carcinoma is a distinct variant of renal cell carcinoma that shares some overlapping histological and immunohistochemical features of clear cell renal cell carcinoma and papillary renal cell carcinoma. Although the clear cell papillary renal cell carcinoma immunohistochemical profile is well described, clear cell papillary renal cell carcinoma mRNA expression has not been well characterized. We investigated the clear cell papillary renal cell carcinoma gene expression profile using previously identified candidate genes. We selected 17 clear cell papillary renal cell carcinoma, 15 clear cell renal cell carcinoma, and 13 papillary renal cell carcinoma cases for molecular analysis following histological review. cDNA from formalin-fixed paraffin-embedded tissue was prepared. Quantitative real-time PCR targeting alpha-methylacyl coenzyme-A racemase (AMACR), BMP and activin membrane-bound inhibitor homolog (BAMBI), carbonic anhydrase IX (CA9), ceruloplasmin (CP), nicotinamide N-methyltransferase (NNMT), schwannomin-interacting protein 1 (SCHIP1), solute carrier family 34 (sodium phosphate) member 2 (SLC34A2), and vimentin (VIM) was performed. Gene expression data were normalized relative to 28S ribosomal RNA. Clear cell papillary renal cell carcinoma expressed all eight genes at variable levels. Compared with papillary renal cell carcinoma, clear cell papillary renal cell carcinoma expressed more CA9, CP, NNMT, and VIM, less AMACR, BAMBI, and SLC34A2, and similar levels of SCHIP1. Compared with clear cell renal cell carcinoma, clear cell papillary renal cell carcinoma expressed slightly less NNMT, but similar levels of the other seven genes. Although clear cell papillary renal cell carcinoma exhibits a unique molecular signature, it expresses several genes at comparable levels to clear cell renal cell carcinoma relative to papillary renal cell carcinoma. Understanding the molecular pathogenesis of clear cell papillary renal cell carcinoma will have a key role in future sub-classifications of this unique tumor.

Characterization of disease-propagating stem cells responsible for myeloproliferative neoplasm-blast phase.

Chronic myeloproliferative neoplasms (MPN) frequently evolve to a blast phase (BP) that is almost uniformly resistant to induction chemotherapy or hypomethylating agents. We explored the functional properties, genomic architecture, and cell of origin of MPN-BP initiating cells (IC) using a serial NSG mouse xenograft transplantation model. Transplantation of peripheral blood mononuclear cells (MNC) from 7 of 18 patients resulted in a high degree of leukemic cell chimerism and recreated clinical characteristics of human MPN-BP. The function of MPN-BP ICs was not dependent on the presence of JAK2V617F, a driver mutation associated with the initial underlying MPN. By contrast, multiple MPN-BP IC subclones coexisted within MPN-BP MNCs characterized by different myeloid malignancy gene mutations and cytogenetic abnormalities. MPN-BP ICs in 4 patients exhibited extensive proliferative and self-renewal capacity, as demonstrated by their ability to recapitulate human MPN-BP in serial recipients. These MPN-BP IC subclones underwent extensive continuous clonal competition within individual xenografts and across multiple generations, and their subclonal dynamics were consistent with functional evolution of MPN-BP IC. Finally, we show that MPN-BP ICs originate from not only phenotypically identified hematopoietic stem cells, but also lymphoid-myeloid progenitor cells, which were each characterized by differences in MPN-BP initiating activity and self-renewal capacity.

Heterogeneity of primary glioblastoma cells in the expression of caspase-8 and the response to TRAIL-induced apoptosis.

Recent studies suggest that cancer stem cells (CSCs) are responsible for cancer resistance to therapies. We therefore investigated how glioblastoma-derived CSCs respond to the treatment of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL). Neurospheres were generated from glioblastomas, characterized for CSC properties including self-renewal, cell differentiation and xenograft formation capacity, and analyzed for TRAIL-induced apoptosis, CASP8 genomic status, and caspase-8 protein expression. The neurosphere NSC326 was sensitive to TRAIL-induced apoptosis as evidenced by cell death and caspase-8, -3, and -7 enzymatic activities. In contrast, however, the neurosphere NSC189 was TRAIL-resistant. G-banding analysis identified five chromosomally distinguishable cell populations in the neurospheres. Fluorescence in situ hybridization revealed the variation of chromosome 2 copy number in these populations and the loss of CASP8 locus in 2q33-34 region in a small set of cell populations in the neurosphere. Immunohistochemistry of NSC189 cell blocks revealed the lack of caspase-8 protein in a subset of neurosphere cells. Western blotting and immunohistochemistry of human glioblastoma tumors demonstrated the expression of caspase-8 protein in the vast majority of the tumors as compared to normal human brain tissues that lack the caspase-8 expression. This study shows heterogeneity of glioblastomas and derived CSCs in the genomic status of CASP8, expression of caspase-8, and thus responsiveness to TRAIL-induced apoptosis. Clinic trials may consider genomic analysis of the cancer tissue to identify the genomic loss of CASP8 and use it as a genomic marker to predict the resistance of glioblastomas to TRAIL apoptosis pathway-targeted therapies.

An evidence-based approach to establish the functional and clinical significance of copy number variants in intellectual and developmental disabilities.

PURPOSE: Copy number variants have emerged as a major cause of human disease such as autism and intellectual disabilities. Because copy number variants are common in normal individuals, determining the functional and clinical significance of rare copy number variants in patients remains challenging. The adoption of whole-genome chromosomal microarray analysis as a first-tier diagnostic test for individuals with unexplained developmental disabilities provides a unique opportunity to obtain large copy number variant datasets generated through routine patient care. METHODS: A consortium of diagnostic laboratories was established (the International Standards for Cytogenomic Arrays consortium) to share copy number variant and phenotypic data in a central, public database. We present the largest copy number variant case-control study to date comprising 15,749 International Standards for Cytogenomic Arrays cases and 10,118 published controls, focusing our initial analysis on recurrent deletions and duplications involving 14 copy number variant regions. RESULTS: Compared with controls, 14 deletions and seven duplications were significantly overrepresented in cases, providing a clinical diagnosis as pathogenic. CONCLUSION: Given the rapid expansion of clinical chromosomal microarray analysis testing, very large datasets will be available to determine the functional significance of increasingly rare copy number variants. This data will provide an evidence-based guide to clinicians across many disciplines involved in the diagnosis, management, and care of these patients and their families.

Genomic characterization of prenatally detected chromosomal structural abnormalities using oligonucleotide array comparative genomic hybridization.

Detection of chromosomal structural abnormalities using conventional cytogenetic methods poses a challenge for prenatal genetic counseling due to unpredictable clinical outcomes and risk of recurrence. Of the 1,726 prenatal cases in a 3-year period, we performed oligonucleotide array comparative genomic hybridization (aCGH) analysis on 11 cases detected with various structural chromosomal abnormalities. In nine cases, genomic aberrations and gene contents involving a 3p distal deletion, a marker chromosome from chromosome 4, a derivative chromosome 5 from a 5p/7q translocation, a de novo distal 6q deletion, a recombinant chromosome 8 comprised of an 8p duplication and an 8q deletion, an extra derivative chromosome 9 from an 8p/9q translocation, mosaicism for chromosome 12q with added material of initially unknown origin, an unbalanced 13q/15q rearrangement, and a distal 18q duplication and deletion were delineated. An absence of pathogenic copy number changes was noted in one case with a de novo 11q/14q translocation and in another with a familial insertion of 21q into a 19q. Genomic characterization of the structural abnormalities aided in the prediction of clinical outcomes. These results demonstrated the value of aCGH analysis in prenatal cases with subtle or complex chromosomal rearrangements. Furthermore, a retrospective analysis of clinical indications of our prenatal cases showed that approximately 20% of them had abnormal ultrasound findings and should be considered as high risk pregnancies for a combined chromosome and aCGH analysis.

Computerized planning of prostate cryosurgery using variable cryoprobe insertion depth.

The current study presents a computerized planning scheme for prostate cryosurgery using a variable insertion depth strategy. This study is a part of an ongoing effort to develop computerized tools for cryosurgery. Based on typical clinical practices, previous automated planning schemes have required that all cryoprobes be aligned at a single insertion depth. The current study investigates the benefit of removing this constraint, in comparison with results based on uniform insertion depth planning as well as the so-called "pullback procedure". Planning is based on the so-called "bubble-packing method", and its quality is evaluated with bioheat transfer simulations. This study is based on five 3D prostate models, reconstructed from ultrasound imaging, and cryoprobe active length in the range of 15-35 mm. The variable insertion depth technique is found to consistently provide superior results when compared to the other placement methods. Furthermore, it is shown that both the optimal active length and the optimal number of cryoprobes vary among prostate models, based on the size and shape of the target region. Due to its low computational cost, the new scheme can be used to determine the optimal cryoprobe layout for a given prostate model in real time.

Clinical and genomic characterization of distal duplications and deletions of chromosome 4q: study of two cases and review of the literature.

Variable clinical presentations of patients with chromosomally detected deletions in the distal long arm (q) of chromosome 4 have been reported. The lack of molecular characterization of the deletion sizes and deleted genes hinders further genotype-phenotype correlation. Using a validated oligonucleotide array comparative genomic hybridization (oaCGH) analysis, we examined two patients with apparent chromosomal deletions in the distal 4q region. In the first, oaCGH identified a 2.441 megabase (Mb) duplication and a 12.651 Mb deletion at 4q34.1 in a pregnant female who transmitted this aberration to her son. This mother has only learning disabilities while her son had both renal and cardiac anomalies in the newborn period. Unrecognized paternal genetic factors may contribute to the variable expression. The second patient is a 17-year-old female with a history of Pierre Robin sequence, cardiac abnormalities and learning disabilities. She was diagnosed prenatally with a de novo 4q deletion, and oaCGH defined a 16.435 Mb deletion of 4q34.1-4q35.2. Phenotypic comparison and subtractive genomic mapping between these two cases suggested a 4 Mb region possibly harboring a candidate gene for Pierre Robin sequence. Our cases and review of reported cases with genomic findings indicated the presence of familial variants with variable expressivity as well as de novo or inherited pathogenic simple deletion, duplication and complex deletion and duplication in the distal 4q region.

Critical role of ASCT2-mediated amino acid metabolism in promoting leukaemia development and progression.

Amino acid (AA) metabolism is involved in diverse cellular functions, including cell survival and growth, however it remains unclear how it regulates normal hematopoiesis versus leukemogenesis. Here, we report that knockout of Slc1a5 (ASCT2), a transporter of neutral AAs, especially glutamine, results in mild to moderate defects in bone marrow and mature blood cell development under steady state conditions. In contrast, constitutive or induced deletion of Slc1a5 decreases leukemia initiation and maintenance driven by the oncogene MLL-AF9 or Pten deficiency. Survival of leukemic mice is prolonged following Slc1a5 deletion, and pharmacological inhibition of ASCT2 also decreases leukemia development and progression in xenograft models of human acute myeloid leukemia. Mechanistically, loss of ASCT2 generates a global effect on cellular metabolism, disrupts leucine influx and mTOR signaling, and induces apoptosis in leukemic cells. Given the substantial difference in reliance on ASCT2-mediated AA metabolism between normal and malignant blood cells, this in vivo study suggests ASCT2 as a promising therapeutic target for the treatment of leukemia.

Genomic analysis of CD8+ NK/T cell line, 'SRIK-NKL', with array-based CGH (aCGH), SKY/FISH and molecular mapping.

We performed aCGH, SKY/FISH, molecular mapping and expression analyses on a permanent CD8+ NK/T cell line, 'SRIK-NKL' established from a lymphoma (ALL) patient, in attempt to define the fundamental genetic profile of its unique NK phenotypes. aCGH revealed hemizygous deletion of 6p containing genes responsible for hematopoietic functions. The SKY demonstrated that a constitutive reciprocal translocation, rcpt(5;14)(p13.2;q11) is a stable marker. Using somatic hybrids containing der(5) derived from SRIK-NKL, we found that the breakpoint in one homologue of no. 5 is located upstream of IL7R and also that the breakpoint in no. 14 is located within TRA@. The FISH analysis using a BAC which contains TRA@ and its flanking region further revealed a approximately 231kb deletion within 14q11 in the der(5) but not in the normal homologue of no. 14. The RT-PCR analysis detected mRNA for TRA@ transcripts which were extending across, but not including, the deleted region. IL7R was detected at least at mRNA levels. These findings were consistent with the immunological findings that TRA@ and IL7R are both expressed at mRNA levels and TRA@ at cytoplasmic protein levels in SRIK-NKL cells. In addition to rcpt(5;14), aCGH identified novel copy number abnormalities suggesting that the unique phenotype of the SRIK-NKL cell line is not solely due to the TRA@ rearrangement. These findings provide supportive evidence for the notion that SRIK-NKL cells may be useful for studying not only the function of NK cells but also genetic deregulations associated with leukemiogenesis.

Computerized Planning of Cryosurgery Using Bubble Packing: An Experimental Validation on a Phantom Material.

The current study focuses on experimentally validating a planning scheme based on the so-called bubble-packing method. This study is a part of an ongoing effort to develop computerized planning tools for cryosurgery, where bubble packing has been previously developed as a means to find an initial, uniform distribution of cryoprobes within a given domain; the so-called force-field analogy was then used to move cryoprobes to their optimum layout. However, due to the high quality of the cryoprobes' distribution, suggested by bubble packing and its low computational cost, it has been argued that a planning scheme based solely on bubble packing may be more clinically relevant. To test this argument, an experimental validation is performed on a simulated cross-section of the prostate, using gelatin solution as a phantom material, proprietary liquid-nitrogen based cryoprobes, and a cryoheater to simulate urethral warming. Experimental results are compared with numerically simulated temperature histories resulting from planning. Results indicate an average disagreement of 0.8 mm in identifying the freezing front location, which is an acceptable level of uncertainty in the context of prostate cryosurgery imaging.

Cryosurgery planning using bubble packing in 3D.

As part of an ongoing project to develop automated tools for cryosurgery planning, the current study focuses on the development of a 3D bubble packing method. A proof-of-concept for the new method is demonstrated on five prostate models, reconstructed from ultrasound images. The new method is a modification of an established method in 2D. Ellipsoidal bubbles are packed in the volume of the prostate in the current study; such bubbles can be viewed as a first-order approximation of a frozen region around a single cryoprobe. When all cryoprobes are inserted to the same depth, optimum planning was found to occur at about 60% of the length of the prostate (measured from its apex), which leads to cooling of approximately 75% of the prostate volume below a specific temperature threshold of - 22 degrees C. Bubble packing has the potential to dramatically reduce the run time for automated planning.

Colon carcinoma cells harboring PIK3CA mutations display resistance to growth factor deprivation induced apoptosis.

PIK3CA, encoding the p110alpha catalytic subunit of phosphatidylinositol 3-kinase (PI3K), is mutated in a variety of human cancers. We screened the colon cancer cell lines previously established in our laboratory for PIK3CA mutations and found that four of them harbored gain of function mutations. We have now compared a panel of mutant and wild-type cell lines for cell proliferation and survival in response to stress. There was little difference in PI3K activity between mutant PIK3CA-bearing cells (mutant cells) and wild-type PIK3CA-bearing cells (wild-type cells) under optimal growth conditions. However, the mutant cells showed constitutive PI3K activity during growth factor deprivation stress (GFDS), whereas PI3K activity decayed rapidly in the wild-type cells. Importantly, constitutively active PI3K rendered the mutant cells resistant to GFDS-induced apoptosis relative to the wild-type cells, indicating a biological advantage under stress conditions that is imparted by the mutant enzymes. Compared with the wild-type cells, the mutant cells were hypersensitive to the apoptosis induced by the PI3K inhibitor LY294002. In addition, PIK3CA small interfering RNA significantly decreased DNA synthesis and/or induced apoptosis in the mutant cells but not in the wild-type cells. Furthermore, ecotopic expression of a mutant PIK3CA in a nontumorigenic PIK3CA wild-type cell line resulted in resistance to GFDS-induced apoptosis, whereas transfection of wild-type PIK3CA or empty vector had little effect. Taken together, our studies show that mutant PIK3CA increases the capacity for proliferation and survival under environmental stresses, such as GFDS while also imparting greater dependency on the PI3K pathway for proliferation and survival.

Validation of a Customized Bioinformatics Pipeline for a Clinical Next-Generation Sequencing Test Targeting Solid Tumor-Associated Variants.

Bioinformatic analysis is an integral and critical part of clinical next-generation sequencing. It is especially challenging for some pipelines to consistently identify insertions and deletions. We present the validation of an open source tumor amplicon pipeline (OTA-pipeline) for clinical next-generation sequencing targeting solid tumor-associated variants. Raw data generated from 557 TruSight Tumor 26 samples and in silico data were analyzed by the OTA-pipeline and legacy pipeline and compared. Discrepant results were confirmed by orthogonal methods. The OTA-pipeline reported 22 variants that were not detected by the previously validated pipeline, including seven synonymous or intronic single-nucleotide variants, five single-nucleotide variants at frequency <5%, one insertion, and nine deletions. Variant allele frequencies reported by the two pipelines were highly concordant, although a few significant discrepancies were present. Analysis of in silico FASTQ files demonstrated a higher sensitivity of detecting complex insertions and deletions with the OTA-pipeline. The higher sensitivity came at a cost, because false-positive calls were increased in difficult-to-sequence regions. However, these calls were all flagged by our strand bias filter, distinguishing them from true variants. Our validation process provides a model for laboratories that want to establish an in-house bioinformatics pipeline for clinical next-generation sequencing.

Common cancer-driver mutations and their association with abnormally methylated genes in lung adenocarcinoma from never-smokers.

OBJECTIVES: Lung adenocarcinoma in never-smokers accounts for 15-20% of all lung cancer. Although targetable mutations are more prevalent in these tumors, the biological and clinical importance of coexisting and/or mutually exclusive abnormalities is just emerging. This study evaluates the relationships between common genetic and epigenetic aberrations in these tumors. MATERIALS AND METHODS: Next-generation sequencing was employed to screen 20 commonly mutated cancer-driver genes in 112 lung adenocarcinomas from never-smokers. The relationship of these mutations with cancer-related methylation of 59 genes, and geographical/ethnic differences in the prevalence for mutations compared to multiple East Asian never-smoker lung adenocarcinoma cohorts was studied. RESULTS: The most common driver mutation detected in 40% (45/112) of the tumors was EGFR, followed by TP53 (18%), SETD2 (11%), and SMARCA4 (11%). Over 72% (81/112) of the cases have mutation of at least one driver gene. While 30% (34/112) of the tumors have co-mutations of two or more genes, 42% (47/112) have only one driver gene mutation. Differences in the prevalence for some of these mutations were seen between adenocarcinomas in East Asian versus US (mainly Caucasian) never-smokers including a significantly lower rate of EGFR mutation among the US patients. Interestingly, aberrant methylation of multiple cancer-related genes was significantly associated with EGFR wildtype tumors. Among 15 differentially methylated genes by EGFR mutation, 14 were more commonly methylated in EGFR wildtype compared to mutant tumors. These findings were independently validated using publicly available data. CONCLUSION: Most lung adenocarcinomas from never-smokers harbor targetable mutation/co-mutations. In the absence of EGFR mutation that drives 40% of these tumors, EGFR wildtype tumors appear to develop by acquiring aberrant promoter methylation that silences tumor-suppressor genes.

Multi-faceted computational assessment of risk and progression in oligodendroglioma implicates NOTCH and PI3K pathways.

Oligodendrogliomas are diffusely infiltrative gliomas defined by IDH-mutation and co-deletion of 1p/19q. They have highly variable clinical courses, with survivals ranging from 6 months to over 20 years, but little is known regarding the pathways involved with their progression or optimal markers for stratifying risk. We utilized machine-learning approaches with genomic data from The Cancer Genome Atlas to objectively identify molecular factors associated with clinical outcomes of oligodendroglioma and extended these findings to study signaling pathways implicated in oncogenesis and clinical endpoints associated with glioma progression. Our multi-faceted computational approach uncovered key genetic alterations associated with disease progression and shorter survival in oligodendroglioma and specifically identified Notch pathway inactivation and PI3K pathway activation as the most strongly associated with MRI and pathology findings of advanced disease and poor clinical outcome. Our findings that Notch pathway inactivation and PI3K pathway activation are associated with advanced disease and survival risk will pave the way for clinically relevant markers of disease progression and therapeutic targets to improve clinical outcomes. Furthermore, our approach demonstrates the strength of machine learning and computational methods for identifying genetic events critical to disease progression in the era of big data and precision medicine.