Search

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

Abstract Background Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of 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.

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

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

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.

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

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

Abstract 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).

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.

Abstract 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

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

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

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

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.

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

Summary: 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.

Abstract Multiplex assays of variant effect (MAVEs), such as deep mutational scans and massively parallel reporter assays, test thousands of sequence variants in a single experiment. Despite the importance of MAVE data for basic and clinical research, there is no standard resource for their discovery and distribution. Here, we present MaveDB (https://www.mavedb.org), a public repository for large-scale measurements of sequence variant impact, designed for interoperability with applications to interpret these datasets. We also describe the first such application, MaveVis, which retrieves, visualizes, and contextualizes variant effect maps. Together, the database and applications will empower the community to mine these powerful datasets.

A number of computational methods have been developed for calling structural variants (SVs) using short read sequencing data. Here, the authors perform a comprehensive benchmarking analysis comparing 10 general-purpose callers and provide recommendations for both users and methods developers.

There is increasing evidence that epigenetic mechanisms contribute to therapeutic resistance in cancer. Here the authors study AML patient samples and a mouse model of non-genetic resistance and find that transcriptional plasticity drives stable epigenetic resistance, and identify regulators of enhancer function as important modulators of resistance.

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.

The inflammasome is normally activated by pathogens to induce tissue inflammation. Here the authors show that, in mouse experimental colitis models, Nlrp1 inflammasome sensor activates IL-18 to reduce beneficial colonic Clostridiales species, thereby decreasing microbial butyrate and its protective effects on colitis.

Abstract Background Scabies is worldwide one of the most common, yet neglected, parasitic skin infections, affecting a wide range of mammals including humans. Limited treatment options and evidence of emerging mite resistance against the currently used drugs drive our research to explore new therapeutic candidates. Previously, we discovered a multicopy family of genes encoding cysteine proteases with their catalytic sites inactivated by mutation (SMIPP-Cs). This protein family is unique in parasitic scabies mites and is absent in related non-burrowing mites. We postulated that the SMIPP-Cs have evolved as an adaptation to the parasitic lifestyle of the scabies mite. To formulate testable hypotheses for their functions and to propose possible strategies for translational research we investigated whether the SMIPP-Cs are common to all scabies mite varieties and where within the mite body as well as when throughout the parasitic life-cycle they are expressed. Results SMIPP-C sequences from human, pig and dog mites were analysed bioinformatically and the phylogenetic relationships between the SMIPP-C multi-copy gene families of human, pig and dog mites were established. Results suggest that amplification of the SMIPP-C genes occurred in a common ancestor and individual genes evolved independently in the different mite varieties. Recombinant human mite SMIPP-C proteins were produced and used for murine polyclonal antibody production. Immunohistology on skin sections from human patients localised the SMIPP-Cs in the mite gut and in mite faeces within in the epidermal skin burrows. SMIPP-C transcription into mRNA in different life stages was assessed in human and pig mites by reverse transcription followed by droplet digital PCR (ddPCR). High transcription levels of SMIPP-C genes were detected in the adult female life stage in comparison to all other life stages. Conclusions The fact that the SMIPP-Cs are unique to three Sarcoptes varieties, present in all burrowing life stages and highly expressed in the digestive system of the infective adult female life stage may highlight an essential role in parasitism. As they are excreted from the gut in scybala they presumably are able to interact or interfere with host proteins present in the epidermis.

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.

STING is an intracellular DNA sensor that can alter response to infection, but in the case of malaria it is unclear how parasite DNA in red blood cells (RBCs) reaches DNA sensors in immune cells. Here the authors show that STING in human monocytes can sense P. falciparum nucleic acids transported from infected RBCs via parasite extracellular vesicles.

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

Abstract Deep mutational scanning is a widely used method for multiplex measurement of functional consequences of protein variants. We developed a new deep mutational scanning statistical model that generates error estimates for each measurement, capturing both sampling error and consistency between replicates. We apply our model to one novel and five published datasets comprising 243,732 variants and demonstrate its superiority in removing noisy variants and conducting hypothesis testing. Simulations show our model applies to scans based on cell growth or binding and handles common experimental errors. We implemented our model in Enrich2, software that can empower researchers analyzing deep mutational scanning data.

Abstract Spiny mice of the genus Acomys display several unique physiological traits, including menstruation and scar-free wound healing; characteristics that are exceedingly rare in mammals, and of considerable interest to the scientific community. These unique attributes, and the potential for spiny mice to accurately model human diseases, are driving increased use of this genus in biomedical research, however little genetic information is accessible for this species. This project aimed to generate a draft transcriptome for the Common spiny mouse (Acomys cahirinus). Illumina sequencing of RNA from 15 organ types (male and female) produced 451 million, 150 bp paired-end reads (92.4Gbp). An extensive survey of de novo transcriptome assembly approaches using Trinity, SOAPdenovo-Trans, and Oases at multiple kmer lengths was conducted, producing 50 single-kmer assemblies from this dataset. Non-redundant transcripts from all assemblies were merged into a meta-assembly using the EvidentialGene tr2aacds pipeline, producing the largest gene catalogue to date for Acomys cahirinus. This study provides the first detailed characterization of the spiny mouse transcriptome. It validates use of the EvidentialGene tr2aacds pipeline in mammals to augment conventional de novo assembly approaches, and provides a valuable scientific resource for further investigation into the unique physiological characteristics inherent in the genus Acomys.

Abstract Background A precise understanding of structural variants (SVs) in DNA is important in the study of cancer and population diversity. Many methods have been designed to identify SVs from DNA sequencing data. However, the problem remains challenging because existing approaches suffer from low sensitivity, precision, and positional accuracy. Furthermore, many existing tools only identify breakpoints, and so not collect related breakpoints and classify them as a particular type of SV. Due to the rapidly increasing usage of high throughput sequencing technologies in this area, there is an urgent need for algorithms that can accurately classify complex genomic rearrangements (involving more than one breakpoint or fusion). Results We present CLOVE, an algorithm for integrating the results of multiple breakpoint or SV callers and classifying the results as a particular SV. CLOVE is based on a graph data structure that is created from the breakpoint information. The algorithm looks for patterns in the graph that are characteristic of more complex rearrangement types. CLOVE is able to integrate the results of multiple callers, producing a consensus call. Conclusions We demonstrate using simulated and real data that re-classified SV calls produced by CLOVE improve on the raw call set of existing SV algorithms, particularly in terms of accuracy. CLOVE is freely available from http://www.github.com/PapenfussLab .

Biomphalaria glabrata is a fresh water snail that acts as a host for trematode Schistosoma mansoni that causes intestinal infection in human. This work describes the genome and transcriptome analyses from 12 different tissues of B glabrata, and identify genes for snail behavior and evolution.

Abstract Background The increasing affordability of DNA sequencing has allowed it to be widely deployed in pathology laboratories. However, this has exposed many issues with the analysis and reporting of variants for clinical diagnostic use. Implementing a high-throughput sequencing (NGS) clinical reporting system requires a diverse combination of capabilities, statistical methods to identify variants, global variant databases, a validated bioinformatics pipeline, an auditable laboratory workflow, reproducible clinical assays and quality control monitoring throughout. These capabilities must be packaged in software that integrates the disparate components into a useable system. Results To meet these needs, we developed a web-based application, PathOS, which takes variant data from a patient sample through to a clinical report. PathOS has been used operationally in the Peter MacCallum Cancer Centre for two years for the analysis, curation and reporting of genetic tests for cancer patients, as well as the curation of large-scale research studies. PathOS has also been deployed in cloud environments allowing multiple institutions to use separate, secure and customisable instances of the system. Increasingly, the bottleneck of variant curation is limiting the adoption of clinical sequencing for molecular diagnostics. PathOS is focused on providing clinical variant curators and pathology laboratories with a decision support system needed for personalised medicine. While the genesis of PathOS has been within cancer molecular diagnostics, the system is applicable to NGS clinical reporting generally. Conclusions The widespread availability of genomic sequencers has highlighted the limited availability of software to support clinical decision-making in molecular pathology. PathOS is a system that has been developed and refined in a hospital laboratory context to meet the needs of clinical diagnostics. The software is available as a set of Docker images and source code at https://github.com/PapenfussLab/PathOS .

Disease monitoring of chronic lymphocytic leukaemia (CLL) is a challenge. Here, the authors show that serial ctDNA analysis in 32 CLL patients allows monitoring of clonal dynamics over time, and identifies the emergence of genomic changes associated with Richter’s syndrome.

The Eμ-Myc lymphoma mouse model has been invaluable in the study of this disease. Here, the authors use multiple sequencing strategies to analyse the tumours in these mice and find recurrent inactivating mutations in Bcor, suggesting that this gene has a negative role in Mycsignalling.

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.

Summary: The ability of cells to stably maintain their fate is governed by specific transcription regulators. Here, we show that the Scalloped (Sd) and Nervous fingers-1 (Nerfin-1) transcription factors physically and functionally interact to maintain medulla neuron fate in the Drosophila melanogaster CNS. Using Targeted DamID, we find that Sd and Nerfin-1 occupy a highly overlapping set of target genes, including regulators of neural stem cell and neuron fate, and signaling pathways that regulate CNS development such as Notch and Hippo. Modulation of either Sd or Nerfin-1 activity causes medulla neurons to dedifferentiate to a stem cell-like state, and this is mediated at least in part by Notch pathway deregulation. Intriguingly, orthologs of Sd and Nerfin-1 have also been implicated in control of neuronal cell fate decisions in both worms and mammals. Our data indicate that this transcription factor pair exhibits remarkable biochemical and functional conservation across metazoans. : Vissers et al. report that the Scalloped and Nerfin-1 transcription factors physically and functionally interact to maintain medulla neuron fate in the Drosophila CNS. Identification of their genome-binding profiles reveal that these transcription factors regulate neuron fate, at least in part, by modulating Notch signaling. Keywords: neuron, cell fate, Drosophila, transcription, brain

Tasmanian devil (Sarcophilus harrisii) pouch young, like other marsupials, are born underdeveloped and immunologically naïve, and are unable to mount an adaptive immune response. The mother’s milk provides nutrients for growth and development as well as providing passive immunity. To better understand immune response in this endangered species, we set out to characterise the genes involved in passive immunity by sequencing and annotating the transcriptome of a devil milk sample collected during mid-lactation. At mid-lactation we expect the young to have heightened immune responses, as they have emerged from the pouch, encountering new pathogens. A total of 233,660 transcripts were identified, including approximately 17,827 unique protein-coding genes and 846 immune genes. The most highly expressed transcripts were dominated by milk protein genes such as those encoding early lactation protein, late lactation proteins, α-lactalbumin, α-casein and β-casein. There were numerous highly expressed immune genes including lysozyme, whey acidic protein, ferritin and major histocompatibility complex I and II. Genes encoding immunoglobulins, antimicrobial peptides, chemokines and immune cell receptors were also identified. The array of immune genes identified in this study reflects the importance of the milk in providing immune protection to Tasmanian devil young and provides the first insight into Tasmanian devil milk.

Nature Communications 8: Article number: 15451 (2017); Published 16 May 2017; Updated 23 August 2017 The original version of this Article contained an error in the spelling of the author Leon Di Stefano, which was incorrectly given as Leon di Stephano. This has now been corrected in both the PDF andHTML versions of the Article.