Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs.

Development of HIV-1 drug resistance mutations (HDRMs) is one of the major reasons for the clinical failure of antiretroviral therapy. Treatment success rates can be improved by applying personalized anti-HIV regimens based on a patient's HDRM profile. However, the sensitivity and specificity of the HDRM profile is limited by the methods used for detection. Sanger-based sequencing technology has traditionally been used for determining HDRM profiles at the single nucleotide variant (SNV) level, but with a sensitivity of only ≥ 20% in the HIV population of a patient. Next Generation Sequencing (NGS) technologies offer greater detection sensitivity (~ 1%) and larger scope (hundreds of samples per run). However, NGS technologies produce reads that are too short to enable the detection of the physical linkages of individual SNVs across the haplotype of each HIV strain present. In this article, we demonstrate that the single-molecule long reads generated using the Third Generation Sequencer (TGS), PacBio RS II, along with the appropriate bioinformatics analysis method, can resolve the HDRM profile at a more advanced quasispecies level. The case studies on patients' HIV samples showed that the quasispecies view produced using the PacBio method offered greater detection sensitivity and was more comprehensive for understanding HDRM situations, which is complement to both Sanger and NGS technologies. In conclusion, the PacBio method, providing a promising new quasispecies level of HDRM profiling, may effect an important change in the field of HIV drug resistance research.

Multiple HIV-1 subtypes and circulating recombinant forms (CRFs) are known to cocirculate in Africa. In West Africa, the high prevalence of CRF02_AG, and cocirculation of subtype A, CRF01_AE, CRF06_cpx, and other complex intersubtype recombinants has been well documented. Mali, situated in the heart of West Africa, is likely to be affected by the spread of recombinant subtypes. However, the dynamics of the spread of HIV-1 recombinant subtypes as well as nonrecombinant HIV-1 group M subtypes in this area have not been systematically assessed. Herein, we undertook genetic analyses on full-length env sequences derived from HIV-1-infected individuals living in the capital city of Mali, Bamako. Of 23 samples we examined, 16 were classified as CRF02_AG and three had a subsubtype A3. Among the remaining HIV-1 strains, CRF06_cpx and CRF09_cpx were each found in two patients. Comparison of phylogenies for six matched pol and full-length env sequences revealed that two strains had discordant subtype/CRF designations between the pol and env regions: one had A3 pol CRF02_AG env and the other had CRF02_AG pol A3 env. Taken together, our study demonstrated the high prevalence of CRF02_AG and complexity of circulating HIV-1 strains in Mali. It also provided evidence of ongoing virus evolution of CRF02_AG, as illustrated by the emergence of more complex CRF02_AG/A3 intersubtype recombinants in this area.

This study evaluated correlation and agreement between version 3 of the Quantiplex human immunodeficiency virus type 1 (HIV-1) RNA assay (v3 branched DNA [bDNA]) and a sensitized Amplicor HIV-1 Monitor assay (reverse transcription [RT]-PCR) for the measurement of HIV RNA. Three hundred eighteen samples from 59 randomly selected, HIV-1-seropositive persons on various drug protocols from the National Institute of Allergy and Infectious Diseases HIV outpatient clinic were studied. The results indicate that v3 bDNA and RT-PCR are highly correlated (r = 0.98) and are in good agreement (mean difference in log(10) copies/ml +/- 2 standard deviations = 0.072 +/- 0.371). The relationship between values obtained by both assays is given by the following equation: log(10)v3 bDNA = -0.0915 + 1.0052.log(10)RT-PCR This represents a 1.026-fold difference between log(10)RT-PCR values and log(10)v3 bDNA values.