Local phylogenetic analysis identifies distinct trends in transmitted HIV drug resistance: implications for public health interventions

Abstract

Background: HIV transmitted drug resistance (TDR) surveillance is usually conducted by sampling from a large population. However, overall TDR prevalence results may be inaccurate for many individual clinical setting. We analyzed HIV genotypes at a tertiary care setting in Ottawa, Ontario in order to evaluate local TDR patterns among sub-populations.

Method: Genotyping reports were digitized from ART naïve patients followed at the Immunodeficiency Clinic at the Ottawa Hospital, between 2008 and 2010. Quality controlled, digitized sequence data were assessed for TDR using the Stanford HIV Database. Patient characteristics were analyzed according to TDR patterns. Finally, a phylogenetic tree was constructed to elucidate the observed pattern of HIV TDR.

Results: Among the 155 clinic patients there was no statistically significantly difference in demographics as compared to the Ontario provincial HIV population. The clinic prevalence of TDR was 12.3%; however, in contrast to the data from Ontario, TDR patterns were inverted with a 21% prevalence among MSM and 5.5% among IDU. Furthermore, nearly 80% of the observed TDR was a D67N/K219Q pattern with 87% of these infections arising from a distinct phylogenetic cluster.

Conclusions: Local patterns of TDR were distinct to what had been observed provincially. Phylogenetic analysis uncovered a cluster of related infections among MSM that appeared more likely to be recent infections. Results support a paradigm of routine local TDR surveillance to identify the sub-populations under care. Furthermore, the routine application of phylogenetic analysis in the TDR surveillance context provides insights into how best to target prevention strategies; and how to correctly measure outcomes.

Figure 1

Maximum likelihood tree of HIV pol sequences generated using GTR model in PHYML with 100 bootstrap replicates. Sequences were trimmed to 986 bp with identified TDR mutation codons removed. 155 study cohort sequences are present in a background of 255 sequences further described in the methods. Orange box shows the D67N/K219Q cluster that was identified. For comparison purposes another cluster identified in the study cohort is identified in blue. ** Indicates the out-group N sequence N.CM.97.YBF106.