Differential RNA splicing as a potentially important driver mechanism in multiple myeloma.

Disruption of the normal splicing patterns of RNA is a major factor in the pathogenesis of a number of diseases. Increasingly research has shown the strong influence that splicing patterns can have on cancer progression. Multiple Myeloma is a molecularly heterogeneous disease classified by the presence of key translocations, gene expression profiles and mutations but the splicing patterns in MM remains largely unexplored. We take a multifaceted approach to define the extent and impact of alternative splicing in MM. We look at the spliceosome component, SF3B1, with hotspot mutations (K700E and K666T/Q) shown to result in an increase in alternative splicing in other cancers. We discovered a number of differentially spliced genes in comparison of the SF3B1 mutant and wild type samples that included, MZB1, DYNLL1, TMEM14C and splicing related genes DHX9, CLASRP, and SNRPE. We identified a broader role for abnormal splicing showing clear differences in the extent of novel splice variants in the different translocation groups. We show that a high number of novel splice loci is associated with adverse survival and an ultra-high risk group. The enumeration of patterns of alternative splicing has the potential to refine MM classification and to aid in the risk stratification of patients.

TarPan: an easily adaptable targeted sequencing panel viewer for research and clinical use.

BACKGROUND: The study of cancer genomics continually matures as the number of patient samples sequenced increases. As more data is generated, oncogenic drivers for specific cancer types are discovered along with their associated risks. This in turn leads to potential treatment strategies that pave the way to precision medicine. However, significant financial and analytical barriers make it infeasible to sequence the entire genome of every patient. In contrast, targeted sequencing panels give reliable information on relevant portions of the genome at a fiscally responsible cost. Therefore, we have created the Targeted Panel (TarPan) Viewer, a software tool, to investigate this type of data. RESULTS: TarPan Viewer helps investigators understand data from targeted sequencing data by displaying the information through a web browser interface. Through this interface, investigators can easily observe copy number changes, mutations, and structural events in cancer samples. The viewer runs in R Shiny with a robust SQLite backend and its input is generated from bioinformatic algorithms reliably described in the literature. Here we show the results from using TarPan Viewer on publicly available follicular lymphoma, breast cancer, and multiple myeloma data. In addition, we have tested and utilized the viewer internally, and this data has been used in high-impact peer-reviewed publications. CONCLUSIONS: We have designed a flexible, simple to setup viewer that is easily adaptable to any type of cancer targeted sequencing, and has already proven its use in a research laboratory environment. Further, we believe with deeper sequencing and/or more targeted application it could be of use in the clinic in conjunction with an appropriate targeted sequencing panel as a cost-effective diagnostic test, especially in cancers such as acute leukemia or diffuse large B-cell lymphoma that require rapid interventions.

Microhomology-mediated end joining drives complex rearrangements and overexpression of MYC and PVT1 in multiple myeloma.

MYC is a widely acting transcription factor and its deregulation is a crucial event in many human cancers. MYC is important biologically and clinically in multiple myeloma, but the mechanisms underlying its dysregulation are poorly understood. We show that MYC rearrangements are present in 36.0% of newly diagnosed myeloma patients, as detected in the largest set of next generation sequencing data to date (n=1,267). Rearrangements were complex and associated with increased expression of MYC and PVT1, but not other genes at 8q24. The highest effect on gene expression was detected in cases where the MYC locus is juxtaposed next to super-enhancers associated with genes such as IGH, IGK, IGL, TXNDC5/BMP6, FAM46C and FOXO3 We identified three hotspots of recombination at 8q24, one of which is enriched for IGH-MYC translocations. Breakpoint analysis indicates primary myeloma rearrangements involving the IGH locus occur through non-homologous end joining, whereas secondary MYC rearrangements occur through microhomology-mediated end joining. This mechanism is different to lymphomas, where non-homologous end joining generates MYC rearrangements. Rearrangements resulted in overexpression of key genes and chromatin immunoprecipitation-sequencing identified that HK2, a member of the glucose metabolism pathway, is directly over-expressed through binding of MYC at its promoter.

