Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

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

Comparative genomic analysis of Acinetobacter baumannii clinical isolates reveals extensive genomic variation and diverse antibiotic resistance determinants

Background

Acinetobacter baumannii is an important nosocomial pathogen that poses a serious health threat to immune-compromised patients. Due to its rapid ability to develop multidrug resistance (MDR), A. baumannii has increasingly become a focus of attention worldwide. To better understand the genetic variation and antibiotic resistance mechanisms of this bacterium at the genomic level, we reported high-quality draft genome sequences of 8 clinical isolates with various sequence types and drug susceptibility profiles.

Results

We sequenced 7 MDR and 1 drug-sensitive clinical A. baumannii isolates and performed comparative genomic analysis of these draft genomes with 16 A. baumannii complete genomes from GenBank. We found a high degree of variation in A. baumannii, including single nucleotide polymorphisms (SNPs) and large DNA fragment variations in the AbaR-like resistance island (RI) regions, the prophage and the type VI secretion system (T6SS). In addition, we found several new AbaR-like RI regions with highly variable structures in our MDR strains. Interestingly, we found a novel genomic island (designated as GI_BJ4) in the drug-sensitive strain BJ4 carrying metal resistance genes instead of antibiotic resistance genes inserted into the position where AbaR-like RIs commonly reside in other A. baumannii strains. Furthermore, we showed that diverse antibiotic resistance determinants are present outside the RIs in A. baumannii, including antibiotic resistance-gene bearing integrons, the bla_OXA-23-containing transposon Tn2009, and chromosomal intrinsic antibiotic resistance genes.

Conclusions

Our comparative genomic analysis revealed that extensive genomic variation exists in the A. baumannii genome. Transposons, genomic islands and point mutations are the main contributors to the plasticity of the A. baumannii genome and play critical roles in facilitating the development of antibiotic resistance in the clinical isolates.

Electronic supplementary material

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

Genomes of Alteromonas australica,a world apart

Background

Alteromonas is a genus of marine bacteria that is very easy to isolate and grow in the laboratory. There are genomes available of the species Alteromonas macleodii from different locations around the world and an Alteromonas sp. isolated from a sediment in Korea. We have analyzed the genomes of two strains classified by 16S rRNA (>99% similarity) as the recently described species Alteromonas australica, and isolated from opposite ends of the world; A. australica DE170 was isolated in the South Adriatic (Mediterranean) at 1000 m depth while A. australica H17^T was isolated from a sea water sample collected in St Kilda Beach, Tasman Sea.

Results

Although these two strains belong to a clearly different species from A. macleodii, the overall synteny is well preserved and the flexible genomic islands seem to code for equivalent functions and be located at similar positions. Actually the genomes of all the Alteromonas species known to date seem to preserve synteny quite well with the only exception of the sediment isolate SN2. Among the specific metabolic features found for the A. australica isolates there is the degradation of xylan and production of cellulose as extracellular polymeric substance by DE170 or the potential ethanol/methanol degradation by H17^T.

Conclusions

The genomes of the two A. australica isolates are not more different than those of strains of A. macleodii isolated from the same sample. Actually the recruitment from metagenomes indicates that all the available genomes are found in most tropical-temperate marine samples analyzed and that they live in consortia of several species and multiple clones within each. Overall the hydrolytic activities of the Alteromonas genus as a whole are impressive and fit with its known capabilities to exploit sudden inputs of organic matter in their environment.

Electronic supplementary material

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

Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus

Background

Vibrio parahaemolyticus is a Gram-negative halophilic bacterium. Infections with the bacterium could become systemic and can be life-threatening to immunocompromised individuals. Genome sequences of a few clinical isolates of V. parahaemolyticus are currently available, but the genome dynamics across the species and virulence potential of environmental strains on a genome-scale have not been described before.

Results

