Immunity in Chimpanzees Chronically Infected with Hepatitis C Virus: Role of Minor Quasispecies in Reinfection

We have previously reported that chimpanzees chronically infected with hepatitis C virus (HCV) could be reinfected, even with the original infecting strain. In this study we tested the hypothesis that this might reflect the presence of minor quasispecies to which there was little or no immunity. To evaluate this hypothesis, we sequenced multiple clones taken at intervals after primary infection and rechallenge from four chronically infected chimpanzees. The inoculum used in these studies (HCV-H, genotype 1a) revealed 17 separate variants among 46 clones sequenced. Following challenge, each of the four challenged animals showed marked alterations of their quasispecies distribution. The new variants, which appeared 1 to 6 weeks after challenge, were either identical to or closely resembled variants present in the challenge inoculum. These results, paralleled by an increase in viremia in some of the challenged animals, suggest that quasispecies in the challenge inoculum were responsible for signs of reinfection and that there was little immunity. However, the newly emerged quasispecies completely took over infection in only one animal. In the remaining three chimpanzees the prechallenge quasispecies were able to persist. The natural evolution of infection within chimpanzees resulted in variants able to compete with the inoculum variants. Whether through reexposure or the natural progression of infection, newly emerged quasispecies are likely to play a role in the pathogenesis of chronic HCV infection.Hepatitis C virus (HCV) is estimated to chronically infect about 400 million people worldwide. More than half of these develop chronic active hepatitis, cirrhosis, or hepatocellular carcinoma. The HCV genome consists of a single-stranded RNA molecule approximately 10 kb long which contains a single open reading frame encoding approximately 3,000 amino acids (1, 5). There are at least six genotypes of HCV, and within a given patient the genomes are distributed among quasispecies which show sequence variation, particularly in the variable regions of the genome (4, 9). Hypervariable region 1 (HVR1) is a 27-amino-acid segment in the amino terminus of the second envelope protein which has been identified as the most variable region of the viral genome (11, 20). Sequential changes have been observed during the course of chronic HCV infections in chimpanzees and in humans (4, 11, 12). It has been postulated that these reflect immune system selection of neutralizing epitopes encoded by HVR1 (18, 19) and that persistent infection depends on the ability of the virus to continually evade the effects of neutralizing antibody (7, 10, 15, 17, 20). Due to its variability, HVR1 has been used extensively as an indicator of viral evolution.We have previously reported that chronically infected chimpanzees could seemingly be reinfected, even with the original infecting strain (13). In a recent report a similar phenomenon was observed in patients with posttransfusion hepatitis (6). We postulated that this might reflect the presence of minor quasispecies in the inoculum to which there was little or no immunity (13). Here we test this hypothesis by sequencing multiple clones of HVR1 derived at intervals after initial infection and after rechallenge.     

A Comprehensive Panel of Near-Full-Length Clones and Reference Sequences for Non-Subtype B Isolates of Human Immunodeficiency Virus Type 1