Clonal selection and double-hit events involving tumor suppressor genes underlie relapse in myeloma.

To elucidate the mechanisms underlying relapse from chemotherapy in multiple myeloma, we performed a longitudinal study of 33 patients entered into Total Therapy protocols investigating them using gene expression profiling, high-resolution copy number arrays, and whole-exome sequencing. The study illustrates the mechanistic importance of acquired mutations in known myeloma driver genes and the critical nature of biallelic inactivation events affecting tumor suppressor genes, especially TP53, the end result being resistance to apoptosis and increased proliferation rates, which drive relapse by Darwinian-type clonal evolution. The number of copy number aberration changes and biallelic inactivation of tumor suppressor genes was increased in GEP70 high risk, consistent with genomic instability being a key feature of high risk. In conclusion, the study highlights the impact of acquired genetic events, which enhance the evolutionary fitness level of myeloma-propagating cells to survive multiagent chemotherapy and to result in relapse.

Enhancing cancer clonality analysis with integrative genomics.

INTRODUCTION: It is understood that cancer is a clonal disease initiated by a single cell, and that metastasis, which is the spread of cancer from the primary site, is also initiated by a single cell. The seemingly natural capability of cancer to adapt dynamically in a Darwinian manner is a primary reason for therapeutic failures. Survival advantages may be induced by cancer therapies and also occur as a result of inherent cell and microenvironmental factors. The selected "more fit" clones outmatch their competition and then become dominant in the tumor via propagation of progeny. This clonal expansion leads to relapse, therapeutic resistance and eventually death. The goal of this study is to develop and demonstrate a more detailed clonality approach by utilizing integrative genomics. METHODS: Patient tumor samples were profiled by Whole Exome Sequencing (WES) and RNA-seq on an Illumina HiSeq 2500 and methylation profiling was performed on the Illumina Infinium 450K array. STAR and the Haplotype Caller were used for RNA-seq processing. Custom approaches were used for the integration of the multi-omic datasets. RESULTS: Reported are major enhancements to CloneViz, which now provides capabilities enabling a formal tumor multi-dimensional clonality analysis by integrating: i) DNA mutations, ii) RNA expressed mutations, and iii) DNA methylation data. RNA and DNA methylation integration were not previously possible, by CloneViz (previous version) or any other clonality method to date. This new approach, named iCloneViz (integrated CloneViz) employs visualization and quantitative methods, revealing an integrative genomic mutational dissection and traceability (DNA, RNA, epigenetics) thru the different layers of molecular structures. CONCLUSION: The iCloneViz approach can be used for analysis of clonal evolution and mutational dynamics of multi-omic data sets. Revealing tumor clonal complexity in an integrative and quantitative manner facilitates improved mutational characterization, understanding, and therapeutic assignments.

Revealing the inherent heterogeneity of human malignancies by variant consensus strategies coupled with cancer clonal analysis.

Tumors are heterogeneous in composition. They are composed of cancer cells proper, along with stromal elements that collectively form a microenvironment, all of which are necessary to nurture the malignant process. In addition, many of the stromal cells are modified to support the unique needs of the malignant state. Tumors are composed of a variety of clones or subpopulations of cancer cells, which may differ in karyotype, growth rate, expression of cell surface markers, sensitivity to therapeutics, etc. New tools and methods to provide an improved understanding of tumor clonal architecture are needed to guide therapy.

XLPM: efficient algorithm for the analysis of protein-protein contacts using chemical cross-linking mass spectrometry.

