Search

There is now evidence that gene fusions activating the MAPK pathway are relatively common in pancreatic acinar cell carcinoma with potentially actionable BRAF or RET fusions being found in ~30%. We sought to investigate the incidence of RAF1 fusions in pancreatic malignancies with acinar cell differentiation. FISH testing for RAF1 was undertaken on 30 tumors comprising 25 'pure' acinar cell carcinomas, 2 mixed pancreatic acinar-neuroendocrine carcinomas, 1 mixed acinar cell-low grade neuroendocrine tumor and 2 pancreatoblastomas. RAF1 rearrangements were identified in 5 cases and confirmed by DNA and RNA sequencing to represent oncogenic fusions (GATM-RAF1, GOLGA4-RAF1, PDZRN3-RAF1, HERPUD1-RAF1 and TRIM33-RAF1) and to be mutually exclusive with BRAF and RET fusions, as well as KRAS mutations. Large genome-wide copy number changes were common and included 1q gain and/or 1p loss in all five RAF1 FISH-positive acinar cell carcinomas. RAF1 expression by immunohistochemistry was found in 3 of 5 (60%) of fusion-positive cases and no FISH-negative cases. Phospho-ERK1/2 expression was found in 4 of 5 RAF1-fusion-positive cases. Expression of both RAF1 and phospho-ERK1/2 was heterogeneous and often only detected at the tumor-stroma interface, thus limiting their clinical utility. We conclude that RAF1 gene rearrangements are relatively common in pancreatic acinar cell carcinomas (14.3% to 18.5% of cases) and can be effectively identified by FISH with follow up molecular testing. The combined results of several studies now indicate that BRAF, RET or RAF1 fusions occur in between one third and one-half of these tumors but are extremely rare in other pancreatic malignancies. As these fusions are potentially actionable with currently available therapies, a strong argument can be made to perform FISH or molecular testing on all pancreatic acinar cell carcinomas.

Melanoma is a deadly form of skin cancer that accounts for a disproportionally large proportion of cancer-related deaths in younger people. Compared with most other skin cancers, a feature of melanoma is its high metastatic capacity, although the mechanisms that confer this are not well understood. The Hippo pathway is a key regulator of organ growth and cell fate that is deregulated in many cancers. To analyse the Hippo pathway in cutaneous melanoma, we generated a transcriptional signature of melanoma cells that overexpressed YAP, the key downstream Hippo pathway oncoprotein. YAP-mediated transcriptional activity varied in melanoma cell lines but did not cluster with known genetic drivers of melanomagenesis such as BRAF and NRAS mutations. Instead, it correlated strongly with published gene expression profiles linked to melanoma cell invasiveness and varied throughout the metastatic cascade in melanoma patient tumours. Consistent with this, YAP was both necessary and sufficient for melanoma cell invasion in vitro. In vivo, YAP promoted spontaneous melanoma metastasis, whilst the growth of YAP-expressing primary tumours was impeded. Finally, we identified the YAP target genes AXL, THBS1 and CYR61 as key mediators of YAP-induced melanoma cell invasion. These data suggest that YAP is a critical regulator of melanoma metastasis.

Human immunodeficiency virus 1 (HIV-1) is a retrovirus with a ten-kilobase single-stranded RNA genome. HIV-1 must express all of its gene products from a single primary transcript, which undergoes alternative splicing to produce diverse protein products that include structural proteins and regulatory factors1,2. Despite the critical role of alternative splicing, the mechanisms that drive the choice of splice site are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements3. Here we use dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) to investigate the structure of HIV-1 RNA in cells, and develop an algorithm that we name 'detection of RNA folding ensembles using expectation-maximization' (DREEM), which reveals the alternative conformations that are assumed by the same RNA sequence. Contrary to previous models that have analysed population averages4, our results reveal heterogeneous regions of RNA structure across the entire HIV-1 genome. In addition to confirming that in vitro characterized5 alternative structures for the HIV-1 Rev responsive element also exist in cells, we discover alternative conformations at critical splice sites that influence the ratio of transcript isoforms. Our simultaneous measurement of splicing and intracellular RNA structure provides evidence for the long-standing hypothesis6-8 that heterogeneity in RNA conformation regulates splice-site use and viral gene expression.

Purpose Circulating tumor DNA (ctDNA) allows noninvasive disease monitoring across a range of malignancies. In metastatic melanoma, the extent to which ctDNA reflects changes in metabolic disease burden assessed by 18F-labeled fluorodeoxyglucose positron emission tomography (FDG-PET) is unknown. We assessed the role of ctDNA analysis in combination with FDG-PET to monitor tumor burden and genomic heterogeneity throughout treatment. Patients and Methods We performed a comprehensive analysis of serial ctDNA and FDG-PET in 52 patients who received systemic therapy for metastatic melanoma. Next-generation sequencing and digital polymerase chain reaction were used to analyze plasma samples from the cohort. Results ctDNA levels were monitored across patients with mutant BRAF, NRAS, and BRAF/NRAS wild type disease. Mutant BRAF and NRAS ctDNA levels correlated closely with changes in metabolic disease burden throughout treatment. TERT promoter mutant ctDNA levels also paralleled changes in tumor burden, which provide an alternative marker for disease monitoring. Of note, subcutaneous and cerebral disease sites were not well represented in plasma. Early changes in ctDNA and metabolic disease burden were important indicators of treatment response. Patients with an early decrease in ctDNA post-treatment had improved progression-free survival compared with patients in whom ctDNA levels remained unchanged or increased over time (hazard ratio, 2.6; P = .05). ctDNA analysis contributed key molecular information through the identification of putative resistance mechanisms to targeted therapy. A detailed comparison of the genomic architecture of plasma and multiregional tumor biopsy specimens at autopsy revealed the ability of ctDNA to comprehensively capture genomic heterogeneity across multiple disease sites. Conclusion The findings highlight the powerful role of ctDNA in metastatic melanoma as a complementary modality to functional imaging that allows real-time monitoring of both tumor burden and genomic changes throughout therapy.

