Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy.

Ion channel mutations are an important cause of rare Mendelian disorders affecting brain, heart, and other tissues. We performed parallel exome sequencing of 237 channel genes in a well-characterized human sample, comparing variant profiles of unaffected individuals to those with the most common neuronal excitability disorder, sporadic idiopathic epilepsy. Rare missense variation in known Mendelian disease genes is prevalent in both groups at similar complexity, revealing that even deleterious ion channel mutations confer uncertain risk to an individual depending on the other variants with which they are combined. Our findings indicate that variant discovery via large scale sequencing efforts is only a first step in illuminating the complex allelic architecture underlying personal disease risk. We propose that in silico modeling of channel variation in realistic cell and network models will be crucial to future strategies assessing mutation profile pathogenicity and drug response in individuals with a broad spectrum of excitability disorders.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Dr.Nod: computational framework for discovery of regulatory non-coding drivers in tissue-matched distal regulatory elements.

The discovery of cancer driver mutations is a fundamental goal in cancer research. While many cancer driver mutations have been discovered in the protein-coding genome, research into potential cancer drivers in the non-coding regions showed limited success so far. Here, we present a novel comprehensive framework Dr.Nod for detection of non-coding cis-regulatory candidate driver mutations that are associated with dysregulated gene expression using tissue-matched enhancer-gene annotations. Applying the framework to data from over 1500 tumours across eight tissues revealed a 4.4-fold enrichment of candidate driver mutations in regulatory regions of known cancer driver genes. An overarching conclusion that emerges is that the non-coding driver mutations contribute to cancer by significantly altering transcription factor binding sites, leading to upregulation of tissue-matched oncogenes and down-regulation of tumour-suppressor genes. Interestingly, more than half of the detected cancer-promoting non-coding regulatory driver mutations are over 20 kb distant from the cancer-associated genes they regulate. Our results show the importance of tissue-matched enhancer-gene maps, functional impact of mutations, and complex background mutagenesis model for the prediction of non-coding regulatory drivers. In conclusion, our study demonstrates that non-coding mutations in enhancers play a previously underappreciated role in cancer and dysregulation of clinically relevant target genes.

Sensitive tumour detection and classification using plasma cell-free DNA methylomes.

The use of liquid biopsies for cancer detection and management is rapidly gaining prominence1. Current methods for the detection of circulating tumour DNA involve sequencing somatic mutations using cell-free DNA, but the sensitivity of these methods may be low among patients with early-stage cancer given the limited number of recurrent mutations2-5. By contrast, large-scale epigenetic alterations-which are tissue- and cancer-type specific-are not similarly constrained6 and therefore potentially have greater ability to detect and classify cancers in patients with early-stage disease. Here we develop a sensitive, immunoprecipitation-based protocol to analyse the methylome of small quantities of circulating cell-free DNA, and demonstrate the ability to detect large-scale DNA methylation changes that are enriched for tumour-specific patterns. We also demonstrate robust performance in cancer detection and classification across an extensive collection of plasma samples from several tumour types. This work sets the stage to establish biomarkers for the minimally invasive detection, interception and classification of early-stage cancers based on plasma cell-free DNA methylation patterns.

Prediction of acute myeloid leukaemia risk in healthy individuals.

The incidence of acute myeloid leukaemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 65. Most cases arise without any detectable early symptoms and patients usually present with the acute complications of bone marrow failure1. The onset of such de novo AML cases is typically preceded by the accumulation of somatic mutations in preleukaemic haematopoietic stem and progenitor cells (HSPCs) that undergo clonal expansion2,3. However, recurrent AML mutations also accumulate in HSPCs during ageing of healthy individuals who do not develop AML, a phenomenon referred to as age-related clonal haematopoiesis (ARCH)4-8. Here we use deep sequencing to analyse genes that are recurrently mutated in AML to distinguish between individuals who have a high risk of developing AML and those with benign ARCH. We analysed peripheral blood cells from 95 individuals that were obtained on average 6.3 years before AML diagnosis (pre-AML group), together with 414 unselected age- and gender-matched individuals (control group). Pre-AML cases were distinct from controls and had more mutations per sample, higher variant allele frequencies, indicating greater clonal expansion, and showed enrichment of mutations in specific genes. Genetic parameters were used to derive a model that accurately predicted AML-free survival; this model was validated in an independent cohort of 29 pre-AML cases and 262 controls. Because AML is rare, we also developed an AML predictive model using a large electronic health record database that identified individuals at greater risk. Collectively our findings provide proof-of-concept that it is possible to discriminate ARCH from pre-AML many years before malignant transformation. This could in future enable earlier detection and monitoring, and may help to inform intervention.

