Comparative and bioinformatics analyses of pathogenic bacterial secretomes identified by mass spectrometry in <Emphasis Type="Italic">Burkholderia</Emphasis> species

Secreted proteins (secretomes) play crucial roles during bacterial pathogenesis in both plant and human hosts. The identification and characterization of secretomes in the two plant pathogens Burkholderia glumae BGR1 and B. gladioli BSR3, which cause diseases in rice such as seedling blight, panicle blight, and grain rot, are important steps to not only understand the disease-causing mechanisms but also find remedies for the diseases. Here, we identified two datasets of secretomes in B. glumae BGR1 and B. gladioli BSR3, which consist of 118 and 111 proteins, respectively, using mass spectrometry approach and literature curation. Next, we characterized the functional properties, potential secretion pathways and sequence information properties of secretomes of two plant pathogens in a comparative analysis by various computational approaches. The ratio of potential non-classically secreted proteins (NCSPs) to classically secreted proteins (CSPs) in B. glumae BGR1 was greater than that in B. gladioli BSR3. For CSPs, the putative hydrophobic regions (PHRs) which are essential for secretion process of CSPs were screened in detail at their N-terminal sequences using hidden Markov model (HMM)–based method. Total 31 pairs of homologous proteins in two bacterial secretomes were indicated based on the global alignment (identity ≥ 70%). Our results may facilitate the understanding of the species-specific features of secretomes in two plant pathogenic Burkholderia species.     

An invasive Mimosa in India does not adopt the symbionts of its native relatives

Background and Aims

The large monophyletic genus Mimosa comprises approx. 500 species, most of which are native to the New World, with Central Brazil being the main centre of radiation. All Brazilian Mimosa spp. so far examined are nodulated by rhizobia in the betaproteobacterial genus Burkholderia. Approximately 10 Mya, transoceanic dispersal resulted in the Indian subcontinent hosting up to six endemic Mimosa spp. The nodulation ability and rhizobial symbionts of two of these, M. hamata and M. himalayana, both from north-west India, are here examined, and compared with those of M. pudica, an invasive species.

Methods

Nodules were collected from several locations, and examined by light and electron microscopy. Rhizobia isolated from them were characterized in terms of their abilities to nodulate the three Mimosa hosts. The molecular phylogenetic relationships of the rhizobia were determined by analysis of 16S rRNA, nifH and nodA gene sequences.

Key Results

Both native Indian Mimosa spp. nodulated effectively in their respective rhizosphere soils. Based on 16S rRNA, nifH and nodA sequences, their symbionts were identified as belonging to the alphaproteobacterial genus Ensifer, and were closest to the ‘Old World’ Ensifer saheli, E. kostiensis and E. arboris. In contrast, the invasive M. pudica was predominantly nodulated by Betaproteobacteria in the genera Cupriavidus and Burkholderia. All rhizobial strains tested effectively nodulated their original hosts, but the symbionts of the native species could not nodulate M. pudica.

Conclusions

The native Mimosa spp. in India are not nodulated by the Burkholderia symbionts of their South American relatives, but by a unique group of alpha-rhizobial microsymbionts that are closely related to the ‘local’ Old World Ensifer symbionts of other mimosoid legumes in north-west India. They appear not to share symbionts with the invasive M. pudica, symbionts of which are mostly beta-rhizobial.     

The complete genome of Burkholderia phenoliruptrix strain BR3459a,a symbiont of Mimosa flocculosa: highlighting the coexistence of symbiotic and pathogenic genes

Background

Burkholderia species play an important ecological role related to xenobiosis, the promotion of plant growth, the biocontrol of agricultural diseases, and symbiotic and non-symbiotic biological nitrogen fixation. Here, we highlight our study as providing the first complete genome of a symbiotic strain of B. phenoliruptrix, BR3459a (=CLA1), which was originally isolated in Brazil from nodules of Mimosa flocculosa and is effective in fixing nitrogen in association with this leguminous species.

Results

