Search

Supplementary Figure 2, related to main Figure 3

Supplementary Figure 1, related to main Figure 1

Supplementary Figure 3, related to main Figure 5.

Supplementary Figure 4, related to main Figure 5. Data pertaining to independent validation in the BEAT-AML cohort.

Purpose:RNA splicing is a fundamental biological process that generates protein diversity from a finite set of genes. Recurrent somatic mutations of splicing factor genes are common in some hematologic cancers but are relatively uncommon in acute myeloid leukemia (AML, < 20% of patients). We examined whether RNA splicing differences exist in AML, even in the absence of splicing factor mutations.Experimental Design:We developed a bioinformatics pipeline to study alternative RNA splicing in RNA-sequencing data from large cohorts of patients with AML.Results:We have identified recurrent differential alternative splicing between patients with poor and good prognosis. These splicing events occurred even in patients without any discernible splicing factor mutations. Alternative splicing recurrently occurred in genes with specific molecular functions, primarily related to protein translation. Developing tools to predict the functional impact of alternative splicing on the translated protein, we discovered that approximately 45% of the splicing events directly affected highly conserved protein domains. Several splicing factors were themselves misspliced and the splicing of their target transcripts were altered. Studying differential gene expression in the same patients, we identified that alternative splicing of protein translation genes in ELNAdv patients resulted in the induction of an integrated stress response and upregulation of inflammation-related genes. Finally, using machine learning techniques, we identified a splicing signature of four genes which refine the accuracy of existing risk prognosis schemes and validated it in a completely independent cohort.Conclusions:Our discoveries therefore identify aberrant alternative splicing as a molecular feature of adverse AML with clinical relevance.See related commentary by Bowman, p. 3503

Supplementary Methods - revised, plus Supplementary Tables 1-6

Sterile male Queensland fruit fly, Bactrocera tryoni (Froggatt), fed as immature adults on the plant compound raspberry ketone (RK), show a reduced attraction to cuelure, a synthetic analogue of RK used as an attractant in Male Annihilation Technique. We hypothesized the reduced attraction of RK-fed adult males to cuelure may be a consequence of altered expression of chemoreception genes. A Y-tube olfactometer assay with RK-fed and RK-unfed sterile B. tryoni males tested the subsequent behavioural response to cuelure. Behavioral assays confirmed a significant decrease in attraction of RK-fed sterile males to cuelure. RK-fed, non-responders (to cue-lure) and RK-unfed, responders (to cue-lure) males were sampled and gene expression compared by de novo RNA-seq analysis. A total of 269 genes in fly heads were differentially expressed between replicated groups of RK-fed, cuelure non-responders and RK-unfed, cuelure responders. Among them, 218 genes including 4 chemoreceptor genes were up regulated and 51 genes were down regulated in RK-fed, cuelure non-responders. De novo assembly generated many genes with unknown functions and no significant BLAST hits to homologues in other species. The enriched and suppressed genes reported here, shed light on the transcriptional changes that affect the dynamics of insect responses to chemical stimuli.

Background The origin of most of the Lactobacillus rhamnosus genome sequences lodged in NCBI can be traced to food and faecal isolates followed by blood and tissue sites but with minimal representation from oral and vaginal isolates. However, on the L. rhamnosus phylogenetic tree no apparent clade is linked to the origin of isolation or to the relevant clinical source, except for a distinct clade exclusively shared by L. rhamnosus isolates from early stages of dental pulp infection (LRHMDP2 and LRHMDP3) and from bronchoalveolar lavage (699_LRHA and 708_LRHA) from a critical care patient. These L. rhamnosus strains, LRHMDP2, LRHMDP3, 699_LRHA and 708_LRHA isolated from different continents, display closest genome neighbour gapped identity of 99.95%. The aim of this study was to define a potentially unique complement of genes of clinical relevance shared between these L. rhamnosus clinical isolates in comparison to probiotic L. rhamnosus strains. Results In this analysis we used orthologous protein identification tools such as ProteinOrtho followed by tblastn alignments to identify a novel tyrosine protein phosphatase (wzb)-tyrosine-protein kinase modulator EpsC (wzd)- synteny exopolysaccharide (EPS) cluster. This EPS cluster was specifically conserved in a clade of 5 clinical isolates containing the four L. rhamnosus clinical isolates noted above and Lactobacillus spp. HMSC077C11, a clinical isolate from a neck abscess. The EPS cluster was shared with only two other strains, L. rhamnosus BPL5 and BPL15, which formed a distant clade on the L. rhamnosus phylogenetic tree, with a closest genome neighbour gapped identity of 97.51% with L. rhamnosus LRHMDP2 and LRHMDP3. Exclusivity of this EPS cluster (from those identified before) was defined by five EPS genes, which were specifically conserved between the clade of 5 clinical isolates and L. rhamnosus BPL5 and BPL15 when compared to the remaining L. rhamnosus strains. Comparative genome analysis between the clade of 5 clinical isolates and L. rhamnosus BPL5 and BPL15 showed a set of 58 potentially unique genes characteristic of the clade of 5. Conclusion The potentially unique functional protein orthologs associated with the clade of 5 clinical isolates may provide understanding of fitness under selective pressure.

BackgroundThe enrichment of Gram-negative bacteria of oral origin in the esophageal microbiome has been associated with the development of metaplasia. However, to date, no study has comprehensively assessed the relationships between the esophageal microbiome and the host.MethodsHere, we examine the esophageal microenvironment in gastro-esophageal reflux disease and metaplasia using multi-omics strategies targeting the microbiome and host transcriptome, followed by targeted culture, comparative genomics, and host-microbial interaction studies of bacterial signatures of interest.ResultsProfiling of the host transcriptome from esophageal mucosal biopsies revealed profound changes during metaplasia. Importantly, five biomarkers showed consistent longitudinal changes with disease progression from reflux disease to metaplasia. We showed for the first time that the esophageal microbiome is distinct from the salivary microbiome and the enrichment ofCampylobacterspecies as a consistent signature in disease across two independent cohorts. Shape fitting and matrix correlation identified associations between the microbiome and host transcriptome profiles, with a novel co-exclusion relationship found betweenCampylobacterand napsin B aspartic peptidase. Targeted culture ofCampylobacterspecies from the same cohort revealed a subset of isolates to have a higher capacity to survive within primary human macrophages. Comparative genomic analyses showed these isolates could be differentiated by specific genomic features, one of which was validated to be associated with intracellular fitness. Screening for theseCampylobacterstrain-specific signatures in shotgun metagenomics data from another cohort showed an increase in prevalence with disease progression. Comparative transcriptomic analyses of primary esophageal epithelial cells exposed to theCampylobacterisolates revealed expression changes within those infected with strains with high intracellular fitness that could explain the increased likelihood of disease progression.ConclusionsWe provide a comprehensive assessment of the esophageal microenvironment, identifying bacterial strain-specific signatures with high relevance to progression of metaplasia.

Chemo-resistance in acute myeloid leukemia (AML) patients is driven by leukemic stem cells (LSCs) resulting in high rates of relapse and low overall survival. Here, we demonstrate that upregulation of the splicing factor, RBM17 preferentially marks and sustains LSCs and directly correlates with shorten patient survival. RBM17 knockdown in primary AML cells leads to myeloid differentiation and impaired colony formation and in vivo engraftment. Integrative multi-omics analyses show that RBM17 repression leads to inclusion of poison exons and production of nonsense-mediated decay (NMD)-sensitive transcripts for pro-leukemic factors and the translation initiation factor, EIF4A2. We show that EIF4A2 is enriched in LSCs and its inhibition impairs primary AML progenitor activity. Proteomic analysis of EIF4A2-depleted AML cells shows recapitulation of the RBM17 knockdown biological effects, including pronounced suppression of proteins involved in ribosome biogenesis. Overall, these results provide a rationale to target RBM17 and/or its downstream NMD-sensitive splicing substrates for AML treatment.

Ocean acidification (OA) can be detrimental to calcifying marine organisms, with stunting of invertebrate larval development one of the most consistent responses. Effects are usually measured by short‐term, within‐generation exposure, an approach that does not consider the potential for adaptation. We examined the genetic response to OA of larvae of the tropical sea urchin Echinometra sp. C. raised on coral reefs that were either influenced by CO2 vents (pH ~ 7.9, future OA condition) or nonvent control reefs (pH 8.2). We assembled a high quality de novo transcriptome of Echinometra embryos (8 hr) and pluteus larvae (48 hr) and identified 68,056 SNPs. We tested for outlier SNPs and functional enrichment in embryos and larvae raised from adults from the control or vent sites. Generally, highest FST values in embryos were observed between sites (intrinsic adaptation, most representative of the gene pool in the spawned populations). This comparison also had the highest number of outlier loci (40). In the other comparisons, classical adaptation (comparing larvae with adults from the control transplanted to either the control or vent conditions) and reverse adaptation (larvae from the vent site returned to the vent or explanted at the control), we only observed modest numbers of outlier SNPs (6–19) and only enrichment in two functional pathways. Most of the outliers detected were silent substitutions without adaptive potential. We conclude that there is little evidence of realized adaptation potential during early development, while some potential (albeit relatively low) exists in the intrinsic gene pool after more than one generation of exposure.

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

A thermostable drug suppresses a pERK-FosB/ΔFosB-VCAM-1 axis, endothelial activation, angiogenesis, and retinal permeability., Vascular permeability and angiogenesis underpin neovascular age-related macular degeneration and diabetic retinopathy. While anti-VEGF therapies are widely used clinically, many patients do not respond optimally, or at all, and small-molecule therapies are lacking. Here, we identified a dibenzoxazepinone BT2 that inhibits endothelial cell proliferation, migration, wound repair in vitro, network formation, and angiogenesis in mice bearing Matrigel plugs. BT2 interacts with MEK1 and inhibits ERK phosphorylation and the expression of FosB/ΔFosB, VCAM-1, and many genes involved in proliferation, migration, angiogenesis, and inflammation. BT2 reduced retinal vascular leakage following rat choroidal laser trauma and rabbit intravitreal VEGF-A165 administration. BT2 suppressed retinal CD31, pERK, VCAM-1, and VEGF-A165 expression. BT2 reduced retinal leakage in rats at least as effectively as aflibercept, a first-line therapy for nAMD/DR. BT2 withstands boiling or autoclaving and several months’ storage at 22°C. BT2 is a new small-molecule inhibitor of vascular permeability and angiogenesis.

Despite advances in algorithms and computational platforms, de-novo genome assembly remains a challenging process. Due to the constant innovation in sequencing technologies (Sanger, SOLiD, Illumina, 454, PacBio and Oxford Nanopore), genome assembly has evolved to respond to the changes in input data type. This paper includes a broad and comparative review of the most recent short-read, long-read and hybrid assembly techniques. In this review, we provide (1) an algorithmic description of the important processes in the workflow that introduces fundamental concepts and improvements; (2) a review of existing software that explains possible options for genome assembly; and (3) a comparison of the accuracy and the performance of existing methods executed on the same computer using the same processing capabilities and using the same set of real and synthetic datasets. Such evaluation allows a fair and precise comparison of accuracy in all aspects. As a result, this paper identifies both the strengths and weaknesses of each method. This comparative review is unique in providing a detailed comparison of a broad spectrum of cutting-edge algorithms and methods.

RNA splicing is a fundamental biological process that generates protein diversity from a finite set of genes. Recurrent somatic mutations of splicing factor genes are relatively uncommon in Acute Myeloid Leukemia (AML, < 20%). We examined whether RNA splicing differences exist in AML even in the absence of splicing factor mutations. Analyzing RNA-seq data from two independent cohorts of AML patients, we identified recurrent differential alternative splicing between patients with poor and good prognosis. These alternative splicing events occurred even in patients without any discernible splicing factor mutations. The alternative splicing events recurrently occurred in genes involved in specific molecular functions, primarily related to protein translation. Developing informatics tools to predict the functional impact of alternative splicing on the translated protein, we discovered that ~45% of the splicing events directly affected highly conserved protein domains. Several splicing factors were themselves misspliced in patients, and the splicing of their target transcripts were also altered. By studying differential gene expression in the same patients, we identified that alternative splicing of protein translation genes in ELNAdv patients resulted in the induction of an integrated stress response and up- regulation of inflammation-related genes. Lastly, using machine learning techniques, we identified a set of four genes whose alternative splicing can refine the accuracy of existing risk prognosis schemes and validated it in a completely independent cohort. Our discoveries therefore identify aberrant alternative splicing as a molecular feature of adverse AML with clinical relevance.

Purpose: RNA splicing is a fundamental biological process that generates protein diversity from a finite set of genes. Recurrent somatic mutations of splicing factor genes are common in some hematologic cancers but are relatively uncommon in acute myeloid leukemia (AML, < 20% of patients). We examined whether RNA splicing differences exist in AML, even in the absence of splicing factor mutations. Experimental Design: We developed a bioinformatics pipeline to study alternative RNA splicing in RNA-sequencing data from large cohorts of patients with AML. Results: We have identified recurrent differential alternative splicing between patients with poor and good prognosis. These splicing events occurred even in patients without any discernible splicing factor mutations. Alternative splicing recurrently occurred in genes with specific molecular functions, primarily related to protein translation. Developing tools to predict the functional impact of alternative splicing on the translated protein, we discovered that approximately 45% of the splicing events directly affected highly conserved protein domains. Several splicing factors were themselves misspliced and the splicing of their target transcripts were altered. Studying differential gene expression in the same patients, we identified that alternative splicing of protein translation genes in ELNAdv patients resulted in the induction of an integrated stress response and upregulation of inflammation-related genes. Finally, using machine learning techniques, we identified a splicing signature of four genes which refine the accuracy of existing risk prognosis schemes and validated it in a completely independent cohort. Conclusions: Our discoveries therefore identify aberrant alternative splicing as a molecular feature of adverse AML with clinical relevance. See related commentary by Bowman, p. 3503

Macroalgae (seaweeds) are essential for the functioning of temperate marine ecosystems, but there is increasing evidence to suggest that their survival is under threat from anthropogenic stressors and disease. Nautella italica R11 is recognized as an aetiological agent of bleaching disease in the red alga, Delisea pulchra. Yet, there is a lack of knowledge surrounding the molecular mechanisms involved in this model host-pathogen interaction. Here we report that mutations in the gene encoding for a LuxR-type quorum sensing transcriptional regulator, RaiR, render N. italica R11 avirulent, suggesting this gene is important for regulating the expression of virulence phenotypes. Using an RNA sequencing approach, we observed a strong transcriptional response of N. italica R11 towards the presence of D. pulchra. In particular, genes involved in oxidative stress resistance, carbohydrate and central metabolism were upregulated in the presence of the host, suggesting a role for these functions in the opportunistic pathogenicity of N. italica R11. Furthermore, we show that RaiR regulates a subset of genes in N. italica R11, including those involved in metabolism and the expression of phage-related proteins. The outcome of this research reveals new functions important for virulence of N. italica R11 and contributes to our greater understanding of the complex factors mitigating microbial diseases in macroalgae.

The human transcriptome consists of various RNA biotypes including multiple types of non-coding RNAs (ncRNAs). Current ncRNA compendia remain incomplete partially because they are almost exclusively derived from the interrogation of small- and polyadenylated RNAs. Here, we present a more comprehensive atlas of the human transcriptome that is derived from matching polyA-, total-, and small-RNA profiles of a heterogeneous collection of nearly 300 human tissues and cell lines. We report on thousands of novel RNA species across all major RNA biotypes, including a hitherto poorly-cataloged class of non-polyadenylated single-exon long non-coding RNAs. In addition, we exploit intron abundance estimates from total RNA-sequencing to test and verify functional regulation by novel non-coding RNAs. Our study represents a substantial expansion of the current catalogue of human ncRNAs and their regulatory interactions. All data, analyses, and results are available in the R2 web portal and serve as a basis to further explore RNA biology and function.

Bacteria capable of dechlorinating the toxic environmental contaminant dichloromethane (DCM, CH2Cl2) are of great interest for potential bioremediation applications. A novel, strictly anaerobic, DCM-fermenting bacterium, “DCMF”, was enriched from organochlorine-contaminated groundwater near Botany Bay, Australia. The enrichment culture was maintained in minimal, mineral salt medium amended with dichloromethane as the sole energy source. PacBio whole genome SMRTTMsequencing of DCMF allowedde novo, gap-free assembly despite the presence of cohabiting organisms in the culture. Illumina sequencing reads were utilised to correct minor indels. The single, circularised 6.44 Mb chromosome was annotated with the IMG pipeline and contains 5,773 predicted protein-coding genes. Based on 16S rRNA gene and predicted proteome phylogeny, the organism appears to be a novel member of thePeptococcaceaefamily. The DCMF genome is large in comparison to known DCM-fermenting bacteria. It includes an abundance of methyltransferases, which may provide clues to the basis of its DCM metabolism, as well as potential to metabolise additional methylated substrates such as quaternary amines. Full annotation has been provided in a custom genome browser and search tool, in addition to multiple sequence alignments and phylogenetic trees for every predicted protein,http://www.slimsuite.unsw.edu.au/research/dcmf/.

TRAM domain proteins present in Archaea and Bacteria have a β-barrel shape with anti-parallel β-sheets that form a nucleic acid binding surface; a structure also present in cold shock proteins (Csps). Aside from protein structures, experimental data defining the function of TRAM domains is lacking. Here, we explore the possible functional properties of a single TRAM domain protein, Ctr3 (cold-responsive TRAM domain protein 3) from the Antarctic archaeon Methanococcoides burtonii that has increased abundance during low temperature growth. Ribonucleic acid (RNA) bound by Ctr3 in vitro was determined using RNA-seq. Ctr3-bound M. burtonii RNA with a preference for transfer (t)RNA and 5S ribosomal RNA, and a potential binding motif was identified. In tRNA, the motif represented the C loop; a region that is conserved in tRNA from all domains of life and appears to be solvent exposed, potentially providing access for Ctr3 to bind. Ctr3 and Csps are structurally similar and are both inferred to function in low temperature translation. The broad representation of single TRAM domain proteins within Archaea compared with their apparent absence in Bacteria, and scarcity of Csps in Archaea but prevalence in Bacteria, suggests they represent distinct evolutionary lineages of functionally equivalent RNA-binding proteins.

The de-novo genome assembly is a challenging computational problem for which several pipelines have been developed. The advent of long-read sequencing technology has resulted in a new set of algorithmic approaches for the assembly process. In this work, we identify that one of these new and fast long-read assembly techniques (using Minimap2 and Miniasm ) can be modified for the short-read assembly process. This possibility motivated us to customize a long-read assembly approach for applications in a short-read assembly scenario. Here, we compare and contrast our proposed de-novo assembly pipeline ( MiniSR ) with three other recently developed programs for the assembly of bacterial and small eukaryotic genomes. We have documented two trade-offs: one between speed and accuracy and the other between contiguity and base-calling errors. Our proposed assembly pipeline shows a good balance in these trade-offs. The resulting pipeline is 6 and 2.2 times faster than the short-read assemblers Spades and SGA, respectively. MiniSR generates assemblies of superior N50 and NGA50 to SGA , although assemblies are less complete and accurate than those from Spades . A third tool, SOAPdenovo2 , is as fast as our proposed pipeline but had poorer assembly quality.

Background & Aims Fecal microbiota transplantation (FMT) can induce remission in patients with ulcerative colitis (UC). In a randomized controlled trial of FMT in patients with active UC, we aimed to identify bacterial taxonomic and functional factors associated with response to therapy. Methods We performed a double-blind trial of 81 patients with active UC randomly assigned to groups that received an initial colonoscopic infusion and then intensive multidonor FMT or placebo enemas, 5 d/wk for 8 weeks. Patients in the FMT group received blended homogenized stool from 3–7 unrelated donors. Patients in the placebo group were eligible to receive open-label FMT after the double-blind study period. We collected 314 fecal samples from the patients at screening, every 4 weeks during treatment, and 8 weeks after the blinded or open-label FMT therapy. We also collected 160 large-bowel biopsy samples from the patients at study entry, at completion of 8 weeks of blinded therapy, and at the end of open-label FMT, if applicable. We analyzed 105 fecal samples from the 14 individual donors (n = 55), who in turn contributed to 21 multidonor batches (n = 50). Bacteria in colonic and fecal samples were analyzed by both 16S ribosomal RNA gene and transcript amplicon sequencing; 285 fecal samples were analyzed by shotgun metagenomics, and 60 fecal samples were analyzed for metabolome features. Results FMT increased microbial diversity and altered composition, based on analyses of colon and fecal samples collected before vs after FMT. Diversity was greater in fecal and colon samples collected before and after FMT treatment from patients who achieved remission compared with patients who did not. Patients in remission after FMT had enrichment of Eubacterium hallii and Roseburia inulivorans compared with patients who did not achieve remission after FMT and had increased levels of short-chain fatty acid biosynthesis and secondary bile acids. Patients who did not achieve remission had enrichment of Fusobacterium gonidiaformans, Sutterella wadsworthensis, and Escherichia species and increased levels of heme and lipopolysaccharide biosynthesis. Bacteroides in donor stool were associated with remission in patients receiving FMT, and Streptococcus species in donor stool was associated with no response to FMT. Conclusions In an analysis of fecal and colonic mucosa samples from patients receiving FMT for active UC and stool samples from donors, we associated specific bacteria and metabolic pathways with induction of remission. These findings may be of value in the design of microbe-based therapies for UC. ClinicalTrials.gov, Number NCT01896635

© 2017 The Author(s) Myelodysplastic syndromes and chronic myelomonocytic leukemia are blood disorders characterized by ineffective hematopoiesis and progressive marrow failure that can transform into acute leukemia. The DNA methyltransferase inhibitor 5-azacytidine (AZA) is the most effective pharmacological option, but only ∼50% of patients respond. A response only manifests after many months of treatment and is transient. The reasons underlying AZA resistance are unknown, and few alternatives exist for non-responders. Here, we show that AZA responders have more hematopoietic progenitor cells (HPCs) in the cell cycle. Non-responder HPC quiescence is mediated by integrin α5 (ITGA5) signaling and their hematopoietic potential improved by combining AZA with an ITGA5 inhibitor. AZA response is associated with the induction of an inflammatory response in HPCs in vivo. By molecular bar coding and tracking individual clones, we found that, although AZA alters the sub-clonal contribution to different lineages, founder clones are not eliminated and continue to drive hematopoiesis even in complete responders.

RNA splicing is a fundamental biological process that generates protein diversity from a finite set of genes. Recurrent somatic mutations of genes involved in RNA splicing are present at high frequency in Myelodysplasia (up to 70%) but less so in Acute Myeloid Leukemia (AML; less than 20%). To investigate whether there were aberrant and recurrent RNA splicing events in the AML transcriptome that were associated with poor prognosis in the absence of splicing factor mutations, we developed a bioinformatics pipeline to systematically annotate and quantify alternative splicing events from RNA-sequencing data (Fig A). We first analysed publicly available RNA-seq data from The Cancer Genome Atlas (TCGA, n=170). We focussed on non-M3 AML patients with no splicing factor mutations (based on reported genomic sequencing and verified by re-analysis of RNA-seq data from all patients) who had received intensive chemotherapy. We segregated these patients based on their European Leukaemia Net (ELN) risk classification and identified 1290 alternatively spliced events impacting 910 genes that were significantly different (FDR In order to prioritise those alternatively spliced events most likely to have a deleterious function, we developed an analytical framework to predict their impact on protein structure (Fig E). 87 alternatively spliced events, 25.81% of the commonly shared splicing events, relating to 78 genes (35.13% of all genes) were predicted to directly alter highly conserved protein domains within the affected genes, leading to either a complete (~25%, Fig E) or a partial loss of a domain (20%, Fig E). These in silico predictions are likely to be an underestimate of the true impact, as splicing alterations mapping to poorly annotated domains or affecting the tertiary structure of proteins would be missed. A number of splicing factors themselves were differentially spliced, with the alternative splicing predicted to have functional consequences. This was exemplified by hnRNPA1, a factor with well-established roles in splicing, is itself alternatively spliced in patients and predicted to be deleterious. Consistent with this, motif scanning analyses indicated that a number of mis-spliced transcripts had hnRNPA1 binding motifs (Fig F). To assess the impact of these alternatively spliced events (that were predicted to also disrupt highly conserved protein domains) on the transcriptome, we simultaneously quantified differential gene expression. IPA analysis of the 602 genes that were differentially expressed between ELNAdv and ELNFav patients and shared between both TCGA and ClinSeq cohorts indicated that they were associated with pathways (Fig G) that were distinct from those associated with aberrantly spliced genes (Fig D). A number of pathways related to inflammation were enriched amongst the genes observed to be upregulated in ELNAdv patients (Fig G). Network analyses integrating the alternatively spliced genes with differentially expressed genes revealed strong interactions (Fig H), indicating functional associations between these biological events. Given these strong network interactions, we investigated the potential prognostic significance of these alternatively spliced events. To this end, we utilised machine-learning methods to derive a "splicing signature" of four mis-spliced genes with a predictive capacity equivalent to the ELN (Fig I). The splicing signature further refined existing risk prediction algorithms to improve the classification of patients (Fig J). Taken together, we report the presence of extensive deregulation of RNA splicing in AML patients even in the absence of splicing factor mutations. Many of these events were shared in patients with adverse outcomes and their impact on the AML transcriptome points towards vulnerabilities that could be targeted. Figure Disclosures Unnikrishnan: Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding. Lehmann:TEVA: Consultancy, Membership on an entity's Board of Directors or advisory committees; Pfizer: Membership on an entity's Board of Directors or advisory committees; Abbive: Membership on an entity's Board of Directors or advisory committees. Pimanda:Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding.

The widespread adoption of high-throughput next-generation sequencing (NGS) technology among the Australian life science research community is highlighting an urgent need to up-skill biologists in tools required for handling and analysing their NGS data. There is currently a shortage of cutting-edge bioinformatics training courses in Australia as a consequence of a scarcity of skilled trainers with time and funding to develop and deliver training courses. To address this, a consortium of Australian research organizations, including Bioplatforms Australia, the Commonwealth Scientific and Industrial Research Organisation and the Australian Bioinformatics Network, have been collaborating with EMBL-EBI training team. A group of Australian bioinformaticians attended the train-the-trainer workshop to improve training skills in developing and delivering bioinformatics workshop curriculum. A 2-day NGS workshop was jointly developed to provide hands-on knowledge and understanding of typical NGS data analysis workflows. The road show–style workshop was successfully delivered at five geographically distant venues in Australia using the newly established Australian NeCTAR Research Cloud. We highlight the challenges we had to overcome at different stages from design to delivery, including the establishment of an Australian bioinformatics training network and the computing infrastructure and resource development. A virtual machine image, workshop materials and scripts for configuring a machine with workshop contents have all been made available under a Creative Commons Attribution 3.0 Unported License. This means participants continue to have convenient access to an environment they had become familiar and bioinformatics trainers are able to access and reuse these resources.

The epithelial response to the opportunistic pathogen Campylobacter concisus is poorly characterised. Here, we assessed the intestinal epithelial responses to two C. concisus strains with different virulence characteristics in Caco-2 cells using RNAseq, and validated a subset of the response using qPCR arrays. C. concisus strains induced distinct response patterns from intestinal epithelial cells, with the toxigenic strain inducing a significantly more amplified response. A range of cellular functions were significantly regulated in a strain-specific manner, including epithelial-to-mesenchymal transition (NOTCH and Hedgehog), cytoskeletal remodeling, tight junctions, inflammatory responses and autophagy. Pattern recognition receptors were regulated, including TLR3 and IFI16, suggesting that nucleic acid sensing was important for epithelial recognition of C. concisus. C. concisus zonula occludens toxin (ZOT) was expressed and purified, and the epithelial response to the toxin was analysed using RNAseq. ZOT upregulated PAR2 expression, as well as processes related to tight junctions and cytoskeletal remodeling. C. concisus ZOT also induced upregulation of TLR3, pro-inflammatory cytokines IL6, IL8 and chemokine CXCL16, as well as the executioner caspase CASP7. Here, we characterise distinct global epithelial responses to C. concisus strains, and the virulence factor ZOT, and provide novel information on mechanisms by which this bacterium may affect the host.

Tilletia indica causes the disease Karnal bunt in wheat. The disease is under international quarantine regulations. Comparative mitochondrial (mt) genome analysis of T. indica (KX394364 and DQ993184) and T. walkeri (EF536375) has found 325 to 328 SNPs, 57 to 60 short InDels (from 1 to 13 nt), two InDels (30 and 61 nt) and five (>200 nt) presence/absence variations (PAVs) between the two species. The mt genomes of both species have identical gene order. The numbers of SNPs and InDels between the mt genomes of the two species are approximately nine times of the corresponding numbers between the two T. indica isolates. There are eight SNPs between T. indica and T. walkeri that resulted in amino acid substitutions in the mt genes of cob, nad2 and nad5. In contrast, there is no amino acid substitution in the mt genes of the T. indica isolates from the SNPs found. The five PAVs present in T. indica (DQ993184) are absent in T. walkeri. Four PAVs are more than 1 kb and are not present in every T. indica isolate. Analysis of their presence and absence separates a collection of T. indica isolates into 11 subgroups. Two PAVs have ORFs for the LAGLIDAG endonuclease and two have ORFs for the GIY-YIG endonuclease family, which are representatives of homing endonuclease genes (HEGs). These intron- encoded HEGs confer intron mobility and account for their fluid distribution in T. indica isolates. The small PAV of 221 bp, present in every T. indica isolate and unique to the species, was used as the genetic fingerprint for the successful development of a rapid, highly sensitive and specific loop mediated isothermal amplification (LAMP) assay. The simple procedure of the LAMP assay and the easy detection formats will enable the assay to be automated for high throughput diagnosis.

With the emergence of glycoproteomics, there is a need to develop bioinformatic tools to identify glycopeptides in protease digests of glycoproteins. GlycoSpectrumScan is a web-based tool that identifies the glycoheterogeneity on a peptide from mass spectrometric data. Two experimental data sets are required as inputs: (1) oligosaccharide compositions of the N- and/or O-linked glycans present in the sample and (2) in silico derived peptide masses of proteolytically digested proteins with a potential number of N- and/or O-glycosylation sites. GlycoSpectrumScan uses MS data, rather than MS/MS data, to identify glycopeptides and determine the relative distribution of N- and O-glycoforms at each site. It is functional for assigning monosaccharide compositions on glycopeptides with single and multiple sites of glycosylation. The algorithm allows the input of raw mass data, including multiply charged ions, making it applicable for both ESI and MALDI data from all mass spectrometer platforms. Manual analysis time for identifying glycosylation heterogeneity at each site on glycoprotein(s) is substantially decreased. The application of this tool to characterize the N- and O-linked glycopeptides from human secretory IgA (sIgA), consisting of secretory component (7 N-linked sites), IgA1 (2 N-linked,or=5 O-linked sites), IgA2 (4 N-linked sites) and J-chain (1 N-linked site) is described. GlycoSpectrumScan is freely available at www.glycospectrumscan.org .

Glycosylation of proteins is important for protein stability, secretion, and localization. In this study, we have investigated the glycan synthesis pathways of 12 filamentous fungi including those of medical/agricultural/industrial importance for which genomes have been recently sequenced. We have adopted a systems biology approach to combine the results from comparative genomics techniques with high confidence information on the enzymes and fungal glycan structures, reported in the literature. From this, we have developed a composite representation of the glycan synthesis pathways in filamentous fungi (both N- and O-linked). The N-glycosylation pathway in the cytoplasm and endoplasmic reticulum was found to be highly conserved evolutionarily across all the filamentous fungi considered in the study. In the final stages of N-glycan synthesis in the Golgi, filamentous fungi follow the high mannose pathway as in Saccharomyces cerevisiae, but the level of glycan mannosylation is reduced. Highly specialized N-glycan structures with galactofuranose residues, phosphodiesters, and other insufficiently trimmed structures have also been identified in the filamentous fungi. O-Linked glycosylation in filamentous fungi was seen to be highly conserved with many mannosyltransferases that are similar to those in S. cerevisiae. However, highly variable and diverse O-linked glycans also exist. We have developed a web resource for presenting the compiled data with user-friendly query options, which can be accessed at www.fungalglycans.org. This resource can assist attempts to remodel glycosylation of recombinant proteins expressed in filamentous fungal hosts.

The B2 phylogenetic group of Escherichia coli contains important pathogens such as extraintestinal pathogenic, adherent-invasive, and uropathogenic strains. In this study, we used comparative genomics and statistical methods to identify genetic variations that define a subset of pathogenic strains belonging to the B2 phylogenetic group. An initial proof of concept analysis indicated that five of the 62 E. coli strains available in the Kyoto Encyclopedia of Genes and Genomes database showed close association with B2 adherent-invasive E. coli, forming a subgroup within the B2 phylogenetic group. The tool, kSNP which uses a k-mer approach, and the statistical phenotype prediction tool PPFS2 were then employed to identify 29 high-resolution SNPs, which reaffirmed this grouping. PPFS2 analysis also provided indications that the clustering of this subgroup was highly consistent, and thus, could have a strong phenotypic basis rather than being only evolutionary. Protein homology analyses identified three proteins to be conserved across this subgrouping, two CRISPR-Cas proteins and a hypothetical protein. Functional analyses of these genetic and protein variations may provide insights into the phenotype of these strains.

We present the first analysis of the human proteome with regard to interactions between proteins. We also compare the human interactome with the available interaction datasets from yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans) and fly (Drosophila melanogaster). Of >70,000 binary interactions, only 42 were common to human, worm and fly, and only 16 were common to all four datasets. An additional 36 interactions were common to fly and worm but were not observed in humans, although a coimmunoprecipitation assay showed that 9 of the interactions do occur in humans. A re-examination of the connectivity of essential genes in yeast and humans indicated that the available data do not support the presumption that the number of interaction partners can accurately predict whether a gene is essential. Finally, we found that proteins encoded by genes mutated in inherited genetic disorders are likely to interact with proteins known to cause similar disorders, suggesting the existence of disease subnetworks. The human interaction map constructed from our analysis should facilitate an integrative systems biology approach to elucidating the cellular networks that contribute to health and disease states.

