Serotype-conversion in Shigella flexneri: identification of a novel bacteriophage,Sf101, from a serotype 7a strain

Background

Shigella flexneri is the major cause of bacillary dysentery in the developing countries. The lipopolysaccharide (LPS) O-antigen of S. flexneri plays an important role in its pathogenesis and also divides S. flexneri into 19 serotypes. All the serotypes with an exception for serotype 6 share a common O-antigen backbone comprising of N-acetylglucosamine and three rhamnose residues. Different serotypes result from modification of the basic backbone conferred by phage-encoded glucosyltransferase and/or acetyltransferase genes, or plasmid-encoded phosphoethanolamine transferase. Recently, a new site for O-acetylation at positions 3 and 4 of Rha^III, in serotypes 1a, 1b, 2a, 5a and Y was shown to be mediated by the oacB gene. Additionally, this gene was shown to be carried by a transposon-like structure inserted upstream of the adrA region on the chromosome.

Results

In this study, a novel bacteriophage Sf101, encoding the oacB gene was isolated and characterised from a serotype 7a strain. The complete sequence of its 38,742 bp genome encoding 66 open reading frames (orfs) was determined. Comparative analysis revealed that phage Sf101 has a mosaic genome, and most of its proteins were >90% identical to the proteins from 12 previously characterised lambdoid phages. In addition, the organisation of Sf101 genes was found to be highly similar to bacteriophage Sf6. Analysis of the Sf101 OacB identified two amino acid substitutions in the protein; however, results obtained by NMR spectroscopy confirmed that Sf101-OacB was functional. Inspection of the chromosomal integration site of Sf101 phage revealed that this phage integrates in the sbcB locus, thus unveiling a new site for integration of serotype-converting phages of S. flexneri, and determining an alternative location of oacB gene in the chromosome. Furthermore, this study identified oacB gene in several serotype 7a isolates from various regions providing evidence of O-acetyl modification in serotype 7a.

Conclusions

This is the first report on the isolation of bacteriophage Sf101 which contains the S. flexneri O-antigen modification gene oacB. Sf101 has a highly mosaic genome and was found to integrate in the sbcB locus. These findings contribute an advance in our current knowledge of serotype converting phages of S. flexneri.

Electronic supplementary material

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

Induced terpene accumulation in Norway spruce inhibits bark beetle colonization in a dose-dependent manner

Background

Tree-killing bark beetles (Coleoptera, Scolytinae) are among the most economically and ecologically important forest pests in the northern hemisphere. Induction of terpenoid-based oleoresin has long been considered important in conifer defense against bark beetles, but it has been difficult to demonstrate a direct correlation between terpene levels and resistance to bark beetle colonization.

Methods

To test for inhibitory effects of induced terpenes on colonization by the spruce bark beetle Ips typographus (L.) we inoculated 20 mature Norway spruce Picea abies (L.) Karsten trees with a virulent fungus associated with the beetle, Ceratocystis polonica (Siem.) C. Moreau, and investigated induced terpene levels and beetle colonization in the bark.

Results

Fungal inoculation induced very strong and highly variable terpene accumulation 35 days after inoculation. Trees with high induced terpene levels (n = 7) had only 4.9% as many beetle attacks (5.1 vs. 103.5 attacks m⁻²) and 2.6% as much gallery length (0.029 m m⁻² vs. 1.11 m m⁻²) as trees with low terpene levels (n = 6). There was a highly significant rank correlation between terpene levels at day 35 and beetle colonization in individual trees. The relationship between induced terpene levels and beetle colonization was not linear but thresholded: above a low threshold concentration of ∼100 mg terpene g⁻¹ dry phloem trees suffered only moderate beetle colonization, and above a high threshold of ∼200 mg terpene g⁻¹ dry phloem trees were virtually unattacked.

Conclusion/Significance

This is the first study demonstrating a dose-dependent relationship between induced terpenes and tree resistance to bark beetle colonization under field conditions, indicating that terpene induction may be instrumental in tree resistance. This knowledge could be useful for developing management strategies that decrease the impact of tree-killing bark beetles.     

Polymorphisms of the Tissue Inhibitor of Metalloproteinase 3 Gene Are Associated with Resistance to High-Altitude Pulmonary Edema (HAPE) in a Japanese Population: A Case Control Study Using Polymorphic Microsatellite Markers

Introduction

High-altitude pulmonary edema (HAPE) is a hypoxia-induced, life-threatening, high permeability type of edema attributable to pulmonary capillary stress failure. Genome-wide association analysis is necessary to better understand how genetics influence the outcome of HAPE.

Materials and Methods

DNA samples were collected from 53 subjects susceptible to HAPE (HAPE-s) and 67 elite Alpinists resistant to HAPE (HAPE-r). The genome scan was carried out using 400 polymorphic microsatellite markers throughout the whole genome in all subjects. In addition, six single nucleotide polymorphisms (SNPs) of the gene encoding the tissue inhibitor of metalloproteinase 3 (TIMP3) were genotyped by Taqman® SNP Genotyping Assays.

Results

The results were analyzed using case-control comparisons. Whole genome scanning revealed that allele frequencies in nine markers were statistically different between HAPE-s and HAPE-r subjects. The SNP genotyping of the TIMP3 gene revealed that the derived allele C of rs130293 was associated with resistance to HAPE [odds ratio (OR) = 0.21, P = 0.0012) and recessive inheritance of the phenotype of HAPE-s (P = 0.0012). A haplotype CAC carrying allele C of rs130293 was associated with resistance to HAPE.

Discussion

This genome-wide association study revealed several novel candidate genes associated with susceptibility or resistance to HAPE in a Japanese population. Among those, the minor allele C of rs130293 (C/T) in the TIMP3 gene was linked to resistance to HAPE; while, the ancestral allele T was associated with susceptibility to HAPE.     

Identification and distribution of the NBS-LRR gene family in the Cassava genome

Background

Plant resistance genes (R genes) exist in large families and usually contain both a nucleotide-binding site domain and a leucine-rich repeat domain, denoted NBS-LRR. The genome sequence of cassava (Manihot esculenta) is a valuable resource for analysing the genomic organization of resistance genes in this crop.

Results

With searches for Pfam domains and manual curation of the cassava gene annotations, we identified 228 NBS-LRR type genes and 99 partial NBS genes. These represent almost 1% of the total predicted genes and show high sequence similarity to proteins from other plant species. Furthermore, 34 contained an N-terminal toll/interleukin (TIR)-like domain, and 128 contained an N-terminal coiled-coil (CC) domain. 63% of the 327 R genes occurred in 39 clusters on the chromosomes. These clusters are mostly homogeneous, containing NBS-LRRs derived from a recent common ancestor.

Conclusions

This study provides insight into the evolution of NBS-LRR genes in the cassava genome; the phylogenetic and mapping information may aid efforts to further characterize the function of these predicted R genes.

Electronic supplementary material

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

Multiple mutations in heterogeneous miltefosine-resistant Leishmania major population as determined by whole genome sequencing

Background

Miltefosine (MF) is the first oral compound used in the chemotherapy against leishmaniasis. Since the mechanism of action of this drug and the targets of MF in Leishmania are unclear, we generated in a step-by-step manner Leishmania major promastigote mutants highly resistant to MF. Two of the mutants were submitted to a short-read whole genome sequencing for identifying potential genes associated with MF resistance.

Methods/Principal Findings

Analysis of the genome assemblies revealed several independent point mutations in a P-type ATPase involved in phospholipid translocation. Mutations in two other proteins—pyridoxal kinase and α-adaptin like protein—were also observed in independent mutants. The role of these proteins in the MF resistance was evaluated by gene transfection and gene disruption and both the P-type ATPase and pyridoxal kinase were implicated in MF susceptibility. The study also highlighted that resistance can be highly heterogeneous at the population level with individual clones derived from this population differing both in terms of genotypes but also susceptibility phenotypes.

Conclusions/Significance

Whole genome sequencing was used to pinpoint known and new resistance markers associated with MF resistance in the protozoan parasite Leishmania. The study also demonstrated the polyclonal nature of a resistant population with individual cells with varying susceptibilities and genotypes.     

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.     

Biogeography and Virulence of Staphylococcus aureus

Background

Staphylococcus aureus is commonly carried asymptomatically in the human anterior nares and occasionally enters the bloodstream to cause invasive disease. Much of the global diversity of S. aureus remains uncharacterised, and is not clear how disease propensity varies between strains, and between host populations.

Methodology

We compared 147 isolates recovered from five kindergartens in Chengdu, China, with 51 isolates contemporaneously recovered from cases of pediatric infection from the main hospital serving this community. The samples were characterised by MLST, the presence/absence of PVL, and antibiotic resistance profiling.

Principal Findings

Genotype frequencies within individual kindergartens differ, but the sample recovered from cases of disease shows a general enrichment of certain MLST genotypes and PVL positive isolates. Genotypes under-represented in the disease sample tend to correspond to a single sequence cluster, and this cluster is more common in China than in other parts of the world.

