The splicing factor RBM17 drives leukemic stem cell maintenance by evading nonsense-mediated decay of pro-leukemic factors

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.

Raspberry ketone diet supplement reduces attraction of sterile male Queensland fruit fly to cuelure by altering expression of chemoreceptor genes

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.

Multi-omics of the esophageal microenvironment identifies signatures associated with progression of Barrett's esophagus

BACKGROUND

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

METHODS

Here, 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.

RESULTS

Profiling 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 of Campylobacter species 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 between Campylobacter and napsin B aspartic peptidase. Targeted culture of Campylobacter species 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 these Campylobacter strain-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 the Campylobacter isolates revealed expression changes within those infected with strains with high intracellular fitness that could explain the increased likelihood of disease progression.

CONCLUSIONS

We provide a comprehensive assessment of the esophageal microenvironment, identifying bacterial strain-specific signatures with high relevance to progression of metaplasia.

The RNA Atlas expands the catalog of human non-coding RNAs

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.

The RNA Atlas expands the catalog of human non-coding RNAs

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.

Thermostable small-molecule inhibitor of angiogenesis and vascular permeability that suppresses a pERK-FosB/ΔFosB-VCAM-1 axis

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 Alterations Induce a Cellular Stress Response Associated with Poor Prognosis in Acute Myeloid Leukemia

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 ELN^Adv 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.

Fast Short Read De-Novo Assembly Using Overlap-Layout-Consensus Approach

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.

Whole genome sequencing of a novel, dichloromethane-fermenting Peptococcaceae from an enrichment culture

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 SMRTTM sequencing of DCMF allowed de 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 the Peptococcaceae family. 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/.

Little evidence of adaptation potential to ocean acidification in sea urchins living in "Future Ocean" conditions at a CO2 vent

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 F ST 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.

Specific Bacteria and Metabolites Associated With Response to Fecal Microbiota Transplantation in Patients With Ulcerative Colitis

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.

Signatures within the esophageal microbiome are associated with host genetics, age, and disease

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.

Integrative Genomics Identifies the Molecular Basis of Resistance to Azacitidine Therapy in Myelodysplastic Syndromes

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.

Small RNA interactome of pathogenic E. coli revealed through crosslinking of RNase E

RNA sequencing studies have identified hundreds of non‐coding RNAs in bacteria, including regulatory small RNA (sRNA). However, our understanding of sRNA function has lagged behind their identification due to a lack of tools for the high‐throughput analysis of RNA–RNA interactions in bacteria. Here we demonstrate that in vivo sRNA–mRNA duplexes can be recovered using UV‐crosslinking, ligation and sequencing of hybrids (CLASH). Many sRNAs recruit the endoribonuclease, RNase E, to facilitate processing of mRNAs. We were able to recover base‐paired sRNA–mRNA duplexes in association with RNase E, allowing proximity‐dependent ligation and sequencing of cognate sRNA–mRNA pairs as chimeric reads. We verified that this approach captures bona fide sRNA–mRNA interactions. Clustering analyses identified novel sRNA seed regions and sets of potentially co‐regulated target mRNAs. We identified multiple mRNA targets for the pathotype‐specific sRNA Esr41, which was shown to regulate colicin sensitivity and iron transport in E. coli. Numerous sRNA interactions were also identified with non‐coding RNAs, including sRNAs and tRNAs, demonstrating the high complexity of the sRNA interactome.

Novel genetic markers define a subgroup of pathogenic Escherichia coli strains belonging to the B2 phylogenetic group

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.

The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Background

The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations.

Results

We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions.

Conclusions

This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.

Tools to covisualize and coanalyze proteomic data with genomes and transcriptomes: validation of genes and alternative mRNA splicing

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 pathogenic potential of Campylobacter concisus strains associated with chronic intestinal diseases

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.

Comparative analyses of Campylobacter concisus strains reveal the genome of the reference strain BAA-1457 is not representative of the species

Background

Several studies have shown that significant genotypic heterogeneity exists among Campylobacter concisus strains. Recently, the genome of C. concisus UNSWCD, isolated from a patient with Crohn's disease, was sequenced.

Results

In this study, comparative analyses were performed between strain UNSWCD and BAA-1457, isolated from a patient with acute gastroenteritis. Searches between C. concisus UNSWCD and BAA-1457 showed that 76% of genes were homologues, whereas those between C. jejuni strains showed 90-91% to be homologues, indicating substantial variation exists within these two C. concisus genomes. More specific bidirectional homology searches identified 1593 genes that are shared between these strains, and 115 and 281 genes unique to UNSWCD and BAA-1457, respectively. Significantly, differences in the type of flagellin glycosylation pathways between the two strains were identified and confirmed by PCR. The protein profiles of UNSWCD, BAA-1457 and a further six strains of C. concisus were compared and analyzed bioinformatically, and this differentiated the strains into four clades. BAA-1457 was found to be highly divergent (average similarity: 56.8%) from the other seven strains (mean average similarity ± standard deviation: 64.7 ± 1.7%). Furthermore, searches for homologues of the 1593 proteins found to be common between UNSWCD and BAA-1457 were conducted against all available bacterial genomes, and 18 proteins were found to be unique to C. concisus, of which 6 were predicted to be secreted, and may represent good markers for detection of this species.

Conclusions

This study has elucidated several features that may be responsible for the heterogeneity that exists among C. concisus strains, and has determined that the strain BAA-1457 is genetically atypical to other C. concisus strains and is not a good candidate reference strain.