HTLV-1 infected T cells cause bone loss via small extracellular vesicles.

Adult T cell leukemia (ATL), caused by infection with human T cell leukemia virus type 1 (HTLV-1), is often complicated by hypercalcemia and osteolytic lesions. Therefore, we studied the communication between patient-derived ATL cells (ATL-PDX) and HTLV-1 immortalized CD4+ T cell lines (HTLV/T) with osteoclasts and their effects on bone mass in mice. Intratibial inoculation of some HTLV/T lead to a profound local decrease in bone mass similar to marrow-replacing ATL-PDX, despite the fact that few HTLV/T cells persisted in the bone. To study the direct effect of HTLV/T and ATL-PDX on osteoclasts, supernatants were added to murine and human osteoclast precursors. ATL-PDX supernatants from hypercalcemic patients promoted formation of mature osteoclasts, while those from HTLV/T were variably stimulatory, but had largely consistent effects between human and murine cultures. Interestingly, this osteoclastic activity did not correlate with expression of osteoclastogenic cytokine RANKL, suggesting an alternative mechanism. HTLV/T and ATL-PDX produce small extracellular vesicles (sEV), known to facilitate HTLV-1 infection. We hypothesized that these sEV also mediate bone loss by targeting osteoclasts. We isolated sEV from both HTLV/T and ATL-PDX, and found they carried most of the activity found in supernatants. In contrast, sEV from uninfected activated T cells had little effect. Analysis of sEV (both active and inactive) by mass spectrometry and electron microscopy confirmed absence of RANKL and intact virus. Viral proteins Tax and Env were only present in sEV from the active, osteoclast-stimulatory group, along with increased representation of proteins involved in osteoclastogenesis and bone resorption. sEV injected over mouse calvaria in the presence of low dose RANKL caused more osteolysis than RANKL alone. Thus, HTLV-1 infection of T cells can cause release of sEV with strong osteolytic potential, providing a mechanism beyond RANKL production that modifies the bone microenvironment, even in the absence of overt leukemia.

PDPR Gene Variants Predisposing to Papillary Thyroid Cancer.

Background: Papillary thyroid cancer (PTC) is the predominant subtype of thyroid cancer (THCA), and it can cluster in families with an autosomal dominant (AD) inheritance pattern. The aim of this study was to identify novel genes and mechanisms underlying PTC susceptibility. Methods: Our previous investigation of 17 AD PTC families led us to conduct a deeper analysis on one family (Family Q) with whole-genome sequencing data from 3 PTC-affected individuals. In addition, 323 sporadic THCA cases from Avatar data and 12 familial adenomatous polyposis (FAP) individuals with secondary THCA were screened for pyruvate dehydrogenase phosphatase regulatory (PDPR) variants. CRISPR-Cas9 was used to create PDPR-deficient THCA (TPC1) and transformed normal thyroid cell lines (N-Thyori3-1) to study the metabolic consequences of PDPR loss. Results: We found truncating PDPR splice donor variants (NM_017990.4:c.361 + 1G>C) in all affected PTC Family Q members, and another PDPR splice donor variant (NM_017990.4:c.443 + 1G>C) in a sporadic PTC case. In addition, an ultra-rare missense variant was found in an FAP-PTC patient. The PDPR-deficient cells presented with elevated phosphorylation of pyruvate dehydrogenase and altered glucose metabolism, implying that PDPR plays an essential part in regulating glucose metabolism in thyroid cells. Conclusions: Our finding of novel truncating germline variants in PDPR in Family Q and additional cohorts suggests a role for PDPR loss in PTC predisposition. Also, somatic and RNA sequencing from the thyroid carcinoma (Firehouse Legacy) data showed that PDPR gene expression is much lower in THCA tumor tissue compared with matching normal tissue. Thus, PDPR appears to have a loss of function effect on THCA tumorigenesis.

multi-GPA-Tree: Statistical approach for pleiotropy informed and functional annotation tree guided prioritization of GWAS results.

