Disruption of metazoan gene regulatory networks in cancer alters the balance of co-expression between genes of unicellular and multicellular origins

BACKGROUND

Metazoans inherited genes from unicellular ancestors that perform essential biological processes such as cell division, metabolism, and protein translation. Multicellularity requires careful control and coordination of these unicellular genes to maintain tissue integrity and homeostasis. Gene regulatory networks (GRNs) that arose during metazoan evolution are frequently altered in cancer, resulting in over-expression of unicellular genes. We propose that an imbalance in co-expression of unicellular (UC) and multicellular (MC) genes is a driving force in cancer.

RESULTS

We combine gene co-expression analysis to infer changes to GRNs in cancer with protein sequence conservation data to distinguish genes with UC and MC origins. Co-expression networks created using RNA sequencing data from 31 tumor types and normal tissue samples are divided into modules enriched for UC genes, MC genes, or mixed UC-MC modules. The greatest differences between tumor and normal tissue co-expression networks occur within mixed UC-MC modules. MC and UC genes not commonly co-expressed in normal tissues form distinct co-expression modules seen only in tumors. The degree of rewiring of genes within mixed UC-MC modules increases with tumor grade and stage. Mixed UC-MC modules are enriched for somatic mutations in cancer genes, particularly amplifications, suggesting an important driver of the rewiring observed in tumors is copy number changes.

CONCLUSIONS

Our study shows the greatest changes to gene co-expression patterns during tumor progression occur between genes of MC and UC origins, implicating the breakdown and rewiring of metazoan gene regulatory networks in cancer development and progression.

A risk-reward examination of sample multiplexing reagents for single cell RNA-Seq

Single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for understanding cellular heterogeneity and function. However the choice of sample multiplexing reagents can impact data quality and experimental outcomes. In this study, we compared various multiplexing reagents, including MULTI-Seq, Hashtag antibody, and CellPlex, across diverse sample types such as human peripheral blood mononuclear cells (PBMCs), mouse embryonic brain and patient-derived xenografts (PDXs). We found that all multiplexing reagents worked well in cell types robust to ex vivo manipulation but suffered from signal-to-noise issues in more delicate sample types. We compared multiple demultiplexing algorithms which differed in performance depending on data quality. We find that minor improvements to laboratory workflows such as titration and rapid processing are critical to optimal performance. We also compared the performance of fixed scRNA-Seq kits and highlight the advantages of the Parse Biosciences kit for fragile samples. Highly multiplexed scRNA-Seq experiments require more sequencing resources, therefore we evaluated CRISPR-based destruction of non-informative genes to enhance sequencing value. Our comprehensive analysis provides insights into the selection of appropriate sample multiplexing reagents and protocols for scRNA-Seq experiments, facilitating more accurate and cost-effective studies.

cellsig plug-in enhances CIBERSORTx signature selection for multidataset transcriptomes with sparse multilevel modelling

MOTIVATION

The precise characterization of cell-type transcriptomes is pivotal to understanding cellular lineages, deconvolution of bulk transcriptomes, and clinical applications. Single-cell RNA sequencing resources like the Human Cell Atlas have revolutionised cell-type profiling. However, challenges persist due to data heterogeneity and discrepancies across different studies. One limitation of prevailing tools such as CIBERSORTx is their inability to address hierarchical data structures and handle nonoverlapping gene sets across samples, relying on filtering or imputation.

RESULTS

Here, we present cellsig, a Bayesian sparse multilevel model designed to improve signature estimation by adjusting data for multilevel effects and modelling for gene-set sparsity. Our model is tailored to large-scale, heterogeneous pseudobulk and bulk RNA sequencing data collections with nonoverlapping gene sets. We tested the performances of cellsig on a novel curated Human Bulk Cell-type Catalogue, which harmonizes 1435 samples across 58 datasets. We show that cellsig significantly enhances cell-type marker gene ranking performance. This approach is valuable for cell-type signature selection, with implications for marker gene validation, single-cell annotation, and deconvolution benchmarks.

AVAILABILITY AND IMPLEMENTATION

Codes and the interactive app are available at https://github.com/stemangiola/cellsig; and the database is available at https://doi.org/10.5281/zenodo.7582421.

The microtubule inhibitor eribulin demonstrates efficacy in platinum-resistant and refractory high-grade serous ovarian cancer patient-derived xenograft models

Background

Despite initial response to platinum-based chemotherapy and PARP inhibitor therapy (PARPi), nearly all recurrent high-grade serous ovarian cancer (HGSC) will acquire lethal drug resistance; indeed, ~15% of individuals have de novo platinum-refractory disease.

Objectives

To determine the potential of anti-microtubule agent (AMA) therapy (paclitaxel, vinorelbine and eribulin) in platinum-resistant or refractory (PRR) HGSC by assessing response in patient-derived xenograft (PDX) models of HGSC.

Design and methods

Of 13 PRR HGSC PDX, six were primary PRR, derived from chemotherapy-naïve samples (one was BRCA2 mutant) and seven were from samples obtained following chemotherapy treatment in the clinic (five were mutant for either BRCA1 or BRCA2 (BRCA1/2), four with prior PARPi exposure), recapitulating the population of individuals with aggressive treatment-resistant HGSC in the clinic. Molecular analyses and in vivo treatment studies were undertaken.