Conclusions/Significance

Virulence propensity likely reflects a synergy between variation in the core genome (MLST) and accessory genome (PVL). By combining evidence form biogeography and virulence we demonstrate the existence of a “native” clade in West China which has lowered virulence, possibility due to acquired host immunity.     

Unusual genome complexity in Lactobacillus salivarius JCM1046

Background

Lactobacillus salivarius strains are increasingly being exploited for their probiotic properties in humans and animals. Dissemination of antibiotic resistance genes among species with food or probiotic-association is undesirable and is often mediated by plasmids or integrative and conjugative elements. L. salivarius strains typically have multireplicon genomes including circular megaplasmids that encode strain-specific traits for intestinal survival and probiotic activity. Linear plasmids are less common in lactobacilli and show a very limited distribution in L. salivarius. Here we present experimental evidence that supports an unusually complex multireplicon genome structure in the porcine isolate L. salivarius JCM1046.

Results

JCM1046 harbours a 1.83 Mb chromosome, and four plasmids which constitute 20% of the genome. In addition to the known 219 kb repA-type megaplasmid pMP1046A, we identified and experimentally validated the topology of three additional replicons, the circular pMP1046B (129 kb), a linear plasmid pLMP1046 (101 kb) and pCTN1046 (33 kb) harbouring a conjugative transposon. pMP1046B harbours both plasmid-associated replication genes and paralogues of chromosomally encoded housekeeping and information-processing related genes, thus qualifying it as a putative chromid. pLMP1046 shares limited sequence homology or gene synteny with other L. salivarius plasmids, and its putative replication-associated protein is homologous to the RepA/E proteins found in the large circular megaplasmids of L. salivarius. Plasmid pCTN1046 harbours a single copy of an integrated conjugative transposon (Tn6224) which appears to be functionally intact and includes the tetracycline resistance gene tetM.

Conclusion

Experimental validation of sequence assemblies and plasmid topology resolved the complex genome architecture of L. salivarius JCM1046. A high-coverage draft genome sequence would not have elucidated the genome complexity in this strain. Given the expanding use of L. salivarius as a probiotic, it is important to determine the genotypic and phenotypic organization of L. salivarius strains. The identification of Tn6224-like elements in this species has implications for strain selection for probiotic applications.

Electronic supplementary material

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

Whole genome sequencing of penicillin-resistant Streptococcus pneumoniae reveals mutations in penicillin-binding proteins and in a putative iron permease

Background

Penicillin resistance in Streptococcus pneumoniae is mediated by a mosaic of genes encoding altered penicillin-binding proteins (PBPs). Nonetheless, S. pneumoniae has also developed non-PBP mechanisms implicated in penicillin resistance. In this study, whole genome sequencing of resistant organisms was used to discover mutations implicated in resistance to penicillin.

Results

We sequenced two S. pneumoniae isolates selected for resistance to penicillin in vitro. The analysis of the genome assemblies revealed that six genes were mutated in both mutants. These included three pbp genes, and three non-pbp genes, including a putative iron permease, spr1178. The nonsense mutation in spr1178 always occurred in the first step of the selection process. Although the mutants had increased resistance to penicillin, the introduction of altered versions of PBPs into a penicillin-susceptible strain by sequential transformation led to strains with a minimal increase in resistance, thus implicating other genes in resistance. The introduction by transformation of the non-PBP recurrent mutations did not increase penicillin resistance, but the introduction of the nonsense mutation in the putative iron permease spr1178 led to a reduced accumulation of reactive oxygen species following exposure to penicillin and to other bactericidal antibiotics as well.

Conclusions

This study indicates that the selection of resistance to penicillin in S. pneumoniae involves the acquisition of mutations conferring tolerance to the antibiotic-induced accumulation of oxidants, which translates into an increased survival that putatively enables the selection of major resistance determinants such as mutations in PBPs.     

High correlation between genotypes and phenotypes of environmental bacteria Comamonas testosteroni strains

Background

Members of Comamonas testosteroni are environmental microorganisms that are usually found in polluted environment samples. They utilize steroids and aromatic compounds but rarely sugars, and show resistance to multiple heavy metals and multiple drugs. However, comprehensive genomic analysis among the C. testosteroni strains is lacked.

Results

