High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations.

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals. Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA). Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies. Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

Recombination Between High-Risk Human Papillomaviruses and Non-Human Primate Papillomaviruses: Evidence of Ancient Host Switching Among Alphapapillomaviruses.

We use all the currently known 405 Papillomavirus (PV) sequences, 343 curated PV sequences from both humans and animals from the PAVE data base, to analyse the recombination dynamics of these viruses at the whole genome levels. After showing some evidence of human and non-human primate PV recombination, we report a comprehensive recombination analysis of all currently known 82 Alphapapillomaviruses (Alpha-PVs). We carried out an exploratory study and found novel recombination events between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 (PpPV1) Papillomaviruses. This is the first evidence of interactions between PVs from different hosts and hence postulates the likelihood of ancient host switching among Alpha-PVs. Notwithstanding these results should be interpreted with caution because the major and minor parents indicated by RDP4 program are simply the sequences in the alignment that most closely resemble the actual parents. We found statistically significant differences between the phylogenies of the PV sequences with recombination regions and PV sequences without recombination regions using the Shimodaira-Hasegawa phylogenetic incongruence testing. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts before speciation of primate species. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. However, a clearer picture of the PVs evolutionary landscape can only be achieved with the incremental discovery of PV sequences, especially from the animal kingdom.

KIR and HLA-C Genetic Polymorphisms Influence Plasma IP-10 Concentration in Antiretroviral Therapy-Naive HIV-Infected Adult Zimbabweans.

Past studies on the relationship between Killer cell Immunoglobulin-like Receptor (KIR) and Human Leukocyte Antigen (HLA) genetic variation and chronic immune activation (CIA) in HIV infection are not uniformly consistent. Moreover, interferon-γ-induced protein 10 (IP-10) is a soluble biomarker of immune activation, with high plasma concentrations predicting accelerated disease progression in HIV infection. Thus, we investigated the association of KIR and HLA-C genetic polymorphisms with plasma IP-10 concentration in 183 treatment-naive chronically HIV-infected adults of Bantu origin from Zimbabwe. KIR genetic variation was determined using allele-specific primer PCR while HLA-C typing was characterized by sequencing. Plasma IP-10 was quantified using enzyme-linked immunosorbent assay. The KIR2DL3 gene was significantly associated with CIA as observed from IP-10 concentrations among KIR2DL3 carriers (265.20 pg/mL, IQR: 179.99-385.19) compared with KIR2DL3 noncarriers (183.56 pg/mL; IQR: 110.98-230.81; p = 0.001) and among KIR2DL3+HLA-C2 carriers (226.23 pg/mL, IQR: 187.96-394.73) compared with KIR2DL3+HLA-C2 noncarriers (212.86 pg/mL, IQR: 160.15-344.99; p = 0.017), respectively. Similarly, IP-10 concentrations were significantly higher (p = 0.030) in the KIR3DS1 carriers (313.86 pg/mL, IQR: 230.05-469.20) compared with KIR3DS1 noncarriers (246.01 pg/mL, IQR: 169.58-373.32). Thus, KIR and HLA-C could be playing important roles in HIV-associated immune activation. The elevation of IP-10 in KIR2DL3 and KIR2DL3+C2 could potentially be explained by increased IFN-γ secretion from activated NK cell activation due to the absence of KIR2DL3's cognate C1 ligand. To the best of our knowledge, this is the first study on a potential link between KIR and HLA-C genetic determinants and plasma IP-10 concentration in this population sample. Future studies are called for in other world populations for biomarkers of disease progression and mechanisms of IP-10 variability in HIV infection.

Human Leukocyte Antigen-A, B, C, DRB1, and DQB1 Allele and Haplotype Frequencies in a Subset of 237 Donors in the South African Bone Marrow Registry.