BACKGROUND: Chemical cross-linking is used for protein-protein contacts mapping and for structural analysis. One of the difficulties in cross-linking studies is the analysis of mass-spectrometry data and the assignment of the site of cross-link incorporation. The difficulties are due to higher charges of fragment ions, and to the overall low-abundance of cross-link species in the background of linear peptides. Cross-linkers non-specific at one end, such as photo-inducible diazirines, may complicate the analysis further. In this report, we design and validate a novel cross-linked peptide mapping algorithm (XLPM) and compare it to StavroX, which is currently one of the best algorithms in this class. RESULTS: We have designed a novel cross-link search algorithm -XLPM - and implemented it both as an online tool and as a downloadable archive of scripts. We designed a filter based on an observation that observation of a b-ion implies observation of a complimentary y-ion with high probability (b-y filter). We validated the b-y filter on the set of linear peptides from NIST library, and demonstrate that it is an effective way to find high-quality mass spectra. Next, we generated cross-linked data from an ssDNA binding protein, Rim1with a specific cross-linker disuccinimidyl suberate, and a semi-specific cross-linker NHS-Diazirine, followed by analysis of the cross-linked products by nanoLC-LTQ-Orbitrap mass spectrometry. The cross-linked data were searched by XLPM and StavroX and the performance of the two algorithms was compared. The cross-links were mapped to the X-ray structure of Rim1 tetramer. Analysis of the mixture of NHS-Diazirine cross-linked ¹5N and ¹4N-labeled Rim1 tetramers yielded ¹5N-labeled to ¹4N-labeled cross-linked peptide pairs, corresponding to C-terminus-to-N-terminus cross-linking, demonstrating interaction between different two Rim1 tetramers. Both XLPM and StavroX were successful in identification of this interaction, with XLPM leading to a better annotation of higher-charged fragments. We also put forward a new method of estimating specificity and sensitivity of identification of a cross-linked residue in the case of a non-specific cross-linker. CONCLUSIONS: The novel cross-link mapping algorithm, XLPM, considerably improves the speed and accuracy of the analysis compared to other methods. The quality selection filter based on b-to-y ions ratio proved to be an effective way to select high quality cross-linked spectra.

Leveraging the new with the old: providing a framework for the integration of historic microarray studies with next generation sequencing.

Next Generation Sequencing (NGS) methods are rapidly providing remarkable advances in our ability to study the molecular profiles of human cancers. However, the scientific discovery offered by NGS also includes challenges concerning the interpretation of large and non-trivial experimental results. This task is potentially further complicated when a multitude of molecular profiling modalities are available, with the goal of a more integrative and comprehensive analysis of the cancer biology.

RNA sequencing of formalin fixed paraffin-embedded heart tissue provides transcriptomic information about chemotherapy-induced cardiotoxicity.

Gene expression of formalin-fixed paraffin-embedded (FFPE) tissue may serve for molecular studies on cardiovascular diseases. Chemotherapeutics, such as doxorubicin (DOX) may cause heart injury, but the mechanisms of these side effects of DOX are not well understood. This study aimed to investigate whether DOX-induced gene expression in archival FFPE heart tissue in experimental rats would correlate with the gene expression in fresh-frozen heart tissue by applying RNA sequencing technology. The results showed RNA from FFPE samples was degraded, resulting in a lower number of uniquely mapped reads. However, DOX-induced differentially expressed genes in FFPE were related to molecular mechanisms of DOX-induced cardiotoxicity, such as inflammation, calcium binding, endothelial dysfunction, senescence, and cardiac hypertrophy signaling. Our data suggest that, despite the limitations, RNA sequencing of archival FFPE heart tissue supports utilizing FFPE tissues from retrospective studies on cardiovascular disorders, including DOX-induced cardiotoxicity.

An app-based ecological momentary assessment of undergraduate student mental Health during the COVID-19 pandemic in Canada (Smart Healthy Campus Version 2.0): Longitudinal study.

This paper presents results from the Smart Healthy Campus 2.0 study/smartphone app, developed and used to collect mental health-related lifestyle data from 86 Canadian undergraduates January-August 2021. Objectives of the study were to 1) address the absence of longitudinal mental health overview and lifestyle-related data from Canadian undergraduate students, and 2) to identify associations between these self-reported mental health overviews (questionnaires) and lifestyle-related measures (from smartphone digital measures). This was a longitudinal repeat measures study conducted over 40 weeks. A 9-item mental health questionnaire was accessible once daily in the app. Two variants of this mental health questionnaire existed; the first was a weekly variant, available each Monday or until a participant responded during the week. The second was a daily variant available after the weekly variant. 6518 digital measure samples and 1722 questionnaire responses were collected. Mixed models were fit for responses to the two questionnaire variants and 12 phone digital measures (e.g. GPS, step counts). The daily questionnaire had positive associations with floors walked, installed apps, and campus proximity, while having negative associations with uptime, and daily calendar events. Daily depression had a positive association with uptime. Daily resilience appeared to have a slight positive association with campus proximity. The weekly questionnaire variant had positive associations with device idling and installed apps, and negative associations with floors walked, calendar events, and campus proximity. Physical activity, weekly, had a negative association with uptime, and a positive association with calendar events and device idling. These lifestyle indicators that associated with student mental health during the COVID-19 pandemic suggest directions for new mental health-related interventions (digital or otherwise) and further efforts in mental health surveillance under comparable circumstances.