To understand the genome bases of the features of C. testosteroni, we sequenced 10 strains of this species and analyzed them together with other related published genome sequences. The results revealed that: 1) the strains of C. testosteroni have genome sizes ranging from 5.1 to 6.0 Mb and G + C contents ranging from 61.1% to 61.8%. The pan-genome contained 10,165 gene families and the core genome contained 3,599 gene families. Heap’s law analysis indicated that the pan-genome of C. testosteroni may be open (α = 0.639); 2) by analyzing 31 phenotypes of 11 available C. testosteroni strains, 99.4% of the genotypes (putative genes) were found to be correlated to the phenotypes, indicating a high correlation between phenotypes and genotypes; 3) gene clusters for nitrate reduction, steroids degradation and metal and multi-drug resistance were found and were highly conserved among all the genomes of this species; 4) the genome similarity of C. testosteroni may be related to the geographical distances.

Conclusions

This work provided an overview on the genomes of C. testosteroni and new genome resources that would accelerate the further investigations of this species. Importantly, this work focused on the analysis of potential genetic determinants for the typical characters and found high correlation between the phenotypes and their corresponding genotypes.

Electronic supplementary material

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

Rapid loss of genetic variation in a founding population of Primula elatior (Primulaceae) after colonization

Background and Aims

Land-use changes and associated extinction/colonization dynamics can have a large impact on population genetic diversity of plant species. The aim of this study was to investigate genetic diversity in a founding population of the self-incompatible forest herb Primula elatior and to elucidate the processes that affect genetic diversity shortly after colonization.

Methods

AFLP markers were used to analyse genetic diversity across three age classes and spatial genetic structure within a founding population of P. elatior in a recently established stand in central Belgium. Parentage analyses were used to assess the amount of gene flow from outside the population and to investigate the contribution of mother plants to future generations.

Results

The genetic diversity of second and third generation plants was significantly reduced compared with that of first generation plants. Significant spatial genetic structure was observed. Parentage analyses showed that <20 % of the youngest individuals originated from parents outside the study population and that >50 % of first and second generation plants did not contribute to seedling recruitment.

Conclusions

These results suggest that a small effective population size and genetic drift can lead to rapid decline of genetic diversity of offspring in founding populations shortly after colonization. This multigenerational study also highlights that considerable amounts of gene flow seem to be required to counterbalance genetic drift and to sustain high levels of genetic diversity after colonization in recently established stands.Key words: AFLP, colonization, forest regeneration, genetic diversity, genetic drift, parentage analysis, spatial genetic structure     

Homoeologous duplicated regions are involved in quantitative resistance of Brassica napus to stem canker

Background

Several major crop species are current or ancient polyploids. To better describe the genetic factors controlling traits of agronomic interest (QTL), it is necessary to understand the structural and functional organisation of these QTL regions in relation to genome duplication. We investigated quantitative resistance to the fungal disease stem canker in Brassica napus, a highly duplicated amphidiploid species, to assess the proportion of resistance QTL located at duplicated positions.

Results

Genome-wide association analysis on a panel of 116 oilseed rape varieties genotyped with 3228 SNP indicated that 321 markers, corresponding to 64 genomic regions, are associated with resistance to stem canker. These genomic regions are relatively equally distributed on the A (53%) and C (47%) genomes of B. napus. Overall, 44% of these regions (28/64) are duplicated homoeologous regions. They are located in duplications of six (E, J, R, T, U and W) of the 24 ancestral blocks that constitute the B. napus genome. Overall, these six ancestral blocks have 34 duplicated copies in the B.napus genome. Almost all of the duplicated copies (82% of the 34 regions) harboured resistance associated markers for stem canker resistance, which suggests structural and functional conservation of genetic factors involved in this trait in B. napus.

Conclusions

Our study provides information on the involvement of duplicated loci in the control of stem canker resistance in B. napus. Further investigation of the similarity/divergence in sequence and gene content of these duplicated regions will provide insight into the conservation and allelic diversity of the underlying genes.

Electronic supplementary material

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

Whole genome characterization of the mechanisms of daptomycin resistance in clinical and laboratory derived isolates of Staphylococcus aureus

Background

Daptomycin remains one of our last-line anti-staphylococcal agents. This study aims to characterize the genetic evolution to daptomycin resistance in S. aureus.

Methods

Whole genome sequencing was performed on a unique collection of isogenic, clinical (21 strains) and laboratory (12 strains) derived strains that had been exposed to daptomycin and developed daptomycin-nonsusceptibility. Electron microscopy (EM) and lipid membrane studies were performed on selected isolates.

Results

