Cholangiocarcinoma (CCA) is a lethal hepatobiliary neoplasm originating from the biliary apparatus. In humans, CCA risk factors include hepatobiliary inflammation and fibrosis. The recently ...identified interleukin (IL)−1 family member, IL‐33, has been shown to be a biliary mitogen which also promotes liver inflammation and fibrosis. Our aim was to generate a mouse model of CCA mimicking the human disease. Ectopic oncogene expression in the biliary tract was accomplished by the Sleeping Beauty transposon transfection system with transduction of constitutively active AKT (myr‐AKT) and Yes‐associated protein. Intrabiliary instillation of the transposon–transposase complex was coupled with lobar bile duct ligation in C57BL/6 mice, followed by administration of IL‐33 for 3 consecutive days. Tumors developed in 72% of the male mice receiving both oncogenes plus IL‐33 by 10 weeks but in only 20% of the male mice transduced with the oncogenes alone. Tumors expressed SOX9 and pancytokeratin (features of CCA) but were negative for HepPar1 (a marker of hepatocellular carcinoma). Substantive overlap with human CCA specimens was revealed by RNA profiling. Not only did IL‐33 induce IL‐6 expression by human cholangiocytes but it likely facilitated tumor development in vivo by an IL‐6–sensitive process as tumor development was significantly attenuated in Il‐6–/– male animals. Furthermore, tumor formation occurred at a similar rate when IL‐6 was substituted for IL‐33 in this model. Conclusion: The transposase‐mediated transduction of constitutively active AKT and Yes‐associated protein in the biliary epithelium coupled with lobar obstruction and IL‐33 administration results in the development of CCA with morphological and biochemical features of the human disease; this model highlights the role of inflammatory cytokines in CCA oncogenesis. (Hepatology 2015;61:1627–1642)
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Cereblon (CRBN) mediates immunomodulatory drug (IMiD) action in multiple myeloma (MM). Using 2 different methodologies, we identified 244 CRBN binding proteins and established relevance to MM biology ...by changes in their abundance after exposure to lenalidomide. Proteins most reproducibly binding CRBN (>fourfold vs controls) included DDB1, CUL4A, IKZF1, KPNA2, LTF, PFKL, PRKAR2A, RANGAP1, and SHMT2. After lenalidomide treatment, the abundance of 46 CRBN binding proteins decreased. We focused attention on 2 of these—IKZF1 and IKZF3. IZKF expression is similar across all MM stages or subtypes; however, IKZF1 is substantially lower in 3 of 5 IMiD-resistant MM cell lines. The cell line (FR4) with the lowest IKZF1 levels also harbors a damaging mutation and a translocation that upregulates IRF4, an IKZF target. Clinical relevance of CRBN-binding proteins was demonstrated in 44 refractory MM patients treated with pomalidomide and dexamethasone therapy in whom low IKZF1 gene expression predicted lack of response (0/11 responses in the lowest expression quartile). CRBN, IKZF1, and KPNA2 levels also correlate with significant differences in overall survival. Our study identifies CRBN-binding proteins and demonstrates that in addition to CRBN, IKZF1, and KPNA2, expression can predict survival outcomes.
•We identified CRBN-binding proteins and confirmed recent work implicating the Ikaros transcription factors as important members.•Low IZKF1 levels predict lack of IMiD responsiveness and shorter OS; IZKF1 expression is a potential predictive biomarker for clinic application.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Tularemia is a geographically widespread, severely debilitating, and occasionally lethal disease in humans. It is caused by infection by a gram-negative bacterium, Francisella tularensis. In order to ...better understand its potency as an etiological agent as well as its potential as a biological weapon, we have completed draft assemblies and report the first complete genomic characterization of five strains belonging to the following different Francisella subspecies (subsp.): the F. tularensis subsp. tularensis FSC033, F. tularensis subsp. holarctica FSC257 and FSC022, and F. tularensis subsp. novicida GA99-3548 and GA99-3549 strains. Here, we report the sequencing of these strains and comparative genomic analysis with recently available public Francisella sequences, including the rare F. tularensis subsp. mediasiatica FSC147 strain isolate from the Central Asian Region. We report evidence for the occurrence of large-scale rearrangement events in strains of the holarctica subspecies, supporting previous proposals that further phylogenetic subdivisions of the Type B clade are likely. We also find a significant enrichment of disrupted or absent ORFs proximal to predicted breakpoints in the FSC022 strain, including a genetic component of the Type I restriction-modification defense system. Many of the pseudogenes identified are also disrupted in the closely related rarely human pathogenic F. tularensis subsp. mediasiatica FSC147 strain, including modulator of drug activity B (mdaB) (FTT0961), which encodes a known NADPH quinone reductase involved in oxidative stress resistance. We have also identified genes exhibiting sequence similarity to effectors of the Type III (T3SS) and components of the Type IV secretion systems (T4SS). One of the genes, msrA2 (FTT1797c), is disrupted in F. tularensis subsp. mediasiatica and has recently been shown to mediate bacterial pathogen survival in host organisms. Our findings suggest that in addition to the duplication of the Francisella Pathogenicity Island, and acquisition of individual loci, adaptation by gene loss in the more recently emerged tularensis, holarctica, and mediasiatica subspecies occurred and was distinct from evolutionary events that differentiated these subspecies, and the novicida subspecies, from a common ancestor. Our findings are applicable to future studies focused on variations in Francisella subspecies pathogenesis, and of broader interest to studies of genomic pathoadaptation in bacteria.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Whole genome comparative studies of many bacterial pathogens have shown an overall high similarity of gene content (>95%) between phylogenetically distinct subspecies. In highly clonal species that ...share the bulk of their genomes subtle changes in gene content and small-scale polymorphisms, especially those that may alter gene expression and protein-protein interactions, are more likely to have a significant effect on the pathogen's biology. In order to better understand molecular attributes that may mediate the adaptation of virulence in infectious bacteria, a comparative study was done to further analyze the evolution of a gene encoding an o-methyltransferase that was previously identified as a candidate virulence factor due to its conservation specifically in highly pathogenic Francisella tularensis subsp. tularensis strains. The o-methyltransferase gene is located in the genomic neighborhood of a known pathogenicity island and predicted site of rearrangement. Distinct o-methyltransferase subtypes are present in different Francisella tularensis subspecies. Related protein families were identified in several host species as well as species of pathogenic bacteria that are otherwise very distant phylogenetically from Francisella, including species of Mycobacterium. A conserved sequence motif profile is present in the mammalian host and pathogen protein sequences, and sites of non-synonymous variation conserved in Francisella subspecies specific o-methyltransferases map proximally to the predicted active site of the orthologous human protein structure. Altogether, evidence suggests a role of the F. t. subsp. tularensis protein in a mechanism of molecular mimicry, similar perhaps to Legionella and Coxiella. These findings therefore provide insights into the evolution of niche-restriction and virulence in Francisella, and have broader implications regarding the molecular mechanisms that mediate host-pathogen relationships.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract Cereblon (CRBN) mediates immunomodulatory drug (IMiD) action in multiple myeloma (MM). We demonstrate here that no patient with very low CRBN expression responded to IMiD plus dexamethasone ...therapy. In 53 refractory MM patients treated with pomalidomide and dexamethasone, CRBN levels predict for decreased response rates and significant differences in PFS (3.0 vs. 8.9 months, p < 0.001) and OS (9.1 vs. 27.2 months, p = 0.01) (lowest quartile vs. highest three quartiles). While higher CRBN levels can serve as a surrogate for low risk disease, our study demonstrates that low CRBN expression can predict resistance to IMiD monotherapy and is a predictive biomarker for survival outcomes.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
The evolution of resistance in Staphylococcus aureus occurs rapidly, and in response to all known antimicrobial treatments. Numerous studies of model species describe compensatory roles of mutations ...in mediating competitive fitness, and there is growing evidence that these mutation types also drive adaptation of S. aureus strains. However, few studies have tracked amino acid changes during the complete evolutionary trajectory of antibiotic adaptation or been able to predict their functional relevance. Here, we have assessed the efficacy of computational methods to predict biological resistance of a collection of clinically known Resistance Associated Mutations (RAMs). We have found that >90% of known RAMs are incorrectly predicted to be functionally neutral by at least one of the prediction methods used. By tracing the evolutionary histories of all of the false negative RAMs, we have discovered that a significant number are reversion mutations to ancestral alleles also carried in the MSSA476 methicillin-sensitive isolate. These genetic reversions are most prevalent in strains following daptomycin treatment and show a tendency to accumulate in biological pathway reactions that are distinct from those accumulating non-reversion mutations. Our studies therefore show that in addition to non-reversion mutations, reversion mutations arise in isolates exposed to new antibiotic treatments. It is possible that acquisition of reversion mutations in the genome may prevent substantial fitness costs during the progression of resistance. Our findings pose an interesting question to be addressed by further clinical studies regarding whether or not these reversion mutations lead to a renewed vulnerability of a vancomycin or daptomycin resistant strain to antibiotics administered at an earlier stage of infection.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Liposarcoma is the most common soft tissue sarcoma, but little is known about the genomic basis of this disease. Given the low cell content of this tumor type, we utilized flow cytometry to isolate ...the diploid normal and aneuploid tumor populations from a well-differentiated liposarcoma prior to array comparative genomic hybridization and whole genome sequencing. This work revealed massive highly focal amplifications throughout the aneuploid tumor genome including MDM2, a gene that has previously been found to be amplified in well-differentiated liposarcoma. Structural analysis revealed massive rearrangement of chromosome 12 and 11 gene fusions, some of which may be part of double minute chromosomes commonly present in well-differentiated liposarcoma. We identified a hotspot of genomic instability localized to a region of chromosome 12 that includes a highly conserved, putative L1 retrotransposon element, LOC100507498 which resides within a gene cluster (NAV3, SYT1, PAWR) where 6 of the 11 fusion events occurred. Interestingly, a potential gene fusion was also identified in amplified DDR2, which is a potential therapeutic target of kinase inhibitors such as dastinib, that are not routinely used in the treatment of patients with liposarcoma. Furthermore, 7 somatic, damaging single nucleotide variants have also been identified, including D125N in the PTPRQ protein. In conclusion, this work is the first to report the entire genome of a well-differentiated liposarcoma with novel chromosomal rearrangements associated with amplification of therapeutically targetable genes such as MDM2 and DDR2.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
To address the need to study frozen clinical specimens using next-generation RNA, DNA, chromatin immunoprecipitation (ChIP) sequencing and protein analyses, we developed a biobank work flow to ...prospectively collect biospecimens from patients with renal cell carcinoma (RCC). We describe our standard operating procedures and work flow to annotate pathologic results and clinical outcomes. We report quality control outcomes and nucleic acid yields of our RCC submissions (N=16) to The Cancer Genome Atlas (TCGA) project, as well as newer discovery platforms, by describing mass spectrometry analysis of albumin oxidation in plasma and 6 ChIP sequencing libraries generated from nephrectomy specimens after histone H3 lysine 36 trimethylation (H3K36me3) immunoprecipitation. From June 1, 2010, through January 1, 2013, we enrolled 328 patients with RCC. Our mean (SD) TCGA RNA integrity numbers (RINs) were 8.1 (0.8) for papillary RCC, with a 12.5% overall rate of sample disqualification for RIN <7. Banked plasma had significantly less albumin oxidation (by mass spectrometry analysis) than plasma kept at 25 °C (P<.001). For ChIP sequencing, the FastQC score for average read quality was at least 30 for 91% to 95% of paired-end reads. In parallel, we analyzed frozen tissue by RNA sequencing; after genome alignment, only 0.2% to 0.4% of total reads failed the default quality check steps of Bowtie2, which was comparable to the disqualification ratio (0.1%) of the 786-O RCC cell line that was prepared under optimal RNA isolation conditions. The overall correlation coefficients for gene expression between Mayo Clinic vs TCGA tissues ranged from 0.75 to 0.82. These data support the generation of high-quality nucleic acids for genomic analyses from banked RCC. Importantly, the protocol does not interfere with routine clinical care. Collections over defined time points during disease treatment further enhance collaborative efforts to integrate genomic information with outcomes.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Rooting phylogenies is critical for understanding evolution, yet the importance, intricacies and difficulties of rooting are often overlooked. For rooting, polymorphic characters among the group of ...interest (ingroup) must be compared to those of a relative (outgroup) that diverged before the last common ancestor (LCA) of the ingroup. Problems arise if an outgroup does not exist, is unknown, or is so distant that few characters are shared, in which case duplicated genes originating before the LCA can be used as proxy outgroups to root diverse phylogenies. Here, we describe a genome-wide expansion of this technique that can be used to solve problems at the other end of the evolutionary scale: where ingroup individuals are all very closely related to each other, but the next closest relative is very distant. We used shared orthologous single nucleotide polymorphisms (SNPs) from 10 whole genome sequences of Coxiella burnetii, the causative agent of Q fever in humans, to create a robust, but unrooted phylogeny. To maximize the number of characters informative about the rooting, we searched entire genomes for polymorphic duplicated regions where orthologs of each paralog could be identified so that the paralogs could be used to root the tree. Recent radiations, such as those of emerging pathogens, often pose rooting challenges due to a lack of ingroup variation and large genomic differences with known outgroups. Using a phylogenomic approach, we created a robust, rooted phylogeny for C. burnetii.
Full text
Available for:
BFBNIB, DOBA, IZUM, KILJ, NMLJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
DNA focused panel sequencing has been rapidly adopted to assess therapeutic targets in advanced/refractory cancer. Integrated Genomic Profiling (IGP) utilising DNA/RNA with tumour/normal comparisons ...in a Clinical Laboratory Improvement Amendments (CLIA) compliant setting enables a single assay to provide: therapeutic target prioritisation, novel target discovery/application and comprehensive germline assessment. A prospective study in 35 advanced/refractory cancer patients was conducted using CLIA-compliant IGP. Feasibility was assessed by estimating time to results (TTR), prioritising/assigning putative therapeutic targets, assessing drug access, ascertaining germline alterations, and assessing patient preferences/perspectives on data use/reporting. Therapeutic targets were identified using biointelligence/pathway analyses and interpreted by a Genomic Tumour Board. Seventy-five percent of cases harboured 1-3 therapeutically targetable mutations/case (median 79 mutations of potential functional significance/case). Median time to CLIA-validated results was 116 days with CLIA-validation of targets achieved in 21/22 patients. IGP directed treatment was instituted in 13 patients utilising on/off label FDA approved drugs (n = 9), clinical trials (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in five patients (two partial response, three stable disease). Although barriers to broader application exist, including the need for wider availability of therapies, IGP in a CLIA-framework is feasible and valuable in selection/prioritisation of anti-cancer therapeutic targets.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK