Custom workflows to improve joint variant calling from multiple related tumour samples: FreeBayesSomatic and Strelka2Pass.

SUMMARY: This work describes two novel workflows for variant calling that extend the widely used algorithms of Strelka2 and FreeBayes to call somatic mutations from multiple related tumour samples and one matched normal sample. We show that these workflows offer higher precision and recall than their single tumour-normal pair equivalents in both simulated and clinical sequencing data. AVAILABILITY AND IMPLEMENTATION: Source code freely available at the following link: https://atlassian.petermac.org.au/bitbucket/projects/DAW/repos/multisamplevariantcalling and executable through Janis (https://github.com/PMCC-BioinformaticsCore/janis) under the GPLv3 licence. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prediction and monitoring of relapse in stage III melanoma using circulating tumor DNA.

BACKGROUND: The advent of effective adjuvant therapies for patients with resected melanoma has highlighted the need to stratify patients based on risk of relapse given the cost and toxicities associated with treatment. Here we assessed circulating tumor DNA (ctDNA) to predict and monitor relapse in resected stage III melanoma. PATIENTS AND METHODS: Somatic mutations were identified in 99/133 (74%) patients through tumor tissue sequencing. Personalized droplet digital PCR (ddPCR) assays were used to detect known mutations in 315 prospectively collected plasma samples from mutation-positive patients. External validation was performed in a prospective independent cohort (n = 29). RESULTS: ctDNA was detected in 37 of 99 (37%) individuals. In 81 patients who did not receive adjuvant therapy, 90% of patients with ctDNA detected at baseline and 100% of patients with ctDNA detected at the postoperative timepoint relapsed at a median follow up of 20 months. ctDNA detection predicted patients at high risk of relapse at baseline [relapse-free survival (RFS) hazard ratio (HR) 2.9; 95% confidence interval (CI) 1.5-5.6; P = 0.002] and postoperatively (HR 10; 95% CI 4.3-24; P < 0.001). ctDNA detection at baseline [HR 2.9; 95% CI 1.3-5.7; P = 0.003 and postoperatively (HR 11; 95% CI 4.3-27; P < 0.001] was also associated with inferior distant metastasis-free survival (DMFS). These findings were validated in the independent cohort. ctDNA detection remained an independent predictor of RFS and DMFS in multivariate analyses after adjustment for disease stage and BRAF mutation status. CONCLUSION: Baseline and postoperative ctDNA detection in two independent prospective cohorts identified stage III melanoma patients at highest risk of relapse and has potential to inform adjuvant therapy decisions.

PD-L1 protein expression in breast cancer is rare, enriched in basal-like tumours and associated with infiltrating lymphocytes.

BACKGROUND: Expression of programmed death ligand 1 (PD-L1) in solid tumours has been shown to predict whether patients are likely to respond to anti-PD-L1 therapies. To estimate the therapeutic potential of PD-L1 inhibition in breast cancer, we evaluated the prevalence and significance of PD-L1 protein expression in a large collection of breast tumours. PATIENTS AND METHODS: Correlations between CD274 (PD-L1) copy number, transcript and protein levels were evaluated in tumours from 418 patients recruited to the METABRIC genomic study. Immunohistochemistry was used to detect PD-L1 protein in breast tumours in tissue microarrays from 5763 patients recruited to the SEARCH population-based study (N = 4079) and the NEAT randomised, controlled trial (N = 1684). RESULTS: PD-L1 protein data was available for 3916 of the possible 5763 tumours from the SEARCH and NEAT studies. PD-L1 expression by immune cells was observed in 6% (235/3916) of tumours and expression by tumour cells was observed in just 1.7% (66/3916). PD-L1 was most frequently expressed in basal-like tumours. This was observed both where tumours were subtyped by combined copy number and expression profiling [39% (17/44) of IntClust 10 i.e. basal-like tumours were PD-L1 immune cell positive; P < 0.001] and where a surrogate IHC-based classifier was used [19% (56/302) of basal-like tumours were PD-L1 immune cell positive; P < 0.001]. Moreover, CD274 (PD-L1) amplification was observed in five tumours of which four were IntClust 10. Expression of PD-L1 by either tumour cells or infiltrating immune cells was positively correlated with infiltration by both cytotoxic and regulatory T cells (P < 0.001). There was a nominally significant association between PD-L1 and improved disease-specific survival (hazard ratio 0.53, 95% confidence interval 0.26-1.07; P = 0.08) in ER-negative disease. CONCLUSIONS: Expression of PD-L1 is rare in breast cancer, markedly enriched in basal-like tumours and is correlated with infiltrating lymphocytes. PD-L1 inhibition may benefit the 19% of patients with basal-like tumours in which the protein is expressed. NEAT CLINICALTRIALSGOV: NCT00003577.

Circulating tumor cells and circulating tumor DNA for precision medicine: dream or reality?

Next-generation sequencing studies have provided further evidence to support the notion that cancer is a disease characterized by Darwinian evolution. Today, we often fail to capture this evolution and treatment decisions, even in the metastatic setting, are often based on analysis of primary tumor diagnosed years ago. Currently, this is considered a major reason for treatment failures in cancer care. Recent technological advances in the detection and characterization of circulating tumor cells and circulating tumor DNA might address this and allow for treatment tailoring based on real-time monitoring of tumor evolution. In this review, we summarize the most important recent findings in the field, focusing on challenges and opportunities in moving these tools forward in clinical practice.

Association between CD8+ T-cell infiltration and breast cancer survival in 12,439 patients.

BACKGROUND: T-cell infiltration in estrogen receptor (ER)-negative breast tumours has been associated with longer survival. To investigate this association and the potential of tumour T-cell infiltration as a prognostic and predictive marker, we have conducted the largest study of T cells in breast cancer to date. PATIENTS AND METHODS: Four studies totalling 12 439 patients were used for this work. Cytotoxic (CD8+) and regulatory (forkhead box protein 3, FOXP3+) T cells were quantified using immunohistochemistry (IHC). IHC for CD8 was conducted using available material from all four studies (8978 samples) and for FOXP3 from three studies (5239 samples)-multiple imputation was used to resolve missing data from the remaining patients. Cox regression was used to test for associations with breast cancer-specific survival. RESULTS: In ER-negative tumours [triple-negative breast cancer and human epidermal growth factor receptor 2 (human epidermal growth factor receptor 2 (HER2) positive)], presence of CD8+ T cells within the tumour was associated with a 28% [95% confidence interval (CI) 16% to 38%] reduction in the hazard of breast cancer-specific mortality, and CD8+ T cells within the stroma with a 21% (95% CI 7% to 33%) reduction in hazard. In ER-positive HER2-positive tumours, CD8+ T cells within the tumour were associated with a 27% (95% CI 4% to 44%) reduction in hazard. In ER-negative disease, there was evidence for greater benefit from anthracyclines in the National Epirubicin Adjuvant Trial in patients with CD8+ tumours [hazard ratio (HR) = 0.54; 95% CI 0.37-0.79] versus CD8-negative tumours (HR = 0.87; 95% CI 0.55-1.38). The difference in effect between these subgroups was significant when limited to cases with complete data (P heterogeneity = 0.04) and approached significance in imputed data (P heterogeneity = 0.1). CONCLUSIONS: The presence of CD8+ T cells in breast cancer is associated with a significant reduction in the relative risk of death from disease in both the ER-negative [supplementary Figure S1, available at Annals of Oncology online] and the ER-positive HER2-positive subtypes. Tumour lymphocytic infiltration may improve risk stratification in breast cancer patients classified into these subtypes. NEAT ClinicalTrials.gov: NCT00003577.

Astronomical algorithms for automated analysis of tissue protein expression in breast cancer.

BACKGROUND: High-throughput evaluation of tissue biomarkers in oncology has been greatly accelerated by the widespread use of tissue microarrays (TMAs) and immunohistochemistry. Although TMAs have the potential to facilitate protein expression profiling on a scale to rival experiments of tumour transcriptomes, the bottleneck and imprecision of manually scoring TMAs has impeded progress. METHODS: We report image analysis algorithms adapted from astronomy for the precise automated analysis of IHC in all subcellular compartments. The power of this technique is demonstrated using over 2000 breast tumours and comparing quantitative automated scores against manual assessment by pathologists. RESULTS: All continuous automated scores showed good correlation with their corresponding ordinal manual scores. For oestrogen receptor (ER), the correlation was 0.82, P<0.0001, for BCL2 0.72, P<0.0001 and for HER2 0.62, P<0.0001. Automated scores showed excellent concordance with manual scores for the unsupervised assignment of cases to 'positive' or 'negative' categories with agreement rates of up to 96%. CONCLUSION: The adaptation of astronomical algorithms coupled with their application to large annotated study cohorts, constitutes a powerful tool for the realisation of the enormous potential of digital pathology.

Aurora kinase A outperforms Ki67 as a prognostic marker in ER-positive breast cancer.

BACKGROUND: Proliferation has emerged as a major prognostic factor in luminal breast cancer. The immunohistochemical (IHC) proliferation marker Ki67 has been most extensively investigated but has not gained widespread clinical acceptance. METHODS: We have conducted a head-to-head comparison of a panel of proliferation markers, including Ki67. Our aim was to establish the marker of the greatest prognostic utility. Tumour samples from 3093 women with breast cancer were constructed as tissue microarrays. We used IHC to detect expression of mini-chromosome maintenance protein 2, Ki67, aurora kinase A (AURKA), polo-like kinase 1, geminin and phospho-histone H3. We used a Cox proportional-hazards model to investigate the association with 10-year breast cancer-specific survival (BCSS). Missing values were resolved using multiple imputation. RESULTS: The prognostic significance of proliferation was limited to oestrogen receptor (ER)-positive breast cancer. Aurora kinase A emerged as the marker of the greatest prognostic significance in a multivariate model adjusted for the standard clinical and molecular covariates (hazard ratio 1.3; 95% confidence interval 1.1-1.5; P=0.005), outperforming all other markers including Ki67. CONCLUSION: Aurora kinase A outperforms other proliferation markers as an independent predictor of BCSS in ER-positive breast cancer. It has the potential for use in routine clinical practice.

Denosumab and invasive cervical root resorption: a case report.

This case report describes an adverse side effect from long-term denosumab therapy in a patient with metastatic breast cancer. The patient presented with extensive invasive cervical root resorption affecting most of her dentition. There were no other identified risk factors. As the treatment of breast cancer evolves leading to improved survival, and the number of cancer survivors increases, it is expected that dental practitioners will see a growing number of patients who are on long-term denosumab. Comprehensive dental examinations including radiographs are warranted to identify unusual or unexpected findings such as invasive cervical root resorption. © 2022 Australian Dental Association.

Comparison of methods for handling missing data on immunohistochemical markers in survival analysis of breast cancer.

BACKGROUND: Tissue micro-arrays (TMAs) are increasingly used to generate data of the molecular phenotype of tumours in clinical epidemiology studies, such as studies of disease prognosis. However, TMA data are particularly prone to missingness. A variety of methods to deal with missing data are available. However, the validity of the various approaches is dependent on the structure of the missing data and there are few empirical studies dealing with missing data from molecular pathology. The purpose of this study was to investigate the results of four commonly used approaches to handling missing data from a large, multi-centre study of the molecular pathological determinants of prognosis in breast cancer. PATIENTS AND METHODS: We pooled data from over 11,000 cases of invasive breast cancer from five studies that collected information on seven prognostic indicators together with survival time data. We compared the results of a multi-variate Cox regression using four approaches to handling missing data - complete case analysis (CCA), mean substitution (MS) and multiple imputation without inclusion of the outcome (MI-) and multiple imputation with inclusion of the outcome (MI+). We also performed an analysis in which missing data were simulated under different assumptions and the results of the four methods were compared. RESULTS: Over half the cases had missing data on at least one of the seven variables and 11 percent had missing data on 4 or more. The multi-variate hazard ratio estimates based on multiple imputation models were very similar to those derived after using MS, with similar standard errors. Hazard ratio estimates based on the CCA were only slightly different, but the estimates were less precise as the standard errors were large. However, in data simulated to be missing completely at random (MCAR) or missing at random (MAR), estimates for MI+ were least biased and most accurate, whereas estimates for CCA were most biased and least accurate. CONCLUSION: In this study, empirical results from analyses using CCA, MS, MI- and MI+ were similar, although results from CCA were less precise. The results from simulations suggest that in general MI+ is likely to be the best. Given the ease of implementing MI in standard statistical software, the results of MI+ and CCA should be compared in any multi-variate analysis where missing data are a problem.

Intraperitoneal distribution imaging in ovarian cancer patients.

BACKGROUND/AIMS: Optimal delivery of intraperitoneal (i.p.) chemotherapy is dependent on adequate drug distribution. An accurate understanding of the limitations of i.p. distribution is critical if we are to improve cytotoxic delivery through this route. METHODS: Using repeated scintigraphic peritoneography we investigated peritoneal fluid distribution in patients receiving i.p. therapy. Both early (1-6 h) and late (24-48 h) images were performed. The peritoneal cavity was divided into six regions; pouch of Douglas, left and right paracolic gutters, left and right subphrenic spaces and the right subhepatic space. The distribution in each region was classified as absent (0), faint (1) or intense (2). A total distribution score was calculated from the summation of distribution values for each of the six regions. Distribution was then graded according to the total distribution score as follows: Grade 0 = 0-3; Grade 1 = 4-6; Grade 2 = 7-9; and Grade 3 = 10-12. RESULTS: Twenty-six participants were included in the study: all 26 patients had early imaging performed and 21 of these also had late imaging. Thirteen (50%) and 15 (71%) patients had grade 1 or 2 i.p. distribution on early and late imaging respectively. The most common abdominal regions to show maldistribution were the subphrenic spaces. CONCLUSIONS: This study highlights the deficiencies of distribution following i.p. drug delivery, with the majority of patients demonstrating multiple regions of faint or absent uptake on scintigraphic peritoneography imaging. Future large clinical studies investigating i.p. chemotherapy should include analyses of i.p. distribution to improve our understanding of optimal drug delivery through this route.

BCL2 in breast cancer: a favourable prognostic marker across molecular subtypes and independent of adjuvant therapy received.

BACKGROUND: Breast cancer is heterogeneous and the existing prognostic classifiers are limited in accuracy, leading to unnecessary treatment of numerous women. B-cell lymphoma 2 (BCL2), an antiapoptotic protein, has been proposed as a prognostic marker, but this effect is considered to relate to oestrogen receptor (ER) status. This study aimed to test the clinical validity of BCL2 as an independent prognostic marker. METHODS: Five studies of 11 212 women with early-stage breast cancer were analysed. Individual patient data included tumour size, grade, lymph node status, endocrine therapy, chemotherapy and mortality. BCL2, ER, progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) levels were determined in all tumours. A Cox model incorporating the time-dependent effects of each variable was used to explore the prognostic significance of BCL2. RESULTS: In univariate analysis, ER, PR and BCL2 positivity was associated with improved survival and HER2 positivity with inferior survival. For ER and PR this effect was time dependent, whereas for BCL2 and HER2 the effect persisted over time. In multivariate analysis, BCL2 positivity retained independent prognostic significance (hazard ratio (HR) 0.76, 95% confidence interval (CI) 0.66-0.88, P<0.001). BCL2 was a powerful prognostic marker in ER- (HR 0.63, 95% CI 0.54-0.74, P<0.001) and ER+ disease (HR 0.56, 95% CI 0.48-0.65, P<0.001), and in HER2- (HR 0.55, 95% CI 0.49-0.61, P<0.001) and HER2+ disease (HR 0.70, 95% CI 0.57-0.85, P<0.001), irrespective of the type of adjuvant therapy received. Addition of BCL2 to the Adjuvant! Online prognostic model, for a subset of cases with a 10-year follow-up, improved the survival prediction (P=0.0039). CONCLUSIONS: BCL2 is an independent indicator of favourable prognosis for all types of early-stage breast cancer. This study establishes the rationale for introduction of BCL2 immunohistochemistry to improve prognostic stratification. Further work is now needed to ascertain the exact way to apply BCL2 testing for risk stratification and to standardise BCL2 immunohistochemistry for this application.