Results

Seven out of thirteen PRR PDX (54%) were sensitive to treatment with the AMA, eribulin (time to progressive disease (PD) ⩾100 days from the start of treatment) and 11 out of 13 PDX (85%) derived significant benefit from eribulin [time to harvest (TTH) for each PDX with p < 0.002]. In 5 out of 10 platinum-refractory HGSC PDX (50%) and one out of three platinum-resistant PDX (33%), eribulin was more efficacious than was cisplatin, with longer time to PD and significantly extended TTH (each PDX p < 0.02). Furthermore, four of these models were extremely sensitive to all three AMA tested, maintaining response until the end of the experiment (120d post-treatment start). Despite harbouring secondary BRCA2 mutations, two BRCA2-mutant PDX models derived from heavily pre-treated individuals were sensitive to AMA. PRR HGSC PDX models showing greater sensitivity to AMA had high proliferative indices and oncogene expression. Two PDX models, both with prior chemotherapy and/or PARPi exposure, were refractory to all AMA, one of which harboured the SLC25A40-ABCB1 fusion, known to upregulate drug efflux via MDR1.

Conclusion

The efficacy observed for eribulin in PRR HGSC PDX was similar to that observed for paclitaxel, which transformed ovarian cancer clinical practice. Eribulin is therefore worthy of further consideration in clinical trials, particularly in ovarian carcinoma with early failure of carboplatin/paclitaxel chemotherapy.

A Matched Molecular and Clinical Analysis of the Epithelioid Haemangioendothelioma Cohort in the Stafford Fox Rare Cancer Program and Contextual Literature Review

BACKGROUND

Epithelioid haemangioendothelioma (EHE) is an ultra-rare malignant vascular tumour with a prevalence of 1 per 1,000,000. It is typically molecularly characterised by a WWTR1::CAMTA1 gene fusion in approximately 90% of cases, or a YAP1::TFE3 gene fusion in approximately 10% of cases. EHE cases are typically refractory to therapies, and no anticancer agents are reimbursed for EHE in Australia.

METHODS

We report a cohort of nine EHE cases with comprehensive histologic and molecular profiling from the Walter and Eliza Hall Institute of Medical Research Stafford Fox Rare Cancer Program (WEHI-SFRCP) collated via nation-wide referral to the Australian Rare Cancer (ARC) Portal. The diagnoses of EHE were confirmed by histopathological and immunohistochemical (IHC) examination. Molecular profiling was performed using the TruSight Oncology 500 assay, the TruSight RNA fusion panel, whole genome sequencing (WGS), or whole exome sequencing (WES).

RESULTS

Molecular analysis of RNA, DNA or both was possible in seven of nine cases. The WWTR1::CAMTA1 fusion was identified in five cases. The YAP1::TFE3 fusion was identified in one case, demonstrating unique morphology compared to cases with the more common WWTR1::CAMTA1 fusion. All tumours expressed typical endothelial markers CD31, ERG, and CD34 and were negative for pan-cytokeratin. Cases with a WWTR1::CAMTA1 fusion displayed high expression of CAMTA1 and the single case with a YAP1::TFE3 fusion displayed high expression of TFE3. Survival was highly variable and unrelated to molecular profile.

CONCLUSIONS

This cohort of EHE cases provides molecular and histopathological characterisation and matching clinical information that emphasises the molecular patterns and variable clinical outcomes and adds to our knowledge of this ultra-rare cancer. Such information from multiple studies will advance our understanding, potentially improving treatment options.

Single-cell RNA sequencing captures patient-level heterogeneity and associated molecular phenotypes in breast cancer pleural effusions

BACKGROUND

Malignant pleural effusions (MPEs) are a common complication of advanced cancers, particularly those adjacent to the pleura, such as lung and breast cancer. The pathophysiology of MPE formation remains poorly understood, and although MPEs are routinely used for the diagnosis of breast cancer patients, their composition and biology are poorly understood. It is difficult to distinguish invading malignant cells from resident mesothelial cells and to identify the directionality of interactions between these populations in the pleura. There is a need to characterize the phenotypic diversity of breast cancer cell populations in the pleural microenvironment, and investigate how this varies across patients.

METHODS

Here, we used single-cell RNA-sequencing to study the heterogeneity of 10 MPEs from seven metastatic breast cancer patients, including three Miltenyi-enriched samples using a negative selection approach. This dataset of almost 65 000 cells was analysed using integrative approaches to compare heterogeneous cell populations and phenotypes.

RESULTS

We identified substantial inter-patient heterogeneity in the composition of cell types (including malignant, mesothelial and immune cell populations), in expression of subtype-specific gene signatures and in copy number aberration patterns, that captured variability across breast cancer cell populations. Within individual MPEs, we distinguished mesothelial cell populations from malignant cells using key markers, the presence of breast cancer subtype expression patterns and copy number aberration patterns. We also identified pleural mesothelial cells expressing a cancer-associated fibroblast-like transcriptomic program that may support cancer growth.

CONCLUSIONS

