Transcriptomic signatures of progression to TB disease among close contacts in Brazil.

BACKGROUND: Approximately 5% of people infected with Mycobacterium tuberculosis progress to tuberculosis (TB) disease without preventive therapy. There is a need for a prognostic test to identify those at highest risk of incident TB, so that therapy can be targeted. We evaluated host blood transcriptomic signatures for progression to TB disease. METHODS: Close contacts (≥4 hours exposure per week) of adult patients with culture-confirmed pulmonary TB were enrolled in Brazil. Investigation for incident, microbiologically-confirmed or clinically-diagnosed pulmonary or extra-pulmonary TB disease through 24 months of follow-up was symptom-triggered. Twenty previously validated blood TB transcriptomic signatures were measured at baseline by real-time quantitative PCR. Prognostic performance for incident TB was tested using receiver operating characteristic curve (ROC) analysis at 6, 9, 12, and 24 months of follow-up. RESULTS: Between June 2015 and June 2019, 1,854 close contacts were enrolled; Twenty-five progressed to incident TB, of whom 13 had microbiologically-confirmed disease. Baseline transcriptomic signature scores were measured in 1,789 close contacts. Prognostic performance for all signatures was best within 6 months of diagnosis. Seven signatures (Gliddon4, Suliman4, Roe3, Roe1, Penn-Nicholson6, Francisco2, and Rajan5) met the minimum World Health Organization target product profile (TPP) for a prognostic test through 6 months; three (Gliddon4, Rajan5, and Duffy9) through 9 months. None met the TPP threshold through 12 or more months of follow-up. CONCLUSIONS: Blood transcriptomic signatures may be useful for predicting TB risk within 9 months of measurement among TB-exposed contacts, to target preventive therapy administration.

Genome-wide association study reveals genetic variants associated with HIV-1C infection in a Botswana study population.

Although there have been many studies of gene variant association with different stages of HIV/AIDS progression in United States and European cohorts, few gene-association studies have assessed genic determinants in sub-Saharan African populations, which have the highest density of HIV infections worldwide. We carried out genome-wide association studies on 766 study participants at risk for HIV-1 subtype C (HIV-1C) infection in Botswana. Three gene associations (AP3B1, PTPRA, and NEO1) were shown to have significant association with HIV-1C acquisition. Each gene association was replicated within Botswana or in the United States-African American or United States-European American AIDS cohorts or in both. Each associated gene has a prior reported influence on HIV/AIDS pathogenesis. Thirteen previously discovered AIDS restriction genes were further replicated in the Botswana cohorts, extending our confidence in these prior AIDS restriction gene reports. This work presents an early step toward the identification of genetic variants associated with and affecting HIV acquisition or AIDS progression in the understudied HIV-1C afflicted Botswana population.

Simulation of African and non-African low and high coverage whole genome sequence data to assess variant calling approaches.

Current variant calling (VC) approaches have been designed to leverage populations of long-range haplotypes and were benchmarked using populations of European descent, whereas most genetic diversity is found in non-European such as Africa populations. Working with these genetically diverse populations, VC tools may produce false positive and false negative results, which may produce misleading conclusions in prioritization of mutations, clinical relevancy and actionability of genes. The most prominent question is which tool or pipeline has a high rate of sensitivity and precision when analysing African data with either low or high sequence coverage, given the high genetic diversity and heterogeneity of this data. Here, a total of 100 synthetic Whole Genome Sequencing (WGS) samples, mimicking the genetics profile of African and European subjects for different specific coverage levels (high/low), have been generated to assess the performance of nine different VC tools on these contrasting datasets. The performances of these tools were assessed in false positive and false negative call rates by comparing the simulated golden variants to the variants identified by each VC tool. Combining our results on sensitivity and positive predictive value (PPV), VarDict [PPV = 0.999 and Matthews correlation coefficient (MCC) = 0.832] and BCFtools (PPV = 0.999 and MCC = 0.813) perform best when using African population data on high and low coverage data. Overall, current VC tools produce high false positive and false negative rates when analysing African compared with European data. This highlights the need for development of VC approaches with high sensitivity and precision tailored for populations characterized by high genetic variations and low linkage disequilibrium.

Dating admixture events is unsolved problem in multi-way admixed populations.