Pathogenic species within the genus Campylobacter are responsible for a considerable burden on global health. Campylobacter concisus is an emergent pathogen that plays a role in acute and chronic gastrointestinal disease. Despite ongoing research on Campylobacter virulence mechanisms, little is known regarding the immunological profile of the host response to Campylobacter infection. In this study, we describe a comprehensive global profile of innate immune responses to C. concisus infection in differentiated THP-1 macrophages infected with an adherent and invasive strain of C. concisus . Using RNA sequencing (RNA-seq), quantitative PCR (qPCR), mass spectrometry, and confocal microscopy, we observed differential expression of pattern recognition receptors and robust upregulation of DNA- and RNA-sensing molecules. In particular, we observed IFI16 inflammasome assembly in C. concisus -infected macrophages. Global profiling of the transcriptome revealed the significant regulation of a total of 8,343 transcripts upon infection with C. concisus , which included the activation of key inflammatory pathways involving CREB1, NF-κB, STAT, and interferon regulatory factor signaling. Thirteen microRNAs and 333 noncoding RNAs were significantly regulated upon infection, including MIR221, which has been associated with colorectal carcinogenesis. This study represents a major advance in our understanding of host recognition and innate immune responses to infection by C. concisus .

Direct links between proteomic and genomic/transcriptomic data are not frequently made, partly because of lack of appropriate bioinformatics tools. To help address this, we have developed the PG Nexus pipeline. The PG Nexus allows users to covisualize peptides in the context of genomes or genomic contigs, along with RNA-seq reads. This is done in the Integrated Genome Viewer (IGV). A Results Analyzer reports the precise base position where LC-MS/MS-derived peptides cover genes or gene isoforms, on the chromosomes or contigs where this occurs. In prokaryotes, the PG Nexus pipeline facilitates the validation of genes, where annotation or gene prediction is available, or the discovery of genes using a "virtual protein"-based unbiased approach. We illustrate this with a comprehensive proteogenomics analysis of two strains of Campylobacter concisus . For higher eukaryotes, the PG Nexus facilitates gene validation and supports the identification of mRNA splice junction boundaries and splice variants that are protein-coding. This is illustrated with an analysis of splice junctions covered by human phosphopeptides, and other examples of relevance to the Chromosome-Centric Human Proteome Project. The PG Nexus is open-source and available from https://github.com/IntersectAustralia/ap11_Samifier . It has been integrated into Galaxy and made available in the Galaxy tool shed.

