MicroRNA expression profile in RAW264.7 cells in response to Brucella melitensis infection

MicroRNA (miRNA) is small non-coding RNA with approximate 22 nt in length. Recent studies indicate that miRNAs play significant roles in pathogen-host interactions. Brucella organisms are Gram-negative facultative intracellular bacteria that cause Brucellosis. Brucella strains infect macrophages and establish chronic infection by altering host life activities including apoptosis and autophagy. Here, we report a comprehensive analysis of miRNA expression profiles in mock- and Brucella-infected RAW264.7 cells using high-throughput sequencing approach. In total, 344 unique miRNAs were co-expressed in the two libraries, in which 57 miRNAs were differentially expressed. Eight differentially expressed miRNAs with high abundance were subjected to further analysis. The GO enrichment analysis suggests that the putative target genes of these differentially expressed miRNAs are involved in apoptosis, autophagy and immune response. In particular, a total of 25 target genes are involved in regulating apoptosis and autophagy, indicating that these miRNAs may play important regulatory roles in the Brucella-host interactions. Furthermore, the interactions of miR-1981 and its target genes, Bcl-2 and Bid, were validated by luciferase assay. The results show that miR-1981 mimic up-regulated the luciferase activity of psiCHECK-2 Bcl-2 3' UTR, but the luciferase activity of psiCHECK-2 Bid 3' UTR was not changed significantly. Taken together, these data provide valuable framework on Brucella induced miRNA expression in RAW264.7 cells, and suggest that Brucella may establish chronic infection by regulating miRNA expression profile.     

Identification of Recombination and Positively Selected Genes in Brucella

Brucella is a facultative intracellular bacterium belongs to the class alpha proteobacteria. It causes zoonotic disease brucellosis to wide range of animals. Brucella species are highly conserved in nucleotide level.Here, we employed a comparative genomics approach to examine the role of homologous recombination and positive selection in the evolution of Brucella. For the analysis, we have selected 19 complete genomes from 8 species of Brucella. Among the 1599 core genome predicted, 24 genes were showing signals of recombination but no significant breakpoint was found. The analysis revealed that recombination events are less frequent and the impact of recombination occurred is negligible on the evolution of Brucella. This leads to the view that Brucella is clonally evolved. On other hand, 56 genes (3.5 % of core genome) were showing signals of positive selection. Results suggest that natural selection plays an important role in the evolution of Brucella. Some of the genes that are responsible for the pathogenesis of Brucella were found positively selected, presumably due to their role in avoidance of the host immune system.     

Comparative genome analysis of pathogenic and non-pathogenic Clavibacter strains reveals adaptations to their lifestyle

Background

The genus Clavibacter harbors economically important plant pathogens infecting agricultural crops such as potato and tomato. Although the vast majority of Clavibacter strains are pathogenic, there is an increasing number of non-pathogenic isolates reported. Non-pathogenic Clavibacter strains isolated from tomato seeds are particularly problematic because they affect the current detection and identification tests for Clavibacter michiganensis subsp. michiganensis (Cmm), which is regulated with a zero tolerance in tomato seed. Their misidentification as pathogenic Cmm hampers a clear judgment on the seed quality and health.

Results

To get more insight in the genetic features linked to the lifestyle of these bacteria, a whole-genome sequence of the tomato seed-borne non-pathogenic Clavibacter LMG 26808 was determined. To gain a better understanding of the molecular determinants of pathogenicity, the genome sequence of LMG 26808 was compared with that of the pathogenic Cmm strain (NCPPB 382). The comparative analysis revealed that LMG 26808 does not contain plasmids pCM1 and pCM2 and also lacks the majority of important virulence factors described so far for pathogenic Cmm. This explains its apparent non-pathogenic nature in tomato plants. Moreover, the genome analysis of LMG 26808 detected sequences from a plasmid originating from a member of Enterobacteriaceae/Klebsiella relative. Genes received that way and coding for antibiotic resistance may provide a competitive advantage for survival of LMG 26808 in its ecological niche. Genetically, LMG 26808 was the most similar to the pathogenic Cmm NCPPB 382 but contained more mobile genetic elements. The genome of this non-pathogenic Clavibacter strain contained also a high number of transporters and regulatory genes.