Towards the integration, annotation and association of historical microarray experiments with RNA-seq.

BACKGROUND: Transcriptome analysis by microarrays has produced important advances in biomedicine. For instance in multiple myeloma (MM), microarray approaches led to the development of an effective disease subtyping via cluster assignment, and a 70 gene risk score. Both enabled an improved molecular understanding of MM, and have provided prognostic information for the purposes of clinical management. Many researchers are now transitioning to Next Generation Sequencing (NGS) approaches and RNA-seq in particular, due to its discovery-based nature, improved sensitivity, and dynamic range. Additionally, RNA-seq allows for the analysis of gene isoforms, splice variants, and novel gene fusions. Given the voluminous amounts of historical microarray data, there is now a need to associate and integrate microarray and RNA-seq data via advanced bioinformatic approaches. METHODS: Custom software was developed following a model-view-controller (MVC) approach to integrate Affymetrix probe set-IDs, and gene annotation information from a variety of sources. The tool/approach employs an assortment of strategies to integrate, cross reference, and associate microarray and RNA-seq datasets. RESULTS: Output from a variety of transcriptome reconstruction and quantitation tools (e.g., Cufflinks) can be directly integrated, and/or associated with Affymetrix probe set data, as well as necessary gene identifiers and/or symbols from a diversity of sources. Strategies are employed to maximize the annotation and cross referencing process. Custom gene sets (e.g., MM 70 risk score (GEP-70)) can be specified, and the tool can be directly assimilated into an RNA-seq pipeline. CONCLUSION: A novel bioinformatic approach to aid in the facilitation of both annotation and association of historic microarray data, in conjunction with richer RNA-seq data, is now assisting with the study of MM cancer biology.

Usability survey of biomedical question answering systems.

We live in an age of access to more information than ever before. This can be a double-edged sword. Increased access to information allows for more informed and empowered researchers, while information overload becomes an increasingly serious risk. Thus, there is a need for intelligent information retrieval systems that can summarize relevant and reliable textual sources to satisfy a user's query. Question answering is a specialized type of information retrieval with the aim of returning precise short answers to queries posed as natural language questions. We present a review and comparison of three biomedical question answering systems: askHERMES (http://www.askhermes.org/), EAGLi (http://eagl.unige.ch/EAGLi/), and HONQA (http://services.hon.ch/cgi-bin/QA10/qa.pl).

Community organization perspectives on COVID-19 vaccine hesitancy and how they increased COVID-19 vaccine confidence: a Canadian Immunization Research Network, social sciences and humanities network study.

Background: COVID-19 vaccines play a critical role in reducing the morbidity and mortality associated with SARS-CoV-2 infection and despite vaccine availability, disparities in COVID-19 vaccine uptake among Canadian subgroups exist. Community organizations are uniquely situated to relay important vaccine messaging around all vaccines, understand components of vaccine hesitancy, and facilitate vaccine uptake within the communities they serve. The objective of this research was to solicit community organizations perspectives specific to COVID-19 vaccines and explore strategies of increasing vaccine uptake within their communities. Methods: A qualitative focus group study was held in the spring of 2021 with 40 community organizations from across the country. Discussions focused on COVID-19 vaccine communication and awareness within their communities, vaccine misinformation, and strategies to increase vaccine acceptance and access. Data were analyzed utilizing thematic and inductive techniques. Results: Vaccine hesitancy was identified among staff and clients. Vaccine confidence, complacency, convenience, and mistrust in government and authorities were identified as contributors to vaccine hesitancy. Community organizations utilized innovative and novel methods to encourage vaccine uptake and increase vaccine confidence. Leveraging established trusting relationships was key to successful messaging within communities. Conclusion: Community organizations used innovative methods, built on established trust, to increase vaccine confidence within their communities and among their staff. Community agencies played an important role in COVID-19 vaccine uptake within subgroups of the Canadian population. Community organizations are key public health partners and play a critical role in increasing COVID-19 vaccine confidence.

