Substrate entropy in enzyme enantioselectivity: an experimental and molecular modeling study of a lipase

The temperature dependence of the enantioselectivity of Candida antarctica lipase B for 3-hexanol, 2-butanol, 3-methyl-2-butanol, 3,3-dimethyl-2-butanol, and 1-bromo-2-butanol revealed that the differential activation entropy, deltaR-SdeltaS, was as significant as the differential activation enthalpy, DeltaR-SdeltaH, to the enantiomeric ratio, E. 1-Bromo-2-butanol, with isosteric substituents, displayed the largest deltaR-SdeltaS. 3-Hexanol displayed, contrary to other sec-alcohols, a positive deltaR-SdeltaS. In other words, for 3-hexanol the preferred R-enantiomer is not only favored by enthalpy but also by entropy. Molecular dynamics (MD) simulations and systematic search calculations of the substrate accessible volume within the active site revealed that the (R)-3-hexanol transition state (TS) accessed a larger volume within the active site than the (S)-3-hexanol TS. This correlates well with the higher TS entropy of (R)-3-hexanol. In addition, this enantiomer did also yield a higher number of allowed conformations, N, from the systematic search routines, than did the S-enantiomer. The substrate accessible volume was greater for the enantiomer preferred by entropy also for 2-butanol. For 3,3-dimethyl-2-butanol, however, neither MD-simulations nor systematic search calculations yielded substrate accessible volumes that correlate to TS entropy. Ambiguous results were achieved for 3-methyl-2-butanol.     

Anthocyanidin synthase from Gerbera hybrida catalyzes the conversion of (+)-catechin to cyanidin and a novel procyanidin

Anthocyanidins were proposed to derive from (+)-naringenin via (2R,3R)-dihydroflavonol(s) and (2R,3S,4S)-leucocyanidin(s) which are eventually oxidized by anthocyanidin synthase (ANS). Recently, the role of ANS has been put into question, because the recombinant enzyme from Arabidopsis exhibited primarily flavonol synthase (FLS) activity with negligible ANS activity. This and other studies led to the proposal that ANS as well as FLS may select for dihydroflavonoid substrates carrying a "beta-face" C-3 hydroxyl group and initially form the 3-geminal diol by "alpha-face" hydroxylation. Assays with recombinant ANS from Gerbera hybrida fully supported the proposal and were extended to catechin and epicatechin isomers as potential substrates to delineate the enzyme specificity. Gerbera ANS converted (+)-catechin to two major and one minor product, whereas ent(-)-catechin (2S,3R-trans-catechin), (-)-epicatechin, ent(+)-epicatechin (2S,3S-cis-epicatechin) and (-)-gallocatechin were not accepted. The K(m) value for (+)-catechin was determined at 175 microM, and the products were identified by LC-MS(n) and NMR as the 4,4-dimer of oxidized (+)-catechin (93%), cyanidin (7%) and quercetin (trace). When these incubations were repeated in the presence of UDP-glucose:flavonoid 3-O-glucosyltransferase from Fragariaxananassa (FaGT1), the product ratio shifted to cyanidin 3-O-glucoside (60%), cyanidin (14%) and dimeric oxidized (+)-catechin (26%) at an overall equivalent rate of conversion. The data appear to identify (+)-catechin as another substrate of ANS in vivo and shed new light on the mechanism of its catalysis. Moreover, the enzymatic dimerization of catechin monomers is reported for the first time suggesting a role for ANS beyond the oxidation of leucocyanidins.     

The photosystem II-associated Cah3 in Chlamydomonas enhances the O2 evolution rate by proton removal

Water oxidation in photosystem II (PSII) is still insufficiently understood and is assumed to involve HCO(3)(-). A Chlamydomonas mutant lacking a carbonic anhydrase associated with the PSII donor side shows impaired O(2) evolution in the absence of HCO(3)(-). The O(2) evolution for saturating, continuous illumination (R(O2)) was slower than in the wild type, but was elevated by HCO(3)(-) and increased further by Cah3. The R(O2) limitation in the absence of Cah3/HCO(3)(-) was amplified by H(2)O/D(2)O exchange, but relieved by an amphiphilic proton carrier, suggesting a role of Cah3/HCO(3)(-) in proton translocation. Chlorophyll fluorescence indicates a Cah3/HCO(3)(-) effect at the donor side of PSII. Time-resolved delayed fluorescence and O(2)-release measurements suggest specific effects on proton-release steps but not on electron transfer. We propose that Cah3 promotes proton removal from the Mn complex by locally providing HCO(3)(-), which may function as proton carrier. Without Cah3, proton removal could become rate limiting during O(2) formation and thus, limit water oxidation under high light. Our results underlie the general importance of proton release at the donor side of PSII during water oxidation.     