Our dataset presents the first unbiased assessment of breast cancer-associated MPEs at a single cell resolution, providing the community with a valuable resource for the study of MPEs. Our work highlights the molecular and cellular diversity captured in MPEs and motivates the potential use of these clinically relevant biopsies in the development of targeted therapeutics for patients with advanced breast cancer.

sccomp: Robust differential composition and variability analysis for single-cell data

Cellular omics such as single-cell genomics, proteomics, and microbiomics allow the characterization of tissue and microbial community composition, which can be compared between conditions to identify biological drivers. This strategy has been critical to revealing markers of disease progression, such as cancer and pathogen infection. A dedicated statistical method for differential variability analysis is lacking for cellular omics data, and existing methods for differential composition analysis do not model some compositional data properties, suggesting there is room to improve model performance. Here, we introduce sccomp, a method for differential composition and variability analyses that jointly models data count distribution, compositionality, group-specific variability, and proportion mean-variability association, being aware of outliers. sccomp provides a comprehensive analysis framework that offers realistic data simulation and cross-study knowledge transfer. Here, we demonstrate that mean-variability association is ubiquitous across technologies, highlighting the inadequacy of the very popular Dirichlet-multinomial distribution. We show that sccomp accurately fits experimental data, significantly improving performance over state-of-the-art algorithms. Using sccomp, we identified differential constraints and composition in the microenvironment of primary breast cancer.

Taurine deficiency as a driver of aging

Aging is associated with changes in circulating levels of various molecules, some of which remain undefined. We find that concentrations of circulating taurine decline with aging in mice, monkeys, and humans. A reversal of this decline through taurine supplementation increased the health span (the period of healthy living) and life span in mice and health span in monkeys. Mechanistically, taurine reduced cellular senescence, protected against telomerase deficiency, suppressed mitochondrial dysfunction, decreased DNA damage, and attenuated inflammaging. In humans, lower taurine concentrations correlated with several age-related diseases and taurine concentrations increased after acute endurance exercise. Thus, taurine deficiency may be a driver of aging because its reversal increases health span in worms, rodents, and primates and life span in worms and rodents. Clinical trials in humans seem warranted to test whether taurine deficiency might drive aging in humans.

A systematic review of computational methods for designing efficient guides for CRISPR DNA base editor systems

In only a few years, as a breakthrough technology, clustered regularly interspaced short palindromic repeats/CRISPR-associated protein (CRISPR/Cas) gene-editing systems have ushered in the era of genome engineering with a plethora of applications. One of the most promising CRISPR tools, so-called base editors, opened an exciting avenue for exploring new therapeutic approaches through controlled mutagenesis. However, the efficiency of a base editor guide varies depending on several biological determinants, such as chromatin accessibility, DNA repair proteins, transcriptional activity, factors related to local sequence context and so on. Thus, the success of genetic perturbation directed by CRISPR/Cas base-editing systems relies on an optimal single guide RNA (sgRNA) design, taking those determinants into account. Although there is 11 commonly used software to design guides specifically for base editors, only three of them investigated and implemented those biological determinants into their models. This review presents the key features, capabilities and limitations of all currently available software with a particular focus on predictive model-based algorithms. Here, we summarize existing software for sgRNA design and provide a base for improving the efficiency of existing available software suites for precise target base editing.

Synthetic Epigenetic Reprogramming of Mesenchymal to Epithelial States Using the CRISPR/dCas9 Platform in Triple Negative Breast Cancer

Epithelial-mesenchymal transition (EMT) is a reversible transcriptional program invoked by cancer cells to drive cancer progression. Transcription factor ZEB1 is a master regulator of EMT, driving disease recurrence in poor-outcome triple negative breast cancers (TNBCs). Here, this work silences ZEB1 in TNBC models by CRISPR/dCas9-mediated epigenetic editing, resulting in highly-specific and nearly complete suppression of ZEB1 in vivo, accompanied by long-lasting tumor inhibition. Integrated "omic" changes promoted by dCas9 linked to the KRAB domain (dCas9-KRAB) enabled the discovery of a ZEB1-dependent-signature of 26 genes differentially-expressed and -methylated, including the reactivation and enhanced chromatin accessibility in cell adhesion loci, outlining epigenetic reprogramming toward a more epithelial state. In the ZEB1 locus transcriptional silencing is associated with induction of locally-spread heterochromatin, significant changes in DNA methylation at specific CpGs, gain of H3K9me3, and a near complete erasure of H3K4me3 in the ZEB1 promoter. Epigenetic shifts induced by ZEB1-silencing are enriched in a subset of human breast tumors, illuminating a clinically-relevant hybrid-like state. Thus, the synthetic epi-silencing of ZEB1 induces stable "lock-in" epigenetic reprogramming of mesenchymal tumors associated with a distinct and stable epigenetic landscape. This work outlines epigenome-engineering approaches for reversing EMT and customizable precision molecular oncology approaches for targeting poor outcome breast cancers.

Targeting homologous recombination deficiency in uterine leiomyosarcoma

BACKGROUND

Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting.

METHODS

A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting.

RESULTS

All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs.

CONCLUSIONS

Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

