Exome Sequencing of the t(4;14) and t(11;14) Translocation Specific Subgroups of MM

Blood ◽  
2011 ◽  
Vol 118 (21) ◽  
pp. 1817-1817
Author(s):  
Christopher P Wardell ◽  
Brian A Walker ◽  
David Johnson ◽  
Iwanka Kozarewa ◽  
Kerry Fenwick ◽  
...  

Abstract Abstract 1817 The two most frequent etiological translocations in multiple myeloma (MM) are t(4;14), which deregulates FGFR3 and MMSET and has a poor outcome and t(11;14) which directly deregulates cyclin D1 and has an indolent course. The t(11;14) is present at 10–15% in both MGUS and MM but the t(4;14) is in only 3–4% of MGUS compared to 11% in MM. Consequently it is thought that patients with a t(4;14) have a less stable disease which progresses more quickly to myeloma than other subtypes. In order to address the hypothesis that cases with the t(4;14) are more prone to acquire mutations and so progress, we have compared the number and mutational spectrum of cases with these two variants. DNA was extracted from CD138-selected bone marrow cells from 10 t(4;14) and 12 t(11;14) cases of newly presenting MM. 50 ng of genomic DNA was used to capture the exome using the SureSelect Human All Exon 50Mb target enrichment set (Agilent). We have previously validated this approach and shown it to have parity with approaches using larger starting amounts of DNA. Libraries were prepared from tumor and non-tumor DNA from the same patient and sequenced using 76 bp paired end reads on a GAIIx (Illumina). Samples were sequenced to a median depth of 61x, with 99% >1x and 85% >20x exomic coverage. Following base calling and quality control metrics the raw fastq reads were aligned to the reference human genome. The Genome Analysis Tool Kit was used to call indels and single nucleotide variants (SNVs), with BreakDancer used to detect structural variants. These variant calls were recalibrated and soft filters applied to remove potential false-positives using dbSNP, HapMap and the thousand genomes project as truth sets. Variants that occurred in both the normal and tumor samples were filtered out and the tumor-specific variants were annotated using Annovar. As well as the identification of commonly affected genes, functional annotation enrichment analysis was used to identify commonly affected pathways. The group of 22 cases sequenced at the exome level showed a mutation spectrum that comprised 32,000 SNVs and 1,800 indels per patient, with 1,600 SNVs and 500 indels in the tumor sample only. Structural and copy number variants inferred from this data were also identified and validated previous results using other technologies. We identified 250 SNVs and indels, per patient, that were not in dbSNP and constitute tumor-acquired mutations. We were able to validate some of these mutations that we had previously analyzed using other platforms (98% concordance). The Ti/Tv ratio of mutations was not consistent with any specific exposure or mechanism. The distribution of indels was biased toward insertions rather than deletions, with both predominantly being multiples of three to produce in-frame mutants. In total sequencing data from 60 exomes is available and pathway analysis of the SNVs newly identified confirmed the deregulation of pathways previously identified as being mutated in myeloma, in addition we also identify novel deregulated genes and pathways. We note a consistent increase in the number of variants in the t(4;14). Each tumor had on average 60 non-synonymous SNVs per sample with a range of 29 to 101, some patients being clear outliers. There was a bias to an excess number of mutations within the t(4;14) group which did not reach statistical significance. Importantly, the overlap between the SNVs in individual patients was limited with few consistently mutated genes across the sample set as a whole. In contrast, pathway analysis of the genes mutated in these two different entities shows marked similarities, with more frequent involvement of genes mediating cell adhesion in the t(4;14)s. Although the t(4;14) group had a greater number of mutations, a larger number of genes were affected in the t(11;14) group with the number of mutated genes in two or more samples being 111 versus 237, respectively. This observation implies a more consistent group of genes are deregulated in the t(4;14) group, suggesting that they are under greater selective pressure than in the t(11;14) group. In this work we show a higher mutation frequency but with more limited numbers of genes affected in the t(4;14) group compared to the t(11;14) group. Overall, the data are consistent within the two etiologically distinct groups of MM having a similar spectrum of mutations driving disease progression, with a focus on pathway deregulation rather than any single gene. Disclosures: No relevant conflicts of interest to declare.