Genomic comparisons with other pathogenic and non-pathogenic Burkholderia strains grouped B. phenoliruptrix BR3459a with plant-associated beneficial and environmental species, although it shares a high percentage of its gene repertoire with species of the B. cepacia complex (Bcc) and "pseudomallei" group. The genomic analyses showed that the bce genes involved in exopolysaccharide production are clustered together in the same genomic region, constituting part of the Group III cluster of non-pathogenic bacteria. Regarding environmental stresses, we highlight genes that might be relevant in responses to osmotic, heat, cold and general stresses. Furthermore, a number of particularly interesting genes involved in the machinery of the T1SS, T2SS, T3SS, T4ASS and T6SS secretion systems were identified. The xenobiotic properties of strain BR3459a were also investigated, and some enzymes involved in the degradation of styrene, nitrotoluene, dioxin, chlorocyclohexane, chlorobenzene and caprolactam were identified. The genomic analyses also revealed a large number of antibiotic-related genes, the most important of which were correlated with streptomycin and novobiocin. The symbiotic plasmid showed high sequence identity with the symbiotic plasmid of B. phymatum. Additionally, comparative analysis of 545 housekeeping genes among pathogenic and non-pathogenic Burkholderia species strongly supports the definition of a new genus for the second branch, which would include BR3459a.

Conclusions

The analyses of B. phenoliruptrix BR3459a showed key property of fixing nitrogen that together with genes for high tolerance to environmental stresses might explain a successful strategy of symbiosis in the tropics. The strain also harbours interesting sets of genes with biotechnological potential. However, the resemblance of certain genes to those of pathogenic Burkholderia raise concerns about large-scale applications in agriculture or for bioremediation.

Electronic supplementary material

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

Genome analysis of Flaviramulus ichthyoenteri Th78T in the family Flavobacteriaceae: insights into its quorum quenching property and potential roles in fish intestine

Background

Intestinal microbes play significant roles in fish and can be possibly used as probiotics in aquaculture. In our previous study, Flaviramulus ichthyoenteri Th78^T, a novel species in the family Flavobacteriaceae, was isolated from fish intestine and showed strong quorum quenching (QQ) ability. To identify the QQ enzymes in Th78^T and explore the potential roles of Th78^T in fish intestine, we sequenced the genome of Th78^T and performed extensive genomic analysis.

Results

An N-acyl homoserine lactonase FiaL belonging to the metallo-β-lactamase superfamily was identified and the QQ activity of heterologously expressed FiaL was confirmed in vitro. FiaL has relatively little similarity to the known lactonases (25.2 ~ 27.9% identity in amino acid sequence). Various digestive enzymes including alginate lyases and lipases can be produced by Th78^T, and enzymes essential for production of B vitamins such as biotin, riboflavin and folate are predicted. Genes encoding sialic acid lyases, sialidases, sulfatases and fucosidases, which contribute to utilization of mucus, are present in the genome. In addition, genes related to response to different stresses and gliding motility were also identified. Comparative genome analysis shows that Th78^T has more specific genes involved in carbohydrate transport and metabolism compared to other two isolates in Flavobacteriaceae, both isolated from sediments.

Conclusions

The genome of Th78^T exhibits evident advantages for this bacterium to survive in the fish intestine, including production of QQ enzyme, utilization of various nutrients available in the intestine as well as the ability to produce digestive enzymes and vitamins, which also provides an application prospect of Th78^T to be used as a probiotic in aquaculture.

Electronic supplementary material

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

Genome based analysis of type-I polyketide synthase and nonribosomal peptide synthetase gene clusters in seven strains of five representative Nocardia species

Background

Actinobacteria of the genus Nocardia usually live in soil or water and play saprophytic roles, but they also opportunistically infect the respiratory system, skin, and other organs of humans and animals. Primarily because of the clinical importance of the strains, some Nocardia genomes have been sequenced, and genome sequences have accumulated. Genome sizes of Nocardia strains are similar to those of Streptomyces strains, the producers of most antibiotics. In the present work, we compared secondary metabolite biosynthesis gene clusters of type-I polyketide synthase (PKS-I) and nonribosomal peptide synthetase (NRPS) among genomes of representative Nocardia species/strains based on domain organization and amino acid sequence homology.