ahctf1 and kras mutations combine to amplify oncogenic stress and restrict liver overgrowth in a zebrafish model of hepatocellular carcinoma

The nucleoporin (NUP) ELYS, encoded by AHCTF1, is a large multifunctional protein with essential roles in nuclear pore assembly and mitosis. Using both larval and adult zebrafish models of hepatocellular carcinoma (HCC), in which the expression of an inducible mutant kras transgene (kras^G12V) drives hepatocyte-specific hyperplasia and liver enlargement, we show that reducing ahctf1 gene dosage by 50% markedly decreases liver volume, while non-hyperplastic tissues are unaffected. We demonstrate that in the context of cancer, ahctf1 heterozygosity impairs nuclear pore formation, mitotic spindle assembly, and chromosome segregation, leading to DNA damage and activation of a Tp53-dependent transcriptional programme that induces cell death and cell cycle arrest. Heterozygous expression of both ahctf1 and ranbp2 (encoding a second nucleoporin), or treatment of heterozygous ahctf1 larvae with the nucleocytoplasmic transport inhibitor, Selinexor, completely blocks kras^G12V-driven hepatocyte hyperplasia. Gene expression analysis of patient samples in the liver hepatocellular carcinoma (LIHC) dataset in The Cancer Genome Atlas shows that high expression of one or more of the transcripts encoding the 10 components of the NUP107-160 subcomplex, which includes AHCTF1, is positively correlated with worse overall survival. These results provide a strong and feasible rationale for the development of novel cancer therapeutics that target ELYS function and suggest potential avenues for effective combinatorial treatments.

Removing unwanted variation from large-scale RNA sequencing data with PRPS

Accurate identification and effective removal of unwanted variation is essential to derive meaningful biological results from RNA sequencing (RNA-seq) data, especially when the data come from large and complex studies. Using RNA-seq data from The Cancer Genome Atlas (TCGA), we examined several sources of unwanted variation and demonstrate here how these can significantly compromise various downstream analyses, including cancer subtype identification, association between gene expression and survival outcomes and gene co-expression analysis. We propose a strategy, called pseudo-replicates of pseudo-samples (PRPS), for deploying our recently developed normalization method, called removing unwanted variation III (RUV-III), to remove the variation caused by library size, tumor purity and batch effects in TCGA RNA-seq data. We illustrate the value of our approach by comparing it to the standard TCGA normalizations on several TCGA RNA-seq datasets. RUV-III with PRPS can be used to integrate and normalize other large transcriptomic datasets coming from multiple laboratories or platforms.

Integrating deep mutational scanning and low-throughput mutagenesis data to predict the impact of amino acid variants

BACKGROUND

Evaluating the impact of amino acid variants has been a critical challenge for studying protein function and interpreting genomic data. High-throughput experimental methods like deep mutational scanning (DMS) can measure the effect of large numbers of variants in a target protein, but because DMS studies have not been performed on all proteins, researchers also model DMS data computationally to estimate variant impacts by predictors.

RESULTS

In this study, we extended a linear regression-based predictor to explore whether incorporating data from alanine scanning (AS), a widely used low-throughput mutagenesis method, would improve prediction results. To evaluate our model, we collected 146 AS datasets, mapping to 54 DMS datasets across 22 distinct proteins.

CONCLUSIONS

We show that improved model performance depends on the compatibility of the DMS and AS assays, and the scale of improvement is closely related to the correlation between DMS and AS results.

Epithelial-to-Mesenchymal Transition Supports Ovarian Carcinosarcoma Tumorigenesis and Confers Sensitivity to Microtubule Targeting with Eribulin

Ovarian carcinosarcoma (OCS) is an aggressive and rare tumor type with limited treatment options. OCS is hypothesized to develop via the combination theory, with 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). In this study, we analyzed DNA variants from isolated carcinoma and sarcoma components to show that OCS from 18 women is monoclonal. RNA sequencing indicated that the carcinoma components were more mesenchymal when compared with pure epithelial ovarian carcinomas, supporting the conversion theory and suggesting that EMT is important in the formation of these tumors. Preclinical OCS models were used to test the efficacy of microtubule-targeting drugs, including eribulin, which has previously been shown to reverse EMT characteristics in breast cancers and induce differentiation in sarcomas. Vinorelbine and eribulin more effectively inhibited OCS growth than standard-of-care platinum-based chemotherapy, and treatment with eribulin reduced mesenchymal characteristics and N-MYC expression in OCS patient-derived xenografts. Eribulin treatment resulted in an accumulation of intracellular cholesterol in OCS cells, which triggered a downregulation of the mevalonate pathway and prevented further cholesterol biosynthesis. Finally, eribulin increased expression of genes related to immune activation and increased the intratumoral accumulation of CD8+ T cells, supporting exploration of immunotherapy combinations in the clinic. Together, these data indicate that EMT plays a key role in OCS tumorigenesis and support the conversion theory for OCS histogenesis. Targeting EMT using eribulin could help improve OCS patient outcomes.

SIGNIFICANCE

Genomic analyses and preclinical models of ovarian carcinosarcoma support the conversion theory for disease development and indicate that microtubule inhibitors could be used to suppress EMT and stimulate antitumor immunity.