Genome-wide association studies (GWAS) have successfully identified over two hundred thousand genotype-trait associations. Yet some challenges remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), most with small or moderate effect sizes, making them difficult to detect. Second, many complex traits share a common genetic basis due to 'pleiotropy' and and though few methods consider it, leveraging pleiotropy can improve statistical power to detect genotype-trait associations with weaker effect sizes. Third, currently available statistical methods are limited in explaining the functional mechanisms through which genetic variants are associated with specific or multiple traits. We propose multi-GPA-Tree to address these challenges. The multi-GPA-Tree approach can identify risk SNPs associated with single as well as multiple traits while also identifying the combinations of functional annotations that can explain the mechanisms through which risk-associated SNPs are linked with the traits. First, we implemented simulation studies to evaluate the proposed multi-GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that multi-GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs for multiple traits. Second, we applied multi-GPA-Tree to a systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), and to a Crohn's disease (CD) and ulcertive colitis (UC) GWAS, and functional annotation data including GenoSkyline and GenoSkylinePlus. Our results demonstrate that multi-GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits.

graph-GPA 2.0: improving multi-disease genetic analysis with integration of functional annotation data.

Genome-wide association studies (GWAS) have successfully identified a large number of genetic variants associated with traits and diseases. However, it still remains challenging to fully understand the functional mechanisms underlying many associated variants. This is especially the case when we are interested in variants shared across multiple phenotypes. To address this challenge, we propose graph-GPA 2.0 (GGPA 2.0), a statistical framework to integrate GWAS datasets for multiple phenotypes and incorporate functional annotations within a unified framework. Our simulation studies showed that incorporating functional annotation data using GGPA 2.0 not only improves the detection of disease-associated variants, but also provides a more accurate estimation of relationships among diseases. Next, we analyzed five autoimmune diseases and five psychiatric disorders with the functional annotations derived from GenoSkyline and GenoSkyline-Plus, along with the prior disease graph generated by biomedical literature mining. For autoimmune diseases, GGPA 2.0 identified enrichment for blood-related epigenetic marks, especially B cells and regulatory T cells, across multiple diseases. Psychiatric disorders were enriched for brain-related epigenetic marks, especially the prefrontal cortex and the inferior temporal lobe for bipolar disorder and schizophrenia, respectively. In addition, the pleiotropy between bipolar disorder and schizophrenia was also detected. Finally, we found that GGPA 2.0 is robust to the use of irrelevant and/or incorrect functional annotations. These results demonstrate that GGPA 2.0 can be a powerful tool to identify genetic variants associated with each phenotype or those shared across multiple phenotypes, while also promoting an understanding of functional mechanisms underlying the associated variants.

Genetic analysis of a malignant meningioma and associated metastases.

To identify genes altered in a highly aggressive metastatic meningioma primary as well as its metastases. Exome sequencing of a primary anaplastic meningioma and metastatic lesions in which DNA could be extracted and compared to germline DNA. Genetic analysis of the metastatic sites found 31 common mutations among the primary tumor and two metastatic sites. Additionally, genetic mutations were identified which were either infrequently (MUC3A, ALDH1A3, HOXA1) or not at all previously described in meningiomas (CASS4, CMKLR1). Exome sequencing of a metastatic meningioma and its distant metastases outside the CNS identified mutations that were not previously well described.

GPA-Tree: statistical approach for functional-annotation-tree-guided prioritization of GWAS results.

MOTIVATION: In spite of great success of genome-wide association studies (GWAS), multiple challenges still remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), each with small or moderate effect sizes. Second, our understanding of the functional mechanisms through which genetic variants are associated with complex traits is still limited. To address these challenges, we propose GPA-Tree and it simultaneously implements association mapping and identifies key combinations of functional annotations related to risk-associated SNPs by combining a decision tree algorithm with a hierarchical modeling framework. RESULTS: First, we implemented simulation studies to evaluate the proposed GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs and identifying the true combinations of functional annotations with high accuracy. Second, we applied GPA-Tree to a systemic lupus erythematosus (SLE) GWAS and functional annotation data including GenoSkyline and GenoSkylinePlus. The results from GPA-Tree highlight the dysregulation of blood immune cells, including but not limited to primary B, memory helper T, regulatory T, neutrophils and CD8+ memory T cells in SLE. These results demonstrate that GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits. AVAILABILITY AND IMPLEMENTATION: The GPATree software is available at https://dongjunchung.github.io/GPATree/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Thyroid Carcinomas That Occur in Familial Adenomatous Polyposis Patients Recurrently Harbor Somatic Variants in APC, BRAF, and KTM2D.