The chloroform-respiring bacterium Dehalobacter UNSWDHB was isolated from subsurface soil contaminated with a mixture of organohalides, including chloroform. Here, we present its 3.2-Mb genome .

Background In spite of its association with gastroenteritis and inflammatory bowel diseases, the isolation of Campylobacter concisus from both diseased and healthy individuals has led to controversy regarding its role as an intestinal pathogen. One proposed reason for this is the presence of high genetic diversity among the genomes of C. concisus strains. Results In this study the genomes of six C. concisus strains were sequenced, assembled and annotated including two strains isolated from Crohn’s disease patients (UNSW2 and UNSW3), three from gastroenteritis patients (UNSW1, UNSWCS and ATCC 51562) and one from a healthy individual (ATCC 51561). The genomes of C. concisus BAA-1457 and UNSWCD, available from NCBI, were included in subsequent comparative genomic analyses. The Pan and Core genomes for the sequenced C. concisus strains consisted of 3254 and 1556 protein coding genes, respectively. Conclusion Genes were identified with specific conservation in C. concisus strains grouped by phenotypes such as invasiveness, adherence, motility and diseased states. Phylogenetic trees based on ribosomal RNA sequences and concatenated host-related pathways for the eight C. concisus strains were generated using the neighbor-joining method, of which the 16S rRNA gene and peptidoglycan biosynthesis grouped the C. concisus strains according to their pathogenic phenotypes. Furthermore, 25 non-synonymous amino acid changes with 14 affecting functional domains, were identified within proteins of conserved host-related pathways, which had possible associations with the pathogenic potential of C. concisus strains. Finally, the genomes of the eight C. concisus strains were compared to the nine available genomes of the well-established pathogen Campylobacter jejuni, which identified several important differences in the respiration pathways of these two species. Our findings indicate that C. concisus strains are genetically diverse, and suggest the genomes of this bacterium contain respiration pathways and modifications in the peptidoglycan layer that may play an important role in its virulence.