Small-molecule inhibitors of PfATP4, a Plasmodium falciparum protein that is believed to pump Na+ out of the parasite while importing H+, are on track to become much-needed new antimalarial drugs. The spiroindolone cipargamin is poised to become the first PfATP4 inhibitor to reach the field, having performed strongly in Phase 1 and 2 clinical trials. Previous attempts to generate cipargamin-resistant parasites in the laboratory have yielded parasites with reduced susceptibility to the drug; however, the highest 50% inhibitory concentration reported to date is 24 nM. Here, we show that P. falciparum parasites can acquire a clinically-significant level of resistance to cipargamin that enables them to withstand micromolar concentrations of the drug. Independent experiments to generate high-level cipargamin resistance using different protocols and strains led to the same change each time – a G358S mutation in PfATP4. Parasites with this mutation showed high-level resistance not only to cipargamin, but also to the dihydroisoquinolone (+)-SJ733. However, for certain other (less clinically advanced) PfATP4-associated compounds the G358S mutation in PfATP4 conferred only moderate resistance or no resistance. The G358S mutation in PfATP4 did not affect parasite susceptibility to antimalarials that do not target PfATP4. The G358S mutation in PfATP4, and the equivalent mutation in the Toxoplasma gondii ATP4 homologue (G419S), decreased the sensitivity of the Na+-ATPase activity of ATP4 to inhibition by cipargamin and (+)-SJ733, and decreased the sensitivity of parasites expressing these ATP4 mutations to disruption of parasite Na+ regulation by cipargamin- and (+)-SJ733. The G358S mutation in PfATP4 reduced the affinity of the protein for Na+ and was associated with an increase in the parasite’s resting cytosolic Na+ concentration; however, no significant defect in parasite growth rate was observed. Our findings suggest that codon 358 in pfatp4 should be monitored closely in the field as a molecular marker for cipargamin resistance, and that PfATP4 inhibitors in clinical development should be tested for their activity against PfATP4G358S parasites.

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), a diverse family of multi-domain proteins expressed on the surface of malaria-infected erythrocytes, is an important target of protective immunity against malaria. Our group recently studied transcription of the var genes encoding PfEMP1 in individuals from Papua, Indonesia with severe or uncomplicated malaria. We cloned and expressed domains from 32 PfEMP1s including 22 that were upregulated in severe malaria and 10 that were upregulated in uncomplicated malaria, using a wheat germ cell-free expression system. We used Luminex technology to measure IgG antibodies to these 32 domains and control proteins in 63 individuals (11 children). At presentation to hospital, levels of antibodies to PfEMP1 domains were either higher in uncomplicated malaria or were not significantly different between groups. Using principal components analysis, antibodies to three of 32 domains were highly discriminatory between groups. These included two domains upregulated in severe malaria, a DBLβ13 domain and a CIDRα1.6 domain (which has been previously implicated in severe malaria pathogenesis), and a DBLδ domain that was upregulated in uncomplicated malaria. Antibody to control non-PfEMP1 antigens did not differ with disease severity. Antibodies to PfEMP1 domains differ with malaria severity. Lack of antibodies to locally expressed PfEMP1 types, including both domains previously associated with severe malaria and newly identified targets, may in part explain malaria severity in Papuan adults.

The accurate identification and effective removal of unwanted variation are essential to derive meaningful biological results from RNA-seq data, especially when the data come from large and complex studies. We have used The Cancer Genome Atlas (TCGA) RNA-seq data to show that library size, batch effects, and tumor purity are major sources of unwanted variation across all TCGA RNA-seq datasets and that existing gold standard approaches to normalizations fail to remove this unwanted variation. Additionally, we illustrate how different sources of unwanted variation can compromise downstream analyses, including gene co-expression, association between gene expression and survival outcomes, and cancer subtype identifications. Here, we propose the use of a novel strategy, pseudo-replicates of pseudo-samples (PRPS), to deploy the Removing Unwanted Variation III (RUV-III) method to remove different sources of unwanted variation from large and complex gene expression studies. Our approach requires at least one roughly known biologically homogenous subclass of samples shared across sources of unwanted variation. To create PRPS, we first need to identify the sources of unwanted variation, which we will call batches in the data. Then the gene expression measurements of biologically homogeneous sets of samples are averaged within batches, and the results called pseudo-samples. Pseudo-samples with the same biology and different batches are then defined to be pseudo-replicates and used in RUV-III as replicates. The variation between pseudo-samples of a set pseudo-replicates is mainly unwanted variation. We illustrate the value of our approach by comparing it to the TCGA normalizations on several TCGA RNA-seq datasets. RUV-III with PRPS can be used for any large genomics project involving multiple labs, technicians, or platforms.

Studies of the gut microbiome have focused on its bacterial composition. We aimed to characterize the gut fungal microbiome (mycobiome) across pregnancy in women with and without type 1 diabetes.Faecal samples (n = 162) were collected from 70 pregnant women (45 with and 25 without type 1 diabetes) across all trimesters. Fungi were analysed by internal transcribed spacer 1 amplicon sequencing. Markers of intestinal inflammation (faecal calprotectin) and intestinal epithelial integrity (serum intestinal fatty acid binding protein; I-FABP), and serum antibodies to Saccharomyces cerevisiae (ASCA) were measured.Women with type 1 diabetes had decreased fungal alpha diversity by the third trimester, associated with an increased abundance of Saccharomyces cerevisiae that was inversely related to the abundance of the anti-inflammatory butyrate-producing bacterium Faecalibacterium prausnitzii. Women with type 1 diabetes had higher concentrations of calprotectin, I-FABP and ASCA.Women with type 1 diabetes exhibit a shift in the gut mycobiome across pregnancy associated with evidence of gut inflammation and impaired intestinal barrier function. The relevance of these findings to the higher rate of pregnancy complications in type 1 diabetes warrants further study.

Summary StructuralVariantAnnotation is an R/Bioconductor package that provides a framework for decoupling downstream analysis of structural variant breakpoints from upstream variant calling methods. It standardizes the representational format from BEDPE, or any of the three different notations supported by VCF into a breakpoint GRanges data structure suitable for use by the wider Bioconductor ecosystem. It handles both transitive breakpoints and duplication/insertion notational differences of identical variants—both common scenarios when comparing short/long read-based call sets that confound downstream analysis. StructuralVariantAnnotation provides the caller-agnostic foundation needed for a R/Bioconductor ecosystem of structural variant annotation, classification and interpretation tools able to handle both simple and complex genomic rearrangements. Availability and implementation StructuralVariantAnnotation is implemented in R and available for download as the Bioconductor StructuralVariantAnnotation package. Details can be found at https://www.bioconductor.org/packages/release/bioc/html/StructuralVariantAnnotation.html. It has been released under a GPL license.

The nucleoporin ELYS, encoded by AHCTF1, is a large multifunctional protein with essential roles in nuclear pore assembly and mitosis. Using a zebrafish model of hepatocellular carcinoma, in which the expression of an inducible mutant kras transgene (krasG12V) drives hepatocyte-specific hyperplasia and liver enlargement, we show that reducing ahctf1 gene dosage by 50% markedly shrinks tumour burden, while non-hyperplastic tissues are unaffected. We demonstrate that ahctf1 heterozygosity impairs nuclear pore formation, mitotic spindle assembly and chromosome segregation, leading to DNA damage and activation of TP53-dependent and independent mechanisms of cell death and cell cycle arrest. This selective vulnerability of cancer cells to mild disruption of Elys function uncovers a novel synthetic lethal interaction between ahctf1 and kras mutations that could be exploited therapeutically. Heterozygous expression of both ahctf1 and ranbp2, or treatment of heterozygous ahctf1 larvae with the nucleocytoplasmic transport inhibitor, Selinexor, completely blocked krasG12V-driven hepatocyte hyperplasia, revealing promising avenues for combinatorial treatments.

BackgroundThe biological significance of structural variation is now more widely recognized. However, due to the lack of available tools for downstream analysis, including processing and annotating, interpretation of structural variant calls remains a challenge.FindingsHere we present svaRetro and svaNUMT, R packages that provide functions for annotating novel genomic events such as non-reference retro-copied transcripts and nuclear integration of mitochondrial DNA. We evaluate the performance of these packages to detect events using simulations and public benchmarking datasets, and annotate processed transcripts in a public structural variant database.ConclusionssvaRetro and svaNUMT provide efficient, modular tools for downstream identification and annotation of structural variant calls.

Background The gut microbiome changes in response to a range of environmental conditions, life events and disease states. Pregnancy is a natural life event that involves major physiological adaptation yet studies of the microbiome in pregnancy are limited and their findings inconsistent. Pregnancy with type 1 diabetes (T1D) is associated with increased maternal and fetal risks but the gut microbiome in this context has not been characterized. By whole metagenome sequencing (WMS), we defined the taxonomic composition and function of the gut bacterial microbiome across 70 pregnancies, 36 in women with T1D. Results Women with and without T1D exhibited compositional and functional changes in the gut microbiome across pregnancy. Profiles in women with T1D were distinct, with an increase in bacteria that produce lipopolysaccharides and a decrease in those that produce short-chain fatty acids, especially in the third trimester. In addition, women with T1D had elevated concentrations of fecal calprotectin, a marker of intestinal inflammation, and serum intestinal fatty acid-binding protein (I-FABP), a marker of intestinal epithelial damage. Conclusions Women with T1D exhibit a shift towards a more pro-inflammatory gut microbiome during pregnancy, associated with evidence of intestinal inflammation. These changes could contribute to the increased risk of pregnancy complications in women with T1D and are potentially modifiable by dietary means.