Association of DNA methylation signatures with cognitive performance among smokers and ex-smokers.

INTRODUCTION: Alterations in DNA methylation profiles have been associated with cancer, and can be influenced by environmental factors such as smoking. A small but growing literature indicates there are reproducible and robust differences in methylation levels among smokers, never smokers, and ex-smokers. Here, we compared differences in salivary DNA methylation levels among current and ex-smokers (at least 2 years abstinent). METHODS: Smokers (n=26) and ex-smokers (n=30) provided detailed smoking histories, completed the Paced Auditory Serial Addition Test (PASAT), and submitted a saliva sample. Whole-genome DNA methylation from saliva was performed, and ANCOVA models and a receiver operating characteristic (ROC) curve were used for the differences between groups and the performance of significant CpG sites. RESULTS: After controlling for race, age, and gender, smokers had significantly lower methylation levels than ex-smokers in two CpG sites: cg05575921 (AHRR) and cg21566642 (ALPPL2). Based on the ROC analyses, both CpGs had strong classification potentials (cg05575921 AUC=0.97 and cg21566642 AUC=0.93) in differentiating smoking status. Across all subjects, the percent methylation of cg05575921 (AHRR) and cg21566642 (ALPPL2) positively correlated with the length of the last quit attempt (r=0.65 and 0.64, respectively, p<0.001) and PASAT accuracy (r=0.29 and 0.30, respectively, p<0.05). CONCLUSIONS: In spite of the small sample size and preliminary research, our results replicate previously reported differences in AHRR hypomethylation among smokers. Furthermore, we show that the duration of smoking abstinence is associated with a recovery of methylation in ex-smokers, which may be linked to a reduced risk of smoking-associated diseases. The association with cognitive performance suggests that the hypomethylation of AHRR in saliva may reflect systemic exposure to cigarette-related toxicants that negatively affect cognitive performance, and should be validated in larger studies.

Increased Utilization of Low-Dose CT for Lung Cancer Screening at an Arkansas Community Oncology Clinic.

BACKGROUND: Low-dose CT (LDCT) is underused in Arkansas for lung cancer screening, a rural state with a high incidence of lung cancer. The objective was to determine whether offering free LDCT increased the number of high-risk individuals screened in a rural catchment area. METHODS: There were 5,402 patients enrolled in screening at Highlands Oncology, a community oncology clinic in Northwest Arkansas, from 2013 to 2020. Screenings were separated into time periods: period 1 (10 months for-fee), period 2 (10 months free with targeted advertisements and primary care outreach), and period 3 (62 months free with only primary care outreach). In all, 5,035 high-risk participants were eligible for analysis based on National Comprehensive Cancer Network Clinical Practice Guidelines in Oncology. Enrollment rates, incidence densities (IDs), Cox proportional hazard models, and Kaplan-Meier curves were performed to investigate differences between enrollment periods and high-risk groups. RESULTS: Patient volume increased drastically once screenings were offered free of charge (period 1 = 4.6 versus period 2 = 66.0 and period 3 = 69.8 average patients per month). Incidence density per 1,000 person-years increased through each period (IDPeriod 1 = 17.2; IDPeriod 2 = 20.8; IDPeriod 3 = 25.5 cases). Cox models revealed significant differences in lung cancer risk between high-risk groups (P = .012) but not enrollment periods (P = .19). Kaplan-Meier lung cancer-free probabilities differed significantly between high-risk groups (log-rank P = .00068) but not enrollment periods (log-rank P = .18). CONCLUSIONS: This study suggests that eligible patients are more receptive to free LDCT screening, despite most insurances not having a required copay for eligible patients.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma.

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