Non-subtype B viruses cause the vast majority of new human immunodeficiency virus type 1 (HIV-1) infections worldwide and are thus the major focus of international vaccine efforts. Although their geographic dissemination is carefully monitored, their immunogenic and biological properties remain largely unknown, in part because well-characterized virological reference reagents are lacking. In particular, full-length clones and sequences are rare, since subtype classification is frequently based on small PCR-derived viral fragments. There are only five proviral clones available for viruses other than subtype B, and these represent only 3 of the 10 proposed (group M) sequence subtypes. This lack of reference sequences also confounds the identification and analysis of mosaic (recombinant) genomes, which appear to be arising with increasing frequency in areas where multiple sequence subtypes cocirculate. To generate a more representative panel of non-subtype B reference reagents, we have cloned (by long PCR or lambda phage techniques) and sequenced 10 near-full-length HIV-1 genomes (lacking less than 80 bp of long terminal repeat sequences) from primary isolates collected at major epicenters of the global AIDS pandemic. Detailed phylogenetic analyses identified six that represented nonrecombinant members of HIV-1 subtypes A (92UG037.1), C (92BR025.8), D (84ZR085.1 and 94UG114.1), F (93BR020.1), and H (90CF056.1), the last two comprising the first full-length examples of these subtypes. Four others were found to be complex mosaics of subtypes A and C (92RW009.6), A and G (92NG083.2 and 92NG003.1), and B and F (93BR029.4), again emphasizing the impact of intersubtype recombination on global HIV-1 diversification. Although a number of clones had frameshift mutations or translational stop codons in major open reading frames, all the genomes contained a complete set of genes and three had intact genomic organizations without inactivating mutations. Reconstruction of one of these (94UG114.1) yielded replication-competent virus that grew to high titers in normal donor peripheral blood mononuclear cell cultures. This panel of non-subtype B reference genomes should prove valuable for structure-function studies of genetically diverse viral gene products, the generation of subtype-specific immunological reagents, and the production of DNA- and protein-based subunit vaccines directed against a broader spectrum of viruses.One critical question facing current AIDS vaccine development efforts is to what extent human immunodeficiency virus type 1 (HIV-1) genetic variation has to be considered in the design of candidate vaccines (11, 21, 41, 72). Phylogenetic analyses of globally circulating viral strains have identified two distinct groups of HIV-1 (M and O) (33, 45, 61, 62), and 10 sequence subtypes (A to J) have been proposed within the major group (M) (29, 30, 45, 72). Sequence variation among viruses belonging to these different lineages is extensive, with envelope amino acid sequence variation ranging from 24% between different subtypes to 47% between the two different groups. Given this extent of diversity, the question has been raised whether immunogens based on a single virus strain can be expected to elicit immune responses effective against a broad spectrum of viruses or whether vaccine preparations should include mixtures of genetically divergent antigens and/or be tailored toward locally circulating strains (11, 21, 41, 72). This is of particular concern in developing countries, where multiple subtypes of HIV-1 are known to cocirculate and where subtype B viruses (which have been the source of most current candidate vaccine preparations [10, 21]) are rare or nonexistent (5, 24, 40, 72).Although the extent of global HIV-1 variation is well defined, little is known about the biological consequences of this genetic diversity and its impact on cellular and humoral immune responses in the infected host. In particular, it remains unknown whether subtype-specific differences in virus biology exist that have to be considered for vaccine design. Thus far, such differences have not been identified. For example, several studies have shown that there is no correlation between HIV-1 genetic subtypes and neutralization serotypes (38, 42, 46, 68). Some viruses are readily neutralized, while most are relatively neutralization resistant (42). Although the reasons for these different susceptibilities remain unknown, it is clear that neutralization is not a function of the viral genotype (38, 42, 46, 68). Similarly, recent studies have identified vigorous cross-clade cytotoxic T-lymphocyte (CTL) reactivities in individuals infected with viruses from several different clades (3, 6), as well as in recipients of a clade B vaccine (15). These results are very encouraging, since they suggest that CTL cross-recognition among HIV-1 clades is much more prevalent than previously anticipated and that immunogens based on a limited number of variants may be able to elicit a broad CTL response (6). Nevertheless, it would be premature to conclude that HIV-1 variation poses no problem for AIDS vaccine design. Only a comprehensive analysis of genetically defined representatives of the various groups and subtypes will allow us to judge whether certain variants differ in fundamental viral properties and whether such differences will have to be incorporated into vaccine strategies. Obviously, such studies require well-characterized reference reagents, in particular full-length and replication-competent molecular clones that can be used for functional and biological studies.Full-length reference sequences representing the various subtypes are also urgently needed for phylogenetic comparisons. Recent analyses of subgenomic (23, 52, 54, 58) as well as full-length (7, 18, 53, 60) HIV-1 sequences identified a surprising number of HIV-1 strains which clustered in different subtypes in different parts of their genome. All of these originated from geographic regions where multiple subtypes cocirculated and are the results of coinfections with highly divergent viruses (52, 60, 62). Detailed phylogenetic characterization revealed that most of them have a complex genome structure with multiple points of crossover (7, 18, 53, 60). Some recombinants, like the “subtype E” viruses, which are in fact A/E recombinants (7, 18), have a widespread geographic dissemination and are responsible for much of the Asian HIV-1 epidemic (69, 70). In other areas, recombinants appear to be generated with increasing frequencies since many randomly chosen isolates exhibit evidence of mosaicism (4, 8, 31, 66, 71). Since recombination provides the opportunity for evolutionary leaps with genetic consequences that are far greater than those of the steady accumulation of individual mutations, the impact of recombination on viral properties must be monitored. We therefore need full-length nonrecombinant reference sequences for all major HIV-1 groups and subtypes before we can map and characterize the extent of intersubtype recombination.The number of molecular reagents for non-subtype B viruses is very limited. There are currently only five full-length, nonrecombinant molecular clones available for viruses other than subtype B (45), and these represent only three of the proposed (group M) subtypes (A, C, and D). Moreover, only three clones (all derived from subtype D viruses) are replication competent and thus useful for studies requiring functional gene products (45, 48, 65). Given the unknown impact of genetic variation on correlates of immune protection, subtype-specific reagents are critically needed for phylogenetic, immunological, and biological studies. In this paper, we report the cloning (by long PCR and lambda techniques) of 10 near-full-length HIV-1 genomes from isolates previously classified as non-subtype B viruses. Detailed phylogenetic analysis showed that six comprise nonmosaic representatives of five major subtypes, including two for which full-length representatives have not been reported. Four others were identified as complex intersubtype recombinants, again emphasizing the prevalence of hybrid genomes among globally circulating HIV-1 strains. We also describe a strategy for the biological evaluation of long-PCR-derived genomes and report the generation of a replication-competent provirus by this approach. The effect of these reagents on vaccine development is discussed.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