Background Prostate cancer is caused by genomic aberrations in normal epithelial cells, however clinical translation of findings from analyses of cancer cells alone has been very limited. A deeper understanding of the tumour microenvironment is needed to identify the key drivers of disease progression and reveal novel therapeutic opportunities. Results In this study, the experimental enrichment of selected cell-types, the development of a Bayesian inference model for continuous differential transcript abundance, and multiplex immunohistochemistry permitted us to define the transcriptional landscape of the prostate cancer microenvironment along the disease progression axis. An important role of monocytes and macrophages in prostate cancer progression and disease recurrence was uncovered, supported by both transcriptional landscape findings and by differential tissue composition analyses. These findings were corroborated and validated by spatial analyses at the single-cell level using multiplex immunohistochemistry. Conclusions This study advances our knowledge concerning the role of monocyte-derived recruitment in primary prostate cancer, and supports their key role in disease progression, patient survival and prostate microenvironment immune modulation.

Background There is a clear need for novel approaches to malaria vaccine development. We aimed to develop a genetically attenuated blood-stage vaccine and test its safety, infectivity, and immunogenicity in healthy volunteers. Our approach was to target the gene encoding the knob-associated histidine-rich protein (KAHRP), which is responsible for the assembly of knob structures at the infected erythrocyte surface. Knobs are required for correct display of the polymorphic adhesion ligand P. falciparum erythrocyte membrane protein 1 (PfEMP1), a key virulence determinant encoded by a repertoire of var genes. Methods The gene encoding KAHRP was deleted from P. falciparum 3D7 and a master cell bank was produced in accordance with Good Manufacturing Practice. Eight malaria naïve males were intravenously inoculated (day 0) with 1800 (2 subjects), 1.8 × 105 (2 subjects), or 3 × 106 viable parasites (4 subjects). Parasitemia was measured using qPCR; immunogenicity was determined using standard assays. Parasites were rescued into culture for in vitro analyses (genome sequencing, cytoadhesion assays, scanning electron microscopy, var gene expression). Results None of the subjects who were administered with 1800 or 1.8 × 105 parasites developed parasitemia; 3/4 subjects administered 3× 106 parasites developed significant parasitemia, first detected on days 13, 18, and 22. One of these three subjects developed symptoms of malaria simultaneously with influenza B (day 17; 14,022 parasites/mL); one subject developed mild symptoms on day 28 (19,956 parasites/mL); and one subject remained asymptomatic up to day 35 (5046 parasites/mL). Parasitemia rapidly cleared with artemether/lumefantrine. Parasitemia induced a parasite-specific antibody and cell-mediated immune response. Parasites cultured ex vivo exhibited genotypic and phenotypic properties similar to inoculated parasites, although the var gene expression profile changed during growth in vivo. Conclusions This study represents the first clinical investigation of a genetically attenuated blood-stage human malaria vaccine. A P. falciparum 3D7 kahrp– strain was tested in vivo and found to be immunogenic but can lead to patent parasitemia at high doses. Trial registration Australian New Zealand Clinical Trials Registry (number: ACTRN12617000824369; date: 06 June 2017).

Optical barcoding reveals enhanced TNFα activation and polyclonality in lung, but not liver, metastases from breast cancer., Intratumoral heterogeneity is a driver of breast cancer progression, but the nature of the clonal interactive network involved in this process remains unclear. Here, we optimized the use of optical barcoding to visualize and characterize 31 cancer subclones in vivo. By mapping the clonal composition of thousands of metastases in two clinically relevant sites, the lungs and liver, we found that metastases were highly polyclonal in lungs but not in the liver. Furthermore, the transcriptome of the subclones varied according to their metastatic niche. We also identified a reversible niche-driven signature that was conserved in lung and liver metastases collected during patient autopsies. Among this signature, we found that the tumor necrosis factor–α pathway was up-regulated in lung compared to liver metastases, and inhibition of this pathway affected metastasis diversity. These results highlight that the cellular and molecular heterogeneity observed in metastases is largely dictated by the tumor microenvironment.

Background The human microbiome plays an important role in cancer. Accumulating evidence indicates that commensal microbiome-derived DNA may be represented in minute quantities in the cell-free DNA of human blood and could possibly be harnessed as a new cancer biomarker. However, there has been limited use of rigorous experimental controls to account for contamination, which invariably affects low-biomass microbiome studies. Results We apply a combination of 16S-rRNA-gene sequencing and droplet digital PCR to determine if the specific detection of cell-free microbial DNA (cfmDNA) is possible in metastatic melanoma patients. Compared to matched stool and saliva samples, the absolute concentration of cfmDNA is low but significantly above the levels detected from negative controls. The microbial community of plasma is strongly influenced by laboratory and reagent contaminants introduced during the DNA extraction and sequencing processes. Through the application of an in silico decontamination strategy including the filtering of amplicon sequence variants (ASVs) with batch dependent abundances and those with a higher prevalence in negative controls, we identify known gut commensal bacteria, such as Faecalibacterium, Bacteroides and Ruminococcus, and also other uncharacterised ASVs. We analyse additional plasma samples, highlighting the potential of this framework to identify differences in cfmDNA between healthy and cancer patients. Conclusions Together, these observations indicate that plasma can harbour a low yet detectable level of cfmDNA. The results highlight the importance of accounting for contamination and provide an analytical decontamination framework to allow the accurate detection of cfmDNA for future biomarker studies in cancer and other diseases.

Significance Targeting CDK4/6 shows great promise in the treatment of many solid cancers, including melanoma. This study has uncovered a mechanism of action of CDK4/6 inhibitors in regulating the MDM4 oncogene and the tumor suppressor, p53. In melanoma, palbociclib activates p53 via modulating PRMT5-dependent alternate MDM4 pre-mRNA splicing, which results in a robust decrease in MDM4 protein expression. Loss of palbociclib regulation of the PRMT5–MDM4 axis leads to drug resistance. Dual inhibition of CDK4/6 and PRMT5 potently suppresses melanoma tumor growth and is well tolerated in vivo. Our findings not only have immediate implications for the advancement of CDK4/6 inhibition for treating melanoma, but also more generally offer insights into CDK4/6 mechanism of action., Cyclin-dependent kinase 4/6 (CDK4/6) inhibitors are an established treatment in estrogen receptor-positive breast cancer and are currently in clinical development in melanoma, a tumor that exhibits high rates of CDK4 activation. We analyzed melanoma cells with acquired resistance to the CDK4/6 inhibitor palbociclib and demonstrate that the activity of PRMT5, a protein arginine methyltransferase and indirect target of CDK4, is essential for CDK4/6 inhibitor sensitivity. By indirectly suppressing PRMT5 activity, palbociclib alters the pre-mRNA splicing of MDM4, a negative regulator of p53, leading to decreased MDM4 protein expression and subsequent p53 activation. In turn, p53 induces p21, leading to inhibition of CDK2, the main kinase substituting for CDK4/6 and a key driver of resistance to palbociclib. Loss of the ability of palbociclib to regulate the PRMT5–MDM4 axis leads to resistance. Importantly, combining palbociclib with the PRMT5 inhibitor GSK3326595 enhances the efficacy of palbociclib in treating naive and resistant models and also delays the emergence of resistance. Our studies have uncovered a mechanism of action of CDK4/6 inhibitors in regulating the MDM4 oncogene and the tumor suppressor, p53. Furthermore, we have established that palbociclib inhibition of the PRMT5–MDM4 axis is essential for robust melanoma cell sensitivity and provide preclinical evidence that coinhibition of CDK4/6 and PRMT5 is an effective and well-tolerated therapeutic strategy. Overall, our data provide a strong rationale for further investigation of novel combinations of CDK4/6 and PRMT5 inhibitors, not only in melanoma but other tumor types, including breast, pancreatic, and esophageal carcinoma.

