Genotypic analysis of human immunodeficiency virus type 1 env V3 loop sequences: bioinformatics prediction of coreceptor usage among 28 infected mother-infant pairs in a drug-naive population

HIV coreceptor tropism determination and mutational pattern identification

In the early stages of infection, Human Immunodeficiency Virus Type 1 (HIV-1) generally selects CCR5 as the primary coreceptor for entering the host cell. As infection progresses, the virus evolves and may exhibit a coreceptor-switch to CXCR4. Accurate determination coreceptor usage and identification key mutational patterns associated tropism switch are essential for selection of appropriate therapies and understanding mechanism of coreceptor change. We developed a classifier composed of two coreceptor-specific weight matrices (CMs) based on a full-scale dataset. For this classifier, we found an AUC of 0.97, an accuracy of 95.21% and an MCC of 0.885 (sensitivity 92.92%; specificity 95.54%) in a ten-fold cross-validation, outperforming all other methods on an independent dataset (13% higher MCC value than geno2pheno and 15% higher MCC value than PSSM). A web server (http://spg.med.tsinghua.edu.cn/CM.html) based on our classifier was provided. Patterns of genetic mutations that occur along with coreceptor transitions were further identified based on the score of each sequence. Six pairs of one-AA mutational patterns and three pairs of two-AA mutational patterns were identified to associate with increasing propensity for X4 tropism. These mutational patterns offered new insights into the mechanism of coreceptor switch and aided in monitoring coreceptor switch.

Evolution of HIV-1 coreceptor usage and coreceptor switching during pregnancy

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. To document the evolution of coreceptor tropism during pregnancy, a longitudinal study of envelope gene sequences was performed in a group of pregnant women infected with HIV-1 of clade B (n=10) or non-B (n=9). Polymerase chain reaction (PCR) amplification of the V1-V3 region was performed on plasma viral RNA, followed by cloning and sequencing. Using geno2pheno and PSSMX4R5, the presence of X4 variants was predicted in nine of 19 subjects (X4 subjects) independent of HIV-1 clade. Six of nine X4 subjects exhibited CD4(+) T cell counts <200 cells/mm(3), and the presence of X4-capable virus was confirmed using a recombinant phenotypic assay in four of seven cases where testing was successful. In five of nine X4 subjects, a statistically significant decline in the geno2pheno false-positive rate was observed during the course of pregnancy, invariably accompanied by progressive increases in the PSSMX4R5 score, the net charge of V3, and the relative representation of X4 sequences. Evolution toward X4 tropism was also echoed in the primary structure of V2, as an accumulation of substitutions associated with CXCR4 tropism was seen in X4 subjects. Results from these experiments provide the first evidence of the ongoing evolution of coreceptor utilization from CCR5 to CXCR4 during pregnancy in a significant fraction of HIV-infected women. These results inform changes in host-pathogen interactions that lead to a directional shaping of viral populations and viral tropism during pregnancy, and provide insights into the biology of HIV transmission from mother to child.

Phylogenetic analysis of human immunodeficiency virus type 1 subtype C env gp120 sequences among four drug-naive families following subsequent heterosexual and vertical transmissions

To characterize phylogenetic relatedness of plasma HIV-1 RNA subtype C env gp120 viral variants capable of establishing an infection following heterosexual and subsequent vertical transmission events a 650-base pair fragment within the C2-V5 subregion was sequenced from four HIV-1-infected families each consisting of biological parent(s), index children (first), and subsequent (second) siblings. None of the family members had received antiretroviral therapy at the time of sample collection. Sequence alignment and analysis were done using Gene Doc, Clustal X, and MEGA software programs. Second siblings' sequences were homogeneous and clustered in a single branch while first siblings' sequences were more heterogeneous, clustering in separate branches, suggestive of more than one donor variants responsible for the infection or evolution from founder variant(s) could have occurred. While the directionality for heterosexual transmission could not be determined, homogeneous viral variants were a unique characteristic of maternal variants as opposed to the more heterogeneous paternal variants. Analysis of families' sequences demonstrated a localized expansion of the subtype C infection. We demonstrated that families' sequences clustered quite closely with other regional HIV-1 subtype C sequences supported by a bootstrap value of 86%, confirming the difficulty of classifying subtype C sequences on a geographic basis. Data are indicative of several mechanisms that may be involved in both vertical and heterosexual transmission. Larger studies are warranted to address the caveats of this study and build on the strengths. Our study could be the beginning of family-based HIV-1 intervention research in Zimbabwe.