Background: Familial adenomatous polyposis (FAP) is a condition typically caused by pathogenic germline mutations in the APC gene. In addition to colon polyps, individuals with FAP have a substantially increased risk of developing papillary thyroid cancer (PTC). Little is known about the events underlying this association, and the prevalence of somatic "second-hit" mutations in APC is controversial. Methods: Whole-genome sequencing was performed on paired thyroid tumor and normal DNA from 12 FAP patients who developed PTC. Somatic mutation profiles were compared with clinical characteristics and previously sequenced sporadic PTC cases. Germline variant profiling was performed to assess the prevalence of variants in genes previously shown to have a role in PTC predisposition. Results: All 12 patients harbored germline mutations in APC, consistent with FAP. Seven patients also had somatic mutations in APC, and seven patients harbored somatic mutations in KMT2D, which encodes a lysine methyl transferase. Mutation of these genes is extremely rare in sporadic PTCs. Notably, only two of the tumors harbored the somatic BRAF p.V600E mutation, which is the most common driver mutation found in sporadic PTCs. Six tumors displayed a cribriform-morular variant of PTC (PTC-CMV) histology, and all six had somatic mutations in APC. Additionally, nine FAP-PTC patients had rare germline variants in genes that were previously associated with thyroid carcinoma. Conclusions: Our data indicate that FAP-associated PTCs typically have distinct mutations compared with sporadic PTCs. Roughly half of the thyroid cancers that arise in FAP patients have somatic "second-hits" in APC, which is associated with PTC-CMV histology. Somatic BRAF p.V600E variants also occur in some FAP patients, a novel finding. We speculate that in carriers of heterozygous pathogenic mutations of tumor suppressor genes such as APC, a cooperating second-hit somatic variant may occur in a different gene such as KTM2D or BRAF, leading to differences in phenotypes. The role of germline variance in genes other than APC (9 of the 12 patients in this series) needs further research.

hsa-miR-9-3p and hsa-miR-9-5p as Post-Transcriptional Modulators of DNA Topoisomerase IIα in Human Leukemia K562 Cells with Acquired Resistance to Etoposide.

DNA topoisomerase IIα protein (TOP2α) 170 kDa (TOP2α/170) is an important target for anticancer agents whose efficacy is often attenuated by chemoresistance. Our laboratory has characterized acquired resistance to etoposide in human leukemia K562 cells. The clonal resistant subline K/VP.5 contains reduced TOP2α/170 mRNA and protein levels compared with parental K562 cells. The aim of this study was to determine whether microRNA (miRNA)-mediated mechanisms play a role in drug resistance via decreased expression of TOP2α/170. miRNA-sequencing revealed that human miR-9-3p and miR-9-5p were among the top six of those overexpressed in K/VP.5 compared with K562 cells; validation by quantitative polymerase chain reaction demonstrated overexpression of both miRNAs. miRNA recognition elements (MREs) for both miRNAs are present in the 3'-untranslated region (UTR) of TOP2α/170. Transfecting K562 cells with a reporter plasmid harboring the TOP2α/170 3'-UTR together with either miR-9-3p or miR-9-5p mimics resulted in a statistically significant decrease in luciferase expression. Mutating the miR-9-3p or miR-9-5p MREs prevented this decrease, demonstrating direct interaction between these miRNAs and TOP2α/170 mRNA. Transfection of K562 cells with miR-9-3p or miR-9-5p mimics led to decreased TOP2α/170 protein levels without a change in TOP2α/170 mRNA and resulted in attenuated etoposide-induced DNA damage (gain-of-miRNA-inhibitory function). Conversely, transfection of miR-9-3p or miR-9-5p inhibitors in K/VP.5 cells (overexpressed miR-9 and low TOP2α/170) led to increased TOP2α/170 protein expression without a change in TOP2α/170 mRNA levels and resulted in enhancement of etoposide-induced DNA damage (loss-of-miRNA-inhibitory function). Taken together, these results strongly suggest that these miRNAs play a role in and are potential targets for circumvention of acquired resistance to etoposide. SIGNIFICANCE STATEMENT: Results presented here indicate that miR-9-3p and miR-9-5p decrease DNA topoisomerase IIα protein 170 kDa expression levels in acquired resistance to etoposide. These findings contribute new information about and potential strategies for circumvention of drug resistance by modulation of microRNA levels. Furthermore, increased expression of miR-9-3p and miR-9-5p in chemoresistant cancer cells may support their validation as biomarkers of responsiveness to DNA topoisomerase II-targeted therapy.