Molecular characteristics of screen-detected vs symptomatic breast cancers and their impact on survival.

BACKGROUND: Several recent studies have shown that screen detection remains an independent prognostic factor after adjusting for disease stage at presentation. This study compares the molecular characteristics of screen-detected with symptomatic breast cancers to identify if differences in tumour biology may explain some of the survival benefit conferred by screen detection. METHODS: A total of 1379 women (aged 50-70 years) with invasive breast cancer from a large population-based case-control study were included in the analysis. Individual patient data included tumour size, grade, lymph node status, adjuvant therapy, mammographic screening status and mortality. Immunohistochemistry was performed on tumour samples using 11 primary antibodies to define five molecular subtypes. The effect of screen detection compared with symptomatic diagnosis on survival was estimated after adjustment for grade, nodal status, Nottingham Prognostic Index (NPI) and the molecular markers. RESULTS: Fifty-six per cent of the survival benefit associated with screen-detected breast cancer was accounted for by a shift in the NPI, a further 3-10% was explained by the biological variables and more than 30% of the effect remained unexplained. CONCLUSION: Currently known biomarkers remain limited in their ability to explain the heterogeneity of breast cancer fully. A more complete understanding of the biological profile of breast tumours will be necessary to assess the true impact of tumour biology on the improvement in survival seen with screen detection.

