Deriving enzymatic and taxonomic signatures of metagenomes from short read data

Background  

We propose a method for deriving enzymatic signatures from short read metagenomic data of unknown species. The short read data are converted to six pseudo-peptide candidates. We search for occurrences of Specific Peptides (SPs) on the latter. SPs are peptides that are indicative of enzymatic function as defined by the Enzyme Commission (EC) nomenclature. The number of SP hits on an ensemble of short reads is counted and then converted to estimates of numbers of enzymatic genes associated with different EC categories in the studied metagenome. Relative amounts of different EC categories define the enzymatic spectrum, without the need to perform genomic assemblies of short reads.     

Microbial degradation of alkylbenzenesulphonates. Metabolism of homologues of short alkyl-chain length by an Alcaligenes sp

1. An organism isolated from sewage and identified as an Alcaligenes sp. utilized benzenesulphonate, toluene-p-sulphonate or phenylethane-p-sulphonate as sole source of carbon and energy for growth. Higher alkylbenzenesulphonate homologues and the hydrocarbons, benzene, toluene, phenylethane and 1-phenyldodecane were not utilized. 2. 2-Phenylpropanesulphonate was metabolized to 4-isopropylcatechol. 3. 1-Phenylpropanesulphonate was metabolized to an ortho-diol, which was tentatively identified, in the absence of an authentic specimen, as 4-n-propylcatechol. 4. In the presence of 4-isopropylcatechol, which inhibited catechol 2,3-dioxygenase, 4-ethylcatechol accumulated in cultures growing on phenylethane-p-sulphonate. 5. Authentic samples of catechol, 3-methylcatechol, 4-methylcatechol, 4-ethylcatechol and 3-isopropylcatechol were oxidized by heat-treated extracts to the corresponding 2-hydroxyalkylmuconic semialdehydes. Ring cleavage occurred between C-2 and C-3. 6. The catechol derived from 1-phenylpropanesulphonate was oxygenated by catechol 2,3-dioxygenase to a compound with all the properties of a 2-hydroxyalkylmuconic semialdehyde, but it was not rigorously identified. 7. The catechol 2,3-dioxygenase induced by growth on benzenesulphonate, toluene-p-sulphonate or phenylethane-p-sulphonate showed a constant ratio of specific activities with catechol, 3-methylcatechol, 4-methylcatechol and 4-ethylcatechol that was independent of the growth substrate. At 60°C, activity towards these substrates declined at an identical first-order rate. 8. Enzymes of the `ortho' pathway of catechol metabolism were present in small amounts in cells grown on benzenesulphonate, toluene-p-sulphonate or phenylethane-p-sulphonate. 9. The catechol 1,2-dioxygenase oxidized the alkylcatechols, but the rates and the total extents of oxidation were less than for catechol itself. The oxidation products of these alkylcatechols were not further metabolized.     

Synthesis and hydrolysis by cysteine and serine proteases of short internally quenched fluorogenic peptides

We developed sensitive substrates for cysteine proteases and specific substrates for serine proteases based on short internally quenched fluorescent peptides, Abz-F-R-X-EDDnp, where Abz (ortho-aminobenzoic acid) is the fluorescent donor, EDDnp [N-(ethylenediamine)-2,4-dinitrophenyl amide] is the fluorescent quencher, and X are natural amino acids. This series of peptides is compared to the commercially available Z-F-R-MCA, where Abz and X replace carbobenzoxy (Z) and methyl-7-aminocoumarin amide (MCA), respectively; and EDDnp can be considered a P(2)' residue. Whereas MCA is the fluorescent probe and cannot be modified, in the series Abz-F-R-X-EDDnp the amino acids X give the choice of matching the specificity of the S(1)' enzyme subsite, increasing the substrate specificity for a particular protease. All Abz-F-R-X-EDDnp synthesized peptides (for X = Phe, Leu, Ile, Ala, Pro, Gln, Ser, Lys, and Arg) were assayed with papain, human cathepsin L and B, trypsin, human plasma, and tissue kallikrein. Abz-F-R-L-EDDnp was the best substrate for papain and Abz-F-R-R-EDDnp or Abz-F-R-A-EDDnp was the more susceptible to cathepsin L. Abz-F-R-L-EDDnp was able to detect papain in the range of 1 to 15 pM. Human plasma kallikrein hydrolyzed Abz-F-R-R-EDDnp with significant efficiency (k(cat)/K(m) = 1833 mM(-1) s(-1)) and tissue kallikrein was very selective, hydrolyzing only the peptides Abz-F-R-A-EDDnp (k(cat)/K(m) = 2852 mM(-1) s(-1)) and Abz-F-R-S-EDDnp (k(cat)/K(m) = 4643 mM(-1) s(-1)). All Abz-F-R-X-EDDnp peptides were resistant to hydrolysis by thrombin and activated factor X.     

A new approach to stabilize complementary complexes formed by short oligonucleotides and their derivatives.

A proposed approach to stabilize complementary complexes is based on the stabilizing effect of oligonucleotide derivatives bearing polyaromatic residues on the adjacent short oligonucleotides if their recognition sites make a tandem on the target. This approach has been used to enhance efficiency and specificity of site-specific modification of nucleic acids.     

Probing the initiation complex formation on E coli ribosomes using short complementary DNA oligomers.

Interactions between Escherichia coli 16S rRNA sequences (as components of 30S ribosomal subunits or tight-couple 70S ribosomes) with the ligands poly(U), poly(AGU), tRNAPhe, tRNAfMet, and the initiation factors have been studied. The ligands were employed as competitors for selected sites on 16S rRNA known to be accessible for hybridization to cDNA oligomers, regions 517-528, 1397-1404, and 1534-1542. The binding of cDNAs 1534-1541 and 1398-1403 decreased in the presence of the ligand pair poly(U)/tRNAPhe. Only the binding of cDNA 1534-1541 was affected by poly(AGU), while none of the complementary DNA oligomer binding was affected by tRNAPhe or tRNAfMet alone. The poly(AGU)/tRNAfMet ligand pair caused an additional decline in the binding of cDNA 1534-1541, relative to that caused by poly(AGU) alone, but the ligand pair did not affect the binding of the cDNA oligomers 517-528 or 1398-1403. The inclusion of the initiation factors did not significantly alter the binding level decreases observed for cDNA 1534-1541 in the presence of mRNAs or tRNA. At the 517-528 and 1398-1403 regions, the inclusion of the initiation factors, in either the presence or absence of the other ligands, caused a large decrease in the binding of the cDNA oligomers. The oligomers complementary to 16S bases 517-528 and 1398-1403 did not bind to tight-couple or reassociated 70S ribosomes. The data are discussed in terms of the decoding site hypothesis, and in terms of the mRNA alignment mechanism proposed by Trifonov [1].     

Simian homologues of Epstein-Barr virus

Gamma-herpesviruses closely related to the Epstein-Barr virus (EBV) are known to naturally infect Old World non-human primates and are classified in the same lymphocryptovirus (LCV) genera. LCV infecting humans and Old World primates share similar biology, and recent studies have demonstrated that these viruses share a similar repertoire of viral genes. Surprisingly, the latent infection genes associated with cell growth transformation demonstrate the most striking sequence divergence, but the functional mechanisms for these genes are generally well conserved. The recent discovery of LCVs naturally infecting New World primates has rewritten the old paradigm of LCV host range restriction to humans and Old World non-human primates, so that these viruses are more widespread than previously believed. However, the New World LCV genome has significant and interesting differences from EBV and other Old World LCVs despite similar biological properties. Thus, the simian homologues of EBV can provide an important animal model for studying LCV pathogenesis, and the similarities and differences that have evolved among these related viruses can provide a unique perspective towards a better understanding of EBV.     

Novel simian homologues of Epstein-Barr virus

Thirty different lymphocryptoviruses (LCV), 26 of them novel, were detected in primates by a panherpesvirus PCR assay. Nineteen LCV from chimpanzees, bonobos, gorillas, and other Old World primates were closely related to Epstein-Barr virus (EBV), the type species of the genus lymphocryptovirus. Seven LCV originating from New World primates were related to callitrichine herpesvirus 3 (CalHV-3), the first recognized New World LCV. Importantly, a second LCV from gorillas and three LCV from orangutans and gibbons were only distantly related to EBV and CalHV-3. They were tentatively assigned to a novel genogroup of Old World primate LCV. The work described in the paper may also help identify an as yet unknown human LCV.     

Simian homologues of human herpesvirus 8

Gamma-herpesviruses can be found in most primates including Old World an New World monkeys. The gamma-herpesvirinae are grouped into two classes: lymphocryptoviruses (gamma1) and rhadinoviruses (gamma2). The lymphocryptoviruses include Epstein-Barr virus, lymphocryptovirus of rhesus monkeys, and Herpesvirus papio of baboons. Rhadinoviruses that infect New World monkeys include Herpesvirus saimiri, whose natural host is the squirrel monkey, and Herpesvirus ateles, which infects spider monkeys. Rhadinoviruses that infect hominoids and Old World monkeys include Kaposi's sarcoma-associated herpesvirus, also known as HHV-8, and rhesus monkey rhadinovirus.     

Semisynthetic engineering of proteinase inhibitor homologues

A semisynthetic approach to modulate the inhibitory specificity of aprotinin, the Kunitz trypsin inhibitor from bovine mast cells, is described. By the use of peptide-chemical procedures a single amino acid of its reactive site can be replaced by any other coded or non-coded amino acid. Thus, a series of aprotinin homologues have been prepared which demonstrate the individual contribution of a single side chain to the inhibition of a particular target proteinase and enable specific inhibitors to be designed.     

Cyclin-dependent kinase homologues of Plasmodium falciparum

The intraerythrocytic asexual cycle of the malarial parasite is complex and atypical: during schizogony the parasite undergoes multiple rounds of DNA replication and asynchronous nuclear division without cytokinesis. This cell cycle deviates from the classical eukaryotic cell cycle model where, 'DNA replicates only once per cell cycle'. A clear understanding of the molecular switches that control this unusual developmental cycle would be of great interest, both in terms of fundamental Plasmodium biology and in terms of novel potential drug target identification. In recent years considerable effort has been made to identify the malarial orthologues of the cyclin-dependent kinases, which are key regulators of the orderly progression of the eukaryotic cell cycle. This review focuses on the current state-of-knowledge of Plasmodium falciparum cyclin-dependent kinase-like kinases and their regulators.     

Current understanding of mammalian TRP homologues

Calcium influx into the cell from the extracellular medium is crucial for important processes including muscle contraction, secretion and gene expression. This calcium influx is mainly mediated through calcium influx channels, which on the basis of their activation mechanism can be subdivided in voltage-gated calcium channels, which have already been thoroughly characterized and non-voltage-gated calcium permeable channels. This latter group includes ion channels activated by binding of extra and intracellular messengers, mechanical stress or depletion of intracellular calcium stores. Currently little molecular data is available concerning this class of calcium influx channels. However, recent studies have indicated that members of the transient receptor potential (TRP) family of ion channels can function as calcium influx channels both in excitable and non-excitable tissues. On the basis of structural information the TRP family is subdivided in three main subfamilies: the TRPC (canonical) group, the TRPV (vanilloid) group and the TRPM (melastatin) group. The cloning and characterization of members of this cation channel family has exploded during recent years, leading to a plethora of data concerning TRPs in a variety of tissues and species, including mammals, insects and yeast. This review summarizes the currently available information concerning members of the TRP family expressed in mammalian tissues.     

Phylogenetic analysis of archaeal PCNA homologues

Proliferating cell nuclear antigen (PCNA) is an essential component of the DNA replication and repair machinery in the domain Eucarya. Eukaryotes and euryarchaeotes, which belong to one subdomain of Archaea, possess a single PCNA homologue, whereas two distinct PCNA homologues have been identified from Sulfolobus solfataricus, which belongs to the other archaeal subdomain, Crenarchaeota. We have cloned and sequenced two genes of PCNA homologues from the thermoacidophilic crenarchaeon Sulfurisphaera ohwakuensis. These genes, referred to as the Soh PCNA A gene and the Soh PCNA B gene, were found to encode 245 amino acids (aa) (27 kDa) and 248 aa (27 kDa), respectively. In deduced amino acid sequences of both PCNA homologues, the motif L/I-A-P-K/R, implicated in binding of PCNA with replication factor C (RFC), was identified. Phylogenetic analysis of all available archaeal PCNA homologues suggests that crenarchaeal homologues are divided into two groups. Group A consists of Soh PCNA A, one of the S. solfataricus PCNA homologues, and one of the Aeropyrum pernix PCNA homologues. The other crenarchaeal homologues form group B. Crenarchaeal PCNA homologues constitute a monophyletic subfamily. These results suggest that the evolution of crenarchaeal PCNA homologues has been characterized by one or two gene duplication events, which are assumed to have occurred after the split of the crenarchaeal and euryarchaeal lineages. Received: July 10, 2000 / Accepted: September 26, 2000     

Furanosidase superfamily: search of homologues

The furanosidase superfamily contains GH32, GH43, GH62, GH68, GH117, DUF377, and DUF1861 families of glycoside hydrolases and their homologues. Catalytic domains of these families have five-bladed beta-propeller tertiary structure. Iterative screening of the protein database allowed to support their relationship as well as evolutionary connections with domains from GH33 and GH93 families of glycoside hydrolases. The latter two have structure of the six-bladed beta-propeller. Among revealed homologues we found 441 unclassified proteins. These proteins are combined into 39 groups based on homology: FURAN1-FURAN39. FURAN8 and FURAN36 can be considered as separate subfamilies within GH43 and GH32 families of glycoside hydrolases, respectively. The remaining 37 groups are new families of hypothetical glycoside hydrolases.     

Furanosidase superfamily: Search of homologues

The furanosidase superfamily contains the GH32, GH43, GH62, GH68, GH117, DUF377 (GH130), and DUF1861 families of glycoside hydrolases and their homologues. Catalytic domains of these families have five-bladed β-propeller tertiary structure. Iterative screening of the protein database supports of their relationship as well as evolutionary connections with domains from GH33 and GH93 families of glycoside hydrolases. The latter two have the structure of the six-bladed β-propeller. Among detected homologues we found 441 unclassified proteins. These proteins are combined into 39 groups based on homology: FURAN1-FURAN39. FURAN8 and FURAN36 can be considered as separate subfamilies within the GH43 and GH32 families of glycoside hydrolases, respectively. The remaining 37 groups are new families of hypothetical glycoside hydrolases.     

Database mining of plant peptide homologues

In plant-pathogen interactions, pathogens employ secreted molecules, known as effectors to overcome physical barriers, modulate plant immunity, and facilitate colonization. Among these diverse effectors, some are found to mimic the plant peptides, to target host’s peptide receptors, and intervene in the peptide-regulated defense pathways and/or plant development. To better understand how pathogens have co-evolved with their plant hosts in order to improve disease management, we explored the presence of plant peptide mimics in microbes by bioinformatic analysis. In total, 36 novel peptide mimics belong to five plant peptide families were detected in bacterial and fungal kingdoms. Among them, phytosulfokine homologues were widely distributed in 22 phytopathogens and one bacterium, thereby constituted the largest proportion of the identified mimics. The putative functional peptide region is well conserved between plant and microbes, while the existence of a putative signal peptide varies between species. Our findings will increase understanding of plant-pathogen interactions, and provide new ideas for future studies of pathogenic mechanisms and disease management.     

Biomarker signatures of aging

Because people age differently, age is not a sufficient marker of susceptibility to disabilities, morbidities, and mortality. We measured nineteen blood biomarkers that include constituents of standard hematological measures, lipid biomarkers, and markers of inflammation and frailty in 4704 participants of the Long Life Family Study (LLFS), age range 30–110 years, and used an agglomerative algorithm to group LLFS participants into clusters thus yielding 26 different biomarker signatures. To test whether these signatures were associated with differences in biological aging, we correlated them with longitudinal changes in physiological functions and incident risk of cancer, cardiovascular disease, type 2 diabetes, and mortality using longitudinal data collected in the LLFS. Signature 2 was associated with significantly lower mortality, morbidity, and better physical function relative to the most common biomarker signature in LLFS, while nine other signatures were associated with less successful aging, characterized by higher risks for frailty, morbidity, and mortality. The predictive values of seven signatures were replicated in an independent data set from the Framingham Heart Study with comparable significant effects, and an additional three signatures showed consistent effects. This analysis shows that various biomarker signatures exist, and their significant associations with physical function, morbidity, and mortality suggest that these patterns represent differences in biological aging. The signatures show that dysregulation of a single biomarker can change with patterns of other biomarkers, and age‐related changes of individual biomarkers alone do not necessarily indicate disease or functional decline.     

Distinctive signatures of recursion

Although recursion has been hypothesized to be a necessary capacity for the evolution of language, the multiplicity of definitions being used has undermined the broader interpretation of empirical results. I propose that only a definition focused on representational abilities allows the prediction of specific behavioural traits that enable us to distinguish recursion from non-recursive iteration and from hierarchical embedding: only subjects able to represent recursion, i.e. to represent different hierarchical dependencies (related by parenthood) with the same set of rules, are able to generalize and produce new levels of embedding beyond those specified a priori (in the algorithm or in the input). The ability to use such representations may be advantageous in several domains: action sequencing, problem-solving, spatial navigation, social navigation and for the emergence of conventionalized communication systems. The ability to represent contiguous hierarchical levels with the same rules may lead subjects to expect unknown levels and constituents to behave similarly, and this prior knowledge may bias learning positively. Finally, a new paradigm to test for recursion is presented. Preliminary results suggest that the ability to represent recursion in the spatial domain recruits both visual and verbal resources. Implications regarding language evolution are discussed.     

Interaction of the TNF homologues BLyS and APRIL with the TNF receptor homologues BCMA and TACI

BLyS (also called TALL-1, THANK, or BAFF) [1] [2] [3] [4] is a member of the tumor necrosis factor (TNF) gene family that stimulates proliferation and immunoglobulin production by B cells. BLyS interacts with the TNF receptor (TNFR) homologue TACI (transmembrane activator and CAML-interactor) [5], and treatment of mice with a TACI-Fc fusion protein abolishes germinal center formation after antigenic challenge [6]. Here we report a novel interaction between BLyS and another TNFR homologue, BCMA (B cell maturation antigen) [7] [8]. Further, the TNF homologue APRIL [9], a close relative of BLyS, also bound to BCMA and TACI. BLyS or APRIL activated nuclear factor-kappaB (NF-kappaB) through TACI and BCMA, and each ligand stimulated immunoglobulin M (IgM) production by peripheral blood B cells. These results define a dual ligand-receptor system that may play an important role in humoral immunity.