N-Terminal-oriented Proteogenomics of the Marine Bacterium Roseobacter Denitrificans Och114 using N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP) Labeling and Diagonal Chromatography

Given the ease of whole genome sequencing with next-generation sequencers, structural and functional gene annotation is now purely based on automated prediction. However, errors in gene structure are frequent, the correct determination of start codons being one of the main concerns. Here, we combine protein N termini derivatization using (N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP Ac-OSu) as a labeling reagent with the COmbined FRActional DIagonal Chromatography (COFRADIC) sorting method to enrich labeled N-terminal peptides for mass spectrometry detection. Protein digestion was performed in parallel with three proteases to obtain a reliable automatic validation of protein N termini. The analysis of these N-terminal enriched fractions by high-resolution tandem mass spectrometry allowed the annotation refinement of 534 proteins of the model marine bacterium Roseobacter denitrificans OCh114. This study is especially efficient regarding mass spectrometry analytical time. From the 534 validated N termini, 480 confirmed existing gene annotations, 41 highlighted erroneous start codon annotations, five revealed totally new mis-annotated genes; the mass spectrometry data also suggested the existence of multiple start sites for eight different genes, a result that challenges the current view of protein translation initiation. Finally, we identified several proteins for which classical genome homology-driven annotation was inconsistent, questioning the validity of automatic annotation pipelines and emphasizing the need for complementary proteomic data. All data have been deposited to the ProteomeXchange with identifier PXD000337.Recent developments in mass spectrometry and bioinformatics have established proteomics as a common and powerful technique for identifying and quantifying proteins at a very broad scale, but also for characterizing their post-translational modifications and interaction networks (1, 2). In addition to the avalanche of proteomic data currently being reported, many genome sequences are established using next-generation sequencing, fostering proteomic investigations of new cellular models. Proteogenomics is a relatively recent field in which high-throughput proteomic data is used to verify coding regions within model genomes to refine the annotation of their sequences (2–8). Because genome annotation is now fully automated, the need for accurate annotation for model organisms with experimental data is crucial. Many projects related to genome re-annotation of microorganisms with the help of proteomics have been recently reported, such as for Mycoplasma pneumoniae (9), Rhodopseudomonas palustris (10), Shewanella oneidensis (11), Thermococcus gammatolerans (12), Deinococcus deserti (13), Salmonella thyphimurium (14), Mycobacterium tuberculosis (15, 16), Shigella flexneri (17), Ruegeria pomeroyi (18), and Candida glabrata (19), as well as for higher organisms such as Anopheles gambiae (20) and Arabidopsis thaliana (4, 5).The most frequently reported problem in automatic annotation systems is the correct identification of the translational start codon (21–23). The error rate depends on the primary annotation system, but also on the organism, as reported for Halobacterium salinarum and Natromonas pharaonis (24), Deinococcus deserti (21), and Ruegeria pomeroyi (18), where the error rate is estimated above 10%. Identification of a correct translational start site is essential for the genetic and biochemical analysis of a protein because errors can seriously impact subsequent biological studies. If the N terminus is not correctly identified, the protein will be considered in either a truncated or extended form, leading to errors in bioinformatic analyses (e.g. during the prediction of its molecular weight, isoelectric point, cellular localization) and major difficulties during its experimental characterization. For example, a truncated protein may be heterologously produced as an unfolded polypeptide recalcitrant to structure determination (25). Moreover, N-terminal modifications, which are poorly documented in annotation databases, may occur (26, 27).Unfortunately, the poor polypeptide sequence coverage obtained for the numerous low abundance proteins in current shotgun MS/MS proteomic studies implies that the overall detection of N-terminal peptides obtained in proteogenomic studies is relatively low. Different methods for establishing the most extensive list of protein N termini, grouped under the so-called “N-terminomics” theme, have been proposed to selectively enrich or improve the detection of these peptides (2, 28, 29). Large N-terminome studies have recently been reported based on resin-assisted enrichment of N-terminal peptides (30) or terminal amine isotopic labeling of substrates (TAILS) coupled to depletion of internal peptides with a water-soluble aldehyde-functionalized polymer (31–35). Among the numerous N-terminal-oriented methods (2), specific labeling of the N terminus of intact proteins with N-tris(2,4,6-trimethoxyphenyl)phosphonium acetyl succinamide (TMPP-Ac-OSu)¹ has proven reliable (21, 36–39). TMPP-derivatized N-terminal peptides have interesting properties for further LC-MS/MS mass spectrometry: (1) an increase in hydrophobicity because of the trimethoxyphenyl moiety added to the peptides, increasing their retention times in reverse phase chromatography, (2) improvement of their ionization because of the introduction of a positively charged group, and (3) a much simpler fragmentation pattern in tandem mass spectrometry. Other reported approaches rely on acetylation, followed by trypsin digestion, and then biotinylation of free amino groups (40); guanidination of lysine lateral chains followed by N-biotinylation of the N termini and trypsin digestion (41); or reductive amination of all free amino groups with formaldehyde preceeding trypsin digestion (42). Recently, we applied the TMPP method to the proteome of the Deinococcus deserti bacterium isolated from upper sand layers of the Sahara desert (13). This method enabled the detection of N-terminal peptides allowing the confirmation of 278 translation initiation codons, the correction of 73 translation starts, and the identification of non-canonical translation initiation codons (21). However, most TMPP-labeled N-terminal peptides are hidden among the more abundant internal peptides generated after proteolysis of a complex proteome, precluding their detection. This results in disproportionately fewer N-terminal validations, that is, 5 and 8% of total polypeptides coded in the theoretical proteomes of Mycobacterium smegmatis (37) and Deinococcus deserti (21) with a total of 342 and 278 validations, respectively.An interesting chromatographic method to fractionate peptide mixtures for gel-free high-throughput proteome analysis has been developed over the last years and applied to various topics (43, 44). This technique, known as COmbined FRActional DIagonal Chromatography (COFRADIC), uses a double chromatographic separation with a chemical reaction in between to change the physico-chemical properties of the extraneous peptides to be resolved from the peptides of interest. Its previous applications include the separation of methionine-containing peptides (43), N-terminal peptide enrichment (45, 46), sulfur amino acid-containing peptides (47), and phosphorylated peptides (48). COFRADIC was identified as the best method for identification of N-terminal peptides of two archaea, resulting in the identification of 240 polypeptides (9% of the theoretical proteome) for Halobacterium salinarum and 220 (8%) for Natronomonas pharaonis (24).Taking advantage of both the specificity of TMPP labeling, the resolving power of COFRADIC for enrichment, and the increase in information through the use of multiple proteases, we performed the proteogenomic analysis of a marine bacterium from the Roseobacter clade, namely Roseobacter denitrificans OCh114. This novel approach allowed us to validate and correct 534 unique proteins (13% of the theoretical proteome) with TMPP-labeled N-terminal signatures obtained using high-resolution tandem mass spectrometry. We corrected 41 annotations and detected five new open reading frames in the R. denitrificans genome. We further identified eight distinct proteins showing direct evidence for multiple start sites.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]