The Brn-3a transcription factor induces neuronal process outgrowth and the coordinate expression of genes encoding synaptic proteins.

The Brn-3a POU family transcription factor is expressed only in posmitotic neurons in the central nervous system and identifies the first differentiated neurons to appear in the midbrain, hindbrain, and spinal cord during development. This factor is also induced when undifferentiated proliferating ND7 cells cease dividing and differentiate to a mature neuronal-like phenotype bearing numerous neurite processes. We show that overexpression of Brn-3a in undifferentiated ND7 cells induces a mature neuronal phenotype characterized by process outgrowth and the induction of genes encoding synaptic proteins, although the cells continue to proliferate. In contrast, the closely related factors Brn-3b and Brn-3c do not have this effect. Although the N-terminal activation domain of Brn-3a is required for maximum induction of neurite outgrowth and gene expression, these effects are primarily dependent on the DNA binding POU domain, which also acts as an activation domain. Overexpression of the isolated POU domain of Brn-3a is sufficient to induce neurite outgrowth, while the ability of full-length Brn-3a to do so is abolished by mutating a single amino acid in the Brn-3a POU homeodomain to its equivalent in Brn-3b. Thus, Brn-3a appears to play a critical role in the specification of the mature neuronal phenotype, acting by stimulating the expression of genes whose products are required for process outgrowth and synapse formation.