Tracing the origins of relapse in acute myeloid leukaemia to stem cells.

In acute myeloid leukaemia, long-term survival is poor as most patients relapse despite achieving remission. Historically, the failure of therapy has been thought to be due to mutations that produce drug resistance, possibly arising as a consequence of the mutagenic properties of chemotherapy drugs. However, other lines of evidence have pointed to the pre-existence of drug-resistant cells. For example, deep sequencing of paired diagnosis and relapse acute myeloid leukaemia samples has provided direct evidence that relapse in some cases is generated from minor genetic subclones present at diagnosis that survive chemotherapy, suggesting that resistant cells are generated by evolutionary processes before treatment and are selected by therapy. Nevertheless, the mechanisms of therapy failure and capacity for leukaemic regeneration remain obscure, as sequence analysis alone does not provide insight into the cell types that are fated to drive relapse. Although leukaemia stem cells have been linked to relapse owing to their dormancy and self-renewal properties, and leukaemia stem cell gene expression signatures are highly predictive of therapy failure, experimental studies have been primarily correlative and a role for leukaemia stem cells in acute myeloid leukaemia relapse has not been directly proved. Here, through combined genetic and functional analysis of purified subpopulations and xenografts from paired diagnosis/relapse samples, we identify therapy-resistant cells already present at diagnosis and two major patterns of relapse. In some cases, relapse originated from rare leukaemia stem cells with a haematopoietic stem/progenitor cell phenotype, while in other instances relapse developed from larger subclones of immunophenotypically committed leukaemia cells that retained strong stemness transcriptional signatures. The identification of distinct patterns of relapse should lead to improved methods for disease management and monitoring in acute myeloid leukaemia. Moreover, the shared functional and transcriptional stemness properties that underlie both cellular origins of relapse emphasize the importance of developing new therapeutic approaches that target stemness to prevent relapse.

SCALP syndrome with a germline heterozygous DOCK6 mutation and somatic mosaic NRAS Q61R mutation.

We present a case of SCALP syndrome, which was diagnosed in a male infant with the characteristic findings of sebaceous nevi, central nervous system malformations, aplasia cutis congenita, limbal dermoid, and giant congenital melanocytic nevi, or pigmented nevi. We identified a germline compound heterozygous DOCK6 mutation and a somatic mosaic NRAS Q61R mutation in the giant congenital melanocytic nevus. This report will increase clinician awareness of SCALP syndrome and augment the literature in characterizing this rare syndrome, including its genetic background.

Environmental pesticide exposure and non-Hodgkin lymphoma survival: a population-based study.

BACKGROUND: There is evidence indicating that pesticide exposure is a risk factor for non-Hodgkin lymphoma (NHL) development. However, the association between pesticide exposure and NHL survival is not well-established. METHODS: Using the California Cancer Registry, we identified patients with a first primary diagnosis of NHL from 2010 to 2016 and linked these patients with CalEnviroScreen 3.0 to obtain production agriculture pesticide exposure to 70 chemicals from the state-mandated Pesticide Use Reporting (PUR) by census tract from 2012 to 2014. In addition, data from PUR was integrated into a geographic information system that employs land-use data to estimate cumulative exposure to specific pesticides previously associated with NHL (glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4-dimethylamine salt) between 10 years prior up to 1 year after NHL diagnosis. Multivariable Cox proportional hazards regression models were used to evaluate the association between total pesticide exposure from CalEnviroScreen 3.0 and individual pesticide exposure from geographic land use data and lymphoma-specific and overall survival. RESULTS: Among 35,808 NHL patients identified, 44.2% were exposed to pesticide in their census tract of residence. Glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4-dimethylamine salt exposure was observed in 34.1%, 26.0%, 10.6%, 14.0%, and 12.8% of NHL patients, respectively. Total pesticide exposure at the time of diagnosis was not associated with lymphoma-specific or overall survival. In addition, no association was consistently found between glyphosate, organophosphorus, carbamate, phenoxyherbicide, and 2,4 dimethylamine salt exposure and lymphoma-specific or overall survival. CONCLUSIONS: Although we found no consistent associations between agricultural pesticide exposure at the neighborhood level and worse survival, these results provide a platform for designing future studies to determine the association between pesticide and NHL.