Novel regulation of calcium inhibition of the inositol 1,4,5-trisphosphate receptor calcium-release channel

The inositol 1,4,5-trisphosphate (InsP3) receptor (InsP3R), a Ca2+-release channel localized to the endoplasmic reticulum, plays a critical role in generating complex cytoplasmic Ca2+ signals in many cell types. Three InsP3R isoforms are expressed in different subcellular locations, at variable relative levels with heteromultimer formation in different cell types. A proposed reason for this diversity of InsP3R expression is that the isoforms are differentially inhibited by high cytoplasmic free Ca2+ concentrations ([Ca2+]i), possibly due to their different interactions with calmodulin. Here, we have investigated the possible roles of calmodulin and bath [Ca2+] in mediating high [Ca2+]i inhibition of InsP3R gating by studying single endogenous type 1 InsP3R channels through patch clamp electrophysiology of the outer membrane of isolated Xenopus oocyte nuclei. Neither high concentrations of a calmodulin antagonist nor overexpression of a dominant-negative Ca2+-insensitive mutant calmodulin affected inhibition of gating by high [Ca2+]i. However, a novel, calmodulin-independent regulation of [Ca2+]i inhibition of gating was revealed: whereas channels recorded from nuclei kept in the regular bathing solution with [Ca2+] approximately 400 nM were inhibited by 290 muM [Ca2+]i, exposure of the isolated nuclei to a bath solution with ultra-low [Ca2+] (<5 nM, for approximately 300 s) before the patch-clamp experiments reversibly relieved Ca2+ inhibition, with channel activities observed in [Ca2+]i up to 1.5 mM. Although InsP3 activates gating by relieving high [Ca2+]i inhibition, it was nevertheless still required to activate channels that lacked high [Ca2+]i inhibition. Our observations suggest that high [Ca2+]i inhibition of InsP3R channel gating is not regulated by calmodulin, whereas it can be disrupted by environmental conditions experienced by the channel, raising the possibility that presence or absence of high [Ca2+]i inhibition may not be an immutable property of different InsP3R isoforms. Furthermore, these observations support an allosteric model in which Ca2+ inhibition of the InsP3R is mediated by two Ca2+ binding sites, only one of which is sensitive to InsP3.     

Genetics and Mechanism of Resistance to Meloidogyne arenaria in Peanut Germplasm

Segregation of resistance to Meloidogyne arenaria in six BC₅F₂ peanut breeding populations was examined in greenhouse tests. Chi-square analysis indicated that segregation of resistance was consistent with resistance being conditioned by a single gene in three breeding populations (TP259-3, TP262-3, and TP271-2), whereas two resistance genes may be present in the breeding populations TP259-2, TP263-2, and TP268-3. Nematode development in clonally propagated lines of resistant individuals of TP262-3 and TP263-2 was compared to that of the susceptible cultivar Florunner. Juvenile nematodes readily penetrated roots of all peanut genotypes, but rate of development was slower (P = 0.05) in the resistant genotypes than in Florunner. Host cell necrosis indicative of a hypersensitive response was not consistently observed in resistant genotypes of either population. Three RFLP loci linked to resistance at distances of 4.2 to 11.0 centiMorgans were identified. Resistant and susceptible alleles for RFLP loci R2430E and R2545E were quite distinct and are useful for identifying individuals homozygous for resistance in segregating populations.     

Endothelial P2X7 receptors’ expression is reduced by schistosomiasis

Endothelial cells control vascular tone, permeability and leukocyte transmigration and are modulated by pro-inflammatory mediators. Schistosomiasis is an intravascular disease associated with inflammation, therefore altering endothelial cells’ phenotype. Purinergic P2X7 receptors (P2X7R) play an important role in inflammation; however, the impact of the disease upon endothelial P2X7R function or expression has not been explored. Using ethidium bromide uptake to investigate P2X7R function, we observed that the effects of ATP (3 mM) and the P2X7R agonist 3′-O-(4-benzoyl)-ATP (BzATP) were smaller in mesenteric endothelial cells from the Schistosoma mansoni-infected group than in the control group. In the control group, BzATP induced endothelial nitric oxide production, which was blocked by the P2X7R antagonists KN-62 and A740003. However, in the infected group, we observed a reduced effect of BzATP and no effect of both P2X7R antagonists, suggesting a downregulation of endothelial P2X7R in schistosomiasis. We observed similar results in both infected and P2X7R^−/− groups, which were also comparable to data obtained with KN-62- or A740004-treated control cells. Data from Western blot and immunocytochemistry assays confirmed the reduced expression of P2X7R in the infected group. In conclusion, our data show a downregulation of P2X7R in schistosomiasis infection, which likely limits the infection-related endothelial damage.     