Repression of a herpes simplex virus immediate-early promoter by the Oct-2 transcription factor is dependent on an inhibitory region at the N terminus of the protein.

The B-cell form of the Oct-2 transcription factor Oct 2.1 can activate the herpes simplex virus immediate-early gene 3 (IE3) promoter, whereas the neuronally expressed Oct 2.4 and 2.5 forms of the protein, which contain a different C terminus, can repress this promoter. Here we show that partial or full deletion of the C terminus of Oct 2.1 in the presence of an intact N terminus results in a protein which can strongly repress the IE3 promoter. In contrast, deletion of the entire N terminus or a short region within it leaving the C terminus intact results in a very strong activator. Deletion of both N and C termini leaving only the isolated POU domain generates only a very weak repressor. The N-terminal region defined in this way can repress a heterologous promoter when linked to the DNA-binding domain of the GAL4 factor, indicating that it can function as an independent inhibitory domain. These results indicate that a specific region within the N terminus common to Oct 2.1, 2.4, and 2.5 plays a critical role in the ability of neuronally expressed forms of Oct-2 to repress the IE3 promoter but can do so only when the C-terminal region of Oct 2.1 is altered or deleted.

Severe hypoglycaemia after long-acting octreotide in a patient with an unrecognized malignant insulinoma.