Advances in human sequencing technologies, coupled with statistical and computational tools, have fostered the development of methods for dating admixture events. These methods have merits and drawbacks in estimating admixture events in multi-way admixed populations. Here, we first provide a comprehensive review and comparison of current methods pertinent to dating admixture events. Second, we assess various admixture dating tools. We do so by performing various simulations. Third, we apply the top two assessed methods to real data of a uniquely admixed population from South Africa. Results reveal that current dating admixture models are not sufficiently equipped to estimate ancient admixtures events and to identify multi-faceted admixture events in complex multi-way admixed populations. We conclude with a discussion of research areas where further work on dating admixture-based methods is needed.

Whole genome sequencing reveals population diversity and variation in HIV-1 specific host genes.

HIV infection continues to be a major global public health issue. The population heterogeneity in susceptibility or resistance to HIV-1 and progression upon infection is attributable to, among other factors, host genetic variation. Therefore, identifying population-specific variation and genetic modifiers of HIV infectivity can catapult the invention of effective strategies against HIV-1 in African populations. Here, we investigated whole genome sequences of 390 unrelated HIV-positive and -negative individuals from Botswana. We report 27.7 million single nucleotide variations (SNVs) in the complete genomes of Botswana nationals, of which 2.8 million were missing in public databases. Our population structure analysis revealed a largely homogenous structure in the Botswana population. Admixture analysis showed elevated components shared between the Botswana population and the Niger-Congo (65.9%), Khoe-San (32.9%), and Europeans (1.1%) ancestries in the population of Botswana. Statistical significance of the mutational burden of deleterious and loss-of-function variants per gene against a null model was estimated. The most deleterious variants were enriched in five genes: ACTRT2 (the Actin Related Protein T2), HOXD12 (homeobox D12), ABCB5 (ATP binding cassette subfamily B member 5), ATP8B4 (ATPase phospholipid transporting 8B4) and ABCC12 (ATP Binding Cassette Subfamily C Member 12). These genes are enriched in the glycolysis and gluconeogenesis (p < 2.84e-6) pathways and therefore, may contribute to the emerging field of immunometabolism in which therapy against HIV-1 infection is being evaluated. Published transcriptomic evidence supports the role of the glycolysis/gluconeogenesis pathways in the regulation of susceptibility to HIV, and that cumulative effects of genetic modifiers in glycolysis/gluconeogenesis pathways may potentially have effects on the expression and clinical variability of HIV-1. Identified genes and pathways provide novel avenues for other interventions, with the potential for informing the design of new therapeutics.

The role of CYP2B6 516G>T polymorphism on efavirenz/nevirapine toxicity. Implications on treatment outcomes: Lessons from Botswana.

ABSTRACT: The two non-nucleoside reverse transcriptase inhibitors (NNRTIs), efavirenz (EFV) and nevirapine (NVP), are currently the core antiretroviral drugs for treatment of HIV in sub-Saharan Africa including Botswana. The drugs are metabolized by Cytochrome P450 2B6 (CYP2B6) liver enzyme. The CYP2B6 gene that encodes for metabolism of these drugs is known to be highly polymorphic. One of the polymorphism in the CYP2B6 gene, 516G>T, particularly the 516T allele, is known to confer poor metabolism of EFV and NVP. This may lead to high levels of plasma drug concentrations and development of treatment toxicities, like central nervous system toxicities, and cutaneous and hepatic toxicities, for EFV and NVP, respectively. The CYP2B6 516G allele on the other hand is associated with an extensive metabolism of the two NNRTIs drugs. We sought to establish association between possible developments of NNRTIs toxicities with CYP2B6 516G>T variation in Botswana.A total of 316 peripheral blood mononuclear cells samples were used in a retrospective view. All the samples were from participants on EFV/NVP-containing regimen with known toxicity output. TaqMan Real-Time PCR approach was applied for assessing CYP2B6 516 allele variation in cases with treatment toxicity and those without. Analysis was performed by chi-square statistics and logistic regression analysis.The rate of poor metabolizers among participants with toxicity and those without toxicity was 18.4% and 15.1%, respectively. The CYP2B6 516 genotype distribution comparisons between the participants with toxicity and those without were not statistically different (chi-square = .326; P = .568).CYP2B6 516 variation was not associated with NNRTI toxicity. No other factors were associated with toxicity when considering age, baseline body mass index, baseline CD4, baseline HIV viral load and adherence. The results were discussed in the context of all the studies done in Botswana to date.

Association of CYP2B6 Genetic Variation with Efavirenz and Nevirapine Drug Resistance in HIV-1 Patients from Botswana.

