RESUMO
Manual interpretation of variants remains rate limiting in precision oncology. The increasing scale and complexity of molecular data generated from comprehensive sequencing of cancer samples requires advanced interpretative platforms as precision oncology expands beyond individual patients to entire populations. To address this unmet need, we introduce a Platform for Oncogenomic Reporting and Interpretation (PORI), comprising an analytic framework that facilitates the interpretation and reporting of somatic variants in cancer. PORI integrates reporting and graph knowledge base tools combined with support for manual curation at the reporting stage. PORI represents an open-source platform alternative to commercial reporting solutions suitable for comprehensive genomic data sets in precision oncology. We demonstrate the utility of PORI by matching 9,961 pan-cancer genome atlas tumours to the graph knowledge base, calculating therapeutically informative alterations, and making available reports describing select individual samples.
Assuntos
Carcinogênese/genética , Neoplasias/genética , Biomarcadores Tumorais , Bases de Dados Genéticas , Variação Genética , Genômica , Humanos , Bases de Conhecimento , Medicina de PrecisãoRESUMO
Advanced and metastatic tumors with complex treatment histories drive cancer mortality. Here we describe the POG570 cohort, a comprehensive whole-genome, transcriptome and clinical dataset, amenable for exploration of the impacts of therapies on genomic landscapes. Previous exposure to DNA-damaging chemotherapies and mutations affecting DNA repair genes, including POLQ and genes encoding Polζ, were associated with genome-wide, therapy-induced mutagenesis. Exposure to platinum therapies coincided with signatures SBS31 and DSB5 and, when combined with DNA synthesis inhibitors, signature SBS17b. Alterations in ESR1, EGFR, CTNNB1, FGFR1, VEGFA and DPYD were consistent with drug resistance and sensitivity. Recurrent noncoding events were found in regulatory region hotspots of genes including TERT, PLEKHS1, AP2A1 and ADGRG6. Mutation burden and immune signatures corresponded with overall survival and response to immunotherapy. Our data offer a rich resource for investigation of advanced cancers and interpretation of whole-genome and transcriptome sequencing in the context of a cancer clinic.
Assuntos
Neoplasias , Humanos , Neoplasias/tratamento farmacológicoRESUMO
We report a case of early-onset pancreatic ductal adenocarcinoma in a patient harboring biallelic MUTYH germline mutations, whose tumor featured somatic mutational signatures consistent with defective MUTYH-mediated base excision repair and the associated driver KRAS transversion mutation p.Gly12Cys. Analysis of an additional 730 advanced cancer cases (N = 731) was undertaken to determine whether the mutational signatures were also present in tumors from germline MUTYH heterozygote carriers or if instead the signatures were only seen in those with biallelic loss of function. We identified two patients with breast cancer each carrying a pathogenic germline MUTYH variant with a somatic MUTYH copy loss leading to the germline variant being homozygous in the tumor and demonstrating the same somatic signatures. Our results suggest that monoallelic inactivation of MUTYH is not sufficient for C:G>A:T transversion signatures previously linked to MUTYH deficiency to arise (N = 9), but that biallelic complete loss of MUTYH function can cause such signatures to arise even in tumors not classically seen in MUTYH-associated polyposis (N = 3). Although defective MUTYH is not the only determinant of these signatures, MUTYH germline variants may be present in a subset of patients with tumors demonstrating elevated somatic signatures possibly suggestive of MUTYH deficiency (e.g., COSMIC Signature 18, SigProfiler SBS18/SBS36, SignatureAnalyzer SBS18/SBS36).
Assuntos
Neoplasias da Mama/genética , Carcinoma Ductal Pancreático/genética , DNA Glicosilases/genética , Mutação , Neoplasias Pancreáticas/genética , Idade de Início , DNA Glicosilases/deficiência , Feminino , Mutação em Linhagem Germinativa , Humanos , Perda de Heterozigosidade , Pessoa de Meia-Idade , Proteínas Proto-Oncogênicas p21(ras)/genéticaRESUMO
Summary: Reliably identifying genomic rearrangements and interpreting their impact is a key step in understanding their role in human cancers and inherited genetic diseases. Many short read algorithmic approaches exist but all have appreciable false negative rates. A common approach is to evaluate the union of multiple tools increasing sensitivity, followed by filtering to retain specificity. Here we describe an application framework for the rapid generation of structural variant consensus, unique in its ability to visualize the genetic impact and context as well as process both genome and transcriptome data. Availability and implementation: http://mavis.bcgsc.ca. Supplementary information: Supplementary data are available at Bioinformatics online.
Assuntos
Genômica , Neoplasias/genética , Software , Biologia Computacional , Humanos , TranscriptomaRESUMO
Eccrine porocarcinomas (EPs) are rare malignant tumours of the intraepidermic sweat gland duct and most often arise from benign eccrine poromas. Some recurrent somatic genomic events have been identified in these malignancies, but very little is known about the complexity of their molecular pathophysiology. We describe the whole genome and whole transcriptome genomic profiling of a metastatic EP in a 66-year-old male patient with a previous history of localized porocarcinoma of the scalp. Whole genome and whole transcriptome genomic profiling was performed on the metastatic EP. Whole genome sequencing was performed on blood-derived DNA in order to allow a comparison between germline and somatic events. We found somatic copy losses of several tumour suppressor genes including APC, PTEN and CDKN2A, CDKN2B and CDKN1A. We identified a somatic hemizygous CDKN2A pathogenic splice site variant. De novo transcriptome assembly revealed abnormal splicing of CDKN2A p14ARF and p16INK4a. Elevated expression of oncogenes EGFR and NOTCH1 was noted and no somatic mutations were found in these genes. Wnt pathway somatic alterations were also observed. In conclusion, our results suggest that the molecular pathophysiology of malignant EP features high complexity and subtle interactions of multiple key genes. Cell cycle dysregulation and CDKN2A loss of function was found to be a new potential driver in EP tumourigenesis. Moreover, the combination of somatic copy number variants and abnormal gene expression perhaps partly related to epigenetic mechanisms, all likely contribute to the development of this rare malignancy in our patient.
RESUMO
X-ray crystallographic data of the carbohydrate mimic MDWNMHAA when bound to an anti-Shigella flexneri Y mAb SYA/J6 indicate the immobilization of water molecules, that is, the presence of "bound" waters, in the active site. Water Ligand Observed via Gradient Spectroscopy (WaterLOGSY) was used in conjunction with saturation transfer difference (STD)-NMR spectroscopy to probe the existence of immobilized water molecules in the complex of MDWNMHAA 1 bound to mAb SYA/J6. Molecular dynamics simulations using the ZymeCAD Molecular Dynamics platform were then used to specify the likely locations of these water molecules. Of note, those waters involved in providing complementarity between the peptide and mAb SYA/J6 remained throughout the course of the simulation. Together, the experimental and computational protocols have been used to identify the bound water molecules present in the antibody-peptide complex.
Assuntos
Anticorpos Monoclonais/química , Anticorpos Monoclonais/metabolismo , Anticorpos/química , Anticorpos/metabolismo , Carboidratos/química , Peptídeos/química , Polissacarídeos Bacterianos/química , Polissacarídeos Bacterianos/imunologia , Shigella flexneri/química , Shigella flexneri/metabolismo , Água/química , Sequência de Carboidratos , Carboidratos/imunologia , Cristalografia por Raios X , Ligantes , Espectroscopia de Ressonância Magnética , Modelos Moleculares , Simulação de Dinâmica Molecular , Estrutura Molecular , Peptídeos/metabolismo , Ligação Proteica , Shigella flexneri/imunologiaRESUMO
Saturation transfer difference (STD)-NMR spectroscopy was used to probe experimentally the bioactive solution conformation of the carbohydrate mimic MDWNMHAA 1 of the O-polysaccharide of Shigella flexneri Y when bound to its complementary antibody, mAb SYA/J6. Molecular dynamics simulations using the ZymeCAD™ Molecular Dynamics platform were also undertaken to give a more accurate picture of the conformational flexibility and the possibilities for bound ligand conformations. The ligand topology, or the dynamic epitope, was mapped with the CORCEMA-ST (COmplete Relaxation and Conformational Exchange Matrix Analysis of Saturation Transfer) program that calculates a total matrix analysis of relaxation and exchange effects to generate predicted STD-NMR intensities from simulation. The comparison of these predicted STD enhancements with experimental data was used to select a representative binding mode. A protocol that employed theoretical STD effects calculated at snapshots during the entire course of a molecular dynamics (MD) trajectory of the peptide bound to the Fv portion of the antibody, and not the averaged atomic positions of receptor-ligand complexes, was also examined. In addition, the R factor was calculated on the basis of STD (fit) to avoid T1 bias, and an effective R factor, R(eff), was defined such that if the calculated STD (fit) for proton k was within error of the experimental STD (fit) for proton k, then that calculated STD (fit) for proton k was not included in the calculation of the R factor. This protocol was effective in deriving the antibody-bound solution conformation of the peptide which also differed from the bound conformation determined by X-ray crystallography; however, several discrepancies between experimental and calculated STD (fit) values were observed. The bound conformation was therefore further refined with a simulated annealing refinement protocol known as STD-NMR intensity-restrained CORCEMA optimization (SICO) to give a more accurate representation of the bound peptide epitope. Further optimization was required in this case, but a satisfactory correlation between experimental and calculated STD values was obtained. Attempts were also made to obtain STD enhancements with a synthetic pentasaccharide hapten, corresponding to the O-polysaccharide, while bound to the antibody. However, unfavorable kinetics of binding in this system prevented sufficient STD build-up, which, in turn, hindered a rigorous analysis via full STD build-up curves.
Assuntos
Anticorpos/metabolismo , Peptídeos/química , Peptídeos/metabolismo , Shigella flexneri/química , Anticorpos/química , Cristalografia por Raios X , Cinética , Ligantes , Espectroscopia de Ressonância Magnética , Modelos Moleculares , Conformação Molecular , Simulação de Dinâmica Molecular , Ligação Proteica , Shigella flexneri/metabolismoRESUMO
UDP-Galactopyranose mutase (UGM) is a flavoenzyme that catalyzes interconversion of UDP-galactopyranose (UDP-Galp) and UDP-galactofuranose (UDP-Galf); its activity depends on FAD redox state. The enzyme is vital to many pathogens, not native to mammals, and is an important drug target. We have probed binding of substrate, UDP-Galp, and UDP to wild type and W160A UGM from K. pneumoniae, and propose that substrate directs recognition loop dynamics by bridging distal FAD and W160 sites; W160 interacts with uracil of the substrate and is functionally essential. Enhanced Trp fluorescence upon substrate binding to UGM indicates conformational changes remote from the binding site because the fluorescence is unchanged upon binding to W70F/W290F UGM where W160 is the sole Trp. MD simulations map these changes to recognition loop closure to coordinate substrate. This requires galactose-FAD interactions as Trp fluorescence is unchanged upon substrate binding to oxidized UGM, or binding of UDP to either form of the enzyme, and MD show heightened recognition loop mobility in complexes with UDP. Consistent with substrate-directed loop closure, UDP binds 10-fold more tightly to oxidized UGM, yet substrate binds tighter to reduced UGM. This requires the W160-U interaction because redox-switched binding affinity of substrate reverses in the W160A mutant where it only binds when oxidized. Without the anchoring W160-U interaction, an alternative binding mode for UDP is detected, and STD-NMR experiments show simultaneous binding of UDP-Galp and UDP to different subsites in oxidized W160A UGM: Substrate no longer directs recognition loop dynamics to coordinate tight binding to the reduced enzyme.
Assuntos
Transferases Intramoleculares/química , Transferases Intramoleculares/metabolismo , Sítios de Ligação , Catálise , Simulação por Computador , Transferases Intramoleculares/genética , Cinética , Mutagênese Sítio-Dirigida , Oxirredução , Conformação Proteica , Espectrometria de Fluorescência , Especificidade por SubstratoRESUMO
UDP-galactopyranose mutase (UGM) is the key enzyme involved in the biosynthesis of Galf. UDP-Galp and UDP-Galf are two natural substrates of UGM. A protocol that combines the use of STD-NMR spectroscopy, molecular modeling, and CORCEMA-ST calculations was applied to the investigation of the binding of UDP-Galf and its C3-fluorinated analogue to UGM from Klebsiella pneumoniae. UDP-Galf and UDP-[3-F]Galf were bound to UGM in a manner similar to that of UDP-Galp. The interconversions of UDP-Galf and UDP-[3-F]Galf to their galactopyranose counterparts were catalyzed by the reduced (active) UGM with different catalytic efficiencies, as observed by NMR spectroscopy. The binding affinities of UDP-Galf and UDP-[3-F]Galf were also compared with those of UDP-Galp and UDP by competition STD-NMR experiments. When UGM was in the oxidized (inactive) state, the binding affinities of UDP-Galf, UDP-Galp, and UDP-[3-F]Galf were of similar magnitudes and were lower than that of UDP. However, when UGM was in the reduced state, UDP-Galp had higher binding affinity compared with UDP. Molecular dynamics (MD) simulations indicated that the "open" mobile loop in UGM "closes" upon binding of the substrates. Combined MD simulations and STD-NMR experiments were used to create models of UGM with UDP-Galf and UDP-[3-F]Galf as bound ligands. Calculated values of saturation-transfer effects with CORCEMA-ST (complete relaxation and conformational exchange matrix analysis of saturation transfer) were compared to the experimental STD effects and permitted differentiation between two main conformational families of the bound ligands. Taken together, these results are used to rationalize the different rates of catalytic turnover of UDP-Galf and UDP-[3-F]Galf and shed light on the mechanism of action of UGM.
Assuntos
Simulação por Computador , Galactose/análogos & derivados , Transferases Intramoleculares/química , Espectroscopia de Ressonância Magnética/métodos , Espectroscopia de Ressonância Magnética/normas , Modelos Químicos , Difosfato de Uridina/análogos & derivados , Galactose/síntese química , Galactose/química , Klebsiella pneumoniae/enzimologia , Conformação Molecular , Ligação Proteica , Padrões de Referência , Difosfato de Uridina/síntese química , Difosfato de Uridina/químicaRESUMO
The effect of mobility on 15N chemical shift/15N-(1)H dipolar coupling (PISEMA) solid state NMR experiments applied to macroscopically oriented beta-barrels is assessed using molecular dynamics simulation data of the NalP autotransporter domain embedded in a DMPC bilayer. In agreement with previous findings for alpha-helices, the fast librational motion of the peptide planes is found to have a considerable effect on the calculated PISEMA spectra. In addition, the dependence of the chemical shift anisotropy (CSA) and dipolar coupling parameters on the calculated spectra is evaluated specifically for the beta-barrel case. It is found that the precise choice of the value of the CSA parameters sigma11, sigma22 and sigma33 has only a minor effect, whereas the choice of the CSA parameter theta shifts the position of the peaks by up to 20 ppm and changes the overall shape of the spectrum significantly. As was found for alpha-helices, the choice of the NH bond distance has a large effect on the dipolar coupling constant used for the calculations. Overall, it is found that the alternating beta-strands in the barrel occupy distinct regions of the PISEMA spectra, forming patterns which may prove useful in peak assignment.