Insulinomas are the most common hormone-producing pancreatic neuroendocrine tumours (NETs), and patients usually present with symptoms secondary to hypoglycaemia. Octreotide has been widely used in the symptomatic treatment of patients with pancreatic NETs, including insulinomas. We describe a case of a patient with a metastatic NET, subsequently identified as a malignant insulinoma, who developed severe hypoglycaemia after treatment with long-acting octreotide.

Listeria outbreak associated with sandwich consumption from a hospital retail shop, United Kingdom.

An outbreak of listeriosis occurred in the Swindon area of the UK in autumn 2003. Five cases were detected in pregnant women. Four of these women were thought to have eaten prepacked sandwiches from a retail outlet in one particular hospital. Sampling at the supplier detected Listeria monocytogenes, which was indistinguishable on molecular testing from the patients isolates. Recent changes in UK food legislation should help diminish the risk of further outbreaks/cases such as ours occurring.

Solid Phase Synthesis of Helically Folded Aromatic Oligoamides.

Aromatic amide foldamers constitute a growing class of oligomers that adopt remarkably stable folded conformations. The folded structures possess largely predictable shapes and open the way toward the design of synthetic mimics of proteins. Important examples of aromatic amide foldamers include oligomers of 7- or 8-amino-2-quinoline carboxylic acid that have been shown to exist predominantly as well-defined helices, including when they are combined with α-amino acids to which they may impose their folding behavior. To rapidly iterate their synthesis, solid phase synthesis (SPS) protocols have been developed and optimized for overcoming synthetic difficulties inherent to these backbones such as low nucleophilicity of amine groups on electron poor aromatic rings and a strong propensity of even short sequences to fold on the solid phase during synthesis. For example, acid chloride activation and the use of microwaves are required to bring coupling at aromatic amines to completion. Here, we report detailed SPS protocols for the rapid production of: (1) oligomers of 8-amino-2-quinolinecarboxylic acid; (2) oligomers containing 7-amino-8-fluoro-2-quinolinecarboxylic acid; and (3) heteromeric oligomers of 8-amino-2-quinolinecarboxylic acid and α-amino acids. SPS brings the advantage to quickly produce sequences having varied main chain or side chain components without having to purify multiple intermediates as in solution phase synthesis. With these protocols, an octamer could easily be synthesized and purified within one to two weeks from Fmoc protected amino acid monomer precursors.