Results

Draft genome sequences of Nocardia asteroides NBRC 15531^T, Nocardia otitidiscaviarum IFM 11049, Nocardia brasiliensis NBRC 14402^T, and N. brasiliensis IFM 10847 were read and compared with published complete genome sequences of Nocardia farcinica IFM 10152, Nocardia cyriacigeorgica GUH-2, and N. brasiliensis HUJEG-1. Genome sizes are as follows: N. farcinica, 6.0 Mb; N. cyriacigeorgica, 6.2 Mb; N. asteroides, 7.0 Mb; N. otitidiscaviarum, 7.8 Mb; and N. brasiliensis, 8.9 - 9.4 Mb. Predicted numbers of PKS-I, NRPS, and PKS-I/NRPS hybrid clusters ranged between 4–11, 7–13, and 1–6, respectively, depending on strains, and tended to increase with increasing genome size. Domain and module structures of representative or unique clusters are discussed in the text.

Conclusion

We conclude the following: 1) genomes of Nocardia strains carry as many PKS-I and NRPS gene clusters as those of Streptomyces strains, 2) the number of PKS-I and NRPS gene clusters in Nocardia strains varies substantially depending on species, and N. brasiliensis strains carry the largest numbers of clusters among the species studied, 3) the seven Nocardia strains studied in the present work have seven common PKS-I and/or NRPS clusters, some of whose products are yet to be studied, and 4) different N. brasiliensis strains have some different gene clusters of PKS-I/NRPS, although the rest of the clusters are common within the N. brasiliensis strains. Genome sequencing suggested that Nocardia strains are highly promising resources in the search of novel secondary metabolites.

Electronic supplementary material

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

Comparative cytogenetic analysis of the genomes of the model grass Brachypodium distachyon and its close relatives

Background and Aims

Brachypodium is a small genus of temperate grasses that comprises 12–15 species. Brachypodium distachyon is now well established as a model species for temperate cereals and forage grasses. In contrast to B. distachyon, other members of the genus have been poorly investigated at the chromosome level or not at all.

Methods

Twenty accessions comprising six species and two subspecies of Brachypodium were analysed cytogenetically. Measurements of nuclear genome size were made by flow cytometry. Chromosomal localization of 18–5·8–25S rDNA and 5S rDNA loci was performed by dual-colour fluorescence in situ hybridization (FISH) on enzymatically digested root-tip meristematic cells. For comparative phylogenetic analyses genomic in situ hybridization (GISH) applied to somatic chromosome preparations was used.

Key Results

All Brachypodium species examined have rather small genomes and chromosomes. Their chromosome numbers and genome sizes vary from 2n = 10 and 0·631 pg/2C in B. distachyon to 2n = 38 and 2·57 pg/2C in B. retusum, respectively. Genotypes with 18 and 28 chromosomes were found among B. pinnatum accessions. GISH analysis revealed that B. pinnatum with 28 chromosomes is most likely an interspecific hybrid between B. distachyon (2n = 10) and B. pinnatum (2n = 18). Two other species, B. phoenicoides and B. retusum, are also allopolyploids and B. distachyon or a close relative seems to be one of their putative ancestral species. In chromosomes of all species examined the 45S rDNA loci are distally distributed whereas loci for 5S rDNA are pericentromeric.

Conclusions

The increasing significance of B. distachyon as a model grass emphasizes the need to understand the evolutionary relationships in the genus Brachypodium and to ensure consistency in the biological nomenclature of its species. Modern molecular cytogenetic techniques such as FISH and GISH are suitable for comparative phylogenetic analyses and may provide informative chromosome- and/or genome-specific landmarks.     

