Evaluation of viral heterogeneity using next-generation sequencing, end-point limiting-dilution and mass spectrometry

Hepatitis C Virus sequence studies mainly focus on the viral amplicon containing the Hypervariable region 1 (HVR1) to obtain a sample of sequences from which several population genetics parameters can be calculated. Recent advances in sequencing methods allow for analyzing an unprecedented number of viral variants from infected patients and present a novel opportunity for understanding viral evolution, drug resistance and immune escape. In the present paper, we compared three recent technologies for amplicon analysis: (i) Next-Generation Sequencing; (ii) Clonal sequencing using End-point Limiting-dilution for isolation of individual sequence variants followed by Real-Time PCR and sequencing; and (iii) Mass spectrometry of base-specific cleavage reactions of a target sequence. These three technologies were used to assess intra-host diversity and inter-host genetic relatedness in HVR1 amplicons obtained from 38 patients (subgenotypes 1a and 1b). Assessments of intra-host diversity varied greatly between sequence-based and mass-spectrometry-based data. However, assessments of inter-host variability by all three technologies were equally accurate in identification of genetic relatedness among viral strains. These results support the application of all three technologies for molecular epidemiology and population genetics studies. Mass spectrometry is especially promising given its high throughput, low cost and comparable results with sequence-based methods.

Coordinated evolution of the hepatitis B virus polymerase

The detection of compensatory mutations that abrogate negative fitness effects of drug-resistance and vaccine-escape mutations indicates the important role of epistatic connectivity in evolution of viruses, especially under the strong selection pressures. Mapping of epistatic connectivity in the form of coordinated substitutions should help to characterize molecular mechanisms shaping viral evolution and provides a tool for the development of novel anti-viral drugs and vaccines. We analyzed coordinated variation among amino acid sites in 370 the hepatitis B virus (HBV) polymerase sequences using Bayesian networks. Among the HBV polymerase domains the spacer domain separating terminal protein from the reverse-transcriptase domain, showed the highest network centrality. Coordinated substitutions preserve the hydrophobicity and charge of Spacer. Maximum likelihood estimates of codon selection showed that Spacer contains the highest number of positively selected sites. Identification of 67% of the domain lacking an ordered structure suggests that Spacer belongs to the class of intrinsically disordered domains and proteins whose crucial functional role in the regulation of transcription, translation and cellular signal transduction has only recently been recognized. Spacer plays a central role in the epistatic network associating substitutions across the HBV genome, including those conferring viral virulence, drug resistance and vaccine escape. The data suggest that Spacer is extensively involved in coordination of HBV evolution.

Characterization of spo0A homologues in diverse Bacillus and Clostridium species identifies a probable DNA-binding domain

Spo0A is a phosphorylation-activated transcription factor of Bacillus subtilis. It is a member of the response regulator superfamily of bacterial signal transduction proteins and controls many of the changes in gene expression that occur during the transition into stationary phase and during the initiation of sporulation. To identify the domains of Spo0A most critical for determining its structural and functional features, presumptive homologues of the spo0A gene were characterized in a collection of eight Bacillus species and six Clostridium species representing phylogenetically diverse members of these genera. An alignment of the partial or complete DNA sequences of these homologues revealed three regions of especially high conservation in the effector domain. We speculate that the most highly conserved of these corresponds to the recognition helix of a putative helix-turn-helix motif, and, therefore, represents the actual DNA-containing surface of the protein. In the case of homologues identified in Bacillus anthracis and Clostridium acetobutylicum and retrieved by polymerase chain reaction amplification, we confirmed by gene-disruption analysis that the homologue actually is required for initiation of sporulation. Apparent homologues of the B. subtilis spoIVB gene were also discovered immediately upstream from the spo0A homologues in all Bacillus and Clostridium species examined. The discovery of homologues of B. subtilis sporulation genes in these diverse species implies that the gene products required for specifying pathways of sporulation-specific gene activation and for determining key morphogenetic changes may be highly conserved and suggests that an approach similar to that undertaken here might be used as a general strategy to retrieve and compare their gene sequences. Exhaustive efforts to detect a spo0A-like gene in non-endospore formers, including close relatives of Bacillus such as Listeria and Staphylococcus, were uniformly unsuccessful, suggesting that regulation of gene activity during the transition into stationary phase mediated by Spo0A-like proteins may be exclusive to the endospore-forming bacteria.