PURPOSE: CYP2B6 liver enzyme metabolizes the two non-nucleoside reverse transcriptase inhibitors Efavirenz (EFV) and Nevirapine (NVP) used in the antiretroviral therapy (ART) regimens for HIV-infected individuals. Polymorphisms of the CYP2B6 gene influence drug levels in plasma and possibly virological outcomes. The aim of this study was to explore the potential impact of CYP2B6 genotype and haplotype variation on the risk of developing EFV/NVP drug resistance mutations (DRMs) in HIV-1 patients receiving EFV-/NVP-containing regimens in Botswana. PATIENTS AND METHODS: Participants were a sub-sample of a larger study (Tshepo study) conducted in Gaborone, Botswana, among HIV-infected individuals taking EFV/NVP containing ART. Study samples were retrieved and assigned to cases (with DRMs) and controls (without DRMs). Four single-nucleotide polymorphisms (SNPs) in the CYP2B6 gene (-82T>C; 516G>T; 785A>G; 983T>C) were genotyped, the haplotypes reconstructed, and the metabolic score assigned. The possible association between drug resistance and several independent factors (baseline characteristics and CYP2B6 genotypes) was assessed by Binary Logistic Regression (BLR) analysis. EFV/NVP resistance status and CYP2B6 haplotypes were also analyzed using Z-test, chi-square and Fisher's exact test statistics. RESULTS: Two hundred and twenty-seven samples were analysed (40 with DRMs, 187 without DRMs). BLR analysis showed an association between EFV/NVP resistance and CYP2B6 516G allele (OR: 2.26; 95% CI: 1.27-4.01; P=0.005). Moreover, haplotype analysis revealed that the proportion of EFV/NVP-resistant infections was higher among CYP2B6 fast than extensive/slow metabolizers (30.8% vs 16.8%; P=0.035), with the 516G allele more represented in the haplotypes of fast than extensive/slow metabolizers (100.0% vs 53.8%; P<0.001). CONCLUSION: We demonstrated that the CYP2B6 516G allele, and even more when combined in fast metabolic haplotypes, is associated with the presence of EFV/NVP resistance, strengthening the need to assess the CYP2B6 genetic profiles in HIV-infected patients in order to improve the virologic outcomes of NNRTI containing ART.

A broad survey of DNA sequence data simulation tools.

In silico DNA sequence generation is a powerful technology to evaluate and validate bioinformatics tools, and accordingly more than 35 DNA sequence simulation tools have been developed. With such a diverse array of tools to choose from, an important question is: Which tool should be used for a desired outcome? This question is largely unanswered as documentation for many of these DNA simulation tools is sparse. To address this, we performed a review of DNA sequence simulation tools developed to date and evaluated 20 state-of-art DNA sequence simulation tools on their ability to produce accurate reads based on their implemented sequence error model. We provide a succinct description of each tool and suggest which tool is most appropriate for the given different scenarios. Given the multitude of similar yet non-identical tools, researchers can use this review as a guide to inform their choice of DNA sequence simulation tool. This paves the way towards assessing existing tools in a unified framework, as well as enabling different simulation scenario analysis within the same framework.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa.

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

Cytomegalovirus Viremia in HIV-1 Subtype C Positive Women at Delivery in Botswana and Adverse Birth/Infant Health Outcomes.

BACKGROUND: We evaluated the association between maternal cytomegalovirus (CMV) viremia during pregnancy and adverse birth and infant health outcomes in HIV-infected mothers and their HIV-exposed uninfected infants. METHODS: HIV-positive women and their infants were followed prospectively from pregnancy through 2 years postpartum in the "Tshipidi" study in Botswana. We analyzed the association between detectable CMV DNA in maternal blood at delivery and adverse birth outcomes (stillbirth, preterm delivery, small for gestational age, or birth defect), as well as infant hospitalization and mortality through 24 months. RESULTS: We measured CMV DNA in blood samples from 350 (77.1%) of 454 HIV-positive women from the Tshipidi study. The median maternal CD4 count was 422 cells/mL, and median HIV-1 RNA at entry was 3.2 log10 copies/mL. Fifty-one (14.6%) women had detectable CMV DNA. In unadjusted analyses, detectable CMV DNA was associated with higher maternal HIV-1 RNA [odds ratio (OR) 1.4, 95% confidence interval (CI): 1.1 to 1.9], presence of a birth defect (OR 9.8, 95% CI: 1.6 to 60.3), and occurrence of any adverse birth outcome (OR 2.0, 95% CI: 1.04 to 3.95). In multivariable analysis, we observed a trend toward association between detectable maternal CMV DNA and occurrence of any adverse birth outcome (adjusted OR 1.9, 95% CI: 0.96 to 3.8). Maternal CMV viremia was not associated with infant hospitalization and/or death by 24 months. CONCLUSIONS: Approximately 1 in 6 HIV-positive women in Botswana had detectable CMV DNA in blood at delivery. The presence of maternal CMV viremia had a borderline association with adverse birth outcomes but not with 24-month morbidity or mortality in HIV-exposed uninfected children.