svaRetro and svaNUMT: modular packages for annotating retrotransposed transcripts and nuclear integration of mitochondrial DNA in genome sequencing data

Nuclear integration of mitochondrial genomes and retrocopied transcript insertion are biologically important but often-overlooked aspects of structural variant (SV) annotation. While tools for their detection exist, these typically rely on reanalysis of primary data using specialised detectors rather than leveraging calls from general purpose structural variant callers. Such reanalysis potentially leads to additional computational expense and does not take advantage of advances in general purpose structural variant calling. Here, we present svaRetro and svaNUMT; R packages that provide functions for annotating novel genomic events, such as nonreference retrocopied transcripts and nuclear integration of mitochondrial DNA. The packages were developed to work within the Bioconductor framework. 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. svaRetro and svaNUMT provide modular, SV-caller agnostic tools for downstream annotation of structural variant calls.

A G358S mutation in the Plasmodium falciparum Na+ pump PfATP4 confers clinically-relevant resistance to cipargamin

Diverse compounds target the Plasmodium falciparum Na⁺ pump PfATP4, with cipargamin and (+)-SJ733 the most clinically-advanced. In a recent clinical trial for cipargamin, recrudescent parasites emerged, with most having a G358S mutation in PfATP4. Here, we show that PfATP4^G358S parasites can withstand micromolar concentrations of cipargamin and (+)-SJ733, while remaining susceptible to antimalarials that do not target PfATP4. The G358S mutation in PfATP4, and the equivalent mutation in Toxoplasma gondii ATP4, decrease the sensitivity of ATP4 to inhibition by cipargamin and (+)-SJ733, thereby protecting parasites from disruption of Na⁺ regulation. The G358S mutation reduces the affinity of PfATP4 for Na⁺ and is associated with an increase in the parasite’s resting cytosolic [Na⁺]. However, no defect in parasite growth or transmissibility is observed. Our findings suggest that PfATP4 inhibitors in clinical development should be tested against PfATP4^G358S parasites, and that their combination with unrelated antimalarials may mitigate against resistance development.

In a recent clinical trial for oral administration of cipargamin in individuals with malaria, there was an emergence of recrudescent parasites with a G358S mutation in PfATP4. In this work, the authors investigate the effect of this mutation on the function of the ATPase, on parasite growth and susceptibility to antimalarial drugs.

Genome-wide interrogation of structural variation reveals novel African-specific prostate cancer oncogenic drivers

BACKGROUND

African ancestry is a significant risk factor for advanced prostate cancer (PCa). Mortality rates in sub-Saharan Africa are 2.5-fold greater than global averages. However, the region has largely been excluded from the benefits of whole genome interrogation studies. Additionally, while structural variation (SV) is highly prevalent, PCa genomic studies are still biased towards small variant interrogation.

METHODS

Using whole genome sequencing and best practice workflows, we performed a comprehensive analysis of SVs for 180 (predominantly Gleason score ≥ 8) prostate tumours derived from 115 African, 61 European and four ancestrally admixed patients. We investigated the landscape and relationship of somatic SVs in driving ethnic disparity (African versus European), with a focus on African men from southern Africa.

RESULTS

Duplication events showed the greatest ethnic disparity, with a 1.6- (relative frequency) to 2.5-fold (count) increase in African-derived tumours. Furthermore, we found duplication events to be associated with CDK12 inactivation and MYC copy number gain, and deletion events associated with SPOP mutation. Overall, African-derived tumours were 2-fold more likely to present with a hyper-SV subtype. In addition to hyper-duplication and deletion subtypes, we describe a new hyper-translocation subtype. While we confirm a lower TMPRSS2-ERG fusion-positive rate in tumours from African cases (10% versus 33%), novel African-specific PCa ETS family member and TMPRSS2 fusion partners were identified, including LINC01525, FBXO7, GTF3C2, NTNG1 and YPEL5. Notably, we found 74 somatic SV hotspots impacting 18 new candidate driver genes, with CADM2, LSAMP, PTPRD, PDE4D and PACRG having therapeutic implications for African patients.

CONCLUSIONS

In this first African-inclusive SV study for high-risk PCa, we demonstrate the power of SV interrogation for the identification of novel subtypes, oncogenic drivers and therapeutic targets. Identifying a novel spectrum of SVs in tumours derived from African patients provides a mechanism that may contribute, at least in part, to the observed ethnic disparity in advanced PCa presentation in men of African ancestry.

Tissue-resident memory T cells from a metastatic vaginal melanoma patient are tumor-responsive T cells and increase after anti-PD-1 treatment

Background

Vaginal 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.

Methods

Using 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.

Results

Primary 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.

Conclusions

In 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.

Unscrambling cancer genomes via integrated analysis of structural variation and copy number