Interrelationships of Rotylenchulus reniformis with Rhizoctonia solani on Cotton

The interrelationships between reniform nematode (Rotylenchulus reniformis) and the cotton (Gossypium hirsutum) seedling blight fungus (Rhizoctonia solani) were studied using three isolates of R. solani, two populations of R. reniformis at multiple inoculum levels, and the cotton cultivars Dehapine 90 (DP 90) and Dehapine 41 (DP 41). Colonization of cotton hypocotyl tissue by R. solani resulted in increases (P ≤ 0.05) in nematode population densities in soil and in eggs recovered from the root systems in both 40- and 90-day-duration experiments. Increases in soil population densities resulted mainly from increases in juveniles. Enhanced reproduction of R. reniformis in the presence of R. solani was consistent across isolates (1, 2, and 3) of R. solani and populations (1 and 2) and inoculum levels (0.5, 2, 4, and 8 individuals/g of soil) of R. reniformis, regardless of cotton cultivar (DP 90 or DP 41). Severity of seedling blight was not influenced by the nematode. Rhizoctonia solani caused reductions (P ≤ 0.05) in cotton growth in 40- and 90-day periods. Rotylenchulus reniformis reduced cotton growth at 90 days. The relationship between nematode inoculum levels and plant growth reductions was linear. At 90 days, the combined effects of these pathogens were antagonistic to plant growth.     

Characterization of Resistance in a Somatic Hybrid of Solanum bulbocastanum and S. tuberosum,to Meloidogyne chitwoodi

A somatic hybrid, CBP-233, between resistant Solanum bulbocastanum (SB-22) and susceptible S. tuberosum (R4) was tested for resistance to Meloidogyne chitwoodi race 1. One week after inoculation, only 0.04-0.4% of the initial inoculum (Pi, 5,000 eggs) as second stage-juveniles infected SB-22 and CBP-233 root systems, compared to 2% in R4. After 8 weeks, the number of M. chitwoodi in SB-22 and CBP-233 roots remained lower (0.3-1.5% of Pi) compared to R4, which increased from 2% to ca. 27%. Development of M. chitwoodi was delayed on SB-22 and CBP-233 by at least 2 weeks, and only half of the infective nematodes established feeding sites and matured in resistant clones compared to 99% in susceptible R4. Necrotic tissue surrounded nematodes that failed to develop in SB-22 and CBP-233. The reproductive factor (ratio of final number of eggs recovered from roots to Pi) was <0.01 for both SB-22 and CBP-233 and 46.8 for R4. Delaying inoculation of CBP-233 from 1 to 3 months after planting did not increase the chance or rate of tuber infection. Only a few M. chitwoodi developed to maturity on CBP-233 tubers and deposited a small number of eggs. SB-22 rarely produced tubers in these experiments, and like CBP-233 were resistant to M. chitwoodi. It appeared that the mechanisms of resistance to M. chitwoodi in roots and tubers of CBP-233 are similar.     

ATP-dependent adenophostin activation of inositol 1,4,5-trisphosphate receptor channel gating: kinetic implications for the durations of calcium puffs in cells

