Publications

The Long-Read Personalized OncoGenomics (POG) dataset comprises a cohort of 189 patient tumors and 41 matched normal samples sequenced using the Oxford Nanopore Technologies PromethION platform. This dataset from the POG program and the Marathon of Hope Cancer Centres Network includes DNA and RNA short-read sequence data, analytics, and clinical information. We show the potential of long-read sequencing for resolving complex cancer-related structural variants, viral integrations, and extrachromosomal circular DNA. Long-range phasing facilitates the discovery of allelically differentially methylated regions (aDMRs) and allele-specific expression, including recurrent aDMRs in the cancer genes RET and CDKN2A. Germline promoter methylation in MLH1 can be directly observed in Lynch syndrome. Promoter methylation in BRCA1 and RAD51C is a likely driver behind homologous recombination deficiency where no coding driver mutation was found. This dataset demonstrates applications for long-read sequencing in precision medicine and is available as a resource for developing analytical approaches using this technology.

The evolutionarily conserved ATG4 cysteine proteases regulate macroautophagy/autophagy through the priming and deconjugation of the Atg8-family proteins. In mammals there are four ATG4 family members (ATG4A, ATG4B, ATG4C, ATG4D) but ATG4D has been relatively understudied. Heightened interest in ATG4D has been stimulated by recent links to human disease. Notably, genetic variations in human were implicated in a heritable neurodevelopmental disorder. Genetic analyses in dogs, along with loss-of-function zebrafish and mouse models, further support a neuroprotective role for ATG4D. Here we discuss the evidence connecting ATG4D to neurological diseases and other pathologies and summarize its roles in both autophagy-dependent and autophagy-independent cellular processes.: ATG: autophagy related; BafA1: bafilomycin A1; BCL2: BCL2 apoptosis regulator; BH3: BCL2 homology region 3; CASP3: caspase 3; EV: extracellular vesicle; GABA: gamma aminobutyric acid; GABARAP: GABA type A receptor-associated protein; GABARAPL1: GABA type A receptor associated protein like 1; GABARAPL2: GABA type A receptor associated protein like 2; GFP: green fluorescent protein; LIR: LC3-interacting region; MAP1LC3: microtubule associated protein 1 light chain 3; MEF: mouse embryonic fibroblast; MYC: MYC proto-oncogene, bHLH transcription factor; PE: phosphatidylethanolamine; PS: phosphatidylserine; QKO: quadruple knockout; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel; SQSTM1: sequestosome 1.

Quality assurance (QA) is essential for precision oncology workflows, in particular in the clinical setting. However, because of numerous variations in laboratory and bioinformatics pipelines, QA practices remain non-standardized, are often ad hoc, and are lacking longitudinal tracking. A selected review of existing software was performed for quality control of Illumina next-generation sequencing data, focusing specifically on generalizable tools that can be integrated into any bioinformatics workflow to easily develop a QA workflow with longitudinal tracking. Although all implementations need to be integrated, validated, and iterated upon to suit individual operations, providing a base suite of options will enable better validation and use of QA in clinical somatic mutation testing for workflows using Illumina next-generation sequencing and beyond.

Somatic hypermutation (SHM) status of the immunoglobulin heavy variable (IGHV) gene plays a crucial role in determining the prognosis and treatment of patients with chronic lymphocytic leukemia (CLL). A common approach for determining SHM status is multiplex polymerase chain reaction and Sanger sequencing of the immunoglobin heavy locus; however, this technique is low throughput, is vulnerable to failure, and does not allow multiplexing with other diagnostic assays.

