De novo transcriptome analyses of host-fungal interactions in oil palm (Elaeis guineensis Jacq.).

BACKGROUND: Basal stem rot (BSR) is a fungal disease in oil palm (Elaeis guineensis Jacq.) which is caused by hemibiotrophic white rot fungi belonging to the Ganoderma genus. Molecular responses of oil palm to these pathogens are not well known although this information is crucial to strategize effective measures to eradicate BSR. In order to elucidate the molecular interactions between oil palm and G. boninense and its biocontrol fungus Trichoderma harzianum, we compared the root transcriptomes of untreated oil palm seedlings with those inoculated with G. boninense and T. harzianum, respectively. RESULTS: Differential gene expression analyses revealed that jasmonate (JA) and salicylate (SA) may act in an antagonistic manner in affecting the hormone biosynthesis, signaling, and downstream defense responses in G. boninense-treated oil palm roots. In addition, G. boninense may compete with the host to control disease symptom through the transcriptional regulation of ethylene (ET) biosynthesis, reactive oxygen species (ROS) production and scavenging. The strengthening of host cell walls and production of pathogenesis-related proteins as well as antifungal secondary metabolites in host plants, are among the important defense mechanisms deployed by oil palm against G. boninense. Meanwhile, endophytic T. harzianum was shown to improve the of nutrition status and nutrient transportation in host plants. CONCLUSION: The findings of this analysis have enhanced our understanding on the molecular interactions of G. boninense and oil palm, and also the biocontrol mechanisms involving T. harzianum, thus contributing to future formulations of better strategies for prevention and treatment of BSR.

ClicO FS: an interactive web-based service of Circos.

UNLABELLED: : We present ClicO Free Service, an online web-service based on Circos, which provides a user-friendly, interactive web-based interface with configurable features to generate Circos circular plots. AVAILABILITY AND IMPLEMENTATION: Online web-service is freely available at http://clicofs.codoncloud.com CONTACT: : soonjoo.yap@codongenomics.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis.

BACKGROUND: Brucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis. RESULTS: Single sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage. CONCLUSION: We showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.

Comparative genomic and phenotypic analyses of pathogenic fungi Neoscytalidium dimidiatum and Bipolaris papendorfii isolated from human skin scraping.

Neoscytalidium dimidiatum and Bipolaris species are fungal plant pathogens that have been reported to cause human diseases. Recently, we have isolated numerous N. dimidiatum and Bipolaris species from the skin scrapings and nails of different patients. In this work, we have sequenced the genome of one strain of N. dimidiatum. The sequenced genome was compared to that of a previously reported Bipolaris papendorfii genome for a better understanding of their complex lifestyle and broad host-range pathogenicity. Both N. dimidiatum UM 880 (~ 43 Mb) and B. papendorfii UM 226 (~ 33 Mb) genomes include 11,015-12,320 putative coding DNA sequences, of which 0.51-2.49% are predicted transposable elements. Analysis of secondary metabolism gene clusters revealed several genes involved in melanin biosynthesis and iron uptake. The arsenal of CAZymes related to plants pathogenicity is comparable between the species, including genes involved in hemicellulose and pectin decomposition. Several important gene encoding keratinolytic peptidases were identified in N. dimidiatum and B. papendorfii, reflecting their potential pathogenic role in causing skin and nail infections. In this study, additional information on the metabolic features of these two species, such as nutritional profiling, pH tolerance, and osmotolerant, are revealed. The genomic characterization of N. dimidiatum and B. papendorfii provides the basis for the future functional studies to gain further insights as to what makes these fungi persist in plants and why they are pathogenic to humans.

Polyhydroxyalkanoate (PHA) synthase genes and PHA-associated gene clusters in Pseudomonas spp. and Janthinobacterium spp. isolated from Antarctica.

