Ziding Zhang

The Bifidobacterium dentium Bd1 Genome Sequence Reflects Its Genetic Adaptation to the Human Oral Cavity

Bifidobacteria, one of the relatively dominant components of the human intestinal microbiota, are considered one of the key groups of beneficial intestinal bacteria (probiotic bacteria). However, in addition to health-promoting taxa, the genus Bifidobacterium also includes Bifidobacterium dentium, an opportunistic cariogenic pathogen. The genetic basis for the ability of B. dentium to survive in the oral cavity and contribute to caries development is not understood. The genome of B. dentium Bd1, a strain isolated from dental caries, was sequenced to completion to uncover a single circular 2,636,368 base pair chromosome with 2,143 predicted open reading frames. Annotation of the genome sequence revealed multiple ways in which B. dentium has adapted to the oral environment through specialized nutrient acquisition, defences against antimicrobials, and gene products that increase fitness and competitiveness within the oral niche. B. dentium Bd1 was shown to metabolize a wide variety of carbohydrates, consistent with genome-based predictions, while colonization and persistence factors implicated in tissue adhesion, acid tolerance, and the metabolism of human saliva-derived compounds were also identified. Global transcriptome analysis demonstrated that many of the genes encoding these predicted traits are highly expressed under relevant physiological conditions. This is the first report to identify, through various genomic approaches, specific genetic adaptations of a Bifidobacterium taxon, Bifidobacterium dentium Bd1, to a lifestyle as a cariogenic microorganism in the oral cavity. In silico analysis and comparative genomic hybridization experiments clearly reveal a high level of genome conservation among various B. dentium strains. The data indicate that the genome of this opportunistic cariogen has evolved through a very limited number of horizontal gene acquisition events, highlighting the narrow boundaries that separate commensals from opportunistic pathogens.

Badh2, Encoding Betaine Aldehyde Dehydrogenase, Inhibits the Biosynthesis of 2-Acetyl-1-Pyrroline, a Major Component in Rice Fragrance

In rice (Oryza sativa), the presence of a dominant Badh2 allele encoding betaine aldehyde dehydrogenase (BADH2) inhibits the synthesis of 2-acetyl-1-pyrroline (2AP), a potent flavor component in rice fragrance. By contrast, its two recessive alleles, badh2-E2 and badh2-E7, induce 2AP formation. Badh2 was found to be transcribed in all tissues tested except for roots, and the transcript was detected at higher abundance in young, healthy leaves than in other tissues. Multiple Badh2 transcript lengths were detected, and the complete, full-length Badh2 transcript was much less abundant than partial Badh2 transcripts. 2AP levels were significantly reduced in cauliflower mosaic virus 35S-driven transgenic lines expressing the complete, but not the partial, Badh2 coding sequences. In accordance, the intact, full-length BADH2 protein (503 residues) appeared exclusively in nonfragrant transgenic lines and rice varieties. These results indicate that the full-length BADH2 protein encoded by Badh2 renders rice nonfragrant by inhibiting 2AP biosynthesis. The BADH2 enzyme was predicted to contain three domains: NAD binding, substrate binding, and oligomerization domains. BADH2 was distributed throughout the cytoplasm, where it is predicted to catalyze the oxidization of betaine aldehyde, 4-aminobutyraldehyde (AB-ald), and 3-aminopropionaldehyde. The presence of null badh2 alleles resulted in AB-ald accumulation and enhanced 2AP biosynthesis. In summary, these data support the hypothesis that BADH2 inhibits 2AP biosynthesis by exhausting AB-ald, a presumed 2AP precursor.

Prediction of mucin-type O-glycosylation sites in mammalian proteins using the composition of k-spaced amino acid pairs

