Tracing hepatitis C virus transmission by using high-throughput sequencing technique

Abstract

An approach for analysis of hepatitis C virus (HCV) quasispecies using Hiseq high-throughput sequencing (hereinafter referred to as Hiseq sequencing) technique was developed and then applied to investigate a possible case of HCV needle sharing transmission. One case of HCV antibody seroconversion (P1) was found in a methadone clinic on January 15, 2015. Four HCV antibody positive injecting drug users (IDUs), P2 to P5, suspected to be involved in needle sharing transmission with P1 during the period (after March 24, 2014) that P1 may be infected with HCV were investigated, and another 28 HCV antibody positive IDUs were selected as controls (C1 to C28). These controls came from the same methadone clinic or lived in the same town with P1. The RNAs were extracted from the plasma specimens and then reverse-transcribed into cDNA. After HCV subtyping, Hiseq sequencing was performed to detect and sequence the HCV quasispecies (263 bp) in the specimens with the same subtype as P1. The frequency of quasispecies was counted and ranked. Intrapersonal and interpersonal genetic distance and phylogenetic tree were calculated. The HCV subtype of specimen P1 was 3b. All the other specimens with the same subtype were P2, C7, C12, C14, C15, C16, C19, C20 and C28. Hiseq sequencing was successfully performed in 9 out of these 10 specimens, and 249 753 to 1 086 333 (average 869 608) cleaned sequences representing 3 to 172 (average 48) unique HCV quasispecies were obtained. The medians (P50) of intrapersonal genetic diversities from the 9 specimens were 0.4% to 12.3%. The P50 (P25, P75) of genetic diversities between P1 and the other 8 specimens were 19.0% (18.4%, 19.8%), 10.4% (2.8%, 18.3%), 19.6% (17.8%, 21.4%),24.9% (23.8%, 26.1%), 19.8% (18.7%, 20.7%), 20.1% (18.9%, 21.2%), 20.6% (20.0%, 21.1%), 23.6% (22.4%, 24.8%). There were no significant difference between the genetic diversities of P1 and P2 and those of P1 and other 7 specimens (H=9.40, P=0.100). The genetic diversities between few HCV quasispecies from P1 and few ones from C7 were 0. Phylogenetic tree analysis indicated that there was no HCV transmission relationship between P1 and P2, but there was HCV transmission relationship between P1 and C7. With the feature of high-throughput, easier operation and lower cost, Hiseq sequencing technique has high practical value in tracing HCV transmission at the quasispecies level.