Immunological non-response and low hemoglobin levels are predictors of incident tuberculosis among HIV-infected individuals on Truvada-based therapy in Botswana.

BACKGROUND: There is a high burden of tuberculosis (TB) in HIV antiretroviral programmes in Africa. However, few studies have looked at predictors of incident TB while on Truvada-based combination antiretroviral therapy (cART) regimens. METHODS: We estimated TB incidence among individuals enrolled into an observational cohort evaluating the efficacy and tolerability of Truvada-based cART in Gaborone, Botswana between 2008 and 2011. We used Cox proportional hazards regressions to determine predictors of incident TB. RESULTS: Of 300 participants enrolled, 45 (15%) had a diagnosis of TB at baseline. During 428 person-years (py) of follow-up, the incidence rate of TB was 3.04/100py (95% CI, 1.69-5.06), with 60% of the cases occurring within 3 months of ART initiation. Incident cases had low baseline CD4+ T cell counts (153cells/mm3 [Q1, Q3: 82, 242]; p = 0.69) and hemoglobin levels (9.2g/dl [Q1, Q3: 8.5,10.1]; p<0.01). In univariate analysis, low BMI (HR = 0.73; 95% CI 0.58-0.91; p = 0.01) and hemoglobin levels <8 g/dl (HR = 10.84; 95%CI: 2.99-40.06; p<0.01) were risk factors for TB. Time to incident TB diagnosis was significantly reduced in patients with poor immunological recovery (p = 0.04). There was no association between baseline viral load and risk of TB (HR = 1.75; 95%CI: 0.70-4.37). CONCLUSION: Low hemoglobin levels prior to initiation of ART are significant predictors of incident tuberculosis. Therefore, there is potential utility of iron biomarkers to identify patients at risk of TB prior to initiation on ART. Furthermore, additional strategies are required for patients with poor immunological recovery to reduce excess risk of TB while on ART.

Correction: Immunological non-response and low hemoglobin levels are predictors of incident tuberculosis among HIV-infected individuals on Truvada-based therapy in Botswana.

[This corrects the article DOI: 10.1371/journal.pone.0192030.].

Chronic and Occult Hepatitis B Virus Infection in Pregnant Women in Botswana.

The hepatitis B virus (HBV) is a global problem; however, the burden of HBV infection in pregnant women in Botswana is unknown. We sought to determine the prevalence of chronic and occult HBV infection in human immunodeficiency virus (HIV)-infected and -uninfected pregnant women in Botswana. Samples from 752 pregnant women were tested for hepatitis B surface antigen (HBsAg), and HBsAg-positive samples were tested for hepatitis B e antigen (HBeAg) and HBV DNA load. Samples that were HBsAg negative were screened for occult HBV infection by determining the HBV DNA load. HBV genotypes were determined based on a 415-base-pair fragment of the surface gene. Among the 752 women tested during pregnancy or early postpartum, 16 (2.1%) (95% confidence interval (CI): 2.0⁻2.2) were HBsAg-positive. The prevalence of chronic HBV infection was higher (3.1%) among HIV-infected (95% CI: 3.0⁻3.2) compared with HIV-uninfected women (1.1%) (95% CI: 1.07⁻1.1, p = 0.057). Among the 622 HBsAg-negative women, the prevalence of occult HBV infection was 6.6% (95% CI: 6.5⁻6.7). Three of thirteen HBsAg-positive participants were HBeAg-positive, and all were HIV-negative. Of the 11 maternal samples successfully genotyped, five (45.5%) were genotype D3, five (45.5%) were genotype A1, and one was genotype E (9%). Low and similar proportions of HIV-infected and -uninfected pregnant women in Botswana had occult or chronic HBV infection. We identified a subset of HIV-negative pregnant women who had high HBV DNA levels and were HBeAg-positive, and thus likely to transmit HBV to their infants.