Complex somatic genomic rearrangements and copy number alterations are hallmarks of nearly all cancers. We have developed an algorithm, LINX, to aid interpretation of structural variant and copy number data derived from short-read, whole-genome sequencing. LINX classifies raw structural variant calls into distinct events and predicts their effect on the local structure of the derivative chromosome and the functional impact on affected genes. Visualizations facilitate further investigation of complex rearrangements. LINX allows insights into a diverse range of structural variation events and can reliably detect pathogenic rearrangements, including gene fusions, immunoglobulin enhancer rearrangements, intragenic deletions, and duplications. Uniquely, LINX also predicts chained fusions that we demonstrate account for 13% of clinically relevant oncogenic fusions. LINX also reports a class of inactivation events that we term homozygous disruptions that may be a driver mutation in up to 9% of tumors and may frequently affect PTEN, TP53, and RB1.

Women with type 1 diabetes exhibit a progressive increase in gut Saccharomyces cerevisiae in pregnancy associated with evidence of gut inflammation

AIMS

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

OUP accepted manuscript

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.

Safety, infectivity and immunogenicity of a genetically attenuated blood-stage malaria vaccine

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 × 10⁵ (2 subjects), or 3 × 10⁶ 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 × 10⁵ parasites developed parasitemia; 3/4 subjects administered 3× 10⁶ 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).

Type 1 diabetes in pregnancy is associated with distinct changes in the composition and function of the gut microbiome

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. Video abstract.

Transcriptome sequencing and multi-plex imaging of prostate cancer microenvironment reveals a dominant role for monocytic cells in progression

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.

GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing

GRIDSS2 is the first structural variant caller to explicitly report single breakends-breakpoints in which only one side can be unambiguously determined. By treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints, GRIDSS2 can explain 47% of somatic centromere copy number changes using single breakends to non-centromere sequence. On a cohort of 3782 deeply sequenced metastatic cancers, GRIDSS2 achieves an unprecedented 3.1% false negative rate and 3.3% false discovery rate and identifies a novel 32-100 bp duplication signature. GRIDSS2 simplifies complex rearrangement interpretation through phasing of structural variants with 16% of somatic calls phasable using paired-end sequencing.

The site of breast cancer metastases dictates their clonal composition and reversible transcriptomic profile

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.

Detection of cell-free microbial DNA using a contaminant-controlled analysis framework

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.

A single-cell RNA expression atlas of normal, preneoplastic and tumorigenic states in the human breast

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.

Interfacing Seurat with the R tidy universe

Motivation

Seurat is one of the most popular software suites for the analysis of single-cell RNA sequencing data. Considering the popularity of the tidyverse ecosystem, which offers a large set of data display, query, manipulation, integration and visualization utilities, a great opportunity exists to interface the Seurat object with the tidyverse. This interface gives the large data science community of tidyverse users the possibility to operate with familiar grammar.

Results

To provide Seurat with a tidyverse-oriented interface without compromising efficiency, we developed tidyseurat, a lightweight adapter to the tidyverse. Tidyseurat displays cell information as a tibble abstraction, allowing intuitively interfacing Seurat with dplyr, tidyr, ggplot2 and plotly packages powering efficient data manipulation, integration and visualization. Iterative analyses on data subsets are enabled by interfacing with the popular nest-map framework.

Availability and implementation

The software is freely available at cran.r-project.org/web/packages/tidyseurat and github.com/stemangiola/tidyseurat.

Supplementary information

Supplementary data are available at Bioinformatics online.

Targeting histone acetylation dynamics and oncogenic transcription by catalytic P300/CBP inhibition

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.

VIRUSBreakend: Viral Integration Recognition Using Single Breakends

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.

Ductal variant prostate carcinoma is associated with a significantly shorter metastasis-free survival

BACKGROUND

Ductal adenocarcinoma is an uncommon prostate cancer variant. Previous studies suggest that ductal variant histology may be associated with worse clinical outcomes, but these are difficult to interpret. To address this, we performed an international, multi-institutional study to describe the characteristics of ductal adenocarcinoma, particularly focussing on the effect of presence of ductal variant cancer on metastasis-free survival.

METHODS

Patients with ductal variant histology from two institutional databases who underwent radical prostatectomies were identified and compared with an independent acinar adenocarcinoma cohort. After propensity score matching, the effect of the presence of ductal adenocarcinoma on time to biochemical recurrence, initiation of salvage therapy and the development of metastatic disease was determined. Deep whole-exome sequencing was performed for selected cases (n = 8).

RESULTS

A total of 202 ductal adenocarcinoma and 2037 acinar adenocarcinoma cases were analysed. Survival analysis after matching demonstrated that patients with ductal variant histology had shorter salvage-free survival (8.1 versus 22.0 months, p = 0.03) and metastasis-free survival (6.7 versus 78.6 months, p < 0.0001). Ductal variant histology was consistently associated with RB1 loss, as well as copy number gains in TAP1, SLC4A2 and EHHADH.

CONCLUSIONS

The presence of any ductal variant adenocarcinoma at the time of prostatectomy portends a worse clinical outcome than pure acinar cancers, with significantly shorter times to initiation of salvage therapies and the onset of metastatic disease. These features appear to be driven by uncoupling of chromosomal duplication from cell division, resulting in widespread copy number aberration with specific gain of genes implicated in treatment resistance.