Downregulation of net phosphorus-uptake capacity is inversely related to leaf phosphorus-resorption proficiency in four species from a phosphorus-impoverished environment

Background and Aims

Previous research has suggested a trade-off between the capacity of plants to downregulate their phosphorus (P) uptake capacity and their efficiency of P resorption from senescent leaves in species from P-impoverished environments.

Methods

To investigate this further, four Australian native species (Banksia attenuata, B. menziesii, Acacia truncata and A. xanthina) were grown in a greenhouse in nutrient solutions at a range of P concentrations [P]. Acacia plants received between 0 and 500 µm P; Banksia plants received between 0 and 10 µm P, to avoid major P-toxicity symptoms in these highly P-sensitive species.

Key Results

For both Acacia species, the net P-uptake rates measured at 10 µm P decreased steadily with increasing P supply during growth. In contrast, in B. attenuata, the net rate of P uptake from a solution with 10 µm P increased linearly with increasing P supply during growth. The P-uptake rate of B. menziesii showed no significant response to P supply in the growing medium. Leaf [P] of the four species supported this finding, with A. truncata and A. xanthina showing an increase up to a saturation value of 19 and 21 mg P g⁻¹ leaf dry mass, respectively (at 500 µm P), whereas B. attenuata and B. menziesii both exhibited a linear increase in leaf [P], reaching 10 and 13 mg P g⁻¹ leaf dry mass, respectively, without approaching a saturation point. The Banksia plants grown at 10 µm P showed mild symptoms of P toxicity, i.e. yellow spots on some leaves and drying and curling of the tips of the leaves. Leaf P-resorption efficiency was 69 % (B. attenuata), 73 % (B. menziesii), 34 % (A. truncata) and 36 % (A. xanthina). The P-resorption proficiency values were 0·08 mg P g⁻¹ leaf dry mass (B. attenuata and B. menziesii), 0·32 mg P g⁻¹ leaf dry mass (A. truncata) and 0·36 mg P g⁻¹ leaf dry mass (A. xanthina). Combining the present results with additional information on P-remobilization efficiency and the capacity to downregulate P-uptake capacity for two other Australian woody species, we found a strong negative correlation between these traits.

Conclusions

It is concluded that species that are adapted to extremely P-impoverished soils, such as many south-western Australian Proteaceae species, have developed extremely high P-resorption efficiencies, but lost their capacity to downregulate their P-uptake mechanisms. The results support the hypothesis that the ability to resorb P from senescing leaves is inversely related to the capacity to downregulate net P uptake, possibly because constitutive synthesis of P transporters is a prerequisite for proficient P remobilization from senescing tissues.     

Involvement of a Quorum-Sensing-Regulated Lipase Secreted by a Clinical Isolate of Burkholderia glumae in Severe Disease Symptoms in Rice

Burkholderia glumae is an emerging rice pathogen in several areas around the world. Closely related Burkholderia species are important opportunistic human pathogens for specific groups of patients, such as patients with cystic fibrosis and patients with chronic granulomatous disease. Here we report that the first clinical isolate of B. glumae, strain AU6208, has retained its capability to be very pathogenic to rice. As previously reported for rice isolate B. glumae BGR1 (and also for the clinical isolate AU6208), TofI or TofR acyl homoserine lactone (AHL) quorum sensing played a pivotal role in rice virulence. We report that AHL quorum sensing in B. glumae AU6208 regulates secreted LipA lipase and toxoflavin, the phytotoxin produced by B. glumae. B. glumae AU6208 lipA mutants were no longer pathogenic to rice, indicating that the lipase is an important virulence factor. It was also established that type strain B. glumae ATCC 33617 did not produce toxoflavin and lipase and was nonpathogenic to rice. It was determined that in strain ATCC 33617 the LuxR family quorum-sensing sensor/regulator TofR was inactive. Introducing the tofR gene of B. glumae AU6208 in strain ATCC 33617 restored its ability to produce toxoflavin and the LipA lipase. This study extends the role of AHL quorum sensing in rice pathogenicity through the regulation of a lipase which was demonstrated to be a virulence factor. It is the first report of a clinical B. glumae isolate retaining strong rice pathogenicity and finally determined that B. glumae can undergo phenotypic conversion through a spontaneous mutation in the tofR regulator.     