The inositol 1,4,5-trisphosphate (InsP(3)) receptor (InsP(3)R) is a ligand-gated intracellular Ca(2+) release channel that plays a central role in modulating cytoplasmic free Ca(2+) concentration ([Ca(2+)](i)). The fungal metabolite adenophostin A (AdA) is a potent agonist of the InsP(3)R that is structurally different from InsP(3) and elicits distinct calcium signals in cells. We have investigated the effects of AdA and its analogues on single-channel activities of the InsP(3)R in the outer membrane of isolated Xenopus laevis oocyte nuclei. InsP(3)R activated by either AdA or InsP(3) have identical channel conductance properties. Furthermore, AdA, like InsP(3), activates the channel by tuning Ca(2+) inhibition of gating. However, gating of the AdA-liganded InsP(3)R has a critical dependence on cytoplasmic ATP free acid concentration not observed for InsP(3)-liganded channels. Channel gating activated by AdA is indistinguishable from that elicited by InsP(3) in the presence of 0.5 mM ATP, although the functional affinity of the channel is 60-fold higher for AdA. However, in the absence of ATP, gating kinetics of AdA-liganded InsP(3)R were very different. Channel open time was reduced by 50%, resulting in substantially lower maximum open probability than channels activated by AdA in the presence of ATP, or by InsP(3) in the presence or absence of ATP. Also, the higher functional affinity of InsP(3)R for AdA than for InsP(3) is nearly abolished in the absence of ATP. Low affinity AdA analogues furanophostin and ribophostin activated InsP(3)R channels with gating properties similar to those of AdA. These results provide novel insights for interpretations of observed effects of AdA on calcium signaling, including the mechanisms that determine the durations of elementary Ca(2+) release events in cells. Comparisons of single-channel gating kinetics of the InsP(3)R activated by InsP(3), AdA, and its analogues also identify molecular elements in InsP(3)R ligands that contribute to binding and activation of channel gating.     

Regio- and stereoselective cyclizations of dianhydro sugar alcohols catalyzed by a chiral (salen)Co(III) complex

The (salen)Co(III)OAc ((R,R)-1 and (S,S)-1) catalyzed cyclizations of the chiral dianhydro sugars, 1,2:5,6-dianhydro-3,4-di-O-methyl-D-glucitol (2), 1,2:5,6-dianhydro-3,4-di-O-methyl-D-mannitol (3), 1,2:5,6-dianhydro-3,4-di-O-methyl-L-iditol (4), and 1,2:4,5-dianhydro-3-O-methyl-L-arabinitol (5), is a facile method for the synthesis of anhydroalditol alcohols. Cyclization of 2 using (R,R)-1 and (S,S)-1 proceeded diastereoselectively to form 2,5-anhydro-3,4-di-O-methyl-D-mannitol (6) and 2,5-anhydro-3,4-di-O-methyl-L-iditol (7), respectively. The cyclization of 3 and 5 is a novel method for obtaining 1,6-anhydro-3,4-di-O-methyl-D-mannitol (11) and a stereoselective route to 1,5-anhydro-3-O-methyl-L-arabinitol (13). It is proposed that the reaction occurs via endo-selective cyclization of an epoxy alcohol produced by the endo-selective ring-opening of one of the two epoxide moieties in the starting material.     

Blazeispirane and protoblazeispirane derivatives from the cultured mycelia of the fungus Agaricus blazei

Two blazeispirane derivatives including blazeispirols G and I were isolated from the cultured mycelia of the fungus Agaricus blazei Murill and were established to be (20S, 22S, 23R, 24S)-14 beta,22: 22,25-diepoxy-5-methoxy-des-A-ergosta-5,7,9-triene-11 alpha,23-diol and (20S, 22S, 23R, 24S)-14 beta,22:22,25-diepoxy-5-methoxy-des-A-ergosta-5,7,9,11-tetraene-23,28-diol by comparison of extensive 1D and 2D NMR spectral data with that of blazeispirol A. Furthermore, four blazeispirol derivatives blazeispirols, U, V, V(1) and Z(1) were isolated form the same source described above. Their structures were determined to be (20S, 22S, 23R, 24S)-14 beta,22:22,25-diepoxy-23-hydroxyergosta-4,6,8,11-tetraen-3-one, (20S, 22S, 23R, 24S)-14 beta,22:22,25-diepoxy-6 alpha,7 alpha,23-trihydroxyergosta-4,8,11-trien-3-one, (20S, 22S, 23R, 24S)-14 beta,22:22,25-diepoxy-6 beta,7 alpha,23-trihydroxyergosta-4,8,11-trien-3-one and (20S, 22S, 23R, 24S)-14 beta,22:22,25-diepoxy-23-hydroxy-4,5-seco-ergosta-6,8-diene-3,5-dione by extensive 1 D and 2D NMR spectral data.     

Comparative proteomic analysis of drought tolerance in the two contrasting Tibetan wild genotypes and cultivated genotype

Background

Drought is one of major abiotic stresses constraining crop productivity worldwide. To adapt to drought stress, plants have evolved sophisticated defence mechanisms. Wild barley germplasm is a treasure trove of useful genes and offers rich sources of genetic variation for crop improvement. In this study, a proteome analysis was performed to identify the genetic resources and to understand the mechanisms of drought tolerance in plants that could result in high levels of tolerance to drought stress.