Human leukocyte antigen- (HLA-) A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 allele and haplotype frequencies were studied in a subset of 237 volunteer bone marrow donors registered at the South African Bone Marrow Registry (SABMR). Hapl-o-Mat software was used to compute allele and haplotype frequencies from individuals typed at various resolutions, with some alleles in multiple allele code (MAC) format. Four hundred and thirty-eight HLA-A, 235 HLA-B, 234 HLA-DRB1, 41 HLA-DQB1, and 29 HLA-C alleles are reported. The most frequent alleles were A∗02:02g (0.096), B∗07:02g (0.082), C∗07:02g (0.180), DQB1∗06:02 (0.157), and DRB1∗15:01 (0.072). The most common haplotype was A∗03:01g~B∗07:02g~C∗07:02g~DQB1∗06:02~DRB1∗15:01 (0.067), which has also been reported in other populations. Deviations from Hardy-Weinberg equilibrium were observed in A, B, and DRB1 loci, with C~DQB1 being the only locus pair in linkage disequilibrium. This study describes allele and haplotype frequencies from a subset of donors registered at SABMR, the only active bone marrow donor registry in Africa. Although the sample size was small, our results form a key resource for future population studies, disease association studies, and donor recruitment strategies.

HLA-C*18:01 and KIR2DL2+C1 genetic variants are associated with low viral load in cART naïve HIV-infected adult Zimbabweans.

INTRODUCTION: Polymorphisms in killer cell immunoglobulin-like receptor (KIR) and human leukocyte antigen (HLA) gene families are implicated in differential outcomes of HIV infection. However, research findings on the influence of KIR and HLA-C polymorphism on HIV disease progression remain inconclusive. We thus investigated the association of KIR and HLA-C gene polymorphisms with plasma HIV load (VL) and CD4+ T lymphocyte (CD4) count in 183 chronically HIV-infected, combination antiretroviral therapy (cART) naïve Zimbabweans of Bantu origin. METHODOLOGY: The presence or absence of 15 KIR genes were determined using sequence specific primer polymerase chain reaction while HLA-C typing was performed using chain termination DNA sequencing. Plasma VL was determined using the Cavidi Exavir viral load version 3 assay while CD4+ T lymphocytes were enumerated using flow cytometry. VLs and CD4 counts were compared between gene/genotype carriers and non-carriers using Mann-Whitney ranksum test. RESULTS: HLA-C*18:01 allele carriers had a significantly lower median log10 VL (2.87copies/mL [IQR;2.3-3.2]) than the non-C*18:01 carriers (3.33copies/mL [IQR; 2.74-3.9]), p = 0.018. Further, median log10 VL was significantly lower in KIR2DL2+C1 carriers (2.745 [IQR; 2.590-2.745]) than non-KIR2DL2+C1 carriers (3.4 [IQR; 2.746-3.412]), p = 0.041. Comparison of CD4 + T lymphocyte counts between C*08:02 allele carriers and non-C*08:02 carriers showed a significantly higher median CD4 count in C*08:02 carriers (548cells/µL [IQR;410-684]) than in non-carriers (428cells/µL [IQR;388-537]), p = 0.034. CONCLUSION: We conclude that the HLA-C*18:01 and KIR2DL2+C1 genetic variants are associated with low VL while the C*08:02 is associated with high CD4+ T lymphocyte count among cART naïve Zimbabwean adults with chronic HIV infection.

Plasma IP-10 Concentrations Correlate Positively with Viraemia and Inversely with CD4 Counts in Untreated HIV Infection.

