Search

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.

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.

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

Background The esophageal microbiome has been proposed to be involved in a range of diseases including the esophageal adenocarcinoma cascade; however, little is currently known about its function and relationship to the host. Here, the esophageal microbiomes of 106 prospectively recruited patients were assessed using 16S rRNA and 18S rRNA amplicon sequencing as well as shotgun sequencing, and associations with age, gender, proton pump inhibitor use, host genetics, and disease were tested. Results The esophageal microbiome was found to cluster into functionally distinct community types (esotypes) defined by the relative abundances of Streptococcus and Prevotella. While age was found to be a significant factor driving microbiome composition, bacterial signatures and functions such as enrichment with Gram-negative oral-associated bacteria and microbial lactic acid production were associated with the early stages of the esophageal adenocarcinoma cascade. Non-bacterial microbes such as archaea, Candida spp., and bacteriophages were also identified in low abundance in the esophageal microbiome. Specific host SNPs in NOTCH2, STEAP2-AS1, and NREP were associated with the composition of the esophageal microbiome in our cohort. Conclusions This study provides the most comprehensive assessment of the esophageal microbiome to date and identifies novel signatures and host markers that can be investigated further in the context of esophageal adenocarcinoma development. Electronic supplementary material The online version of this article (10.1186/s40168-018-0611-4) contains supplementary material, which is available to authorized users.

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.

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.

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.

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.

Background Mycoparasitism, a lifestyle where one fungus is parasitic on another fungus, has special relevance when the prey is a plant pathogen, providing a strategy for biological control of pests for plant protection. Probably, the most studied biocontrol agents are species of the genus Hypocrea/Trichoderma. Results Here we report an analysis of the genome sequences of the two biocontrol species Trichoderma atroviride (teleomorph Hypocrea atroviridis) and Trichoderma virens (formerly Gliocladium virens, teleomorph Hypocrea virens), and a comparison with Trichoderma reesei (teleomorph Hypocrea jecorina). These three Trichoderma species display a remarkable conservation of gene order (78 to 96%), and a lack of active mobile elements probably due to repeat-induced point mutation. Several gene families are expanded in the two mycoparasitic species relative to T. reesei or other ascomycetes, and are overrepresented in non-syntenic genome regions. A phylogenetic analysis shows that T. reesei and T. virens are derived relative to T. atroviride. The mycoparasitism-specific genes thus arose in a common Trichoderma ancestor but were subsequently lost in T. reesei. Conclusions The data offer a better understanding of mycoparasitism, and thus enforce the development of improved biocontrol strains for efficient and environmentally friendly protection of plants.