The polyhydroxyalkanoate (PHA) producing capability of four bacterial strains isolated from Antarctica was reported in a previous study. This study analyzed the PHA synthase genes and the PHA-associated gene clusters from the two antarctic Pseudomonas isolates (UMAB-08 and UMAB-40) and the two antarctic Janthinobacterium isolates (UMAB-56 and UMAB-60) through whole-genome sequence analysis. The Pseudomonas isolates were found to carry PHA synthase genes which fall into two different PHA gene clusters, namely Class I and Class II, which are involved in the biosynthesis of short-chain-length-PHA (SCL-PHA) and medium-chain-length-PHA (MCL-PHA), respectively. On the other hand, the Janthinobacterium isolates carry a Class I and an uncharacterized putative PHA synthase genes. No other gene involved in PHA synthesis was detected in close proximity to the uncharacterized putative PHA synthase gene in the Janthinobacterium isolates, therefore it falls into a separate clade from the ordinary Class I, II, III and IV clades of PHA synthase (PhaC) phylogenetic tree. Multiple sequence alignment showed that the uncharacterized putative PHA synthase gene contains all the highly conserved amino acid residues and the proposed catalytic triad of PHA synthase. PHA biosynthesis and in vitro PhaC enzymatic assay results showed that this uncharacterized putative PHA synthase from Janthinobacterium sp. UMAB-60 is funtional. This report adds new knowledge to the PHA synthase database as we describe scarce information of PHA synthase genes and PHA-associated gene clusters from the antarctic bacterial isolates (extreme and geographically isolated environment) and comparing with those from non-antarctic PHA-producing bacteria.

Geographical distribution of Brucella melitensis inferred from rpoB gene variation.

INTRODUCTION: Currently available tests have limitations for the identification of Brucella species and strains, and their genetic lineage. The genome sequence of the rpoB gene encoding the ß-subunit of DNA-dependent RNA polymerase was investigated for its use in genotyping Brucella melitensis. METHODOLOGY: Complete rpoB gene sequences of globally distributed Brucella melitensis strains were analyzed. Single nucleotides polymorphisms (SNPs) of the rpoB gene sequences were identified and used to type Brucella melitensis strains. RESULTS: Six DNA polymorphisms were identified, of which two (nucleotides 3201 and 558) were novel. Analysis of the geographical distribution of the strains revealed a spatial clustering pattern with rpoB type 1 representing European and American strains, rpoB type 2 representing European, African, and Asian strains, rpoB type 3 representing Mediterranean strains, and rpoB type 4 representing African (C3201T) and European (C3201T/T558A) strains. CONCLUSIONS: We report the discovery of two novel SNPs of rpoB gene that can serve as useful markers for epidemiology and geographical tracking of B. melitensis.

Complete Genome Sequence of Pasteurella multocida Serotype A Strain PMTB2.1 Isolated from Buffaloes That Died of Septicemia in Malaysia.

Pasteurella multocida causes pneumonic pasteurellosis and hemorrhagic septicemia (HS) in large ruminants. In this study, we determined the complete genome sequence of P. multocida strain PMTB2.1 capsular serotype A isolated from buffaloes that died of septicemia.

Genomic insight into pathogenicity of dematiaceous fungus Corynespora cassiicola.

Corynespora cassiicola is a common plant pathogen that causes leaf spot disease in a broad range of crop, and it heavily affect rubber trees in Malaysia (Hsueh, 2011; Nghia et al., 2008). The isolation of UM 591 from a patient's contact lens indicates the pathogenic potential of this dematiaceous fungus in human. However, the underlying factors that contribute to the opportunistic cross-infection have not been fully studied. We employed genome sequencing and gene homology annotations in attempt to identify these factors in UM 591 using data obtained from publicly available bioinformatics databases. The assembly size of UM 591 genome is 41.8 Mbp, and a total of 13,531 (≥99 bp) genes have been predicted. UM 591 is enriched with genes that encode for glycoside hydrolases, carbohydrate esterases, auxiliary activity enzymes and cell wall degrading enzymes. Virulent genes comprising of CAZymes, peptidases, and hypervirulence-associated cutinases were found to be present in the fungal genome. Comparative analysis result shows that UM 591 possesses higher number of carbohydrate esterases family 10 (CE10) CAZymes compared to other species of fungi in this study, and these enzymes hydrolyses wide range of carbohydrate and non-carbohydrate substrates. Putative melanin, siderophore, ent-kaurene, and lycopene biosynthesis gene clusters are predicted, and these gene clusters denote that UM 591 are capable of protecting itself from the UV and chemical stresses, allowing it to adapt to different environment. Putative sterigmatocystin, HC-toxin, cercosporin, and gliotoxin biosynthesis gene cluster are predicted. This finding have highlighted the necrotrophic and invasive nature of UM 591.

Oil Palm Defensin: A Thermal Stable Peptide that Restricts the Mycelial Growth of Ganoderma boninense.