T-Cell Receptor Therapy in the Treatment of Ovarian Cancer: A Mini Review

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.

Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling

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.

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

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 parasite Plasmodium falciparum. The complexity of the var multigene 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 of var genes 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 within var genes 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 of Plasmodium Laverania species. 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 of P. 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.

Molecular diagnosis of scabies using a novel probe-based polymerase chain reaction assay targeting high-copy number repetitive sequences in the Sarcoptes scabiei genome

Background

The suboptimal sensitivity and specificity of available diagnostic methods for scabies hampers clinical management, trials of new therapies and epidemiologic studies. Additionally, parasitologic diagnosis by microscopic examination of skin scrapings requires sample collection with a sharp scalpel blade, causing discomfort to patients and difficulty in children. Polymerase chain reaction (PCR)-based diagnostic assays, combined with non-invasive sampling methods, represent an attractive approach. In this study, we aimed to develop a real-time probe-based PCR test for scabies, test a non-invasive sampling method and evaluate its diagnostic performance in two clinical settings.

Methodology/Principal findings

High copy-number repetitive DNA elements were identified in draft Sarcoptes scabiei genome sequences and used as assay targets for diagnostic PCR. Two suitable repetitive DNA sequences, a 375 base pair microsatellite (SSR5) and a 606 base pair long tandem repeat (SSR6), were identified. Diagnostic sensitivity and specificity were tested using relevant positive and negative control materials and compared to a published assay targeting the mitochondrial cox1 gene. Both assays were positive at a 1:100 dilution of DNA from a single mite; no amplification was observed in DNA from samples from 19 patients with other skin conditions nor from house dust, sheep or dog mites, head and body lice or from six common skin bacterial and fungal species. Moderate sensitivity of the assays was achieved in a pilot study, detecting 5/7 (71.4% [95% CI: 29.0% - 96.3%]) of clinically diagnosed untreated scabies patients). Greater sensitivity was observed in samples collected by FLOQ swabs compared to skin scrapings.

Conclusions/Significance

This newly developed qPCR assay, combined with the use of an alternative non-invasive swab sampling technique offers the possibility of enhanced diagnosis of scabies. Further studies will be required to better define the diagnostic performance of these tests.

As scabies control efforts continue to grow, scarcity of diagnostic options hinders success of elimination efforts in endemic areas. Efficiency in large-scale monitoring is further obstructed by invasive sample collection techniques, which are often uncomfortable for patients, and lack sensitivity. We have developed two PCR-based diagnostic assays targeting repetitive DNA elements. These were identified using new data on the S. scabiei genome. Targeting these elements by PCR improved the detection of scabies DNA. Enhanced sensitivity was demonstrated when tested against routine microscopy and a published PCR-based diagnostic assay. When combined with a non-invasive, effective FLOQ swab sampling method, the developed qPCR-based assays may provide a useful complementary tool for diagnosis of scabies, and its application will likely improve scabies control in target populations.

tidybulk: an R tidy framework for modular transcriptomic data analysis

Recently, efforts have been made toward the harmonization of transcriptomic data structures and workflows using the concept of data tidiness, to facilitate modularisation. We present tidybulk, a modular framework for bulk transcriptional analyses that introduces a tidy transcriptomic data structure paradigm and analysis grammar. Tidybulk covers a wide variety of analysis procedures and integrates a large ecosystem of publicly available analysis algorithms under a common framework. Tidybulk decreases coding burden, facilitates reproducibility, increases efficiency for expert users, lowers the learning curve for inexperienced users, and bridges transcriptional data analysis with the tidyverse. Tidybulk is available at R/Bioconductor bioconductor.org/packages/tidybulk.

Combined BRAF, MEK, and CDK4/6 Inhibition Depletes Intratumoral Immune-Potentiating Myeloid Populations in Melanoma

Combined inhibition of BRAF, MEK, and CDK4/6 is currently under evaluation in clinical trials for patients with melanoma harboring a BRAF^V600 mutation. While this triple therapy has potent tumor-intrinsic effects, the impact of this combination on antitumor immunity remains unexplored. Here, using a syngeneic Braf^V600ECdkn2a^-/-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.

Histone modifications associated with gene expression and genome accessibility are dynamically enriched at Plasmodium falciparum regulatory sequences

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.

Author Correction: Determination of RNA structural diversity and its role in HIV-1 RNA splicing

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

Unifying package managers, workflow engines, and containers: Computational reproducibility with BioNix

MOTIVATION

A challenge for computational biologists is to make our analyses reproducible-i.e. to rerun, combine, and share, with the assurance that equivalent runs will generate identical results. Current best practice aims at this using a combination of package managers, workflow engines, and containers.

RESULTS

We present BioNix, a lightweight library built on the Nix deployment system. BioNix manages software dependencies, computational environments, and workflow stages together using a single abstraction: pure functions. This lets users specify workflows in a clean, uniform way, with strong reproducibility guarantees.

AVAILABILITY AND IMPLEMENTATION

BioNix is implemented in the Nix expression language and is released on GitHub under the 3-clause BSD license: https://github.com/PapenfussLab/bionix (biotools:BioNix) (BioNix, RRID:SCR_017662).

High-quality nuclear genome for Sarcoptes scabiei-A critical resource for a neglected parasite

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.