BackgroundAs one of the most common protein post-translational modifications, glycosylation is involved in a variety of important biological processes. Computational identification of glycosylation sites in protein sequences becomes increasingly important in the post-genomic era. A new encoding scheme was employed to improve the prediction of mucin-type O-glycosylation sites in mammalian proteins.ResultsA new protein bioinformatics tool, CKSAAP_OGlySite, was developed to predict mucin-type O-glycosylation serine/threonine (S/T) sites in mammalian proteins. Using the composition of k-spaced amino acid pairs (CKSAAP) based encoding scheme, the proposed method was trained and tested in a new and stringent O-glycosylation dataset with the assistance of Support Vector Machine (SVM). When the ratio of O-glycosylation to non-glycosylation sites in training datasets was set as 1:1, 10-fold cross-validation tests showed that the proposed method yielded a high accuracy of 83.1% and 81.4% in predicting O-glycosylated S and T sites, respectively. Based on the same datasets, CKSAAP_OGlySite resulted in a higher accuracy than the conventional binary encoding based method (about +5.0%). When trained and tested in 1:5 datasets, the CKSAAP encoding showed a more significant improvement than the binary encoding. We also merged the training datasets of S and T sites and integrated the prediction of S and T sites into one single predictor (i.e. S+T predictor). Either in 1:1 or 1:5 datasets, the performance of this S+T predictor was always slightly better than those predictors where S and T sites were independently predicted, suggesting that the molecular recognition of O-glycosylated S/T sites seems to be similar and the increase of the S+T predictors accuracy may be a result of expanded training datasets. Moreover, CKSAAP_OGlySite was also shown to have better performance when benchmarked against two existing predictors.ConclusionBecause of CKSAAP encodings ability of reflecting characteristics of the sequences surrounding mucin-type O-glycosylation sites, CKSAAP_ OGlySite has been proved more powerful than the conventional binary encoding based method. This suggests that it can be used as a competitive mucin-type O-glycosylation site predictor to the biological community. CKSAAP_OGlySite is now available at http://bioinformatics.cau.edu.cn/zzd_lab/CKSAAP_OGlySite/.

Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae.

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice worldwide. The fungal pathogen is notorious for its ability to overcome host resistance. To better understand its genetic variation in nature, we sequenced the genomes of two field isolates, Y34 and P131. In comparison with the previously sequenced laboratory strain 70-15, both field isolates had a similar genome size but slightly more genes. Sequences from the field isolates were used to improve genome assembly and gene prediction of 70-15. Although the overall genome structure is similar, a number of gene families that are likely involved in plant-fungal interactions are expanded in the field isolates. Genome-wide analysis on asynonymous to synonymous nucleotide substitution rates revealed that many infection-related genes underwent diversifying selection. The field isolates also have hundreds of isolate-specific genes and a number of isolate-specific gene duplication events. Functional characterization of randomly selected isolate-specific genes revealed that they play diverse roles, some of which affect virulence. Furthermore, each genome contains thousands of loci of transposon-like elements, but less than 30% of them are conserved among different isolates, suggesting active transposition events in M. oryzae. A total of approximately 200 genes were disrupted in these three strains by transposable elements. Interestingly, transposon-like elements tend to be associated with isolate-specific or duplicated sequences. Overall, our results indicate that gain or loss of unique genes, DNA duplication, gene family expansion, and frequent translocation of transposon-like elements are important factors in genome variation of the rice blast fungus.

Genetic analysis and morphological identification of pilus-like structures in members of the genus Bifidobacterium

BackgroundCell surface pili in Gram positive bacteria have been reported to orchestrate the colonization of host tissues, evasion of immunity and the development of biofilms. So far, little if any information is available on the presence of pilus-like structures in human gut commensals like bifidobacteria.Results and discussionIn this report, Atomic Force Microscopy (AFM) of various bifidobacterial strains belonging to Bifidobacterium bifidum, Bifidobacterium longum subsp. longum, Bifidobacterium dentium, Bifidobacterium adolescentis and Bifidobacterium animalis subsp. lactis revealed the existence of appendages resembling pilus-like structures. Interestingly, these microorganisms harbour two to six predicted pilus gene clusters in their genome, with each organized in an operon encompassing the major pilin subunit-encoding gene (designated fim A or fim P) together with one or two minor pilin subunit-encoding genes (designated as fim B and/or fim Q), and a gene encoding a sortase enzyme (str A). Quantitative Real Time (qRT)-PCR analysis and RT-PCR experiments revealed a polycistronic mRNA, encompassing the fim A/P and fim B/Q genes, which are differentially expressed upon cultivation of bifidobacteria on various glycans.

