ABSTRACT
BACKGROUND: 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. RESULTS: We present a clustering-based analysis pipeline integrated with PacBio sequencing data for profiling highly similar gene isoforms. This approach was first evaluated in comparison to de novo assembly of 454 reads using a benchmark admixture containing 10 known, cloned msg genes encoding the major surface glycoprotein of Pneumocystis jirovecii. All 10 msg isoforms were successfully reconstructed with the expected length (~1.5 kb) and correct sequence by the new approach, while 454 reads could not be correctly assembled using various assembly programs. When using an additional benchmark admixture containing 22 known P. jirovecii msg isoforms, this approach accurately reconstructed all but 4 these isoforms in their full-length (~3 kb); these 4 isoforms were present in low concentrations in the admixture. Finally, when applied to the original clinical sample from which the 22 known msg isoforms were cloned, this approach successfully identified not only all known isoforms accurately (~3 kb each) but also 48 novel isoforms. CONCLUSIONS: PacBio sequencing integrated with the clustering-based analysis pipeline achieves high-throughput and high-resolution discrimination of highly similar sequences, and can serve as a new approach for genome-wide characterization of gene isoforms and other highly repetitive sequences.
ABSTRACT
Abstract 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.