BISR-RNAseq: an efficient and scalable RNAseq analysis workflow with interactive report generation.

BACKGROUND: RNA sequencing has become an increasingly affordable way to profile gene expression patterns. Here we introduce a workflow implementing several open-source softwares that can be run on a high performance computing environment. RESULTS: Developed as a tool by the Bioinformatics Shared Resource Group (BISR) at the Ohio State University, we have applied the pipeline to a few publicly available RNAseq datasets downloaded from GEO in order to demonstrate the feasibility of this workflow. Source code is available here: workflow: https://code.bmi.osumc.edu/gadepalli.3/BISR-RNAseq-ICIBM2019 and shiny: https://code.bmi.osumc.edu/gadepalli.3/BISR_RNASeq_ICIBM19. Example dataset is demonstrated here: https://dataportal.bmi.osumc.edu/RNA_Seq/. CONCLUSION: The workflow allows for the analysis (alignment, QC, gene-wise counts generation) of raw RNAseq data and seamless integration of quality analysis and differential expression results into a configurable R shiny web application.

MiR-34a regulates the invasive capacity of canine osteosarcoma cell lines.

BACKGROUND: Osteosarcoma (OSA) is the most common bone tumor in children and dogs; however, no substantial improvement in clinical outcome has occurred in either species over the past 30 years. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression and play a fundamental role in cancer. The purpose of this study was to investigate the potential contribution of miR-34a loss to the biology of canine OSA, a well-established spontaneous model of the human disease. METHODOLOGY AND PRINCIPAL FINDINGS: RT-qPCR demonstrated that miR-34a expression levels were significantly reduced in primary canine OSA tumors and canine OSA cell lines as compared to normal canine osteoblasts. In canine OSA cell lines stably transduced with empty vector or pre-miR-34a lentiviral constructs, overexpression of miR-34a inhibited cellular invasion and migration but had no effect on cell proliferation or cell cycle distribution. Transcriptional profiling of canine OSA8 cells possessing enforced miR-34a expression demonstrated dysregulation of numerous genes, including significant down-regulation of multiple putative targets of miR-34a. Moreover, gene ontology analysis of down-regulated miR-34a target genes showed enrichment of several biological processes related to cell invasion and motility. Lastly, we validated changes in miR-34a putative target gene expression, including decreased expression of KLF4, SEM3A, and VEGFA transcripts in canine OSA cells overexpressing miR-34a and identified KLF4 and VEGFA as direct target genes of miR-34a. Concordant with these data, primary canine OSA tumor tissues demonstrated increased expression levels of putative miR-34a target genes. CONCLUSIONS: These data demonstrate that miR-34a contributes to invasion and migration in canine OSA cells and suggest that loss of miR-34a may promote a pattern of gene expression contributing to the metastatic phenotype in canine OSA.

Discovery and characterization of the feline miRNAome.