Plant defensins are plant defence peptides that have many different biological activities, including antifungal, antimicrobial, and insecticidal activities. A cDNA (EgDFS) encoding defensin was isolated from Elaeis guineensis. The open reading frame of EgDFS contained 231 nucleotides encoding a 71-amino acid protein with a predicted molecular weight at 8.69 kDa, and a potential signal peptide. The eight highly conserved cysteine sites in plant defensins were also conserved in EgDFS. The EgDFS sequence lacking 30 amino acid residues at its N-terminus (EgDFSm) was cloned into Escherichia coli BL21 (DE3) pLysS and successfully expressed as a soluble recombinant protein. The recombinant EgDFSm was found to be a thermal stable peptide which demonstrated inhibitory activity against the growth of G. boninense possibly by inhibiting starch assimilation. The role of EgDFSm in oil palm defence system against the infection of pathogen G. boninense was discussed.

Candidate Essential Genes in Burkholderia cenocepacia J2315 Identified by Genome-Wide TraDIS.

Burkholderia cenocepacia infection often leads to fatal cepacia syndrome in cystic fibrosis patients. However, antibiotic therapy rarely results in complete eradication of the pathogen due to its intrinsic resistance to many clinically available antibiotics. Recent attention has turned to the identification of essential genes as the proteins encoded by these genes may serve as potential targets for development of novel antimicrobials. In this study, we utilized TraDIS (Transposon Directed Insertion-site Sequencing) as a genome-wide screening tool to facilitate the identification of B. cenocepacia genes essential for its growth and viability. A transposon mutant pool consisting of approximately 500,000 mutants was successfully constructed, with more than 400,000 unique transposon insertion sites identified by computational analysis of TraDIS datasets. The saturated library allowed for the identification of 383 genes that were predicted to be essential in B. cenocepacia. We extended the application of TraDIS to identify conditionally essential genes required for in vitro growth and revealed an additional repertoire of 439 genes to be crucial for B. cenocepacia growth under nutrient-depleted conditions. The library of B. cenocepacia mutants can subsequently be subjected to various biologically related conditions to facilitate the discovery of genes involved in niche adaptation as well as pathogenicity and virulence.

Genomic Analyses of Cladophialophora bantiana, a Major Cause of Cerebral Phaeohyphomycosis Provides Insight into Its Lifestyle, Virulence and Adaption in Host.

Cladophialophora bantiana is a dematiaceous fungus with a predilection for causing central nervous system (CNS) infection manifesting as brain abscess in both immunocompetent and immunocompromised patients. In this paper, we report comprehensive genomic analyses of C. bantiana isolated from the brain abscess of an immunocompetent man, the first reported case in Malaysia and Southeast Asia. The identity of the fungus was determined using combined morphological analysis and multilocus phylogeny. The draft genome sequence of a neurotrophic fungus, C. bantiana UM 956 was generated using Illumina sequencing technology to dissect its genetic fundamental and basic biology. The assembled 37.1 Mb genome encodes 12,155 putative coding genes, of which, 1.01% are predicted transposable elements. Its genomic features support its saprophytic lifestyle, renowned for its versatility in decomposing hemicellulose and pectin components. The C. bantiana UM 956 was also found to carry some important putative genes that engaged in pathogenicity, iron uptake and homeostasis as well as adaptation to various stresses to enable the organism to survive in hostile microenvironment. This wealth of resource will further catalyse more downstream functional studies to provide better understanding on how this fungus can be a successful and persistent pathogen in human.

Expression profiles of defence related cDNAs in oil palm (Elaeis guineensis Jacq.) inoculated with mycorrhizae and Trichoderma harzianum Rifai T32.

Basal stem rot is one of the major diseases of oil palm (Elaies guineensis Jacq.) caused by pathogenic Ganoderma species. Trichoderma and mycorrhizae were proposed to be able to reduce the disease severity. However, their roles in improving oil palm defence system by possibly inducing defence-related genes in the host are not well characterized. To better understand that, transcript profiles of eleven putative defence-related cDNAs in the roots of oil palm inoculated with Trichoderma harzianum T32 and mycorrhizae at different time points were studied. Transcripts encoding putative Bowman-Birk protease inhibitor (EgBBI2) and defensin (EgDFS) increased more than 2 fold in mycorrhizae-treated roots at 6 weeks post inoculation (wpi) compared to those in controls. Transcripts encoding putative dehydrin (EgDHN), glycine-rich RNA binding protein (EgGRRBP), isoflavone reductase (EgIFR), type 2 ribosome inactivating protein (EgT2RIP), and EgDFS increased in the oil palm roots treated with T. harzianum at 6 and/or 12 wpi compared to those in the controls. Some of these genes were also expressed in oil palm roots treated with Ganoderma boninense. This study provides an insight of some defence-related genes induced by Trichoderma and mycorrhizae, and their roles as potential agents to boost the plant defence system.