Differential regulation of toxoflavin production and its role in the enhanced virulence of Burkholderia gladioli

Burkholderia gladioli is a causal agent of bacterial panicle blight and sheath/grain browning in rice in many countries. Many strains produce the yellow pigment toxoflavin, which is highly toxic to plants, fungi, animals and microorganisms. Although there have been several studies on the toxoflavin biosynthesis system of B. glumae, it is still unclear how B. gladioli activates toxoflavin biosynthesis. In this study, we explored the genomic organization of the toxoflavin system of B. gladioli and its biological functions using comparative genomic analysis between toxoflavin‐producing strains (B. glumae BGR1 and B. gladioli BSR3) and a strain not producing toxoflavin (B. gladioli KACC11889). The latter exhibits normal physiological characteristics similar to other B. gladioli strains. Burkholderia gladioli KACC11889 possesses all the genes involved in toxoflavin biosynthesis, but lacks the quorum‐sensing (QS) system that functions as an on/off switch for toxoflavin biosynthesis. These data suggest that B. gladioli has evolved to use the QS signalling cascade of toxoflavin production (TofI/TofR of QS → ToxJ or ToxR → tox operons) similar to that in B. glumae. However, some strains may have evolved to eliminate toxoflavin production through deletion of the QS genes. In addition, we demonstrate that the toxoflavin biosynthetic system enhances the virulence of B. gladioli. These findings provide another line of evidence supporting the differential regulation of the toxoflavin system in Burkholderia strains.     

Comparative analysis of mitochondrial genomes between the hau cytoplasmic male sterility (CMS) line and its iso-nuclear maintainer line in Brassica juncea to reveal the origin of the CMS-associated gene orf288

Background

Cytoplasmic male sterility (CMS) is not only important for exploiting heterosis in crop plants, but also as a model for investigating nuclear-cytoplasmic interaction. CMS may be caused by mutations, rearrangement or recombination in the mitochondrial genome. Understanding the mitochondrial genome is often the first and key step in unraveling the molecular and genetic basis of CMS in plants. Comparative analysis of the mitochondrial genome of the hau CMS line and its maintainer line in B. juneca (Brassica juncea) may help show the origin of the CMS-associated gene orf288.

Results

Through next-generation sequencing, the B. juncea hau CMS mitochondrial genome was assembled into a single, circular-mapping molecule that is 247,903 bp in size and 45.08% in GC content. In addition to the CMS associated gene orf288, the genome contains 35 protein-encoding genes, 3 rRNAs, 25 tRNA genes and 29 ORFs of unknown function. The mitochondrial genome sizes of the maintainer line and another normal type line “J163-4” are both 219,863 bp and with GC content at 45.23%. The maintainer line has 36 genes with protein products, 3 rRNAs, 22 tRNA genes and 31 unidentified ORFs. Comparative analysis the mitochondrial genomes of the hau CMS line and its maintainer line allowed us to develop specific markers to separate the two lines at the seedling stage. We also confirmed that different mitotypes coexist substoichiometrically in hau CMS lines and its maintainer lines in B. juncea. The number of repeats larger than 100 bp in the hau CMS line (16 repeats) are nearly twice of those found in the maintainer line (9 repeats). Phylogenetic analysis of the CMS-associated gene orf288 and four other homologous sequences in Brassicaceae show that orf288 was clearly different from orf263 in Brassica tournefortii despite of strong similarity.

Conclusion

The hau CMS mitochondrial genome was highly rearranged when compared with its iso-nuclear maintainer line mitochondrial genome. This study may be useful for studying the mechanism of natural CMS in B. juncea, performing comparative analysis on sequenced mitochondrial genomes in Brassicas, and uncovering the origin of the hau CMS mitotype and structural and evolutionary differences between different mitotypes.