MTOR signaling orchestrates stress-induced mutagenesis, facilitating adaptive evolution in cancer

In microorganisms, evolutionarily conserved mechanisms facilitate adaptation to harsh conditions through stress-induced mutagenesis (SIM). Analogous processes may underpin progression and therapeutic failure in human cancer. We describe SIM in multiple in vitro and in vivo models of human cancers under nongenotoxic drug selection, paradoxically enhancing adaptation at a competing intrinsic fitness cost. A genome-wide approach identified the mechanistic target of rapamycin (MTOR) as a stress-sensing rheostat mediating SIM across multiple cancer types and conditions. These observations are consistent with a two-phase model for drug resistance, in which an initially rapid expansion of genetic diversity is counterbalanced by an intrinsic fitness penalty, subsequently normalizing to complete adaptation under the new conditions. This model suggests synthetic lethal strategies to minimize resistance to anticancer therapy.

Determination of RNA structural diversity and its role in HIV-1 RNA splicing

Human immunodeficiency virus-1 (HIV-1) is a retrovirus with a 10-kb single-stranded RNA genome. HIV-1 must express all of its gene products from the same primary transcript, which undergoes alternative splicing to produce diverse protein products, including structural proteins and regulatory factors^1,2. Despite the critical role of alternative splicing, the mechanisms driving splice-site choice are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements³. We use DMS-MaPseq to probe the structure of HIV-1 RNA in cells and develop an algorithm called Detection of RNA folding Ensembles using Expectation-Maximization (DREEM), which reveals alternative conformations assumed by the same RNA sequence. Contrary to previous models, which analyzed population averages⁴, our results reveal the widespread heterogeneous nature of HIV-1 RNA structure. In addition to confirming that in vitro characterized alternative structures for the HIV-1 Rev Responsive Element (RRE) 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 hypothesis^5–7 that RNA conformation heterogeneity regulates splice site usage and viral gene expression.

Two of a kind: transmissible Schwann cell cancers in the endangered Tasmanian devil (Sarcophilus harrisii)

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.

Tentacle Transcriptomes of the Speckled Anemone (Actiniaria: Actiniidae: Oulactis sp.): Venom-Related Components and Their Domain Structure

Cnidarians are one of the oldest known animal lineages (ca. 700 million years), with a unique envenomation apparatus to deliver a potent mixture of peptides and proteins. Some peptide toxins from cnidarian venom have proven therapeutic potential. Here, we use a transcriptomic/proteomic strategy to identify sequences with similarity to known venom protein families in the tentacles of the endemic Australian 'speckled anemone' (Oulactis sp.). Illumina RNASeq data were assembled de novo. Annotated sequences in the library were verified by cross-referencing individuals' transcriptomes or protein expression evidence from LC-MS/MS data. Sequences include pore-forming toxins, phospholipases, peptidases, neurotoxins (sodium and potassium channel modulators), cysteine-rich secretory proteins and defensins (antimicrobial peptides). Fewer than 4% of the sequences in the library occurred across the three individuals examined, demonstrating high sequence variability of an individual's arsenal. We searched for actinoporins in Oulactis sp. to assess sequence similarity to the only described toxins (OR-A and -G) for this genus and examined the domain architecture of venom-related peptides and proteins. The novel putative actinoporin of Oulactis sp. has a greater similarity to other species in the Actiniidae family than to O. orientalis. Venom-related sequences have an architecture that occurs in single, repeat or multi-domain combinations of venom-related (e.g. ShK-like) and non-venom (e.g. whey acid protein) domains. This study has produced the first transcriptomes for an endemic Australian sea anemone species and the genus Oulactis, while identifying nearly 400 novel venom-related peptides and proteins for future structural and functional analyses and venom evolution studies.

Higher frequency of vertebrate-infecting viruses in the gut of infants born to mothers with type 1 diabetes

BACKGROUND

Microbial exposures in utero and early life shape the infant microbiome, which can profoundly impact on health. Compared to the bacterial microbiome, very little is known about the virome. We set out to characterize longitudinal changes in the gut virome of healthy infants born to mothers with or without type 1 diabetes using comprehensive virome capture sequencing.

METHODS

Healthy infants were selected from Environmental Determinants of Islet Autoimmunity (ENDIA), a prospective cohort of Australian children with a first-degree relative with type 1 diabetes, followed from pregnancy. Fecal specimens were collected three-monthly in the first year of life.

RESULTS

Among 25 infants (44% born to mothers with type 1 diabetes) at least one virus was detected in 65% (65/100) of samples and 96% (24/25) of infants during the first year of life. In total, 26 genera of viruses were identified and >150 viruses were differentially abundant between the gut of infants with a mother with type 1 diabetes vs without. Positivity for any virus was associated with maternal type 1 diabetes and older infant age. Enterovirus was associated with older infant age and maternal smoking.

CONCLUSIONS

We demonstrate a distinct gut virome profile in infants of mothers with type 1 diabetes, which may influence health outcomes later in life. Higher prevalence and greater number of viruses observed compared to previous studies suggests significant underrepresentation in existing virome datasets, arising most likely from less sensitive techniques used in data acquisition.