2019 ◽  
Vol 21 (Supplement_6) ◽  
pp. vi111-vi111
Author(s):  
Matija Snuderl ◽  
Kasthuri Kannan ◽  
Jean-Pierre Gagner ◽  
Elad Mashiach ◽  
Matthias Karajannis ◽  
...  

Abstract BACKGROUND While most hemangioblastomas (~70%) are sporadic and occur predominantly in the cerebellum, they may present as well as familial form associated with von Hippel-Lindau (VHL) syndrome, an autosomal dominant disorder caused by germline mutations of the VHL gene that trigger nuclear translocation of hypoxia-inducible factor (HIF)-1α and angiogenesis. Although inactivation of VHL, a tumor suppressor gene, has been observed in hemangioblastomas, the underlying pathogenic mechanisms responsible for familial and sporadic hemangioblastomas remain incompletely understood. METHODS Whole exome sequencing of cerebellar hemangioblastoma tumors and matched blood leukocytes from 24 patients, age 24–63, was performed. After preparation and amplification of barcoded libraries, exomes were captured using Kapa Biosystems methodology and paired-end sequenced on Illumina HiSeq 2500 to an average 100-fold coverage. Following read alignment to hg19 genome, tumor and germline (leukocyte) sequences were compared, and pathogenic single nucleotide variants (SNVs) identified and validated by re-sequencing followed by pathway analysis. Additionally, tumor RNA isolated using Maxwell Promega was sequenced on Illumina instrument and the expression counts determined and normalized. RESULTS We found 314 pathogenic and/or highly deleterious mutations (both germline and somatic) with a median of 13 mutations per patient. Five patients had VHL syndrome (germline VHL mutation) and 4 carried somatic VHL mutations. Among the VHL tumors, 82 mutations were identified, including HNF1B, NOTCH1 and TCF7L1, suggesting a potential contribution of altered RNA metabolism based upon pathway analysis. Among all hemangioblastomas, germline growth factor receptor variants (FGFR4 p.G388R (14/23 (61%) patients), IGF1R, PDGFRA and TYK2) known to activate STAT3 signaling and induce HIF-1α and angiogenesis, were identified. Non-hierarchical clustering of RNA sequencing data revealed two transcriptionally-distinct subtypes of hemangioblastomas. CONCLUSIONS Our findings indicate that hemangioblastomas can also occur by germline mutations known to activate STAT3 signaling, which may have significant implication in genetic testing and counseling of patients with hemangioblastomas.


Blood ◽  
2013 ◽  
Vol 122 (21) ◽  
pp. 3118-3118
Author(s):  
Rakesh Bam ◽  
Sathisha Upparahalli Venkateshaiah ◽  
Xin Li ◽  
Sharmin Khan ◽  
Wen Ling ◽  
...  