Conclusions

The genome sequence of the non-pathogenic Clavibacter strain LMG 26808 and the comparative analyses with other pathogenic Clavibacter strains provided a better understanding of the genetic bases of virulence and adaptation mechanisms present in the genus Clavibacter.

Electronic supplementary material

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

Genotyping of Giardia in Dutch patients and animals: a phylogenetic analysis of human and animal isolates

Giardia duodenalis (syn. Giardia lamblia, Giardia intestinalis) is a protozoan organism that can infect the intestinal tract of many animal species including mammals. Genetic heterogeneity of G. duodenalis is well described but the zoonotic potential is still not clear. In this study, we analysed 100 Giardia DNA samples directly isolated from human stool specimens, to get more insight in the different G. duodenalis assemblages present in the Dutch human population. Results showed that these human isolates could be divided into two main Assemblages A and B within the G. duodenalis group on the basis of PCR assays specific for the Assemblages A and B and the DNA sequences of 18S ribosomal RNA and the glutamate dehydrogenase (gdh) genes. Genotyping results showed that G. duodenalis isolates originating from Dutch human patients belonged in 35% of the cases to Assemblage A (34/98) and in 65% of the cases to Assemblage B (64/98) whereas two human cases remained negative in all assays tested. In addition, we compared these human samples with animal samples from the Netherlands and human and animal samples from other countries. A phylogenetic analysis was carried out on the DNA sequences obtained from these Giardia and those available in GenBank. Using gdh DNA sequence analysis, human and animal Assemblage A and B Giardia isolates could be identified. However, phylogenetic analysis revealed different sub-clustering for human and animal isolates where host-species-specific assemblages (C, D, E, F and G) could be identified. The geographic origin of the human and animal samples was not a discriminating factor.     

A Brucella Virulence Factor Targets Macrophages to Trigger B-cell Proliferation

Brucella spp. and Trypanosoma cruzi are two intracellular pathogens that have no evolutionary common origins but share a similar lifestyle as they establish chronic infections for which they have to circumvent the host immune response. Both pathogens have a virulence factor (prpA in Brucella and tcPrac in T. cruzi) that induces B-cell proliferation and promotes the establishment of the chronic phase of the infectious process. We show here that, even though PrpA promotes B-cell proliferation, it targets macrophages in vitro and is translocated to the cytoplasm during the intracellular replication phase. We observed that PrpA-treated macrophages induce the secretion of a soluble factor responsible for B-cell proliferation and identified nonmuscular myosin IIA (NMM-IIA) as a receptor required for binding and function of this virulence factor. Finally, we show that the Trypanosoma cruzi homologue of PrpA also targets macrophages to induce B-cell proliferation through the same receptor, indicating that this virulence strategy is conserved between a bacterial and a protozoan pathogen.     

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

Background

In recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.

Results

Here, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.

Conclusion

The significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1190-9) contains supplementary material, which is available to authorized users.     

Clostridium sordellii genome analysis reveals plasmid localized toxin genes encoded within pathogenicity loci

Background

Clostridium sordellii can cause severe infections in animals and humans, the latter associated with trauma, toxic shock and often-fatal gynaecological infections. Strains can produce two large clostridial cytotoxins (LCCs), TcsL and TcsH, related to those produced by Clostridium difficile, Clostridium novyi and Clostridium perfringens, but the genetic basis of toxin production remains uncharacterised.

Results