Activation of the iNOS gene promoter by Brn-3 POU family transcription factors is dependent upon the octamer motif in the promoter.

The promoter of the gene encoding the inducible nitric oxide synthase (iNOS) contains an octamer motif which is of importance for its activation by specific stimuli. We show that in contrast to the promoter of the neuronal nitric oxide synthase gene (nNOS) which is strongly activated by the Oct-2 octamer-binding POU family transcription factor, the iNOS gene is only weakly activated by Oct-2 via its octamer motif. Unlike the nNOS promoter, however, the iNOS promoter is strongly activated by the POU family transcription factors Brn-3a and Brn-3b. This activation is dependent upon the octamer motif in the iNOS promoter and requires the activation domain located within the POU domain of Brn-3a or Brn-3b but not the N-terminal activation domain of Brn-3a. Thus different but related POU proteins play important roles in the regulation of the genes encoding different forms of nitric oxide synthase.

Two pathways of iron-catalyzed oxidation of bilirubin: effect of desferrioxamine and trolox, and comparison with microsomal oxidation.

The bilirubin-degrading activity of liver microsomes from rats induced with 3-methylcholanthrene has been shown to be markedly stimulated by addition of 3,3',4,4',5,5'-hexabromobiphenyl, a polyhalogenated chemical which resembles in size and shape the most effective inducers of cytochrome P450IA1, but lacks the structural features necessary for it to be metabolised. The degradation of bilirubin by this microsomal system has been compared to oxidation by a chemical model system involving H2O2 and Fe-EDTA (ethylenediaminetetraacetic acid). In both systems bilirubin disappearance was accompanied by bleaching. However, when either desferrioxamine or Trolox were present in the chemical model system, the rate of bilirubin oxidation was greatly enhanced and, at the same time, bilirubin was largely or entirely converted to biliverdin, a pathway of oxidation which proceeds by dehydrogenation. In the presence of desferrioxamine, biliverdin was also further oxidised to an unidentified red pigment.

The different inhibitory domains of the Oct-2 transcription factor have distinct functional activities.

The Oct-2 POU family transcription factor contains three distinct regions whose deletion reduces its ability to inhibit transcription via its octamer binding site. Here we show that only one of these inhibitory domains is capable of also inhibiting the activity of activating molecules bound at adjacent sites upstream of a TATA box-containing promoter whereas the other two regions are inactive in this assay. None of the three regions is able to achieve this effect when located upstream of the same promoter containing an initiator motif. The mechanisms of action of these domains and their role in the functioning of the Oct-2 factor are discussed.