Results

A greenhouse pot experiment was performed to compare proteomic characteristics of two contrasting Tibetan wild barley genotypes (drought-tolerant XZ5 and drought-sensitive XZ54) and cv. ZAU3, in response to drought stress at soil moisture content 10 % (SMC10) and 4 % (SMC4) and subsequently 2 days (R1) and 5 days (R2) of recovery. More than 1700 protein spots were identified that are involved in each gel, wherein 132, 92, 86, 242 spots in XZ5 and 261, 137, 156, 187 in XZ54 from SMC10, SMC4, R1 and R2 samples were differentially expressed by drought over the control, respectively. Thirty-eight drought-tolerance-associated proteins were identified using mass spectrometry and data bank analysis. These proteins were categorized mainly into photosynthesis, stress response, metabolic process, energy and amino-acid biosynthesis. Among them, 6 protein spots were exclusively expressed or up-regulated under drought stress in XZ5 but not in XZ54, including melanoma-associated antigen p97, type I chlorophyll a/b-binding protein b, glutathione S-transferase 1, ribulosebisphosphate carboxylase large chain. Moreover, type I chlorophyll a/b-binding protein b was specifically expressed in XZ5 (Spots A4, B1 and C3) but not in both of XZ54 and ZAU3. These proteins may play crucial roles in drought-tolerance in XZ5. Coding Sequences (CDS) of rbcL and Trx-M genes from XZ5, XZ54 and ZAU3 were cloned and sequenced. CDS length of rbcL and Trx-M was 1401 bp (the partial-length CDS region) and 528 bp (full-length CDS region), respectively, encoding 467 and 176 amino acids. Comparison of gene sequences among XZ5, XZ54 and ZAU3 revealed 5 and 2 SNPs for rbcL and Trx-M, respectively, with two 2 SNPs of missense mutation in the both genes.

Conclusions

Our findings highlight the significance of specific-proteins associated with drought tolerance, and verified the potential value of Tibetan wild barley in improving drought tolerance of barley as well as other cereal crops.

Electronic supplementary material

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

Isoform-specific regulation of the inositol 1,4,5-trisphosphate receptor by O-linked glycosylation

The inositol 1,4,5-trisphosphate receptor (InsP(3)R), an intracellular calcium channel, has three isoforms with >65% sequence homology, yet the isoforms differ in their function and regulation by post-translational modifications. We showed previously that InsP(3)R-1 is functionally modified by O-linked β-N-acetylglucosamine glycosylation (O-GlcNAcylation) (Rengifo, J., Gibson, C. J., Winkler, E., Collin, T., and Ehrlich, B. E. (2007) J. Neurosci. 27, 13813-13821). We now report the effect of O-GlcNAcylation on InsP(3)R-2 and InsP(3)R-3. Analysis of AR4-2J cells, a rat pancreatoma cell line expressing predominantly InsP(3)R-2, showed no detectable O-GlcNAcylation of InsP(3)R-2 and no significant functional changes despite the presence of the enzymes for addition (O-β-N-acetylglucosaminyltransferase) and removal (O-β-N-acetylglucosaminidase) of the monosaccharide. In contrast, InsP(3)R-3 in Mz-ChA-1 cells, a human cholangiocarcinoma cell line expressing predominantly InsP(3)R-3, was functionally modified by O-GlcNAcylation. Interestingly, the functional impact of O-GlcNAcylation on the InsP(3)R-3 channel was opposite the effect measured with InsP(3)R-1. Addition of O-GlcNAc by O-β-N-acetylglucosaminyltransferase increased InsP(3)R-3 single channel open probability. Incubation of Mz-ChA-1 cells in hyperglycemic medium caused an increase in the InsP(3)-dependent calcium release from the endoplasmic reticulum. The dynamic and inducible nature of O-GlcNAcylation and the InsP(3)R isoform specificity suggest that this form of modification of InsP(3)R and subsequent changes in intracellular calcium transients are important in physiological and pathophysiological processes.     

Identification and characterization of two critical sequences in SV40PolyA that activate the green fluorescent protein reporter gene