Electronic supplementary material

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

A triallelic genetic male sterility locus in Brassica napus: an integrative strategy for its physical mapping and possible local chromosome evolution around it

Background and Aims

Spontaneous male sterility is an advantageous trait for both constructing efficient pollination control systems and for understanding the developmental process of the male reproductive unit in many crops. A triallelic genetic male-sterile locus (BnMs5) has been identified in Brassica napus; however, its complicated genome structure has greatly hampered the isolation of this locus. The aim of this study was to physically map BnMs5 through an integrated map-based cloning strategy and analyse the local chromosomal evolution around BnMs5.

Methods

A large F₂ population was used to integrate the existing genetic maps around BnMs5. A map-based cloning strategy in combination with comparative mapping among B. napus, Arabidopsis, Brassica rapa and Brassica oleracea was employed to facilitate the identification of a target bacterial artificial chromosome (BAC) clone covering the BnMs5 locus. The genomic sequences from the Brassica species were analysed to reveal the regional chromosomal evolution around BnMs5.

Key Results

BnMs5 was finally delimited to a 0·3-cM genetic fragment from an integrated local genetic map, and was anchored on the B. napus A8 chromosome. Screening of a B. napus BAC clone library and identification of the positive clones validated that JBnB034L06 was the target BAC clone. The closest flanking markers restrict BnMs5 to a 21-kb region on JBnB034L06 containing six predicted functional genes. Good collinearity relationship around BnMs5 between several Brassica species was observed, while violent chromosomal evolutionary events including insertions/deletions, duplications and single nucleotide mutations were also found to have extensively occurred during their divergence.

Conclusions

This work represents major progress towards the molecular cloning of BnMs5, as well as presenting a powerful, integrative method to mapping loci in plants with complex genomic architecture, such as the amphidiploid B. napus.     

Genomic regions involved in yield potential detected by genome-wide association analysis in Japanese high-yielding rice cultivars

Background

High-yielding cultivars of rice (Oryza sativa L.) have been developed in Japan from crosses between overseas indica and domestic japonica cultivars. Recently, next-generation sequencing technology and high-throughput genotyping systems have shown many single-nucleotide polymorphisms (SNPs) that are proving useful for detailed analysis of genome composition. These SNPs can be used in genome-wide association studies to detect candidate genome regions associated with economically important traits. In this study, we used a custom SNP set to identify introgressed chromosomal regions in a set of high-yielding Japanese rice cultivars, and we performed an association study to identify genome regions associated with yield.

Results

An informative set of 1152 SNPs was established by screening 14 high-yielding or primary ancestral cultivars for 5760 validated SNPs. Analysis of the population structure of high-yielding cultivars showed three genome types: japonica-type, indica-type and a mixture of the two. SNP allele frequencies showed several regions derived predominantly from one of the two parental genome types. Distinct regions skewed for the presence of parental alleles were observed on chromosomes 1, 2, 7, 8, 11 and 12 (indica) and on chromosomes 1, 2 and 6 (japonica). A possible relationship between these introgressed regions and six yield traits (blast susceptibility, heading date, length of unhusked seeds, number of panicles, surface area of unhusked seeds and 1000-grain weight) was detected in eight genome regions dominated by alleles of one parental origin. Two of these regions were near Ghd7, a heading date locus, and Pi-ta, a blast resistance locus. The allele types (i.e., japonica or indica) of significant SNPs coincided with those previously reported for candidate genes Ghd7 and Pi-ta.

Conclusions

Introgression breeding is an established strategy for the accumulation of QTLs and genes controlling high yield. Our custom SNP set is an effective tool for the identification of introgressed genome regions from a particular genetic background. This study demonstrates that changes in genome structure occurred during artificial selection for high yield, and provides information on several genomic regions associated with yield performance.

Electronic supplementary material

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