Pancreatic ductal adenocarcinoma (PDAC) exhibits elevated levels of autophagy, which promote tumor progression and treatment resistance. ATG4B is an autophagy-related cysteine protease under consideration as a potential therapeutic target, but it is largely unexplored in PDAC. Here, we investigated the clinical and functional relevance of ATG4B expression in PDAC. Using two PDAC patient cohorts, we found that low ATG4B mRNA or protein expression is associated with worse patient survival outcomes, poorly differentiated PDAC tumors and a lack of survival benefit from adjuvant chemotherapy. In PDAC cell lines, ATG4B knockout reduced proliferation, abolished processing of LC3B (also known as MAP1LC3B), and reduced GABARAP and GABARAPL1 levels, but increased ATG4A levels. ATG4B and ATG4A double knockout lines displayed a further reduction in proliferation, characterized by delays in G1-S phase transition and mitosis. Pro-LC3B accumulated aberrantly at the centrosome with a concomitant increase in centrosomal proteins PCM1 and CEP131, which was rescued by exogenous ATG4B. The two-stage cell cycle defects following ATG4B and ATG4A loss have important therapeutic implications for PDAC.

Hematopoietic stem and progenitor cells (HSPCs) originate from an endothelial-to-hematopoietic transition (EHT) during embryogenesis. Characterization of early hemogenic endothelial (HE) cells is required to understand what drives hemogenic specification and to accurately define cells capable of undergoing EHT. Using Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-seq), we define the early subpopulation of pre-HE cells based on both surface markers and transcriptomes. We identify the transcription factor Meis1 as an essential regulator of hemogenic cell specification in the embryo prior to Runx1 expression. Meis1 is expressed at the earliest stages of EHT and distinguishes pre-HE cells primed towards the hemogenic trajectory from the arterial endothelial cells that continue towards a vascular fate. Endothelial-specific deletion of Meis1 impairs the formation of functional Runx1-expressing HE which significantly impedes the emergence of pre-HSPC via EHT. Our findings implicate Meis1 in a critical fate-determining step for establishing EHT potential in endothelial cells.

Germline structural variants (SVs) are challenging to resolve by conventional genetic testing assays. Long-read sequencing has improved the global characterization of SVs, but its sensitivity at cancer susceptibility loci has not been reported. Nanopore long-read genome sequencing was performed for nineteen individuals with pathogenic copy number alterations in BRCA1, BRCA2, CHEK2 and PALB2 identified by prior clinical testing. Fourteen variants, which spanned single exons to whole genes and included a tandem duplication, were accurately represented. Defining the precise breakpoints of SVs in BRCA1 and CHEK2 revealed unforeseen allelic heterogeneity and informed the mechanisms underlying the formation of recurrent deletions. Integrating read-based and statistical phasing further helped define extended haplotypes associated with founder alleles. Long-read sequencing is a sensitive method for characterizing private, recurrent and founder SVs underlying breast cancer susceptibility. Our findings demonstrate the potential for nanopore sequencing as a powerful genetic testing assay in the hereditary cancer setting.

We recently described a 16-gene expression signature for improved risk stratification of acute myeloid leukemia (AML) patients called the AML Prognostic Score (APS). A subset of APS-high-risk AML patients showed increased levels of focal adhesion kinase (FAK), encoded by the Protein Tyrosine Kinase 2 (PTK2) gene, which was correlated with RUNX1 mutations. RUNX1 mutant cells are more sensitive to PTK2 inhibitors. As we were not able to detect RUNX1-binding sites in the PTK2 promoter, we hypothesized that RUNX1 might regulate micro(mi)RNAs that repress PTK2, such that loss-of-function RUNX1 mutations would result in reduced miRNA expression and derepression of PTK2. Examination of paired RNA-seq and miRNA-seq data from 301 AML cases revealed two miRNAs that positively correlated with RUNX1 expression, contained RUNX1-binding sites in their promoters and were predicted to target PTK2. We show that the hsa-let7a-2-3p and hsa-miR-135a-5p promoters are regulated by RUNX1, and that PTK2 is a direct target of both miRNAs. Even in the absence of RUNX1 mutations, hsa-let7a-2-3p and hsa-miR-135a-5p regulate PTK2 expression, and reduced expression of these two miRNAs sensitizes AML cells to PTK2 inhibition. These data explain how RUNX1 regulates PTK2, and identify potential miRNA biomarkers for targeting AML with PTK2 inhibitors.

We present a genome assembly of (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.