Abstract Primary human myeloma (MM) cells do not survive in culture while current in vitro and in vivo systems for growing these cells are limited to coculture with specific bone marrow (BM) cell type or growth in immunodeficient animal model. The aim of the study was to determine long-term survival and interaction of primary MM plasma cells with a healthy adult human BM that include immune cells capable of functional activation. This system is different from the autologous BM culture that is already affected by the disease. Whole BM cells from healthy donors were cultured in αMEM medium supplemented with 10% FBS and 10% serum pooled from MM patients. Following 7-9 days the cultures were composed of adherent and nonadherent cellular compartments. The nonadherent compartment contained typical BM hematopoietic cells such as monocytes, B and T lymphocytes and NK and normal plasma cells as assessed by flow cytometry, while the adherent compartment contained cells that morphologically resemble macrophages, osteoclasts, megakaryocytes and fibroblast-like cells. At this culture stage, CD138-selected MM cells from 20 patients were added to the BM cultures (4:1 BM:MM cell ratio) and survival and growth of MM cells were determined after 7 days by assessing proportion of CD45low/intermediate/CD38high MM plasma cells among total number of cells. MM and BM cell viability was constantly high (>90%) in cocultures. Subsets of primary MM plasma cells, regardless of molecular risk or subtype, were survived and detected in all cases while in 14 of 20 experiments, number of MM plasma cells was increased by 58±12% (p<0.0005, n=14). MM cell proliferation following long-term coculture was evident by the loss of cell membrane PKH26 dye or by BudR uptake in dividing cocultured MM cells. Growth of primary MM was superior in cocultures supplemented with patient serum compared to healthy donor serum. In additional study, we stably infected IL6- or stroma-dependent MM lines, or two primary MM cell cases capable of passaging in SCID-hu mice with EGFP/luciferase construct and demonstrated increased MM cell growth in all experiments in coculture using bioluminescence analysis (statistical significance range: p<0.04 to p<0.0003). Growth of OPM2 MM line was also enhanced in coculture compared to culture alone. The coculture conditions protected OPM2 cells from dexamethasone but not bortezomib while proportion of MM cell killing by lenalidomide was enhanced compared to culture of OPM2 cells alone. To assess the effect of MM cells on BM cells in coculture, global gene expression profile was performed on BM cells cultured alone or plasma cell-depleted BM after coculture with MM cells from 4 patients. Among the top underexpressed genes we identified immunoglobulin genes related to polyclonal plasma cells, extracellular factors associated with osteoblastogenesis (e.g. MGP, IGFBP2), WNT signaling (e.g. SOX4, LRP1, LRP6) and TGFb bioavailability (e.g. FBN1, LTBP1). Top upregulated genes include immuneregulatory factors such as PROK2, LRG1, OLFM4 and IL16, and cellular markers (e.g. ARG1 expressed by MDSCs). This culture system demonstrates the ability of primary MM cells to interact with and to survive in coculture with healthy adult BM that was first cultivated by patients' serum and is appropriate for studying MM-microenvironment interaction, characterization of MM cell subpopulations capable of long term survival and targeted therapy. Disclosures: No relevant conflicts of interest to declare.


2019 ◽  
Author(s):  
Navid Ahmadinejad ◽  
Shayna Troftgruben ◽  
Carlo Maley ◽  
Junwen Wang ◽  
Li Liu