Prediction of Ubiquitination Sites by Using the Composition of k-Spaced Amino Acid Pairs

As one of the most important reversible protein post-translation modifications, ubiquitination has been reported to be involved in lots of biological processes and closely implicated with various diseases. To fully decipher the molecular mechanisms of ubiquitination-related biological processes, an initial but crucial step is the recognition of ubiquitylated substrates and the corresponding ubiquitination sites. Here, a new bioinformatics tool named CKSAAP_UbSite was developed to predict ubiquitination sites from protein sequences. With the assistance of Support Vector Machine (SVM), the highlight of CKSAAP_UbSite is to employ the composition of k-spaced amino acid pairs surrounding a query site (i.e. any lysine in a query sequence) as input. When trained and tested in the dataset of yeast ubiquitination sites (Radivojac et al, Proteins, 2010, 78: 365–380), a 100-fold cross-validation on a 1∶1 ratio of positive and negative samples revealed that the accuracy and MCC of CKSAAP_UbSite reached 73.40% and 0.4694, respectively. The proposed CKSAAP_UbSite has also been intensively benchmarked to exhibit better performance than some existing predictors, suggesting that it can be served as a useful tool to the community. Currently, CKSAAP_UbSite is freely accessible at http://protein.cau.edu.cn/cksaap_ubsite/. Moreover, we also found that the sequence patterns around ubiquitination sites are not conserved across different species. To ensure a reasonable prediction performance, the application of the current CKSAAP_UbSite should be limited to the proteome of yeast.

MDP25, A Novel Calcium Regulatory Protein, Mediates Hypocotyl Cell Elongation by Destabilizing Cortical Microtubules in Arabidopsis

This study reveals that the Arabidopsis thaliana microtubule-destabilizing protein MDP25 regulates hypocotyl elongation. The microtubule-destabilizing activity of MDP25 is regulated by calcium, which causes MDP25 to disassociate from the plasma membrane and move into the cytosol. We identified a novel microtubule mechanism that regulates hypocotyl cell elongation in Arabidopsis. The regulation of hypocotyl elongation is important for plant growth. Microtubules play a crucial role during hypocotyl cell elongation. However, the molecular mechanism underlying this process is not well understood. In this study, we describe a novel Arabidopsis thaliana microtubule-destabilizing protein 25 (MDP25) as a negative regulator of hypocotyl cell elongation. We found that MDP25 directly bound to and destabilized microtubules to enhance microtubule depolymerization in vitro. The seedlings of mdp25 mutant Arabidopsis lines had longer etiolated hypocotyls. In addition, MDP25 overexpression resulted in significant overall shortening of hypocotyl cells, which exhibited destabilized cortical microtubules and abnormal cortical microtubule orientation, suggesting that MDP25 plays a crucial role in the negative regulation of hypocotyl cell elongation. Although MDP25 localized to the plasma membrane under normal conditions, increased calcium levels in cells caused MDP25 to partially dissociate from the plasma membrane and move into the cytosol. Cellular MDP25 bound to and destabilized cortical microtubules, resulting in their reorientation, and subsequently inhibited hypocotyl cell elongation. Our results suggest that MDP25 exerts its function on cortical microtubules by responding to cytoplasmic calcium levels to mediate hypocotyl cell elongation.

INTERACTION OF POSITIONAL ISOMERS OF QUERCETIN GLUCURONIDES WITH THE TRANSPORTER ABCC2 (CMOAT, MRP2)