Coming of age: ten years of next-generation sequencing technologies.

Since the completion of the human genome project in 2003, extraordinary progress has been made in genome sequencing technologies, which has led to a decreased cost per megabase and an increase in the number and diversity of sequenced genomes. An astonishing complexity of genome architecture has been revealed, bringing these sequencing technologies to even greater advancements. Some approaches maximize the number of bases sequenced in the least amount of time, generating a wealth of data that can be used to understand increasingly complex phenotypes. Alternatively, other approaches now aim to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions. These and other strategies are providing researchers and clinicians a variety of tools to probe genomes in greater depth, leading to an enhanced understanding of how genome sequence variants underlie phenotype and disease.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns.

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Wastewater surveillance for SARS-CoV-2 in a small coastal community: Effects of tourism on viral presence and variant identification among low prevalence populations.

Wastewater-based epidemiology has been used to measure SARS-CoV-2 prevalence in cities worldwide as an indicator of community health, however, few longitudinal studies have followed SARS-CoV-2 in wastewater in small communities from the start of the pandemic or evaluated the influence of tourism on viral loads. Therefore the objective of this study was to use measurements of SARS-CoV-2 in wastewater to monitor viral trends and variants in a small island community over a twelve-month period beginning May 1, 2020, before the community re-opened to tourists. Wastewater samples were collected weekly and analyzed to detect and quantify SARS-CoV-2 genome copies. Sanger sequencing was used to determine genome sequences from total RNA extracted from wastewater samples positive for SARS-CoV-2. Visitor data was collected from the local Chamber of Commerce. We performed Poisson and linear regression to determine if visitors to the Cedar Key Chamber of Commerce were positively associated with SARS-CoV-2-positive wastewater samples and the concentration of SARS-CoV-2 RNA. Results indicated that weekly wastewater samples were negative for SARS-CoV-2 until mid-July when positive samples were recorded in four of five consecutive weeks. Additional positive results were recorded in November and December 2020, as well as January, March, and April 2021. Tourism data revealed that the SARS-CoV-2 RNA concentration in wastewater increased by 1.06 Log10 genomic copies/L per 100 tourists weekly. Sequencing from six positive wastewater samples yielded two complete sequences of SARS-CoV-2, two overlapping sequences, and two low yield sequences. They show arrival of a new variant SARS-CoV-2 in January 2021. Our results demonstrate the utility of wastewater surveillance for SARS-CoV-2 in a small community. Wastewater surveillance and viral genome sequencing suggest that population mobility likely plays an important role in the introduction and circulation of SARS-CoV-2 variants among communities experiencing high tourism and who have a small population size.

Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line.

The SK-BR-3 cell line is one of the most important models for HER2+ breast cancers, which affect one in five breast cancer patients. SK-BR-3 is known to be highly rearranged, although much of the variation is in complex and repetitive regions that may be underreported. Addressing this, we sequenced SK-BR-3 using long-read single molecule sequencing from Pacific Biosciences and develop one of the most detailed maps of structural variations (SVs) in a cancer genome available, with nearly 20,000 variants present, most of which were missed by short-read sequencing. Surrounding the important ERBB2 oncogene (also known as HER2), we discover a complex sequence of nested duplications and translocations, suggesting a punctuated progression. Full-length transcriptome sequencing further revealed several novel gene fusions within the nested genomic variants. Combining long-read genome and transcriptome sequencing enables an in-depth analysis of how SVs disrupt the genome and sheds new light on the complex mechanisms involved in cancer genome evolution.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia.

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

Medical Procedure Services in Internal Medicine Residencies in the US: a Systematic Review and Meta-Analysis.