The domestic cat is an important human companion animal that can also serve as a relevant model for ~250 genetic diseases, many metabolic and degenerative conditions, and forms of cancer that are analogous to human disorders. MicroRNAs (miRNAs) play a crucial role in many biological processes and their dysregulation has a significant impact on important cellular pathways and is linked to a variety of diseases. While many species already have a well-defined and characterized miRNAome, miRNAs have not been carefully studied in cats. As a result, there are no feline miRNAs present in the reference miRNA databases, diminishing the usefulness of medical research on spontaneous disease in cats for applicability to both feline and human disease. This study was undertaken to define and characterize the cat miRNAome in normal feline tissues. High-throughput sequencing was performed on 12 different normal cat tissues. 271 candidate feline miRNA precursors, encoding a total of 475 mature sequences, were identified, including several novel cat-specific miRNAs. Several analyses were performed to characterize the discovered miRNAs, including tissue distribution of the precursors and mature sequences, genomic distribution of miRNA genes and identification of clusters, and isomiR characterization. Many of the miRNAs were regulated in a tissue/organ-specific manner.

Differential mutation frequencies in metastatic cutaneous squamous cell carcinomas versus primary tumors.

BACKGROUND: Exome and targeted sequencing studies have identified potential driver mutations for a variety of tumor types. Cutaneous squamous cell carcinoma (cSCC) is one of the most highly mutated cancers but typically is associated with low rates of metastasis and high survival rates. Nevertheless, metastatic cSCC is a significant health threat; up to 8800 individuals die each year of this disease. METHODS: Because it is difficult to predict which cSCCs are more likely to metastasize, and because to the best of the authors' knowledge there are no targeted therapies specifically designated for patients with metastatic cSCC, exome and/or targeted sequencing of 18 metastatic and 10 primary cSCCs was performed to identify mutations that were more frequent in metastatic tumors and might be targeted for therapeutic benefit. The authors compared their results with published sequencing results of an additional 223 primary tumors and 68 metastatic cSCCs. RESULTS: The authors identified genes demonstrating higher mutation frequencies in metastatic cSCC compared with primary tumors, including the chromatin remodeling gene lysine methyltransferase 2D (KMT2D) and the classic skin tumor suppressor tumor protein p53 (TP53), which was found to be mutated in 54% of primary tumors compared with 85% of metastatic tumors (P<.0001). CONCLUSIONS: These studies appear to uncover potential pathways that are important in metastatic cSCC and that broaden understanding of the biology contributing to aggressive tumor behavior. These results may lead to new therapeutic strategies. Cancer 2017;123:1184-1193. © 2016 American Cancer Society.

Impact of Pre-Analytical Variables on Cancer Targeted Gene Sequencing Efficiency.

Tumor specimens are often preserved as formalin-fixed paraffin-embedded (FFPE) tissue blocks, the most common clinical source for DNA sequencing. Herein, we evaluated the effect of pre-sequencing parameters to guide proper sample selection for targeted gene sequencing. Data from 113 FFPE lung tumor specimens were collected, and targeted gene sequencing was performed. Libraries were constructed using custom probes and were paired-end sequenced on a next generation sequencing platform. A PCR-based quality control (QC) assay was utilized to determine DNA quality, and a ratio was generated in comparison to control DNA. We observed that FFPE storage time, PCR/QC ratio, and DNA input in the library preparation were significantly correlated to most parameters of sequencing efficiency including depth of coverage, alignment rate, insert size, and read quality. A combined score using the three parameters was generated and proved highly accurate to predict sequencing metrics. We also showed wide read count variability within the genome, with worse coverage in regions of low GC content like in KRAS. Sample quality and GC content had independent effects on sequencing depth, and the worst results were observed in regions of low GC content in samples with poor quality. Our data confirm that FFPE samples are a reliable source for targeted gene sequencing in cancer, provided adequate sample quality controls are exercised. Tissue quality should be routinely assessed for pre-analytical factors, and sequencing depth may be limited in genomic regions of low GC content if suboptimal samples are utilized.

Somatic Mutation Spectrum of Non-Small-Cell Lung Cancer in African Americans: A Pooled Analysis.