Devil facial tumour disease (DFTD) comprises two genetically distinct transmissible cancers (DFT1 and DFT2) endangering the survival of the Tasmanian devil (Sarcophilus harrisii) in the wild. DFT1 first arose from a cell of the Schwann cell lineage; however, the tissue-of-origin of the recently discovered DFT2 cancer is unknown. In this study, we compared the transcriptome and proteome of DFT2 tumours to DFT1 and normal Tasmanian devil tissues to determine the tissue-of-origin of the DFT2 cancer. Our findings demonstrate that DFT2 expresses a range of Schwann cell markers and exhibits expression patterns consistent with a similar origin to the DFT1 cancer. Furthermore, DFT2 cells express genes associated with the repair response to peripheral nerve damage. These findings suggest that devils may be predisposed to transmissible cancers of Schwann cell origin. The combined effect of factors such as frequent nerve damage from biting, Schwann cell plasticity and low genetic diversity may allow these cancers to develop on rare occasions. The emergence of two independent transmissible cancers from the same tissue in the Tasmanian devil presents an unprecedented opportunity to gain insight into cancer development, evolution and immune evasion in mammalian species.

Melanoma is usually driven by mutations in BRAF or NRAS, which trigger hyperactivation of MAPK signaling. However, MAPK-targeted therapies are not sustainably effective in most patients. Accordingly, characterizing mechanisms that co-operatively drive melanoma progression is key to improving patient outcomes. One possible mechanism is the Hippo signaling pathway, which regulates cancer progression via its central oncoproteins YAP and TAZ, although is thought to be only rarely affected by direct mutation. As YAP hyperactivation occurs in uveal melanoma, we investigated this oncogene in cutaneous melanoma. YAP protein expression was elevated in most benign nevi and primary cutaneous melanomas but present at only very low levels in normal melanocytes. In patient-derived xenografts and melanoma cell lines, we observed variable reliance of cell viability on Hippo pathway signaling that was independent of TAZ activity and also of classical melanoma driver mutations such as BRAF and NRAS. Finally, in genotyping studies of melanoma, we observed the first ever hyperactivating YAP mutations in a human cancer, manifest as seven distinct missense point mutations that caused serine to alanine transpositions. Strikingly, these mutate four serine residues known to be targeted by the Hippo pathway and we show that they lead to hyperactivation of YAP. Implications: Our studies highlight the YAP oncoprotein as a potential therapeutic target in select subgroups of melanoma patients, although successful treatment with anti-YAP therapies will depend on identification of biomarkers additional to YAP protein expression.

Prostate cancer is a leading cause of morbidity and cancer-related death worldwide. Androgen deprivation therapy (ADT) is the cornerstone of management for advanced disease. The use of these therapies is associated with multiple side effects, including metabolic syndrome and truncal obesity. At the same time, obesity has been associated with both prostate cancer development and disease progression, linked to its effects on chronic inflammation at a tissue level. The connection between ADT, obesity, inflammation and prostate cancer progression is well established in clinical settings; however, an understanding of the changes in adipose tissue at the molecular level induced by castration therapies is missing. Here, we investigated the transcriptional changes in periprostatic fat tissue induced by profound ADT in a group of patients with high-risk tumours compared to a matching untreated cohort. We find that the deprivation of androgen is associated with a pro-inflammatory and obesity-like adipose tissue microenvironment. This study suggests that the beneficial effect of therapies based on androgen deprivation may be partially counteracted by metabolic and inflammatory side effects in the adipose tissue surrounding the prostate.

Many of the small, cysteine-rich ion-channel modulatory peptides found in Cnidaria are distantly related to vertebrate defensins (of the trans-defensin superfamily). Transcriptomic and proteomic studies of the endemic Australian speckled sea anemone (Oulactis sp.) yielded homologous peptides to known defensin sequences. We extended these data using existing and custom-built hidden Markov models to extract defensin-like families from the transcriptomes of seven endemic Australian cnidarian species. Newly sequenced transcriptomes include three species of Actiniaria (true sea anemones); the speckled anemone (Oulactis sp.), Oulactis muscosa, Dofleinia cf. armata and a species of Corallimorpharia, Rhodactis sp. We analyzed these novel defensin-like sequences along with published homologues to study the evolution of their physico-chemical properties in vertebrate and invertebrate fauna. The cnidarian trans-defensins form a distinct cluster within the chemical space of the superfamily, with a unique set of motifs and biophysical properties. This cluster contains identifiable subgroups, whose distribution in chemical space also correlates with the divergent evolution of their structures. These sequences, currently restricted to cnidarians, form an evolutionarily distinct clade within the trans-defensin superfamily.

The Tasmanian devil (Sarcophilus harrisii) is threatened by a contagious cancer, known as Devil Facial Tumour Disease (DFTD). A highly diverse T-cell receptor (TCR) repertoire is crucial for successful host defence against cancers. By investigating TCR beta chain diversity in devils of different ages, we show that the T-cell repertoire in devils constricts in their second year of life, which may explain the higher DFTD prevalence in older devils. Unexpectedly, we also observed a pronounced decline in TCR diversity and T cell clonal expansion in devils after DFTD infection. These findings overturned the previous assumption that DFTD did not directly impact host immunity., Yuanyuan Cheng et al. showed that the T-cell repertoire diversity of Tasmanian devils diminishes during their second year of life which may explain the prevalence of devil facial tumor disease in older devils. Infection with this disease also impacts T-cell diversity highlighting a previously unknown effect of the devil facial tumor disease on host immunity.