On average, six coding region mutations were observed across the genome in the clinical daptomycin exposed strains, whereas only two mutations on average were seen in the laboratory exposed pairs. All daptomycin-nonsusceptible strains had a mutation in a phospholipid biosynthesis gene. This included mutations in the previously described mprF gene, but also in other phospholipid biosynthesis genes, including cardiolipin synthase (cls2) and CDP-diacylglycerol-glycerol-3-phosphate 3-phosphatidyltransferase (pgsA). EM and lipid membrane composition analyses on two clinical pairs showed that the daptomycin-nonsusceptible strains had a thicker cell wall and an increase in membrane lysyl-phosphatidylglycerol.

Conclusion

Point mutations in genes coding for membrane phospholipids are associated with the development of reduced susceptibility to daptomycin in S. aureus. Mutations in cls2 and pgsA appear to be new genetic mechanisms affecting daptomycin susceptibility in S. aureus.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

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

Single-molecule sequencing reveals the molecular basis of multidrug-resistance in ST772 methicillin-resistant Staphylococcus aureus

Background

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of hospital-associated infection, but there is growing awareness of the emergence of multidrug-resistant lineages in community settings around the world. One such lineage is ST772-MRSA-V, which has disseminated globally and is increasingly prevalent in India. Here, we present the complete genome sequence of DAR4145, a strain of the ST772-MRSA-V lineage from India, and investigate its genomic characteristics in regards to antibiotic resistance and virulence factors.

Results

Sequencing using single-molecule real-time technology resulted in the assembly of a single continuous chromosomal sequence, which was error-corrected, annotated and compared to nine draft genome assemblies of ST772-MRSA-V from Australia, Malaysia and India. We discovered numerous and redundant resistance genes associated with mobile genetic elements (MGEs) and known core genome mutations that explain the highly antibiotic resistant phenotype of DAR4145. Staphylococcal toxins and superantigens, including the leukotoxin Panton-Valentinin Leukocidin, were predominantly associated with genomic islands and the phage φ-IND772PVL. Some of these mobile resistance and virulence factors were variably present in other strains of the ST772-MRSA-V lineage.

Conclusions

The genomic characteristics presented here emphasize the contribution of MGEs to the emergence of multidrug-resistant and highly virulent strains of community-associated MRSA. Antibiotic resistance was further augmented by chromosomal mutations and redundancy of resistance genes. The complete genome of DAR4145 provides a valuable resource for future investigations into the global dissemination and phylogeography of ST772-MRSA-V.

Electronic supplementary material

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

Mono-allelic retrotransposon insertion addresses epigenetic transcriptional repression in human genome

Background

Retrotransposons have been extensively studied in plants and animals and have been shown to have an impact on human genome dynamics and evolution. Their ability to move within genomes gives retrotransposons to affect genome instability.

Methods

we examined the polymorphic inserted AluYa5, evolutionary young Alu, in the progesterone receptor gene to determine the effects of Alu insertion on molecular environment. We used mono-allelic inserted cell lines which carry both Alu-present and Alu-absent alleles. To determine the epigenetic change and gene expression, we performed restriction enzyme digestion, Pyrosequencing, and Chromatin Immunoprecipitation.

Results

We observed that the polymorphic insertion of evolutionally young Alu causes increasing levels of DNA methylation in the surrounding genomic area and generates inactive histone tail modifications. Consequently the Alu insertion deleteriously inactivates the neighboring gene expression.

Conclusion

The mono-allelic Alu insertion cell line clearly showed that polymorphic inserted repetitive elements cause the inactivation of neighboring gene expression, bringing aberrant epigenetic changes.     

Gene Amplification and Point Mutations in Pyrimidine Metabolic Genes in 5-Fluorouracil Resistant Leishmania infantum

Background

The human protozoan parasites Leishmania are prototrophic for pyrimidines with the ability of both de novo biosynthesis and uptake of pyrimidines.

Methodology/Principal Findings

Five independent L. infantum mutants were selected for resistance to the pyrimidine analogue 5-fluorouracil (5-FU) in the hope to better understand the metabolism of pyrimidine in Leishmania. Analysis of the 5-FU mutants by comparative genomic hybridization and whole genome sequencing revealed in selected mutants the amplification of DHFR-TS and a deletion of part of chromosome 10. Point mutations in uracil phosphorybosyl transferase (UPRT), thymidine kinase (TK) and uridine phosphorylase (UP) were also observed in three individual resistant mutants. Transfection experiments confirmed that these point mutations were responsible for 5-FU resistance. Transport studies revealed that one resistant mutant was defective for uracil and 5-FU import.

Conclusion/Significance

This study provided further insights in pyrimidine metabolism in Leishmania and confirmed that multiple mutations can co-exist and lead to resistance in Leishmania.     

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.