Here we present genome sequences of four V. parahaemolyticus clinical strains from stool samples of patients and five environmental strains in Hong Kong. Phylogenomics analysis based on single nucleotide polymorphisms revealed a clear distinction between the clinical and environmental isolates. A new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus was found in clincial strains. In addition, a novel small genomic island frequently found among clinical isolates was reported. A few environmental strains were found harboring virulence genes and prophage elements, indicating their virulence potential. A unique biphenyl degradation pathway was also reported. A database for V. parahaemolyticus (http://kwanlab.bio.cuhk.edu.hk/vp) was constructed here as a platform to access and analyze genome sequences and annotations of the bacterium.

Conclusions

We have performed a comparative genomics analysis of clinical and environmental strains of V. parahaemolyticus. Our analyses could facilitate understanding of the phylogenetic diversity and niche adaptation of this bacterium.

Electronic supplementary material

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

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

Background

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

Results

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409–1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

Conclusions

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

Electronic supplementary material

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

Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human

Background

The community composition of the human microbiome is known to vary at distinct anatomical niches. But little is known about the nature of variations, if any, at the genome/sub-genome levels of a specific microbial community across different niches. The present report aims to explore, as a case study, the variations in gene repertoire of 28 Prevotella reference genomes derived from different body-sites of human, as reported earlier by the Human Microbiome Consortium.

Results

The pan-genome for Prevotella remains “open”. On an average, 17% of predicted protein-coding genes of any particular Prevotella genome represent the conserved core genes, while the remaining 83% contribute to the flexible and singletons. The study reveals exclusive presence of 11798, 3673, 3348 and 934 gene families and exclusive absence of 17, 221, 115 and 645 gene families in Prevotella genomes derived from human oral cavity, gastro-intestinal tracts (GIT), urogenital tract (UGT) and skin, respectively. Distribution of various functional COG categories differs significantly among the habitat-specific genes. No niche-specific variations could be observed in distribution of KEGG pathways.

Conclusions

Prevotella genomes derived from different body sites differ appreciably in gene repertoire, suggesting that these microbiome components might have developed distinct genetic strategies for niche adaptation within the host. Each individual microbe might also have a component of its own genetic machinery for host adaptation, as appeared from the huge number of singletons.

Electronic supplementary material

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

High resolution assembly and characterization of genomes of Canadian isolates of Salmonella Enteritidis

Background

There is a need to characterize genomes of the foodborne pathogen, Salmonella enterica serovar Enteritidis (SE) and identify genetic information that could be ultimately deployed for differentiating strains of the organism, a need that is yet to be addressed mainly because of the high degree of clonality of the organism. In an effort to achieve the first characterization of the genomes of SE of Canadian origin, we carried out massively parallel sequencing of the nucleotide sequence of 11 SE isolates obtained from poultry production environments (n = 9), a clam and a chicken, assembled finished genomes and investigated diversity of the SE genome.

Results

The median genome size was 4,678,683 bp. A total of 4,833 chromosomal genes defined the pan genome of our field SE isolates consisting of 4,600 genes present in all the genomes, i.e., core genome, and 233 genes absent in at least one genome (accessory genome). Genome diversity was demonstrable by the presence of 1,360 loci showing single nucleotide polymorphism (SNP) in the core genome which was used to portray the genetic distances by means of a phylogenetic tree for the SE isolates. The accessory genome consisted mostly of previously identified SE prophage sequences as well as two, apparently full- sized, novel prophages namely a 28 kb sequence provisionally designated as SE-OLF-10058 (3) prophage and a 43 kb sequence provisionally designated as SE-OLF-10012 prophage.

Conclusions

The number of SNPs identified in the relatively large core genome of SE is a reflection of substantial diversity that could be exploited for strain differentiation as shown by the development of an informative phylogenetic tree. Prophage sequences can also be exploited for SE strain differentiation and lineage tracking. This work has laid the ground work for further studies to develop a readily adoptable laboratory test for the subtyping of SE.

Electronic supplementary material

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

Mutations upstream of fabI in triclosan resistant Staphylococcus aureus strains are associated with elevated fabI gene expression

Background

The enoyl-acyl carrier protein (ACP) reductase enzyme (FabI) is the target for a series of antimicrobial agents including novel compounds in clinical trial and the biocide triclosan. Mutations in fabI and heterodiploidy for fabI have been shown to confer resistance in S. aureus strains in a previous study. Here we further determined the fabI upstream sequence of a selection of these strains and the gene expression levels in strains with promoter region mutations.

Results

Mutations in the fabI promoter were found in 18% of triclosan resistant clinical isolates, regardless the previously identified molecular mechanism conferring resistance. Although not significant, a higher rate of promoter mutations were found in strains without previously described mechanisms of resistance. Some of the mutations identified in the clinical isolates were also detected in a series of laboratory mutants. Microarray analysis of selected laboratory mutants with fabI promoter region mutations, grown in the absence of triclosan, revealed increased fabI expression in three out of four tested strains. In two of these strains, only few genes other than fabI were upregulated. Consistently with these data, whole genome sequencing of in vitro selected mutants identified only few mutations except the upstream and coding regions of fabI, with the promoter mutation as the most probable cause of fabI overexpression. Importantly the gene expression profiling of clinical isolates containing similar mutations in the fabI promoter also showed, when compared to unrelated non-mutated isolates, a significant up-regulation of fabI.

Conclusions

In conclusion, we have demonstrated the presence of C34T, T109G, and A101C mutations in the fabI promoter region of strains with fabI up-regulation, both in clinical isolates and/or laboratory mutants. These data provide further observations linking mutations upstream fabI with up-regulated expression of the fabI gene.

Electronic supplementary material

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

Identification of Nocardia species using matrix-assisted laser desorption/ionization–time-of-flight mass spectrometry

Background

The MALDI (matrix-assisted laser desorption/ionization) Biotyper system for bacterial identification has already been utilized in clinical microbiology laboratories as a successful clinical application of protoemics. However, in cases of Nocardia, mass spectra suitable for MALDI Biotyper identification are often not obtained if such specimens are processed like general bacteria. This problem is related to the insufficiencies in bacterial spectrum databases that preclude accurate specimen identification. Here, we developed a bacterial processing method to improve mass spectra from specimens of the genus Nocardia. In addition, with the new processing method, we constructed a novel in-house bacterial database that combines a commercial database and mass spectra of Nocardia strains from the Department of Clinical Laboratory at Chiba University Hospital (DCLC) and the Medical Mycology Research Center at Chiba University (MMRC).

Results

The newly developed method (Nocardia Extraction Method at DCLC [NECLC]) based on ethanol-formic acid extraction (EFAE) improved mass spectra obtained from Nocardia specimens. The Nocardia in-house database at Chiba University Hospital (NDCUH) was then successfully validated. In brief, prior to introduction of the NECLC and NDCUH, 10 of 64 (15.6%) clinical isolates were identified at the species level and 16 isolates (25.0%) could only be identified at the genus level. In contrast, after the introduction, 58 isolates (90.6%) were identified at the species level and 6 isolates (9.4%) were identified at the genus level.

Conclusions

The results of this study suggest that MALDI-TOF (time-of-flight) Biotyper system can identify Nocardia accurately in a short time in combination with a simple processing method and an in-house database.

Electronic supplementary material

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

Genomic analysis of the emergence of 20th century epidemic dysentery

Background

Shigella dysenteriae type 1 (Sd1) causes recurrent epidemics of dysentery associated with high mortality in many regions of the world. Sd1 infects humans at very low infectious doses (10 CFU), and treatment is complicated by the rapid emergence of antibiotic resistant Sd1 strains. Sd1 is only detected in the context of human infections, and the circumstances under which epidemics emerge and regress remain unknown.

Results

Phylogenomic analyses of 56 isolates collected worldwide over the past 60 years indicate that the Sd1 clone responsible for the recent pandemics emerged at the turn of the 20^th century, and that the two world wars likely played a pivotal role for its dissemination. Several lineages remain ubiquitous and their phylogeny indicates several recent intercontinental transfers. Our comparative genomics analysis reveals that isolates responsible for separate outbreaks, though closely related to one another, have independently accumulated antibiotic resistance genes, suggesting that there is little or no selection to retain these genes in-between outbreaks. The genomes appear to be subjected to genetic drift that affects a number of functions currently used by diagnostic tools to identify Sd1, which could lead to the potential failure of such tools.

Conclusions

Taken together, the Sd1 population structure and pattern of evolution suggest a recent emergence and a possible human carrier state that could play an important role in the epidemic pattern of infections of this human-specific pathogen. This analysis highlights the important role of whole-genome sequencing in studying pathogens for which epidemiological or laboratory investigations are particularly challenging.

Electronic supplementary material

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

The isolation and characterization of two Stenotrophomonas maltophilia bacteriophages capable of cross-taxonomic order infectivity

Background

A rapid worldwide increase in the number of human infections caused by the extremely antibiotic resistant bacterium Stenotrophomonas maltophilia is prompting alarm. One potential treatment solution to the current antibiotic resistance dilemma is “phage therapy”, the clinical application of bacteriophages to selectively kill bacteria.

Results

Towards that end, phages DLP1 and DLP2 (vB_SmaS-DLP_1 and vB_SmaS-DLP_2, respectively) were isolated against S. maltophilia strain D1585. Host range analysis for each phage was conducted using 27 clinical S. maltophilia isolates and 11 Pseudomonas aeruginosa strains. Both phages exhibit unusually broad host ranges capable of infecting bacteria across taxonomic orders. Transmission electron microscopy of the phage DLP1 and DLP2 morphology reveals that they belong to the Siphoviridae family of bacteriophages. Restriction fragment length polymorphism analysis and complete genome sequencing and analysis indicates that phages DLP1 and DLP2 are closely related but different phages, sharing 96.7 % identity over 97.2 % of their genomes. These two phages are also related to P. aeruginosa phages vB_Pae-Kakheti_25 (PA25), PA73, and vB_PaeS_SCH_Ab26 (Ab26) and more distantly related to Burkholderia cepacia complex phage KL1, which together make up a taxonomic sub-family. Phages DLP1 and DLP2 exhibited significant differences in host ranges and growth kinetics.

Conclusions

The isolation and characterization of phages able to infect two completely different species of bacteria is an exciting discovery, as phages typically can only infect related bacterial species, and rarely infect bacteria across taxonomic families, let alone across taxonomic orders.

Electronic supplementary material

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

RiTE database: a resource database for genus-wide rice genomics and evolutionary biology

Background

Comparative evolutionary analysis of whole genomes requires not only accurate annotation of gene space, but also proper annotation of the repetitive fraction which is often the largest component of most if not all genomes larger than 50 kb in size.

Results

Here we present the Rice TE database (RiTE-db) - a genus-wide collection of transposable elements and repeated sequences across 11 diploid species of the genus Oryza and the closely-related out-group Leersia perrieri. The database consists of more than 170,000 entries divided into three main types: (i) a classified and curated set of publicly-available repeated sequences, (ii) a set of consensus assemblies of highly-repetitive sequences obtained from genome sequencing surveys of 12 species; and (iii) a set of full-length TEs, identified and extracted from 12 whole genome assemblies.

Conclusions

This is the first report of a repeat dataset that spans the majority of repeat variability within an entire genus, and one that includes complete elements as well as unassembled repeats. The database allows sequence browsing, downloading, and similarity searches. Because of the strategy adopted, the RiTE-db opens a new path to unprecedented direct comparative studies that span the entire nuclear repeat content of 15 million years of Oryza diversity.

Electronic supplementary material

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

Comparative genome sequencing reveals chemotype-specific gene clusters in the toxigenic black mold Stachybotrys

Background

The fungal genus Stachybotrys produces several diverse toxins that affect human health. Its strains comprise two mutually-exclusive toxin chemotypes, one producing satratoxins, which are a subclass of trichothecenes, and the other producing the less-toxic atranones. To determine the genetic basis for chemotype-specific differences in toxin production, the genomes of four Stachybotrys strains were sequenced and assembled de novo. Two of these strains produce atranones and two produce satratoxins.

Results

Comparative analysis of these four 35-Mbp genomes revealed several chemotype-specific gene clusters that are predicted to make secondary metabolites. The largest, which was named the core atranone cluster, encodes 14 proteins that may suffice to produce all observed atranone compounds via reactions that include an unusual Baeyer-Villiger oxidation. Satratoxins are suggested to be made by products of multiple gene clusters that encode 21 proteins in all, including polyketide synthases, acetyltransferases, and other enzymes expected to modify the trichothecene skeleton. One such satratoxin chemotype-specific cluster is adjacent to the core trichothecene cluster, which has diverged from those of other trichothecene producers to contain a unique polyketide synthase.

Conclusions

The results suggest that chemotype-specific gene clusters are likely the genetic basis for the mutually-exclusive toxin chemotypes of Stachybotrys. A unified biochemical model for Stachybotrys toxin production is presented. Overall, the four genomes described here will be useful for ongoing studies of this mold’s diverse toxicity mechanisms.

Electronic supplementary material

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

Comprehensive molecular,genomic and phenotypic analysis of a major clone of Enterococcus faecalis MLST ST40

Background

Enterococcus faecalis is a multifaceted microorganism known to act as a beneficial intestinal commensal bacterium. It is also a dreaded nosocomial pathogen causing life-threatening infections in hospitalised patients. Isolates of a distinct MLST type ST40 represent the most frequent strain type of this species, distributed worldwide and originating from various sources (animal, human, environmental) and different conditions (colonisation/infection). Since enterococci are known to be highly recombinogenic we determined to analyse the microevolution and niche adaptation of this highly distributed clonal type.

Results

We compared a set of 42 ST40 isolates by assessing key molecular determinants, performing whole genome sequencing (WGS) and a number of phenotypic assays including resistance profiling, formation of biofilm and utilisation of carbon sources. We generated the first circular closed reference genome of an E. faecalis isolate D32 of animal origin and compared it with the genomes of other reference strains. D32 was used as a template for detailed WGS comparisons of high-quality draft genomes of 14 ST40 isolates. Genomic and phylogenetic analyses suggest a high level of similarity regarding the core genome, also demonstrated by similar carbon utilisation patterns. Distribution of known and putative virulence-associated genes did not differentiate between ST40 strains from a commensal and clinical background or an animal or human source. Further analyses of mobile genetic elements (MGE) revealed genomic diversity owed to: (1) a modularly structured pathogenicity island; (2) a site-specifically integrated and previously unknown genomic island of 138 kb in two strains putatively involved in exopolysaccharide synthesis; and (3) isolate-specific plasmid and phage patterns. Moreover, we used different cell-biological and animal experiments to compare the isolate D32 with a closely related ST40 endocarditis isolate whose draft genome sequence was also generated. D32 generally showed a greater capacity of adherence to human cell lines and an increased pathogenic potential in various animal models in combination with an even faster growth in vivo (not in vitro).

Conclusion

Molecular, genomic and phenotypic analysis of representative isolates of a major clone of E. faecalis MLST ST40 revealed new insights into the microbiology of a commensal bacterium which can turn into a conditional pathogen.

Electronic supplementary material

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

Complete sequencing of Novosphingobium sp. PP1Y reveals a biotechnologically meaningful metabolic pattern

Background

Novosphingobium sp. strain PP1Y is a marine α-proteobacterium adapted to grow at the water/fuel oil interface. It exploits the aromatic fraction of fuel oils as a carbon and energy source. PP1Y is able to grow on a wide range of mono-, poly- and heterocyclic aromatic hydrocarbons. Here, we report the complete functional annotation of the whole Novosphingobium genome.

Results

PP1Y genome analysis and its comparison with other Sphingomonadal genomes has yielded novel insights into the molecular basis of PP1Y’s phenotypic traits, such as its peculiar ability to encapsulate and degrade the aromatic fraction of fuel oils. In particular, we have identified and dissected several highly specialized metabolic pathways involved in: (i) aromatic hydrocarbon degradation; (ii) resistance to toxic compounds; and (iii) the quorum sensing mechanism.

Conclusions

In summary, the unraveling of the entire PP1Y genome sequence has provided important insight into PP1Y metabolism and, most importantly, has opened new perspectives about the possibility of its manipulation for bioremediation purposes.

Electronic supplementary material

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

Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity

Background

So-called 936-type phages are among the most frequently isolated phages in dairy facilities utilising Lactococcus lactis starter cultures. Despite extensive efforts to control phage proliferation and decades of research, these phages continue to negatively impact cheese production in terms of the final product quality and consequently, monetary return.

Results

Whole genome sequencing and in silico analysis of three 936-type phage genomes identified several putative (orphan) methyltransferase (MTase)-encoding genes located within the packaging and replication regions of the genome. Utilising SMRT sequencing, methylome analysis was performed on all three phages, allowing the identification of adenine modifications consistent with N-6 methyladenine sequence methylation, which in some cases could be attributed to these phage-encoded MTases. Heterologous gene expression revealed that M.Phi145I/M.Phi93I and M.Phi93DAM, encoded by genes located within the packaging module, provide protection against the restriction enzymes HphI and DpnII, respectively, representing the first functional MTases identified in members of 936-type phages.

Conclusions

SMRT sequencing technology enabled the identification of the target motifs of MTases encoded by the genomes of three lytic 936-type phages and these MTases represent the first functional MTases identified in this species of phage. The presence of these MTase-encoding genes on 936-type phage genomes is assumed to represent an adaptive response to circumvent host encoded restriction-modification systems thereby increasing the fitness of the phages in a dynamic dairy environment.

Electronic supplementary material

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

Haemophilus influenzae: using comparative genomics to accurately identify a highly recombinogenic human pathogen

Background

Haemophilus influenzae is an opportunistic bacterial pathogen that exclusively colonises humans and is associated with both acute and chronic disease. Despite its clinical significance, accurate identification of H. influenzae is a non-trivial endeavour. H. haemolyticus can be misidentified as H. influenzae from clinical specimens using selective culturing methods, reflecting both the shared environmental niche and phenotypic similarities of these species. On the molecular level, frequent genetic exchange amongst Haemophilus spp. has confounded accurate identification of H. influenzae, leading to both false-positive and false-negative results with existing speciation assays.

Results

Whole-genome single-nucleotide polymorphism data from 246 closely related global Haemophilus isolates, including 107 Australian isolate genomes generated in this study, were used to construct a whole-genome phylogeny. Based on this phylogeny, H. influenzae could be differentiated from closely related species. Next, a H. influenzae-specific locus, fucP, was identified, and a novel TaqMan real-time PCR assay targeting fucP was designed. PCR specificity screening across a panel of clinically relevant species, coupled with in silico analysis of all species within the order Pasteurellales, demonstrated that the fucP assay was 100 % specific for H. influenzae; all other examined species failed to amplify.

Conclusions

This study is the first of its kind to use large-scale comparative genomic analysis of Haemophilus spp. to accurately delineate H. influenzae and to identify a species-specific molecular signature for this species. The fucP assay outperforms existing H. influenzae targets, most of which were identified prior to the next-generation genomics era and thus lack validation across a large number of Haemophilus spp. We recommend use of the fucP assay in clinical and research laboratories for the most accurate detection and diagnosis of H. influenzae infection and colonisation.

Electronic supplementary material

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