Alu repeats or Line-1-ORF2 (ORF2) inhibit expression of the green fluorescent protein (GFP) gene when inserted downstream of this gene in the vector pEGFP-C1. In this work, we studied cis-acting elements that eliminated the repression of GFP gene expression induced by Alu and ORF2 and sequence characteristics of these elements. We found that sense and antisense PolyA of simian virus 40 (SV40PolyA, 240 bp) eliminated the repression of GFP gene expression when inserted between the GFP gene and the Alu (283 bp) repeats or ORF2 (3825 bp) in pAlu14 (14 tandem Alu repeats were inserted downstream of the GFP gene in the vector pEGFP-C1) or pORF2. Antisense SV40PolyA (PolyAas) induced stronger gene expression than its sense orientation (PolyA). Of four 60-bp segments of PolyAas (1F1R, 2F2R, 3F3R and 4F4R) inserted independently into pAlu14, only two (2F2R and 3F3R) eliminated the inhibition of GFP gene expression induced by Alu repeats. Deletion analysis revealed that a 17 nucleotide AT repeat (17ntAT; 5'-AAAAAAATGCTTTATTT-3') in 2F2R and the fragment 3F38d9 (5'-ATAAACAAGTTAACAACA ACAATTGCATT-3') in 3F3R were critical sequences for activating the GFP gene. Sequence and structural analyses showed that 17ntAT and 3F38d9 included imperfect palindromes and may form a variety of unstable stem-loops. We suggest that the presence of imperfect palindromes and unstable stem-loops in DNA enhancer elements plays an important role in GFP gene activation.     

Synthesis of anellated carbasugars from (--)-quinic acid

(3R,4R,5R)-3-[(tert-Butyl-dimethylsilyl)oxy]-4,5-(isopropylidenedioxy)-1-cyclohexanone (2) reacted with carbon disulfide and methyl iodide in the presence of sodium hydride to furnish (3R,4R,5R)-5-[(tert-butyl-dimethylsilyl)oxy]-3,4-(isopropylidenedioxy)-2-[bis(methylthio)methylene]-1-cyclohexanone (3). 2 and N,N-dimethylformamide dimethyl acetal afforded (2E,3R,4R,5R)-5-[(tert-butyl-dimethylsilyl)oxy]-2-(dimethylaminomethylene)-3,4-(isopropylidenedioxy)-1-cyclohexanone (4). These push-pull activated methylenecyclohexanones 3 and 4 underwent a ring closure reaction with hydrazine hydrate and methylhydrazine, respectively, to give pyrazoloanellated carbasugars. Treatment of 3 with formamidinium, acetamidinium and benzamidinium salts, respectively, in the presence of sodium methanolate yielded three (5R,6R,7R)-7-[(tert-butyl-dimethylsilyl)oxy]-5,6,7,8-tetrahydro-5,6-(isopropylidenedioxy)benzo[d]pyrimidines.     

Metagenomic analysis of the Rhinopithecus bieti fecal microbiome reveals a broad diversity of bacterial and glycoside hydrolase profiles related to lignocellulose degradation

Background

The animal gastrointestinal tract contains a complex community of microbes, whose composition ultimately reflects the co-evolution of microorganisms with their animal host and the diet adopted by the host. Although the importance of gut microbiota of humans has been well demonstrated, there is a paucity of research regarding non-human primates (NHPs), especially herbivorous NHPs.

Results

In this study, an analysis of 97,942 pyrosequencing reads generated from Rhinopithecus bieti fecal DNA extracts was performed to help better understanding of the microbial diversity and functional capacity of the R. bieti gut microbiome. The taxonomic analysis of the metagenomic reads indicated that R. bieti fecal microbiomes were dominated by Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria phyla. The comparative analysis of taxonomic classification revealed that the metagenome of R. bieti was characterized by an overrepresentation of bacteria of phylum Fibrobacteres and Spirochaetes as compared with other animals. Primary functional categories were associated mainly with protein, carbohydrates, amino acids, DNA and RNA metabolism, cofactors, cell wall and capsule and membrane transport. Comparing glycoside hydrolase profiles of R. bieti with those of other animal revealed that the R. bieti microbiome was most closely related to cow rumen.

Conclusions

Metagenomic and functional analysis demonstrated that R. bieti possesses a broad diversity of bacteria and numerous glycoside hydrolases responsible for lignocellulosic biomass degradation which might reflect the adaptations associated with a diet rich in fibrous matter. These results would contribute to the limited body of NHPs metagenome studies and provide a unique genetic resource of plant cell wall degrading microbial enzymes. However, future studies on the metagenome sequencing of R. bieti regarding the effects of age, genetics, diet and environment on the composition and activity of the metagenomes are required.

Electronic supplementary material

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