BACKGROUND: Chronic immune activation is a feature of HIV infection associated with accelerated HIV disease progression. There is conflicting data on the association of biomarkers of immune activation with traditional markers of HIV disease progression; CD4 counts and viral load (VL). OBJECTIVE: The study aimed to determine the association of biomarkers of immune activation; interferon (IFN)-γ-induced protein 10 (IP-10) and soluble cluster of differentiation 14 (sCD14) in chronic HIV infection with traditional markers of HIV disease progression. METHODS: We collected demographic data, enumerated CD4 counts and quantified VL in 183 antiretroviral therapy (ART)-naive adults with chronic HIV infection. Plasma concentrations of IP-10 and sCD14 were quantified in the ART-naive adults with chronic HIV infection and 75 HIV-uninfected controls. RESULTS: IP-10 concentrations were significantly higher in the HIV-infected group (median; 257.40pg/ml, IQR; 174.08-376.32) than in the HIV-uninfected (median; 86.19pg/ml, IQR; 67.70-116.39) (P<0.001). Similarly, sCD14 concentrations were significantly higher in the HIV-infected (median; 1.45µg/ml, IQR; 1.02-2.16) group than in the controls (median; 0.89µ/ml, IQR; 0.74-1.18) (P<0.001). High log10 IP-10 concentrations were positively correlated with high log10 viral loads (Spearman's correlation coefficient [R]=0.21, P=0.003) and inversely correlated with low CD4 counts (R= -0.19, P=0.011). In contrast, log10 sCD14 was not significantly associated with either log10 viral loads (R=0.03, P=0.707) nor CD4 count (R=-0.04, P=0.568). CONCLUSION: We conclude that plasma sCD14 and IP-10 were elevated in the HIV-infected patients compared to HIV-uninfected individuals possibly due to on-going immune activation. In addition, plasma high concentrations of IP-10 but not sCD14 concentrations are associated with high VL and low CD4 count.

KIR Gene Content Diversity in a Zimbabwean Population: Does KIR2DL2 Have a Role in Protection Against Human Immunodeficiency Virus Infection?

Killer cell immunoglobulin-like receptors (KIRs) mediate natural killer cell function through interaction with their cognate human leukocyte antigen ligands. Thus, KIR gene variants have been implicated in resistance or susceptibility to viral infections. However, research on the role of these variants remains contradictory and inconclusive. In the present study, we investigated KIR gene content diversity and its association with human immunodeficiency virus (HIV) infection in an adult Black Zimbabwean population. Presence or absence of 15 KIR genes was determined in 189 HIV-infected adults and 97 HIV-uninfected blood donors using sequence specific primer polymerase chain reaction. Frequencies of KIR genes, genotypes, and haplotypes were compared between the cases and controls to identify putative associations between KIR gene variants and HIV status. We report in this study the frequencies of 15 KIR genes and 43 KIR genotypes (40 known and 3 novel) among Zimbabweans. Importantly, the frequency of the inhibitory KIR2DL2 gene was significantly higher in the uninfected group (62%) compared to the HIV-infected group (47%) (OR = 0.55, 95% CI: 0.33-0.90, p = 0.019). KIR2DL2/2DL2 homozygosity was also significantly higher in the uninfected group (35%) compared to HIV-infected group (53%) (OR = 0.33, 95% CI: 0.16-0.72, p = 0.005) under a recessive model. We conclude that the KIR2DL2 gene may be involved in protection against HIV infection. It may be possible that inhibitory KIR genes may have an important role to play in HIV acquisition among populations of African origin in whom the activating KIR genes are less frequent compared to among Caucasians.

Human Leukocyte Antigen Diversity: A Southern African Perspective.

Despite the increasingly well-documented evidence of high genetic, ethnic, and linguistic diversity amongst African populations, there is limited data on human leukocyte antigen (HLA) diversity in these populations. HLA is part of the host defense mechanism mediated through antigen presentation to effector cells of the immune system. With the high disease burden in southern Africa, HLA diversity data is increasingly important in the design of population-specific vaccines and the improvement of transplantation therapeutic interventions. This review highlights the paucity of HLA diversity data amongst southern African populations and defines a need for information of this kind. This information will support disease association studies, provide guidance in vaccine design, and improve transplantation outcomes.

Frequency of Betapapillomavirus infections among HIV infected and uninfected Black Zimbabweans with cutaneous lesions.

Human papillomavirus (HPV) types from the Betapapillomavirus (ß-HPV) genus are plentiful in non-melanoma skin cancers and warts among Caucasians, but there is paucity of information among black Africans. To determine the frequency of ß-HPV genotypes in cutaneous infections among Black Zimbabweans, a cross-sectional study was carried out in which blood samples and skin biopsies were collected from patients infected and uninfected with HIV attending a referral hospital. We included 144 participants (72 infected and 72 uninfected with HIV) with clinically apparent cutaneous warts (n = 34), suspected non-melanoma skin cancers (n = 98) and Kaposi sarcoma (KS) (n = 18). The skin biopsies were analyzed for HPV DNA presence and type. ß-HPV DNA was identified among 70% (101/144) and was significantly higher among patients infected with HIV, 79% (57/72) compared to the HIV uninfected 61% (44/72) [OR = 2.42, 95% CI (1.09-5.47), P = 0.018]. All patients with warts, 89% of those with KS and 58% of those with non-melanoma skin cancers were HPV DNA positive and ß-HPV type 14 was identified in nearly half of the study participants 49.3% (71/144). Single HPV infections were observed in 33.7% (34/101) of the participants that were HPV DNA positive, 66.3% (67/101) had multiple HPV types. There was no significant difference between patients infected and uninfected with HIV in terms of multiple HPV infections. The distribution of different HPV types did not reveal any association with age and gender but there was an association between HPV 14 and HIV immune status. ß-HPVs are not uncommon among the Black Zimbabweans with skin lesions.

Prevalence of Extensively Drug Resistant Tuberculosis among Archived Multidrug Resistant Tuberculosis Isolates in Zimbabwe.

We conducted a cross-sectional study of second line drug resistance patterns and genetic diversity of MDR-TB isolates archived at the BRTI-TB Laboratory, Harare, between January 2007 and December 2011. DSTs were performed for second line antituberculosis drugs. XDR-TB strains were defined as MDR-TB strains with resistance to either kanamycin and ofloxacin or capreomycin and ofloxacin. Strain types were identified by spoligotyping. No resistance to any second line drugs was shown in 73% of the isolates, with 23% resistant to one or two drugs but not meeting the definition of XDR-TB. A total of 26 shared types were identified, and 18 (69%) matched preexisting shared types in the current published spoligotype databases. Of the 11 out of 18 clustered SITs, 4 predominant (>6 isolates per shared type) were identified. The most and least abundant types were SIT 1468 (LAM 11-ZWE) with 12 (18%) isolates and SIT 53 (T1) with 6 (9%) isolates, respectively. XDR-TB strains are rare in Zimbabwe, but the high proportion of "pre-XDR-TB" strains and treatment failure cases is of concern. The genetic diversity of the MDR-TB strains showed no significant association between SITs and drug resistance.

Surveillance of transmitted antiretroviral drug resistance among HIV-1 infected women attending antenatal clinics in Chitungwiza, Zimbabwe.

The rapid scale-up of highly active antiretroviral therapy (HAART) and use of single dose Nevirapine (SD NVP) for prevention of mother-to-child transmission (pMTCT) have raised fears about the emergence of resistance to the first line antiretroviral drug regimens. A cross-sectional study was conducted to determine the prevalence of primary drug resistance (PDR) in a cohort of young (<25 yrs) HAART-naïve HIV pregnant women attending antenatal clinics in Chitungwiza, Zimbabwe. Whole blood was collected in EDTA for CD4 counts, viral load, serological estimation of duration of infection using the BED Calypte assay and genotyping for drug resistance. Four hundred and seventy-one women, mean age 21 years; SD: 2.1 were enrolled into the study between 2006 and 2007. Their median CD4 count was 371cells/µL; IQR: 255-511 cells/µL. Two hundred and thirty-six samples were genotyped for drug resistance. Based on the BED assay, 27% were recently infected (RI) whilst 73% had long-term infection (LTI). Median CD4 count was higher (p<0.05) in RI than in women with LTI. Only 2 women had drug resistance mutations; protease I85V and reverse transcriptase Y181C. Prevalence of PDR in Chitungwiza, 4 years after commencement of the national ART program remained below WHO threshold limit (5%). Frequency of recent infection BED testing is consistent with high HIV acquisition during pregnancy. With the scale-up of long-term ART programs, maintenance of proper prescribing practices, continuous monitoring of patients and reinforcement of adherence may prevent the acquisition and transmission of PDR.