Rapid phylogenetic and functional classification of short genomic fragments with signature peptides

ABSTRACT: BACKGROUND: Classification is difficult for shotgun metagenomics data from environments such as soils, where the diversity of sequences is high and where reference sequences from close relatives may not exist. Approaches based on sequence-similarity scores must deal with the confounding effects that inheritance and functional pressures exert on the relation between scores and phylogenetic distance, while approaches based on sequence alignment and tree-building are typically limited to a small fraction of gene families. We describe an approach based on finding one or more exact matches between a read and a precomputed set of peptide 10-mers. RESULTS: At even the largest phylogenetic distances, thousands of 10-mer peptide exact matches can be found between pairs of bacterial genomes. Genes that share one or more peptide 10-mers typically have high reciprocal BLAST scores. Among a set of 403 representative bacterial genomes, some 20 million 10-mer peptides were found to be shared. We assign each of these peptides as a signature of a particular node in a phylogenetic reference tree based on the RNA polymerase genes. We classify the phylogeny of a genomic fragment (e.g., read) at the most specific node on the reference tree that is consistent with the phylogeny of observed signature peptides it contains. Using both synthetic data from four newly-sequenced soil-bacterium genomes and ten real soil metagenomics data sets, we demonstrate a sensitivity and specificity comparable to that of the MEGAN metagenomics analysis package using BLASTX against the NR database. Phylogenetic and functional similarity metrics applied to real metagenomics data indicates a signal-to-noise ratio of approximately 400 for distinguishing among environments. Our method assigns ~6.6 Gbp/hr on a single CPU, compared with 25 kbp/hr for methods based on BLASTX against the NR database. CONCLUSIONS: Classification by exact matching against a precomputed list of signature peptides provides comparable results to existing techniques for reads longer than about 300 bp and does not degrade severely with shorter reads. Orders of magnitude faster than existing methods, the approach is suitable now for inclusion in analysis pipelines and appears to be extensible in several different directions.     

Blood groups and protein polymorphisms in five goat breeds (Capra hircus)

Data on allele frequencies at six red cell blood group systems and three blood protein polymorphic loci in five goat breeds are reported. Two blood proteins, albumin and carbonic anhy-drase, were not found to be polymorphic. The B blood group system of goats, like its homologue in cattle and sheep, is highly complex. At least 44 B phenogroups (haplotypes) have been distinguished in this study. Based on the variation in allele frequencies between breeds, genetic distances were calculated. The distances estimated by four different methods were in close agreement with data from the history and geographic origins of the breeds examined.     

Virus inhibition activity of effector memory CD8(+) T cells determines simian immunodeficiency virus load in vaccinated monkeys after vaccine breakthrough infection

The goal of an effective AIDS vaccine is to generate immunity that will prevent human immunodeficiency virus 1 (HIV-1) acquisition. Despite limited progress toward this goal, renewed optimism has followed the recent success of the RV144 vaccine trial in Thailand. However, the lack of complete protection in this trial suggests that breakthroughs, where infection occurs despite adequate vaccination, will be a reality for many vaccine candidates. We previously reported that neutralizing antibodies elicited by DNA prime-recombinant adenovirus serotype 5 (rAd5) boost vaccination with simian immunodeficiency virus strain mac239 (SIVmac239) Gag-Pol and Env provided protection against pathogenic SIVsmE660 acquisition after repeated mucosal challenge. Here, we report that SIV-specific CD8(+) T cells elicited by that vaccine lowered both peak and set-point viral loads in macaques that became infected despite vaccination. These SIV-specific CD8(+) T cells showed strong virus-inhibitory activity (VIA) and displayed an effector memory (EM) phenotype. VIA correlated with high levels of CD107a mobilization and perforin expression in SIV-specific CD8(+) T cells. Remarkably, both the frequency and the number of Gag CM9-specific public clonotypes were strongly correlated with VIA mediated by EM CD8(+) T cells. The ability to elicit such virus-specific EM CD8(+) T cells might contribute substantially to an efficacious HIV/AIDS vaccine, even after breakthrough infection.     

Overlapping epitopes in human immunodeficiency virus type 1 gp120 presented by HLA A, B, and C molecules: effects of viral variation on cytotoxic T-lymphocyte recognition.

Human immunodeficiency virus (HIV)-specific cytotoxic T lymphocytes (CTL) are thought to exert immunologic selection pressure in infected persons, yet few data regarding the effects of this constraint on viral sequence variation in vivo, particularly in the highly variable Env protein, are available. In this study, CD8+ HIV type 1 (HIV-1) envelope-specific CTL clones specific for gp120 were isolated from peripheral blood mononuclear cells of four HIV-infected individuals, all of which recognized the same 25-amino-acid (aa) peptide (aa 371 to 395), which is partially contained in the CD4-binding domain of HIV-1 gp120. Fine mapping studies revealed that two of the clones optimally recognized the 9-aa sequence 375 to 383 (SFNCGGEFF), while the two other clones optimally recognized the epitope contained in the overlapping 9-aa sequence 376 to 384 (FNCGGEFFY). Lysis of target cells by the two clones recognizing aa 375 to 383 was restricted by HLA B15 and Cw4, respectively, whereas both clones recognizing aa 376 to 384 were restricted by HLA A29. Sequence variation, relative to the IIIB strain sequence used to identify CTL clones, was observed in autologous viruses in the epitope-containing region in all four subjects. However, poorly recognized autologous sequence variants were predominantly seen for the A29-restricted clones, whereas the clones specific for SFNCGGEFF continued to recognize the predominant autologous sequences. These results suggest that the HLA profile of an individual may not only be important in determining the specificity of CTL recognition but may also affect the ability to recognize virus variants and suppress escape from CTL recognition. These results also identify overlapping viral CTL epitopes which can be presented by HLA A, B, and C molecules.