Genome Analysis of the First Extensively Drug-Resistant (XDR) Mycobacterium tuberculosis in Malaysia Provides Insights into the Genetic Basis of Its Biology and Drug Resistance.

The outbreak of extensively drug-resistant tuberculosis (XDR-TB) has become an increasing problem in many TB-burdened countries. The underlying drug resistance mechanisms, including the genetic variation favored by selective pressure in the resistant population, are partially understood. Recently, the first case of XDR-TB was reported in Malaysia. However, the detailed genotype family and mechanisms of the formation of multiple drugs resistance are unknown. We sequenced the whole genome of the UM 1072388579 strain with a 2-kb insert-size library and combined with that from previously sequenced 500-bp-insert paired-end reads to produce an improved sequence with maximal sequencing coverage across the genome. In silico spoligotyping and phylogenetic analyses demonstrated that UM 1072388579 strain belongs to an ancestral-like, non-Beijing clade of East Asia lineage. This is supported by the presence of a number of lineage-specific markers, including fadD28, embA, nuoD and pks7. Polymorphism analysis showed that the drug-susceptibility profile is correlated with the pattern of resistance mutations. Mutations in drug-efflux pumps and the cell wall biogenesis pathway such as mmpL, pks and fadD genes may play an important role in survival and adaptation of this strain to its surrounding environment. In this work, fifty-seven putative promoter SNPs were identified. Among them, we identified a novel SNP located at -4 T allele of TetR/acrR promoter as an informative marker to recognize strains of East Asian lineage. Our work indicates that the UM 1072388579 harbors both classical and uncommon SNPs that allow it to escape from inhibition by many antibiotics. This study provides a strong foundation to dissect the biology and underlying resistance mechanisms of the first reported XDR M. tuberculosis in Malaysia.

Draft Genome Sequence of Ochroconis constricta UM 578, Isolated from Human Skin Scraping.

Ochroconis constricta is a soilborne dematiaceous fungus that has never been reported to be associated with human infection. Here we report the first draft genome sequence of strain UM 578, isolated from human skin scraping. The genomic information revealed will contribute to a better understanding of this species.

Expression profiles of putative defence-related proteins in oil palm (Elaeis guineensis) colonized by Ganoderma boninense.

Basal stem rot (BSR) is a major disease of oil palm caused by a pathogenic fungus, Ganoderma boninense. However, the interaction between the host plant and its pathogen is not well characterized. To better understand the response of oil palm to G. boninense, transcript profiles of eleven putative defence-related genes from oil palm were measured by quantitative reverse-transcription (qRT)-PCR in the roots of oil palms treated with G. boninense from 3 to 12 weeks post infection (wpi). These transcripts encode putative Bowman-Birk serine protease inhibitors (EgBBI1 and 2), defensin (EgDFS), dehydrin (EgDHN), early methionine-labeled polypeptides (EgEMLP1 and 2), glycine-rich RNA binding protein (EgGRRBP), isoflavone reductase (EgIFR), metallothionein-like protein (EgMT), pathogenesis-related-1 protein (EgPRP), and type 2 ribosome-inactivating protein (EgT2RIP). The transcript abundance of EgBBI2 increased in G. boninense-treated roots at 3 and 6wpi compared to those of controls; while the transcript abundance of EgBBI1, EgDFS, EgEMLP1, EgMT, and EgT2RIP increased in G. boninense-treated roots at 6 or 12wpi. Meanwhile, the gene expression of EgDHN was up-regulated at all three time points in G. boninense-treated roots. The expression profiles of the eleven transcripts were also studied in leaf samples upon inoculation of G. boninense and Trichoderma harzianum to identify potential biomarkers for early detection of BSR. Two candidate genes (EgEMLP1 and EgMT) that have different profiles in G. boninense-treated leaves compared to those infected by T. harzianum may have the potential to be developed as biomarkers for early detection of G. boninense infection.