The spatio-temporal evolution of multiple myeloma from baseline to relapse-refractory states.

Deciphering Multiple Myeloma evolution in the whole bone marrow is key to inform curative strategies. Here, we perform spatial-longitudinal whole-exome sequencing, including 140 samples collected from 24 Multiple Myeloma patients during up to 14 years. Applying imaging-guided sampling we observe three evolutionary patterns, including relapse driven by a single-cell expansion, competing/co-existing sub-clones, and unique sub-clones at distinct locations. While we do not find the unique relapse sub-clone in the baseline focal lesion(s), we show a close phylogenetic relationship between baseline focal lesions and relapse disease, highlighting focal lesions as hotspots of tumor evolution. In patients with ≥3 focal lesions on positron-emission-tomography at diagnosis, relapse is driven by multiple distinct sub-clones, whereas in other patients, a single-cell expansion is typically seen (p < 0.01). Notably, we observe resistant sub-clones that can be hidden over years, suggesting that a prerequisite for curative therapies would be to overcome not only tumor heterogeneity but also dormancy.

Genome-Wide DNA Methylation Signatures Predict the Early Asymptomatic Doxorubicin-Induced Cardiotoxicity in Breast Cancer.

Chemotherapy with doxorubicin (DOX) may cause unpredictable cardiotoxicity. This study aimed to determine whether the methylation signature of peripheral blood mononuclear cells (PBMCs) prior to and after the first cycle of DOX-based chemotherapy could predict the risk of cardiotoxicity in breast cancer patients. Cardiotoxicity was defined as a decrease in left ventricular ejection fraction (LVEF) by >10%. DNA methylation of PBMCs from 9 patients with abnormal LVEF and 10 patients with normal LVEF were examined using Infinium HumanMethylation450 BeadChip. We have identified 14,883 differentially methylated CpGs at baseline and 18,718 CpGs after the first cycle of chemotherapy, which significantly correlated with LVEF status. Significant differentially methylated regions (DMRs) were found in the promoter and the gene body of SLFN12, IRF6 and RNF39 in patients with abnormal LVEF. The pathway analysis found enrichment for regulation of transcription, mRNA splicing, pathways in cancer and ErbB2/4 signaling. The preliminary results from this study showed that the DNA methylation profile of PBMCs may predict the risk of DOX-induced cardiotoxicity prior to chemotherapy. Further studies with larger cohorts of patients are needed to confirm these findings.

Low­Level Environmental Heavy Metals are Associated with Obesity Among Postmenopausal Women in a Southern State.

Both arsenic and cadmium are reported to be toxic to humans. The use of saliva as a biomarker of low-level exposures to these elements has not been adequately explored, and the putative relationship between exposure and obesity is unclear. This cross-sectional study aims to investigate the relationship between salivary arsenic and cadmium concentrations and their association with obesity. Arsenic and cadmium concentrations were analyzed in human saliva samples by Inductively Coupled Plasma-Mass Spectrometry on 270 randomly selected women who participated in the Arkansas Rural Community Health Study. Multivariable logistic regression was performed to evaluate the association between heavy metal concentrations and obesity. Stratified logistic regression was performed based on menopausal status. Generalized linear models were used to evaluate weight gain velocity. Significant positive associations were observed in postmenopausal women for both arsenic (OR = 4.43, 95% CI 1.91-10.28) and cadmium (OR = 2.72, 95% CI 1.23-5.99) concentrations, as well as significant trends among tertiles (p < 0.01 and p = 0.01, respectively). No relationship with obesity was evident among premenopausal women for either metal. A dose-response relationship was observed between increasing weight gain velocity and increasing metal concentrations. At concentrations well below governmental and industrial standards for acute toxicity, significant associations between obesity and concentration of these heavy metals are evident. The rate at which individuals gain weight is affected by metal concentrations and may play a role in the rapid increase in weight in postmenopausal women. These results might explain, in part, the missing variability in the increasing obesity pandemic in certain population exposed to these environmental toxicants.