Phylogenetic analysis of the genome sequences of 44 strains isolated from human and animal infections in the UK, US and Australia placed the species into four clades. Although all strains originated from animal or clinical disease, only 5 strains contained LCC genes: 4 strains contain tcsL alone and one strain contains tcsL and tcsH. Four toxin-positive strains were found within one clade. Where present, tcsL and tcsH were localised in a pathogenicity locus, similar to but distinct from that present in C. difficile. In contrast to C. difficile, where the LCCs are chromosomally localised, the C. sordellii tcsL and tcsH genes are localised on plasmids. Our data suggest gain and loss of entire toxigenic plasmids in addition to horizontal transfer of the pathogenicity locus. A high quality, annotated sequence of ATCC9714 reveals many putative virulence factors including neuraminidase, phospholipase C and the cholesterol-dependent cytolysin sordellilysin that are highly conserved between all strains studied.

Conclusions

Genome analysis of C. sordellii reveals that the LCCs, the major virulence factors, are localised on plasmids. Many strains do not contain the LCC genes; it is probable that in several of these cases the plasmid has been lost upon laboratory subculture. Our data are consistent with LCCs being the primary virulence factors in the majority of infections, but LCC-negative strains may precipitate certain categories of infection. A high quality genome sequence reveals putative virulence factors whose role in virulence can be investigated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1613-2) contains supplementary material, which is available to authorized users.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

Characterization of the complete chloroplast genome of Hevea brasiliensis reveals genome rearrangement, RNA editing sites and phylogenetic relationships

Rubber tree (Hevea brasiliensis) is an economical plant and widely grown for natural rubber production. However, genomic research of rubber tree has lagged behind other species in the Euphorbiaceae family. We report the complete chloroplast genome sequence of rubber tree as being 161,191 bp in length including a pair of inverted repeats of 26,810 bp separated by a small single copy region of 18,362 bp and a large single copy region of 89,209 bp. The chloroplast genome contains 112 unique genes, 16 of which are duplicated in the inverted repeat. Of the 112 unique genes, 78 are predicted protein-coding genes, 4 are ribosomal RNA genes and 30 are tRNA genes. Relative to other plant chloroplast genomes, we observed a unique rearrangement in the rubber tree chloroplast genome: a 30-kb inversion between the trnE(UUC)-trnS(GCU) and the trnT(GGU)-trnR(UCU). A comparison between the rubber tree chloroplast genes and cDNA sequences revealed 51 RNA editing sites in which most (48 sites) were located in 26 protein coding genes and the other 3 sites were in introns. Phylogenetic analysis based on chloroplast genes demonstrated a close relationship between Hevea and Manihot in Euphorbiaceae and provided a strong support for a monophyletic group of the eurosid I.     

A novel termini analysis theory using HTS data alone for the identification of Enterococcus phage EF4-like genome termini

Background

Enterococcus faecalis and Enterococcus faecium are typical enterococcal bacterial pathogens. Antibiotic resistance means that the identification of novel E. faecalis and E. faecium phages against antibiotic-resistant Enterococcus have an important impact on public health. In this study, the E. faecalis phage IME-EF4, E. faecium phage IME-EFm1, and both their hosts were antibiotic resistant. To characterize the genome termini of these two phages, a termini analysis theory was developed to provide a wealth of terminal sequence information directly, using only high-throughput sequencing (HTS) read frequency statistics.

Results

The complete genome sequences of phages IME-EF4 and IME-EFm1 were determined, and our termini analysis theory was used to determine the genome termini of these two phages. Results showed 9 bp 3′ protruding cohesive ends in both IME-EF4 and IME-EFm1 genomes by analyzing frequencies of HTS reads. For the positive strands of their genomes, the 9 nt 3′ protruding cohesive ends are 5′-TCATCACCG-3′ (IME-EF4) and 5′-GGGTCAGCG-3′ (IME-EFm1). Further experiments confirmed these results. These experiments included mega-primer polymerase chain reaction sequencing, terminal run-off sequencing, and adaptor ligation followed by run-off sequencing.