ABSTRACTUnderstanding intratumor heterogeneity is critical to designing personalized treatments and improving clinical outcomes of cancers. Such investigations require accurate delineation of the subclonal composition of a tumor, which to date can only be reliably inferred from deep-sequencing data (>300x depth). To enable accurate subclonal discovery in tumors sequenced at standard depths (30-50x), we develop a novel computational method that incorporates an adaptive error model into statistical decomposition of mixed populations, which corrects the mean-variance dependency of sequencing data at the subclonal level. Tested on extensive computer simulations and real-world data, this new method, named model-based adaptive grouping of subclones (MAGOS), consistently outperforms existing methods on minimum sequencing depth, decomposition accuracy and computation efficiency. MAGOS supports subclone analysis using single nucleotide variants and copy number variants from one or more samples of an individual tumor. Applications of MAGOS to whole-exome sequencing data of 331 liver cancer samples discovered a significant association between subclonal diversity and patient overall survival. MAGOS is freely available as an R package at github (https://github.com/liliulab/magos).


Cancers ◽  
2019 ◽  
Vol 11 (3) ◽  
pp. 362 ◽  
Author(s):  
Marcos Díaz-Gay ◽  
Sebastià Franch-Expósito ◽  
Coral Arnau-Collell ◽  
Solip Park ◽  
Fran Supek ◽  
...  

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson’s two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.


2021 ◽  
Vol 3 (Supplement_3) ◽  
pp. iii17-iii17
Author(s):  
Lucy Boyce Kennedy ◽  
Amanda E D Van Swearingen ◽  
Jeff Sheng ◽  
Dadong Zhang ◽  
Xiaodi Qin ◽  
...  

Abstract Background MBM have a unique molecular profile compared to ECM. Methods We analyzed a previously published dataset from MD Anderson Cancer Center, including RNA-seq on surgically resected, FFPE MBM and ECM from the same patients. STAR pipeline was used to estimate mRNA abundance. DESeq2 package was used to perform differential gene expression (DGE) analyses. Pathway analysis was performed using Gene Set Enrichment Analysis (GSEA). Paired DGE and GSEA compared MBM vs. lymph node (LN) metastases (n = 16) and MBM vs. skin mets (n = 10). CIBERSORTx estimated relative abundance of immune cell types in MBM and ECM. GATK Mutect2 pipeline was used to call somatic mutations using paired normal tumor samples. Mutations were annotated using the Ensembl Variant Effect Predictor and visualized using the Maftools package in R. RNA-seq was available on 54 human primary cutaneous melanomas (CM). Gene Ontology or KEGG Pathway analysis was performed using goana function of limma package in R. Results Paired GSEA found that autophagy pathways may be up-regulated in MBM vs. LN and MBM vs. skin mets. On a single-gene level, the most strongly up-regulated genes in autophagy pathways were GFAP and HBB. Fold changes in other autophagy-related genes were low and did not reach significance. Comparison between CM which recurred in brain vs. CM which did not recur identified up-regulation of autophagy pathways. CIBERSORTx identified an increased proportion of immune suppressive M2 macrophages compared to tumor suppressive M1 macrophages in MBMs and ECMs. Conclusion Up-regulation of autophagy pathways was observed in patient-matched MBM vs. LN and skin mets. This finding was driven by up-regulation of GFAP and HBB, which could reflect changes in the tumor microenvironment. Higher M2:M1 ratio may contribute to an immune suppressive tumor microenvironment and may be targetable. Validation of our findings in an independent Duke dataset is ongoing.


2020 ◽  
Vol 21 (1) ◽  
Author(s):  
Junhua Rao ◽  
Lihua Peng ◽  
Xinming Liang ◽  
Hui Jiang ◽  
Chunyu Geng ◽  
...  

Abstract Background DNBSEQ™ platforms are new massively parallel sequencing (MPS) platforms that use DNA nanoball technology. Use of data generated from DNBSEQ™ platforms to detect single nucleotide variants (SNVs) and small insertions and deletions (indels) has proven to be quite effective, while the feasibility of copy number variants (CNVs) detection is unclear. Results Here, we first benchmarked different CNV detection tools based on Illumina whole-genome sequencing (WGS) data of NA12878 and then assessed these tools in CNV detection based on DNBSEQ™ sequencing data from the same sample. When the same tool was used, the CNVs detected based on DNBSEQ™ and Illumina data were similar in quantity, length and distribution, while great differences existed within results from different tools and even based on data from a single platform. We further estimated the CNV detection power based on available CNV benchmarks of NA12878 and found similar precision and sensitivity between the DNBSEQ™ and Illumina platforms. We also found higher precision of CNVs shorter than 1 kbp based on DNBSEQ™ platforms than those based on Illumina platforms by using Pindel, DELLY and LUMPY. We carefully compared these two available benchmarks and found a large proportion of specific CNVs between them. Thus, we constructed a more complete CNV benchmark of NA12878 containing 3512 CNV regions. Conclusions We assessed and benchmarked CNV detections based on WGS with DNBSEQ™ platforms and provide guidelines for future studies.


2021 ◽  
Author(s):  
Dou-Dou Ding ◽  
Quan Zhou ◽  
Ze He ◽  
Hong-Xia He ◽  
Man-Zhen Zuo

Abstract Introduction:Epidemiological studies have found that the occurrence of endometrial cancer(EC) is closely related to metabolic diseases, and insulin resistance (IR) plays an important role in the pathogenesis of endometrium, but the specific pathogenesis is still unclear. The purpose of this study is to reveal the relationship between insulin resistance and endothelial cells by gene screening technology. Material and methods:We analyzed one endometrial carcinoma database (GSE106191) and one insulin-resistant database (GSE63992), with Gene Expression Omnibus (GEO) database and Venny online analysis tool, then, we found an add-up to 148 different genes. Functional enrichment analysis of these genes using DAVID showed that they were participated in transcription factor activity,signaling pathways and response to factors, etc. Then used cytoHubba in Cytoscape,we got 25 hub genes.Results: The results showed that the survival time of OGT, IGSF3, TRO, NEURL2 and PIK3C2B was significantly and closely related to EC, and the percentage of gene changes of five central genes ranged from 3% to 10% of a single gene, was also related to the infiltration of seven kinds of immune cells in endometrial carcinoma.Conclusion:The five key genes (OGT,IGSF3, PIK3C2B,TRO and NEURL2) are involved in immune infiltration in the progression of endometrial carcinoma, and there is also a certain mutation probability in gene mutation. This may be the pathogenesis of insulin resistance and endometrial cancer.


2021 ◽  
Author(s):  
Tobias Fehlmann ◽  
Fabian Kern ◽  
Omar Laham ◽  
Christina Backes ◽  
Jeffrey Solomon ◽  
...  

Abstract Analyzing all features of small non-coding RNA sequencing data can be demanding and challenging. To facilitate this process, we developed miRMaster. After the analysis of over 125 000 human samples and 1.5 trillion human small RNA reads over 4 years, we present miRMaster 2 with a wide range of updates and new features. We extended our reference data sets so that miRMaster 2 now supports the analysis of eight species (e.g. human, mouse, chicken, dog, cow) and 10 non-coding RNA classes (e.g. microRNAs, piRNAs, tRNAs, rRNAs, circRNAs). We also incorporated new downstream analysis modules such as batch effect analysis or sample embeddings using UMAP, and updated annotation data bases included by default (miRBase, Ensembl, GtRNAdb). To accommodate the increasing popularity of single cell small-RNA sequencing data, we incorporated a module for unique molecular identifier (UMI) processing. Further, the output tables and graphics have been improved based on user feedback and new output formats that emerged in the community are now supported (e.g. miRGFF3). Finally, we integrated differential expression analysis with the miRNA enrichment analysis tool miEAA. miRMaster is freely available at https://www.ccb.uni-saarland.de/mirmaster2.


2021 ◽  
Vol 12 ◽  
Author(s):  
Tao Feng ◽  
Dechao Wei ◽  
Qiankun Li ◽  
Xiaobing Yang ◽  
Yili Han ◽  
...  

Prostate cancer (PCa) is one of the most common malignancies for males, but very little is known about its pathogenesis. This study aimed to identify novel biomarkers associated with PCa prognosis and elucidate the underlying molecular mechanism. First, The Cancer Genome Atlas (TCGA) RNA-sequencing data were utilized to identify differentially expressed genes (DEGs) between tumor and normal samples. The DEGs were then applied to construct a co-expression and mined using structure network analysis. The magenta module that was highly related to the Gleason score (r = 0.46, p = 3e–26) and tumor stage (r = 0.38, p = 2e–17) was screened. Subsequently, all genes of the magenta module underwent function annotation. From the key module, CCNA2, CKAP2L, NCAPG, and NUSAP1 were chosen as the four candidate genes. Finally, internal (TCGA) and external data sets (GSE32571, GSE70770, and GSE141551) were combined to validate and predict the value of real hub genes. The results show that the above genes are up-regulated in PCa samples, and higher expression levels show significant association with higher Gleason scores and tumor T stage. Moreover, receiver operating characteristic curve and survival analysis validate the excellent value of hub genes in PCa progression and prognosis. In addition, the protein levels of these four genes also remain higher in tumor tissues when compared with normal tissues. Gene set enrichment analysis and gene set variation analysis for a single gene reveal the close relation with cell proliferation. Meanwhile, 11 small molecular drugs that have the potential to treat PCa were also screened. In conclusion, our research identified four potential prognostic genes and several candidate molecular drugs for treating PCa.


2021 ◽  
Vol 4 (1) ◽  
pp. 21
Author(s):  
Austin Bow

The reduction in costs associated with performing RNA-sequencing has driven an increase in the application of this analytical technique; however, restrictive factors associated with this tool have now shifted from budgetary constraints to time required for data processing. The sheer scale of the raw data produced can present a formidable challenge for researchers aiming to glean vital information about samples. Though many of the companies that perform RNA-sequencing provide a basic report for the submitted samples, this may not adequately capture particular pathways of interest for sample comparisons. To further assess these data, it can therefore be necessary to utilize various enrichment and mapping software platforms to highlight specific relations. With the wide array of these software platforms available, this can also present a daunting task. The methodology described herein aims to enable researchers new to handling RNA-sequencing data with a streamlined approach to pathway analysis. Additionally, the implemented software platforms are readily available and free to utilize, making this approach viable, even for restrictive budgets. The resulting tables and nodal networks will provide valuable insight into samples and can be used to generate high-quality graphics for publications and presentations.


Sign in / Sign up

Export Citation Format

Share Document