The exporter ABCC2 (cMOAT, MRP2) is a membrane-bound protein on the apical side of enterocytes and hepatic biliary vessels that transports leukotriene C4, glutathione, some conjugated bile salts, drugs, xenobiotics, and phytonutrients. The latter class includes quercetin, a bioactive flavonoid found in foods such as onions, apples, tea, and wine. There is no available three-dimensional (3D) structure of ABCC2. We have predicted the 3D structure by in silico modeling, showing that 3-[[3-[2-(7-chloroquinolin-2-yl)vinyl]phenyl]-(2-dimethylcarbamoylethylsulfanyl)methylsulfanyl] propionic acid (MK571) binds most tightly to the putative binding site, and then tested the computational prediction experimentally by measuring interaction with all quercetin monoglucuronides occurring in vivo (quercetin substituted with glucuronic acid at the 3-, 3′-, 4′-, and 7-hydroxyl groups). The 4′-O-β-D-glucuronide is predicted in silico to interact most strongly and the 3-O-β-D-glucuronide most weakly, and this prediction is supported experimentally using binding and competition assays on ABCC2-overexpressing baculovirus-infected Sf9 cells. To test the transport in situ, we examined the effect of two ABCC2 inhibitors, MK571 and cyclosporin A, on the transport into the media of quercetin glucuronides produced intracellularly by Caco2 cells. The inhibitors reduced the amount of all quercetin glucuronides in the media. The results show that the molecular model of ABCC2 agrees well with experimentally determined ABCC2-ligand interactions and, importantly, that the interaction of ABCC2 with quercetin glucuronides is dependent on the position and nature of substitution.

Similarity networks of protein binding sites

An increasing attention has been dedicated to the characterization of complex networks within the protein world. This work is reporting how we uncovered networked structures that reflected the structural similarities among protein binding sites. First, a 211 binding sites dataset has been compiled by removing the redundant proteins in the Protein Ligand Database (PLD) (http://www‐mitchell.ch.cam.ac.uk/pld/). Using a clique detection algorithm we have performed all‐against‐all binding site comparisons among the 211 available ones. Within the set of nodes representing each binding site an edge was added whenever a pair of binding sites had a similarity higher than a threshold value. The generated similarity networks revealed that many nodes had few links and only few were highly connected, but due to the limited data available it was not possible to definitively prove a scale‐free architecture. Within the same dataset, the binding site similarity networks were compared with the networks of sequence and fold similarity networks. In the protein world, indications were found that structure is better conserved than sequence, but on its own, sequence was better conserved than the subset of functional residues forming the binding site. Because a binding site is strongly linked with protein function, the identification of protein binding site similarity networks could accelerate the functional annotation of newly identified genes. In view of this we have discussed several potential applications of binding site similarity networks, such as the construction of novel binding site classification databases, as well as the implications for protein molecular design in general and computational chemogenomics in particular. Proteins 2006.

Characterization of the serpin-encoding gene of Bifidobacterium breve 210B

ABSTRACT Members of the serpin (serine protease inhibitor) superfamily have been identified in higher multicellular eukaryotes, as well as in bacteria, although examination of available genome sequences has indicated that homologs of the bacterial serpin-encoding gene (ser) are not widely distributed. In members of the genus Bifidobacterium this gene appears to be present in at least 5, and perhaps up to 9, of the 30 species tested. Moreover, phylogenetic analysis using available bacterial and eukaryotic serpin sequences revealed that bifidobacteria produce serpins that form a separate clade. We characterized the ser210B locus of Bifidobacterium breve 210B, which encompasses a number of genes whose deduced protein products display significant similarity to proteins encoded by corresponding loci found in several other bifidobacteria. Northern hybridization, primer extension, microarray, reverse transcription-PCR (RT-PCR), and quantitative real-time PCR (qRT-PCR) analyses revealed that a 3.5-kb polycistronic mRNA encompassing the ser210B operon with a single transcriptional start site is strongly induced following treatment of B. breve 210B cultures with some proteases. Interestingly, transcription of other bifidobacterial ser homologs appears to be triggered by different proteases.