Conclusion

Using this termini analysis theory, the termini of two newly isolated antibiotic-resistant Enterococcus phages, IME-EF4 and IME-EFm1, were identified as the byproduct of HTS. Molecular biology experiments confirmed the identification. Because it does not require time-consuming wet lab termini analysis experiments, the termini analysis theory is a fast and easy means of identifying phage DNA genome termini using HTS read frequency statistics alone. It may aid understanding of phage DNA packaging.     

Use of rpoB gene sequence analysis for phylogenetic identification of Cronobacter species

Here, we report the RNA polymerase beta-subunit gene (rpoB) as a new molecular marker for the identification of the Cronobacter species. The results indicated that members of the Cronobacter genus are more easily discriminated by rpoB sequencing than 16S rRNA sequencing, and reliable identification could be achieved by rpoB gene sequence comparison.     

Multilocus phylogenetic analysis of the genus Aeromonas

A broad multilocus phylogenetic analysis (MLPA) of the representative diversity of a genus offers the opportunity to incorporate concatenated inter-species phylogenies into bacterial systematics. Recent analyses based on single housekeeping genes have provided coherent phylogenies of Aeromonas. However, to date, a multi-gene phylogenetic analysis has never been tackled. In the present study, the intra- and inter-species phylogenetic relationships of 115 strains representing all Aeromonas species described to date were investigated by MLPA. The study included the independent analysis of seven single gene fragments (gyrB, rpoD, recA, dnaJ, gyrA, dnaX, and atpD), and the tree resulting from the concatenated 4705 bp sequence. The phylogenies obtained were consistent with each other, and clustering agreed with the Aeromonas taxonomy recognized to date. The highest clustering robustness was found for the concatenated tree (i.e. all Aeromonas species split into 100% bootstrap clusters). Both possible chronometric distortions and poor resolution encountered when using single-gene analysis were buffered in the concatenated MLPA tree. However, reliable phylogenetic species delineation required an MLPA including several “bona fide” strains representing all described species.     

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

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

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

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

Campylobacter group II phage CP21 is the prototype of a new subgroup revealing a distinct modular genome organization and host specificity

Background

The application of phages is a promising tool to reduce the number of Campylobacter along the food chain. Besides the efficacy against a broad range of strains, phages have to be safe in terms of their genomes. Thus far, no genes with pathogenic potential (e.g., genes encoding virulence factors) have been detected in Campylobacter phages. However, preliminary studies suggested that the genomes of group II phages may be diverse and prone to genomic rearrangements.

Results

We determined and analysed the genomic sequence (182,761 bp) of group II phage CP21 that is closely related to the already characterized group II phages CP220 and CPt10. The genomes of these phages are comprised of four modules separated by very similar repeat regions, some of which harbouring open reading frames (ORFs). Though, the arrangement of the modules and the location of some ORFs on the genomes are different in CP21 and in CP220/CPt10. In this work, a PCR system was established to study the modular genome organization of other group II phages demonstrating that they belong to different subgroups of the CP220-like virus genus, the prototypes of which are CP21 and CP220. The subgroups revealed different restriction patterns and, interestingly enough, also distinct host specificities, tail fiber proteins and tRNA genes. We additionally analysed the genome of group II phage vB_CcoM-IBB_35 (IBB_35) for which to date only five individual contigs could be determined. We show that the contigs represent modules linked by long repeat regions enclosing some yet not identified ORFs (e.g., for a head completion protein). The data suggest that IBB_35 is a member of the CP220 subgroup.

Conclusion

Campylobacter group II phages are diverse regarding their genome organization. Since all hitherto characterized group II phages contain numerous genes for transposases and homing endonucleases as well as similar repeat regions, it cannot be excluded that these phages are genetically unstable. To answer this question, further experiments and sequencing of more group II phages should be performed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1837-1) contains supplementary material, which is available to authorized users.