BACKGROUND AND AIMS: The number of procedures performed by internal medicine residents in the United States (US) is declining. An increasing proportion of residents do not feel confident performing essential invasive bedside procedures and, upon graduation, desire additional training. Several residency programs have utilized the medical procedure service (MPS) to address this issue. We aim to summarize the current state of evidence by systematically evaluating the effect of the MPS on resident education, comfort, and training, as well as patient safety and procedural outcomes in the US. METHODS: We conducted a systematic review of all studies reporting the use of an MPS with supervision from a board-certified physician in internal medicine residencies in the US. Database search was performed on PubMed, Embase, ERIC, and Cochrane Library from January 2000 to November 2020 for relevant studies. Quality of evidence assessment and random-effects proportion meta-analyses were performed. RESULTS: A total of nine studies reporting on 3879 procedures performed by MPS were identified. Procedures were safely performed, with a pooled complication rate of 2.1% (95% CI: 1.0-3.5) and generally successful, with a pooled success rate of 94.7% (95% CI: 90.8-97.7). The range of procedures performed by residents under MPS was 6.7-72.8 procedures per month (n = 9) compared to 4.3-64.4 procedures (n = 4) without MPS. MPS significantly increased confidence, comfort, and use of appropriate safety measures among residents. CONCLUSION: There are a limited number of published studies on MPS supervised by a board-certified physician in US internal medicine residencies. Procedures performed by MPS are generally successfully completed and safe. MPS benefits internal medicine residents training by improving competency, comfort, and confidence.

Germline and somatic DNA repair gene alterations in prostate cancer.

BACKGROUND: Emerging evidence has suggested that DNA repair gene alterations may be important in prostate cancer pathogenesis. In the current study, the authors sought to characterize alterations in DNA repair pathway genes in both primary and metastatic prostate tumors with attention to tissue distribution as well as specific genomic alterations. METHODS: The authors studied the distribution and type of alterations in 24 genes that are considered important for DNA repair in 944 prostate cancers harvested from localized and metastatic tumors. Tumor DNA underwent hybrid capture for all coding exons of 287 or 395 cancer-related genes plus select introns from 19 or 31 genes frequently rearranged in cancer. Captured libraries were sequenced to a median exon coverage depth of >×500. Specific genomic alterations were characterized and the frequencies of mutations by tissue site (prostate vs metastases) were compared using logistic regression. RESULTS: A total of 152 patients from the cohort of 944 men (16%) harbored a germline or somatic mutation in ≥1 DNA repair genes. The most frequently mutated genes were BRCA2 (11.4%) and ATM (5.8%), followed by MSH6 (2.5%) and MSH2 (2.1%). Mutations were identified in approximately 20.1% of primary prostate tumors compared with 18.8% of bone metastases. When stratified by tissue site, the highest rates of DNA repair mutations were found in solid organ metastases, including brain and visceral metastases, compared with prostate. CONCLUSIONS: DNA repair gene mutations are more common in metastatic than localized prostate tumors. Visceral and other solid organ metastases appear enriched for these mutations compared with localized tumors or bone and lymph node metastases.

Identification of pre-leukaemic haematopoietic stem cells in acute leukaemia.

In acute myeloid leukaemia (AML), the cell of origin, nature and biological consequences of initiating lesions, and order of subsequent mutations remain poorly understood, as AML is typically diagnosed without observation of a pre-leukaemic phase. Here, highly purified haematopoietic stem cells (HSCs), progenitor and mature cell fractions from the blood of AML patients were found to contain recurrent DNMT3A mutations (DNMT3A(mut)) at high allele frequency, but without coincident NPM1 mutations (NPM1c) present in AML blasts. DNMT3A(mut)-bearing HSCs showed a multilineage repopulation advantage over non-mutated HSCs in xenografts, establishing their identity as pre-leukaemic HSCs. Pre-leukaemic HSCs were found in remission samples, indicating that they survive chemotherapy. Therefore DNMT3A(mut) arises early in AML evolution, probably in HSCs, leading to a clonally expanded pool of pre-leukaemic HSCs from which AML evolves. Our findings provide a paradigm for the detection and treatment of pre-leukaemic clones before the acquisition of additional genetic lesions engenders greater therapeutic resistance.

Derivation and Validation of the CREST Model for Very Early Prediction of Circulatory Etiology Death in Patients Without ST-Segment-Elevation Myocardial Infarction After Cardiac Arrest.