Malignant pleural effusions (MPE) are a common complication of advanced cancers, particularly those adjacent to the pleura such as lung and breast cancer and are a frequent complication in metastatic disease. The pathophysiology of MPE formation in advanced breast cancer remains poorly understood, and their composition and biology are understudied. To characterise the phenotypic diversity of malignant pleural effusion, we performed single-cell RNA sequencing on 10 MPEs from 7 metastatic breast cancer patients with diverse molecular subtypes: two triple negative (TNBC) patients, three luminal B patients including one with a rare inflammatory subtype, and two luminal A patients. For all patients, we sequenced cells from the entire MPE, without performing any enrichment or selection, in order to ascertain the cellular composition and molecular phenotypes in an unbiased manner. We identified pronounced differences in the relative proportions of malignant, mesothelial and immune cell populations: both TNBC and two Luminal B patients had extensive malignant cell populations, while the other three patients had few or no detectible malignant cells but extensive immune populations. We also observed heterogeneity in the expression of subtype-specific gene signatures and in copy number aberration patterns that captured variability across malignant cell populations both within and between patients. We observed that most malignant cells retained the molecular subtype diagnosed in the primary tumour, however sub-populations of malignant cells could be identified that mapped to different molecular subtypes, indicative of complexity in the molecular phenotypes of cancer cells disseminated to or metastasising in the pleural cavity. We distinguished mesothelial cell populations from malignant cells using key markers, expression of the PAM50 signatures, and copy number aberration patterns. We found that pleural mesothelial cells expressed a cancer associated fibroblast-like transcriptomic program that may support cancer growth within the pleural cavity through the secretion of growth factors that target malignant cells. Our dataset presents the first unbiased and unselected assessment of breast cancer associated MPEs at single cell resolution, providing the community with a vital resource for the study of MPEs. Our work highlights the molecular and cellular diversity captured in MPEs and advances the use of these clinically relevant biopsies both in monitoring disease progression and in the development of targeted therapeutics for patients with advanced breast cancer. Citation Format: Holly J. Whitfield, Jean Berthelet, Stefano Mangiola, Caroline Bell, Robin L. Anderson, Bhupinder Pal, Anthony T. Papenfuss, Belinda Yeo, Delphine Merino, Melissa J. Davis. Defining the cellular composition and associated molecular phenotypes of malignant and non-malignant cells in breast cancer pleural effusions [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6232.

Ovarian cancer, in particularly high-grade serous ovarian cancer (HGSOC) and ovarian carcinosarcoma (OCS), are highly aggressive and deadly female cancers with limited treatment options. These tumors are generally unresponsive to immune check-point inhibitor (ICI) therapy and are referred to as immunologically “cold” tumors. Cell-based therapy, in particular, adoptive T-cell therapy, is an alternative immunotherapy option that has shown great potential, especially chimeric antigen receptor T cell (CAR-T) therapy in the treatment of hematologic malignancies. However, the efficacy of CAR-T therapy in solid tumors has been modest. This review explores the potential of another cell-based therapy, T-cell receptor therapy (TCR-T) as an alternate treatment option for immunological “cold” OC and OCS tumors.

Epidemiological studies link Sarcoptes scabiei infection and impetigo. Scabies mites can promote Streptococcus pyogenes (Group A Streptococcus) and Staphylococcus aureus infections by breaching the skin barrier and excreting molecules that inhibit host innate immune responses. However, little is known about the composition and the function of the scabies-associated microbiota. Here, high-throughput whole-metagenome sequencing was used to explore the scabies-associated microbiome. Scabies mites including their immediate microenvironments were isolated from two patients with severe scabies in Northern Australia. Two ~45–50 million paired-end reads Illumina libraries were generated of which ~2 (5.1%) and 0.7 million (1.3%) microbial reads were filtered out by mapping to human (hg19) and mite draft genomes. Taxonomic profiling revealed a microbial community dominated by the phylum Firmicutes (A: 79% and B: 59%) and genera that comprise Streptococcus, Staphylococcus, Acinetobacter, and Corynebacterium. Assembly of the metagenome reads resulted in genome bins representing reference genomes of Acinetobacter baumannii, Streptococcus dysgalactiae (Group C/G), Proteus mirablis and Staphylococcus aureus. The contigs contained genes relevant to pathogenicity and antibiotics resistance. Confocal microscopy of a patient skin sample confirmed A. baumannii, Streptococci and S. aureus in scabies mite gut and faeces and the surrounding skin. The study provides fundamental evidence for the association of opportunistic pathogens with scabies infection.

To examine global changes in breast heterogeneity across different states, we determined the single‐cell transcriptomes of > 340,000 cells encompassing normal breast, preneoplastic BRCA1 +/– tissue, the major breast cancer subtypes, and pairs of tumors and involved lymph nodes. Elucidation of the normal breast microenvironment revealed striking changes in the stroma of post‐menopausal women. Single‐cell profiling of 34 treatment‐naive primary tumors, including estrogen receptor (ER)+, HER2+, and triple‐negative breast cancers, revealed comparable diversity among cancer cells and a discrete subset of cycling cells. The transcriptomes of preneoplastic BRCA1 +/– tissue versus tumors highlighted global changes in the immune microenvironment. Within the tumor immune landscape, proliferative CD8+ T cells characterized triple‐negative and HER2+ cancers but not ER+ tumors, while all subtypes comprised cycling tumor‐associated macrophages, thus invoking potentially different immunotherapy targets. Copy number analysis of paired ER+ tumors and lymph nodes indicated seeding by genetically distinct clones or mass migration of primary tumor cells into axillary lymph nodes. This large‐scale integration of patient samples provides a high‐resolution map of cell diversity in normal and cancerous human breast., A large‐scale gene expression resource integrates diverse tissue samples and reveals unexpected heterogeneity of breast cancer subtypes.

Although melanoma is initiated by acquisition of point mutations and limited focal copy number alterations in melanocytes-of-origin, the nature of genetic changes that characterise lethal metastatic disease is poorly understood. Here, we analyze the evolution of human melanoma progressing from early to late disease in 13 patients by sampling their tumours at multiple sites and times. Whole exome and genome sequencing data from 88 tumour samples reveals only limited gain of point mutations generally, with net mutational loss in some metastases. In contrast, melanoma evolution is dominated by whole genome doubling and large-scale aneuploidy, in which widespread loss of heterozygosity sculpts the burden of point mutations, neoantigens and structural variants even in treatment-naïve and primary cutaneous melanomas in some patients. These results imply that dysregulation of genomic integrity is a key driver of selective clonal advantage during melanoma progression., The genetic changes that occur in late stage metastatic melanoma are not well delineated. Here, the authors use rapid autopsy samples from metastatic melanoma patients and show that the late stage in the disease is characterised by whole genome doubling and aneuploidy.

Relative transcript abundance has proven to be a valuable tool for understanding the function of genes in biological systems. For the differential analysis of transcript abundance using RNA sequencing data, the negative binomial model is by far the most frequently adopted. However, common methods that are based on a negative binomial model are not robust to extreme outliers, which we found to be abundant in public datasets. So far, no rigorous and probabilistic methods for detection of outliers have been developed for RNA sequencing data, leaving the identification mostly to visual inspection. Recent advances in Bayesian computation allow large-scale comparison of observed data against its theoretical distribution given in a statistical model. Here we propose ppcseq, a key quality-control tool for identifying transcripts that include outlier data points in differential expression analysis, which do not follow a negative binomial distribution. Applying ppcseq to analyse several publicly available datasets using popular tools, we show that from 3 to 10 percent of differentially abundant transcripts across algorithms and datasets had statistics inflated by the presence of outliers.

Accurate and sensitive somatic structural variant calling is essential to the comprehensive characterisation of genomic rearrangements in cancer. Current methods restrict their reporting of structural variants to breakpoint calls only, discarding critical information about genomic rearrangements involving low mappability sequences. Here we present GRIDSS2, the first structural variant caller to explicitly report single breakends - breakpoints in which only one side can be unambiguously determined. We show that treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints enables GRIDSS2 to identify rearrangements in regions currently considered inaccessible to short read sequencing. GRIDSS2 outperforms recent state-of-the-art tools on cell line, validated patient samples, and large cohort data sets. Using a cohort of 3,782 deeply sequenced metastatic cancers, we show that GRIDSS2 achieves a false negative rate 3.1%, with 7.0% of calls being single breakend variants to low-mappability sequence. GRIDSS2 retains high sensitivity and precision even for small events. Our identification of a novel small (32-100bp) duplication overlaps the 50bp threshold typically used for SV calling. Demonstrating the power single breakend calling has in low mappability regions, we find that 47% of somatic centromeric copy number changes can be explained by breaks to non-centromeric sequence, with chromosome 1 exhibiting a unique centromeric rearrangement signature. Finally, we show that somatic structural variants are highly clustered, with GRIDSS2 able to phase 16% of somatic structural variants in our metastasis cohort from short read sequencing alone.

Motivation Integration of viruses into infected host cell DNA can cause DNA damage and disrupt genes. Recent cost reductions and growth of whole genome sequencing has produced a wealth of data in which viral presence and integration detection is possible. While key research and clinically relevant insights can be uncovered, existing software has not achieved widespread adoption, limited in part due to high computational costs, the inability to detect a wide range of viruses, as well as precision and sensitivity. Results Here, we describe VIRUSBreakend, a high-speed tool that identifies viral DNA presence and genomic integration. It utilizes single breakends, breakpoints in which only one side can be unambiguously placed, in a novel virus-centric variant calling and assembly approach to identify viral integrations with high sensitivity and a near-zero false discovery rate. VIRUSBreakend detects viral integrations anywhere in the host genome including regions such as centromeres and telomeres unable to be called by existing tools. Applying VIRUSBreakend to a large metastatic cancer cohort, we demonstrate that it can reliably detect clinically relevant viral presence and integration including HPV, HBV, MCPyV, EBV and HHV-8. Availability and implementation VIRUSBreakend is part of the Genomic Rearrangement IDentification Software Suite (GRIDSS). It is available under a GPLv3 license from https://github.com/PapenfussLab/VIRUSBreakend. Supplementary information Supplementary data are available at Bioinformatics online.

Complex somatic genomic rearrangement and copy number alterations (CNA) are hallmarks of nearly all cancers. Whilst whole genome sequencing (WGS) in principle allows comprehensive profiling of these events, biological and clinical interpretation remains challenging. We have developed LINX, a novel algorithm which allows interpretation of short-read paired-end WGS derived structural variant and CNA data by clustering raw structural variant calls into distinct events, predicting their impact on the local structure of the derivative chromosome, and annotating their functional impact on affected genes. Novel visualisations facilitate further investigation of complex genomic rearrangements. We show that LINX provides insights into a diverse range of structural variation events including single and double break-junction events, mobile element insertions, complex shattering and high amplification events. We demonstrate that LINX can reliably detect a wide range of pathogenic rearrangements including gene fusions, immunoglobulin enhancer rearrangements, intragenic deletions and duplications. Uniquely, LINX also predicts chained fusions which we demonstrate account for 13% of clinically relevant oncogenic fusions. LINX also reports a class of inactivation events we term homozygous disruptions which may be a driver mutation in up to 8.8% of tumors including frequently affectingPTEN,TP53andRB1, and are likely missed by many standard WGS analysis pipelines.

BackgroundVaginal melanoma (VM) is a rare cancer and has a poor response to immune checkpoint blockade (ICB). CD8+Tissue Resident Memory (TRM) T cells proliferate in response to ICB and correlate with longer survival in metastatic cutaneous melanoma. However, their capacity to respond to VM and their neoantigens is not known.MethodsUsing longitudinal samples, we explored the evolution of VM mutations by whole-exome sequencing and RNAseq, we also defined the immune context using multiplex immunohistochemistry and nanostring pan cancer immune profile. Then using fresh single cell suspensions of the metastatic samples, we explored VM T cells via mass cytometry and single cell RNAseq and T cell receptor sequencing (TCRseq). Finally, we investigated TRM, pre-TRM and exhausted T cell function against melanoma neo-antigens and melanoma differentiation antigens in vitro.ResultsPrimary VM was non-inflamed and devoid of CD8+ TRM cells. In contrast, both metastases showed proliferating CD8+ TRM were clustered at the tumor margin, with increased numbers in the second ICB-refractory metastasis. The first metastasis showed dense infiltration of CD8+ T cells, the second showed immune exclusion with loss of melanoma cell Major histocompatibility complex (MHC)-I expression associated with downregulation of antigen presentation pathway gene expression. CD8+ TRM from both metastases responded to autologous melanoma cells more robustly than all other CD8+ T cell subsets. In addition, CD8+ TRM shared TCR clones across metastases, suggesting a response to common antigens, which was supported by recognition of the same neoantigen by expanded tumor infiltrating lymphocytes.ConclusionsIn this study, we identified TRM clusters in VM metastases from a patient, but not primary disease. We showed TRM location at the tumor margin, and their superior functional response to autologous tumor cells, predicted neoantigens and melanoma differentiation antigens. These CD8+ TRM exhibited the highest tumor-responsive potential and shared their TCR with tumor-infiltrating effector memory T cells. This suggests VM metastases from this patient retain strong antitumor T cell functional responses; however, this response is suppressed in vivo. The loss of VG MHC-I expression is a common immune escape mechanism which was not addressed by anti-PD-1 monotherapy; rather an additional targeted approach to upregulate MHC-I expression is required.

Ovarian carcinosarcoma (OCS) is an aggressive and rare tumour type with limited treatment options. OCS is hypothesised to develop via the combination theory from a single progenitor, resulting in carcinomatous and sarcomatous components, or alternatively via the conversion theory, with the sarcomatous component developing from the carcinomatous component through epithelial-to-mesenchymal transition (EMT). We show OCS from 18 women to be monoclonal through analysis of DNA variants from isolated carcinoma and sarcoma components. RNA sequencing indicated the carcinoma components were more mesenchymal when compared with pure ovarian carcinomas, supporting the conversion theory. We used pre-clinical OCS models to test the efficacy of microtubule-targeting drugs, including eribulin, which has been shown to reverse EMT characteristics. We demonstrated that microtubule inhibitors, vinorelbine and eribulin, were more effective than standard-of-care platinum-based chemotherapy. Eribulin reduced mesenchymal characteristics, N-MYC expression and cholesterol biosynthesis. Finally, eribulin induced a strong immune response, supporting immunotherapy combinations in the clinic.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasitePlasmodium falciparum. The complexity of thevarmultigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region ofvargenes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination withinvargenes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples ofPlasmodium Laveraniaspecies. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics ofP.falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

Combined inhibition of BRAF, MEK, and CDK4/6 is currently under evaluation in clinical trials for patients with melanoma harboring a BRAFV600 mutation. While this triple therapy has potent tumor-intrinsic effects, the impact of this combination on antitumor immunity remains unexplored. Here, using a syngeneic BrafV600ECdkn2a−/−Pten−/− melanoma model, we demonstrated that triple therapy promoted durable tumor control through tumor-intrinsic mechanisms and promoted immunogenic cell death and T-cell infiltration. Despite this, tumors treated with triple therapy were unresponsive to immune checkpoint blockade (ICB). Flow cytometric and single-cell RNA sequencing analyses of tumor-infiltrating immune populations revealed that triple therapy markedly depleted proinflammatory macrophages and cross-priming CD103+ dendritic cells, the absence of which correlated with poor overall survival and clinical responses to ICB in patients with melanoma. Indeed, immune populations isolated from tumors of mice treated with triple therapy failed to stimulate T-cell responses ex vivo. While combined BRAF, MEK, and CDK4/6 inhibition demonstrates favorable tumor-intrinsic activity, these data suggest that collateral effects on tumor-infiltrating myeloid populations may impact antitumor immunity. These findings have important implications for the design of combination strategies and clinical trials that incorporate BRAF, MEK, and CDK4/6 inhibition with immunotherapy for the treatment of patients with melanoma.

Prostate cancer is caused by genomic aberrations in normal epithelial cells, however clinical translation of findings from analyses of cancer cells alone has been very limited. A deeper understanding of the tumour microenvironment is needed to identify the key drivers of disease progression and reveal novel therapeutic opportunities. In this study, the experimental enrichment of selected cell-types and the development of a Bayesian inference model for continuous differential transcript abundance permitted us to define the transcriptional landscape of the prostate cancer microenvironment along the disease progression axis. An important role of monocytes and macrophages in prostate cancer progression and disease recurrence was uncovered, supported by both transcriptional landscape findings and by differential tissue composition analyses. These findings were corroborated and validated by spatial analyses at the single-cell level using multiplex immunohistochemistry. This study advances our knowledge concerning the role of monocyte-derived recruitment in primary prostate cancer, and supports their key role in disease progression, patient survival and prostate microenvironment immune modulation.

Relative transcript abundance has proven to be a valuable tool for inferring the phenotype of biological systems from genetic material. Several methods for the analysis of differential transcript abundance have been developed, and some of the most popular are based on negative binomial models. Although most genes are fitted reasonably well by the negative binomial distribution, the presence of outlier observations that do not fit such models can lead to artifactual identification of significant changes in transcription. Identifying those transcripts for the correct interpretation of results is extremely important. A robust and automated tool for detecting sample/transcript pairs that do not fit a negative binomial regression model is currently lacking. Here we propose ppcseq, a robust statistical framework that models hierarchically sample- and gene-wise features such as sequencing depth bias, the association between mean transcript abundance and its over-dispersion, and provides a theoretical transcript abundance distribution, on which the observed transcript abundance can be tested for outliers. We show using a publicly available data set where nearly 10% of differentially abundant transcripts had fold change inflated by the presence of outliers. This method has broad utility in filtering artifactual results of differential transcript abundance analyses based on a negative binomial framework.

The parasitic mite Sarcoptes scabiei is an economically highly significant parasite of the skin of humans and animals worldwide. In humans, this mite causes a neglected tropical disease (NTD), called scabies. This disease results in major morbidity, disability, stigma and poverty globally and is often associated with secondary bacterial infections. Currently, anti-scabies treatments are not sufficiently effective, resistance to them is emerging and no vaccine is available. Here, we report the first high-quality genome and transcriptomic data for S. scabiei. The genome is 56.6 Mb in size, has a a repeat content of 10.6% and codes for 9,174 proteins. We explored key molecules involved in development, reproduction, host-parasite interactions, immunity and disease. The enhanced ‘omic data sets for S. scabiei represent comprehensive and critical resources for genetic, functional genomic, metabolomic, phylogenetic, ecological and/or epidemiological investigations, and will underpin the design and development of new treatments, vaccines and/or diagnostic tests., Author summary Scabies is a highly significant parasitic disease caused by the mite S. scabiei. This NTD has a major adverse impact in disadvantaged communities around the world, particularly when associated with secondary bacterial infections and clinical complications. Here we report the first high-quality genome and transcriptomic data for S. scabiei and explore molecular aspects of S. scabiei/scabies. This genome (56.6 Mb, encoding ~ 9,200 proteins) provides a solid foundation for fundamental investigations of the molecular biology of the mite, host-parasite interactions and disease processes as well as for translational research to develop new treatments, vaccines and diagnostic tests.

To separate causal effects of histone acetylation on chromatin accessibility and transcriptional output, we used integrated epigenomic and transcriptomic analyses following acute inhibition of major cellular lysine acetyltransferases P300 and CBP in hematological malignancies. We found that catalytic P300/CBP inhibition dynamically perturbs steady-state acetylation kinetics and suppresses oncogenic transcriptional networks in the absence of changes to chromatin accessibility. CRISPR-Cas9 screening identified NCOR1 and HDAC3 transcriptional co-repressors as the principal antagonists of P300/CBP by counteracting acetylation turnover kinetics. Finally, deacetylation of H3K27 provides nucleation sites for reciprocal methylation switching, a feature that can be exploited therapeutically by concomitant KDM6A and P300/CBP inhibition. Overall, this study indicates that the steady-state histone acetylation-methylation equilibrium functions as a molecular rheostat governing cellular transcription that is amenable to therapeutic exploitation as an anti-cancer regimen.

Background The malaria parasite Plasmodium falciparum has an unusually euchromatic genome with poorly conserved positioning of nucleosomes in intergenic sequences and poorly understood mechanisms of gene regulation. Variant histones and histone modifications determine nucleosome stability and recruit trans factors, but their combinatorial contribution to gene regulation is unclear. Results Here, we show that the histone H3 acetylations H3K18ac and H3K27ac and the variant histone Pf H2A.Z are enriched together at regulatory sites upstream of genes. H3K18ac and H3K27ac together dynamically mark regulatory regions of genes expressed during the asexual life cycle. In contrast, H3K4me1 is depleted in intergenic sequence and dynamically depleted upstream of expressed genes. The temporal pattern of H3K27ac and H3K18ac enrichment indicates that they accumulate during S phase and mitosis and are retained at regulatory sequences until at least G1 phase and after cessation of expression of the cognate genes. We integrated our ChIPseq data with existing datasets to show that in schizont stages H3K18ac, H3K27ac and Pf H2A.Z colocalise with the transcription factor PfAP2-I and the bromodomain protein PfBDP1 and are enriched at stably positioned nucleosomes within regions of exposed DNA at active transcriptional start sites. Using transient transfections we showed that sequences enriched with colocalised H3K18ac, H3K27ac and Pf H2A.Z possess promoter activity in schizont stages, but no enhancer-like activity. Conclusions The dynamic H3 acetylations define P. falciparum regulatory sequences and contribute to gene activation. These findings expand the knowledge of the chromatin landscape that regulates gene expression in P. falciparum.

Our understanding of genomic heterogeneity in lung cancer is largely based on the analysis of early-stage surgical specimens. Here we used endoscopic sampling of paired primary and intrathoracic metastatic tumors from 11 lung cancer patients to map genomic heterogeneity inoperable lung cancer with deep whole-genome sequencing. Intra-patient heterogeneity in driver or targetable mutations was predominantly in the form of copy number gain. Private mutation signatures, including patterns consistent with defects in homologous recombination, were highly variable both within and between patients. Irrespective of histotype, we observed a smaller than expected number of private mutations, suggesting that ancestral clones accumulated large mutation burdens immediately prior to metastasis. Single-region whole-genome sequencing of from 20 patients showed that tumors in ever-smokers with the strongest tobacco signatures were associated with germline variants in genes implicated in the repair of cigarette-induced DNA damage. Our results suggest that lung cancer precursors in ever-smokers accumulate large numbers of mutations prior to the formation of frank malignancy followed by rapid metastatic spread. In advanced lung cancer, germline variants in DNA repair genes may interact with the airway environment to influence the pattern of founder mutations, whereas similar interactions with the tumor microenvironment may play a role in the acquisition of mutations following metastasis.

Anti-microbial signaling pathways are normally triggered by innate immune receptors when detecting pathogenic microbes to provide protective immunity. Here we show that the inflammasome sensor Nlrp1 aggravates DSS-induced experimental mouse colitis by limiting beneficial, butyrate-producing Clostridiales in the gut. The colitis-protective effects of Nlrp1 deficiency are thus reversed by vancomycin treatment, but recapitulated with butyrate supplementation in wild-type mice. Moreover, an activating mutation in Nlrp1a increases IL-18 and IFNγ production, and decreases colonic butyrate to exacerbate colitis. We also show that, in patients with ulcerative colitis, increased NLRP1 in inflamed regions of the colon is associated with increased IFN-γ. In this context, NLRP1, IL-18 or IFN-γ expression negatively correlates with the abundance of Clostridiales in human rectal mucosal biopsies. Our data identify the NLRP1 inflammasome to be a key negative regulator of protective, butyrate-producing commensals, which therefore promotes inflammatory bowel disease.

It has long been assumed that p53 suppresses tumor development through induction of apoptosis, possibly with contributions by cell cycle arrest and cell senescence1,2. However, combined deficiency in these three processes does not result in spontaneous tumor formation as observed upon loss of p53, suggesting the existence of additional mechanisms that are critical mediators of p53-dependent tumor suppression function3-5. To define such mechanisms, we performed in vivo shRNA screens targeting p53-regulated genes in sensitized genetic backgrounds. We found that knockdown of Zmat3, Ctsf and Cav1, promoted lymphoma/leukemia development only when PUMA and p21, the critical effectors of p53-driven apoptosis, cell cycle arrest and senescence, were also absent. Notably, loss of the DNA repair gene Mlh1 caused lymphoma in a wild-type background, and its enforced expression was able to delay tumor development driven by loss of p53. Further examination of direct p53 target genes implicated in DNA repair showed that knockdown of Mlh1, Msh2, Rnf144b, Cav1 and Ddit4 accelerated MYC-driven lymphoma development to a similar extent as knockdown of p53. Collectively, these findings demonstrate that extensive functional overlap of several p53-regulated processes safeguards against cancer and that coordination of DNA repair appears to be an important process by which p53 suppresses tumor development.

As a topical cancer immunotherapy, the toll-like receptor 7 ligand imiquimod activates tumor regression via stimulation of immune cell infiltration and cytotoxic responses. Imiquimod also exerts direct pro-apoptotic effects on tumor cells in vitro, but a role for these effects in imiquimod-induced tumor regression remains undefined. We previously demonstrated that cell lines derived from devil facial tumor disease (DFTD), a transmissible cancer threatening the survival of the Tasmanian devil (Sarcophilus harrisii), are sensitive to imiquimod-induced apoptosis. In this study, the pro-apoptotic effects of imiquimod in DFTD have been investigated using RNA-sequencing and label-free quantitative proteomics. This analysis revealed that changes to gene and protein expression in imiquimod treated DFTD cells are consistent with the onset of oxidative and endoplasmic reticulum stress responses, and subsequent activation of the unfolded protein response, autophagy, cell cycle arrest and apoptosis. Imiquimod also regulates the expression of oncogenic pathways, providing a direct mechanism by which this drug may increase tumor susceptibility to immune cytotoxicity in vivo. Our study has provided the first global analysis of imiquimod-induced effects in any tumor cell line. These findings have highlighted the potential of cell stress pathways as therapeutic targets in DFTD, and will allow for improved mechanistic use of imiquimod as a therapy in both the Tasmanian devil and human cancers.

Correction After publication of our article [1] it was brought to our attention that a line of code was missing from our program to combine the within-replicate variance and between-replicate variance. This led to an overestimation of the standard errors calculated using the Enrich2 random-effects model.

Pheochromocytomas (PC) and paragangliomas (PGL) are endocrine tumors for which the genetic and clinicopathological features of metastatic progression remain incompletely understood. As a result, the risk of metastasis from a primary tumor cannot be predicted. Early diagnosis of individuals at high risk of developing metastases is clinically important and the identification of new biomarkers that are predictive of metastatic potential is of high value. Activation of TERT has been associated with a number of malignant tumors, including PC/PGL. However, the mechanism of TERT activation in the majority of PC/PGL remains unclear. As TERT promoter mutations occur rarely in PC/PGL, we hypothesized that other mechanisms – such as structural variations – may underlie TERT activation in these tumors. From 35 PC and four PGL, we identified three primary PCs that developed metastases with elevated TERT expression, each of which lacked TERT promoter mutations and promoter DNA methylation. Using whole genome sequencing, we identified somatic structural alterations proximal to the TERT locus in two of these tumors. In both tumors, the genomic rearrangements led to the positioning of super-enhancers proximal to the TERT promoter, that are likely responsible for the activation of the normally tightly repressed TERT expression in chromaffin cells

Neoplastic growth and many of the hallmark properties of cancer are driven by the disruption of molecular networks established during the emergence of multicellularity. Regulatory pathways and molecules that evolved to impose regulatory constraints upon networks established in earlier unicellular organisms enabled greater communication and coordination between the diverse cell types required for multicellularity, but also created liabilities in the form of points of vulnerability in the network that when mutated or dysregulated facilitate the development of cancer. These factors are usually overlooked in genomic analyses of cancer, but understanding where vulnerabilities to cancer lie in the networks of multicellular species would provide important new insights into how core molecular processes and gene regulation change during tumourigenesis. We describe how the evolutionary origins of genes influence their roles in cancer, and how connections formed between unicellular and multicellular genes that act as key regulatory hubs for normal tissue homeostasis can also contribute to malignant transformation when disrupted. Tumours in general are characterised by increased dependence on unicellular processes for survival, and major dysregulation of the control structures imposed on these processes during the evolution of multicellularity. Mounting molecular evidence suggests altered interactions at the interface between unicellular and multicellular genes play key roles in the initiation and progression of cancer. Furthermore, unicellular network regions activated in cancer show high degrees of robustness and plasticity, conferring increased adaptability to tumour cells by supporting effective responses to environmental pressures such as drug exposure. Examining how the links between multicellular and unicellular regions get disrupted in tumours has great potential to identify novel drivers of cancer, and to guide improvements to cancer treatment by identifying more effective therapeutic strategies. Recent successes in targeting unicellular processes by novel compounds underscore the logic of such approaches. Further gains could come from identifying genes at the interface between unicellular and multicellular processes and manipulating the communication between network regions of different evolutionary ages.

Background High throughput sequencing requires bioinformatics pipelines to process large volumes of data into meaningful variants that can be translated into a clinical report. These pipelines often suffer from a number of shortcomings: they lack robustness and have many components written in multiple languages, each with a variety of resource requirements. Pipeline components must be linked together with a workflow system to achieve the processing of FASTQ files through to a VCF file of variants. Crafting these pipelines requires considerable bioinformatics and IT skills beyond the reach of many clinical laboratories. Results Here we present Canary, a single program that can be run on a laptop, which takes FASTQ files from amplicon assays through to an annotated VCF file ready for clinical analysis. Canary can be installed and run with a single command using Docker containerization or run as a single JAR file on a wide range of platforms. Although it is a single utility, Canary performs all the functions present in more complex and unwieldy pipelines. All variants identified by Canary are 3′ shifted and represented in their most parsimonious form to provide a consistent nomenclature, irrespective of sequencing variation. Further, proximate in-phase variants are represented as a single HGVS ‘delins’ variant. This allows for correct nomenclature and consequences to be ascribed to complex multi-nucleotide polymorphisms (MNPs), which are otherwise difficult to represent and interpret. Variants can also be annotated with hundreds of attributes sourced from MyVariant.info to give up to date details on pathogenicity, population statistics and in-silico predictors. Conclusions Canary has been used at the Peter MacCallum Cancer Centre in Melbourne for the last 2 years for the processing of clinical sequencing data. By encapsulating clinical features in a single, easily installed executable, Canary makes sequencing more accessible to all pathology laboratories. Canary is available for download as source or a Docker image at https://github.com/PapenfussLab/Canary under a GPL-3.0 License. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1950-z) contains supplementary material, which is available to authorized users.