INTRODUCTION: The mutational profile of non-small-cell lung cancer (NSCLC) has become an important tool in tailoring therapy to patients, with clear differences according to the population of origin. African Americans (AAs) have higher lung cancer incidence and mortality than Caucasians, yet discrepant results have been reported regarding the frequency of somatic driver mutations. We hypothesized that NSCLC has a distinct mutational profile in this group. METHODS: We collected NSCLC samples resected from self-reported AAs in five sites from Tennessee, Michigan, and Ohio. Gene mutations were assessed by either SNaPshot or next generation sequencing, and ALK translocations were evaluated by fluorescence in situ hybridization. RESULTS: Two hundred sixty patients were included, mostly males (62.3%) and smokers (86.6%). Eighty-one samples (31.2%) were squamous cell carcinomas. The most frequently mutated genes were KRAS (15.4%), epidermal growth factor receptor (EGFR, 5.0%), PIK3CA (0.8%), BRAF, NRAS, ERBB2, and AKT1 (0.4% each). ALK translocations were detected in two nonsquamous tumors (1.7%), totaling 61 cases (23.5%) with driver oncogenic alterations. Among 179 nonsquamous samples, 54 (30.2%) presented a driver alteration. The frequency of driver alterations altogether was lower than that reported in Caucasians, whereas no difference was detected in either EGFR or KRAS mutations. Overall survival was longer among patients with EGFR mutations. CONCLUSIONS: We demonstrated that NSCLC from AAs has a different pattern of somatic driver mutations than from Caucasians. The majority of driver alterations in this group are yet to be described, which will require more comprehensive panels and assessment of noncanonical alterations.

Resistance mechanisms for the Bruton's tyrosine kinase inhibitor ibrutinib.

BACKGROUND: Ibrutinib is an irreversible inhibitor of Bruton's tyrosine kinase (BTK) and is effective in chronic lymphocytic leukemia (CLL). Resistance to irreversible kinase inhibitors and resistance associated with BTK inhibition have not been characterized. Although only a small proportion of patients have had a relapse during ibrutinib therapy, an understanding of resistance mechanisms is important. We evaluated patients with relapsed disease to identify mutations that may mediate ibrutinib resistance. METHODS: We performed whole-exome sequencing at baseline and the time of relapse on samples from six patients with acquired resistance to ibrutinib therapy. We then performed functional analysis of identified mutations. In addition, we performed Ion Torrent sequencing for identified resistance mutations on samples from nine patients with prolonged lymphocytosis. RESULTS: We identified a cysteine-to-serine mutation in BTK at the binding site of ibrutinib in five patients and identified three distinct mutations in PLCγ2 in two patients. Functional analysis showed that the C481S mutation of BTK results in a protein that is only reversibly inhibited by ibrutinib. The R665W and L845F mutations in PLCγ2 are both potentially gain-of-function mutations that lead to autonomous B-cell-receptor activity. These mutations were not found in any of the patients with prolonged lymphocytosis who were taking ibrutinib. CONCLUSIONS: Resistance to the irreversible BTK inhibitor ibrutinib often involves mutation of a cysteine residue where ibrutinib binding occurs. This finding, combined with two additional mutations in PLCγ2 that are immediately downstream of BTK, underscores the importance of the B-cell-receptor pathway in the mechanism of action of ibrutinib in CLL. (Funded by the National Cancer Institute and others.).

Identification of medium-sized copy number alterations in whole-genome sequencing.

The genome-wide discoveries such as detection of copy number alterations (CNA) from high-throughput whole-genome sequencing data enabled new developments in personalized medicine. The CNAs have been reported to be associated with various diseases and cancers including acute myeloid leukemia. However, there are multiple challenges to the use of current CNA detection tools that lead to high false-positive rates and thus impede widespread use of such tools in cancer research. In this paper, we discuss these issues and propose possible solutions. First, since the entire genome cannot be mapped due to some regions lacking sequence uniqueness, current methods cannot be appropriately adjusted to handle these regions in the analyses. Thus, detection of medium-sized CNAs is also being directly affected by these mappability problems. The requirement for matching control samples is also an important limitation because acquiring matching controls might not be possible or might not be cost efficient. Here we present an approach that addresses these issues and detects medium-sized CNAs in cancer genomes by (1) masking unmappable regions during the initial CNA detection phase, (2) using pool of a few normal samples as control, and (3) employing median filtering to adjust CNA ratios to its surrounding coverage and eliminate false positives.