BACKGROUND: No practical tool quantitates the risk of circulatory-etiology death (CED) immediately after successful cardiopulmonary resuscitation in patients without ST-segment-elevation myocardial infarction. We developed and validated a prediction model to rapidly determine that risk and facilitate triage to individualized treatment pathways. METHODS: With the use of INTCAR (International Cardiac Arrest Registry), an 87-question data set representing 44 centers in the United States and Europe, patients were classified as having had CED or a combined end point of neurological-etiology death or survival. Demographics and clinical factors were modeled in a derivation cohort, and backward stepwise logistic regression was used to identify factors independently associated with CED. We demonstrated model performance using area under the curve and the Hosmer-Lemeshow test in the derivation and validation cohorts, and assigned a simplified point-scoring system. RESULTS: Among 638 patients in the derivation cohort, 121 (18.9%) had CED. The final model included preexisting coronary artery disease (odds ratio [OR], 2.86; confidence interval [CI], 1.83-4.49; P≤0.001), nonshockable rhythm (OR, 1.75; CI, 1.10-2.77; P=0.017), initial ejection fraction<30% (OR, 2.11; CI, 1.32-3.37; P=0.002), shock at presentation (OR, 2.27; CI, 1.42-3.62; P<0.001), and ischemic time >25 minutes (OR, 1.42; CI, 0.90-2.23; P=0.13). The derivation model area under the curve was 0.73, and Hosmer-Lemeshow test P=0.47. Outcomes were similar in the 318-patient validation cohort (area under the curve 0.68, Hosmer-Lemeshow test P=0.41). When assigned a point for each associated factor in the derivation model, the average predicted versus observed probability of CED with a CREST score (coronary artery disease, initial heart rhythm, low ejection fraction, shock at the time of admission, and ischemic time >25 minutes) of 0 to 5 was: 7.1% versus 10.2%, 9.5% versus 11%, 22.5% versus 19.6%, 32.4% versus 29.6%, 38.5% versus 30%, and 55.7% versus 50%. CONCLUSIONS: The CREST model stratified patients immediately after resuscitation according to risk of a circulatory-etiology death. The tool may allow for estimation of circulatory risk and improve the triage of survivors of cardiac arrest without ST-segment-elevation myocardial infarction at the point of care.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

Functional genomics identifies specific vulnerabilities in PTEN-deficient breast cancer.

BACKGROUND: Phosphatase and tensin homolog (PTEN) is one of the most frequently inactivated tumor suppressors in breast cancer. While PTEN itself is not considered a druggable target, PTEN synthetic-sick or synthetic-lethal (PTEN-SSL) genes are potential drug targets in PTEN-deficient breast cancers. Therefore, with the aim of identifying potential targets for precision breast cancer therapy, we sought to discover PTEN-SSL genes present in a broad spectrum of breast cancers. METHODS: To discover broad-spectrum PTEN-SSL genes in breast cancer, we used a multi-step approach that started with (1) a genome-wide short interfering RNA (siRNA) screen of ~ 21,000 genes in a pair of isogenic human mammary epithelial cell lines, followed by (2) a short hairpin RNA (shRNA) screen of ~ 1200 genes focused on hits from the first screen in a panel of 11 breast cancer cell lines; we then determined reproducibility of hits by (3) identification of overlaps between our results and reanalyzed data from 3 independent gene-essentiality screens, and finally, for selected candidate PTEN-SSL genes we (4) confirmed PTEN-SSL activity using either drug sensitivity experiments in a panel of 19 cell lines or mutual exclusivity analysis of publicly available pan-cancer somatic mutation data. RESULTS: The screens (steps 1 and 2) and the reproducibility analysis (step 3) identified six candidate broad-spectrum PTEN-SSL genes (PIK3CB, ADAMTS20, AP1M2, HMMR, STK11, and NUAK1). PIK3CB was previously identified as PTEN-SSL, while the other five genes represent novel PTEN-SSL candidates. Confirmation studies (step 4) provided additional evidence that NUAK1 and STK11 have PTEN-SSL patterns of activity. Consistent with PTEN-SSL status, inhibition of the NUAK1 protein kinase by the small molecule drug HTH-01-015 selectively impaired viability in multiple PTEN-deficient breast cancer cell lines, while mutations affecting STK11 and PTEN were largely mutually exclusive across large pan-cancer data sets. CONCLUSIONS: Six genes showed PTEN-SSL patterns of activity in a large proportion of PTEN-deficient breast cancer cell lines and are potential specific vulnerabilities in PTEN-deficient breast cancer. Furthermore, the NUAK1 PTEN-SSL vulnerability identified by RNA interference techniques can be recapitulated and exploited using the small molecule kinase inhibitor HTH-01-015. Thus, NUAK1 inhibition may be an effective strategy for precision treatment of PTEN-deficient breast tumors.