Campylobacter showae UNSWCD was isolated from a patient with Crohn's disease. Here we present a 2.1 Mb draft assembly of its genome.

Campylobacter concisus has garnered increasing attention due to its association with intestinal disease, thus, the pathogenic potential of strains isolated from different intestinal diseases was investigated. A method to isolate C. concisus was developed and the ability of eight strains from chronic and acute intestinal diseases to adhere to and invade intestinal epithelial cells was determined. Features associated with bacterial invasion were investigated using comparative genomic analyses and the effect of C. concisus on host protein expression was examined using proteomics. Our isolation method from intestinal biopsies resulted in the isolation of three C. concisus strains from children with Crohn's disease or chronic gastroenteritis. Four C. concisus strains from patients with chronic intestinal diseases can attach to and invade host cells using mechanisms such as chemoattraction to mucin, aggregation, flagellum-mediated attachment, "membrane ruffling", cell penetration and damage. C. concisus strains isolated from patients with chronic intestinal diseases have significantly higher invasive potential than those from acute intestinal diseases. Investigation of the cause of this increased pathogenic potential revealed a plasmid to be responsible. 78 and 47 proteins were upregulated and downregulated in cells infected with C. concisus, respectively. Functional analysis of these proteins showed that C. concisus infection regulated processes related to interleukin-12 production, proteasome activation and NF-κB activation. Infection with all eight C. concisus strains resulted in host cells producing high levels of interleukin-12, however, only strains capable of invading host cells resulted in interferon-γ production as confirmed by ELISA. These findings considerably support the emergence of C. concisus as an intestinal pathogen, but more significantly, provide novel insights into the host immune response and an explanation for the heterogeneity observed in the outcome of C. concisus infection. Moreover, response to infection with invasive strains has substantial similarities to that observed in the inflamed mucosa of Crohn's disease patients.

The analysis of expressed sequence tag (EST) datasets offers a rapid and cost-effective approach to elucidate the transcriptome of an organism, but requiring several computational methods for assembly and annotation. ESTExplorer is a comprehensive workflow system for EST data management and analysis. The pipeline uses a ‘distributed control approach’ in which the most appropriate bioinformatics tools are implemented over different dedicated processors. Species-specific repeat masking and conceptual translation are in-built. ESTExplorer accepts a set of ESTs in FASTA format which can be analysed using programs selected by the user. After pre-processing and assembly, the dataset is annotated at the nucleotide and protein levels, following conceptual translation. Users may optionally provide ESTExplorer with assembled contigs for annotation purposes. Functionally annotated contigs/ESTs can be analysed individually. The overall outputs are gene ontologies, protein functional identifications in terms of mapping to protein domains and metabolic pathways. ESTExplorer has been applied successfully to annotate large EST datasets from parasitic nematodes and to identify novel genes as potential targets for parasite intervention. ESTExplorer runs on a Linux cluster and is freely available for the academic community at http://estexplorer.biolinfo.org.