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

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.

Host genetic variation at a locus near CHD1L impacts HIV sequence diversity in a South African population

IMPORTANCE

It has been previously shown that genetic variants near CHD1L on chromosome 1 are associated with reduced HIV VL in African populations. However, the impact of these variants on viral diversity and how they restrict viral replication are unknown. We report on a regional association analysis in a South African population and show evidence of selective pressure by variants near CHD1L on HIV RT and gag. Our findings provide further insight into how genetic variability at this locus contributes to host control of HIV in a South African population.

Direct inference and control of genetic population structure from RNA sequencing data

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

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

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.

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.

Promoter Hypermethylation Analysis of Host Genes in Cervical Cancer Patients With and Without Human Immunodeficiency Virus in Botswana

Background: Epidemics of human immunodeficiency virus (HIV) and cervical cancer are interconnected. DNA hypermethylation of host genes' promoter in cervical lesions has also been recognized as a contributor to cervical cancer progression.

Methods: For this purpose we analyzed promoter methylation of four tumor suppressor genes (RARB, CADM1, DAPK1 and PAX1) and explored their possible association with cervical cancer in Botswana among women of known HIV status. Overall, 228 cervical specimens (128 cervical cancers and 100 non-cancer subjects) were used. Yates-corrected chi-square test and Fisher's exact test were used to explore the association of promoter methylation for each host gene and cancer status. Subsequently, a logistic regression analysis was performed to find which factors, HIV status, high risk-HPV genotypes, patient's age and promoter methylation, were associated with the following dependent variables: cancer status, cervical cancer stage and promoter methylation rate.

Results: In patients with cervical cancer the rate of promoter methylation observed was greater than 64% in all the genes studied. Analysis also showed a higher risk of cervical cancer according to the increased number of methylated promoter genes (OR = 6.20; 95% CI: 3.66–10.51; P < 0.001). RARB methylation showed the strongest association with cervical cancer compared to other genes (OR = 15.25; 95% CI: 6.06–40.0; P < 0.001). Cervical cancer and promoter methylation of RARB and DAPK1 genes were associated with increasing age (OR = 1.12; 95% CI: 1.01-1.26; P = 0.037 and OR = 1.05; 95% CI: 1.00-1.10; P = 0.040). The presence of epigenetic changes at those genes appeared to be independent of HIV status among subjects with cervical cancer. Moreover, we found that cervical cancer stage was influenced by RARB (χ²= 7.32; P = 0.002) and CADM1 (χ²=12.68; P = 0.013) hypermethylation, and HIV status (χ²= 19.93; P = 0.001).

Conclusion: This study confirms the association between invasive cervical cancer and promoter gene methylation of tumor suppressing genes at the site of cancer. HIV infection did not show any association to methylation changes in this group of cervical cancer patients from Botswana. Further studies are needed to better understand the role of HIV in methylation of host genes among cancer subjects leading to cervical cancer progression.

Globally Rare BRCA2 Variants With Founder Haplotypes in the South African Population: Implications for Point-of-Care Testing Based on a Single-Institution BRCA1/2 Next-Generation Sequencing Study

Breast cancer patients historically benefitted from population-based genetic research performed in South Africa, which led to the development of founder-based BRCA1/2 diagnostic tests. With the advent of next-generation sequencing (NGS) technologies, the clinical utility of limited, targeted genetic assays were questioned. The study focused on mining NGS data obtained from an extensive single-institution NGS series (n=763). The aims were to determine (i) the prevalence of the most common recurrent/founder variants in patients referred for NGS directly; and (ii) to explore the data for inferred haplotypes associated with previous and potential new recurrent/founder variants. The identification of additional founder variants was essential for promoting and potentially advancing to rapid founder-based BRCA1/2 point-of-care (POC) technology as a time- and cost-effective alternative. NGS revealed actionable BRCA1/2 variants in 11.1% of patients tested (BRCA1 - 4.7%; BRCA2 - 6.4%), of which 22.4% represented variants currently screened for using first-tier targeted genetic testing. A retrospective investigation into the overall mutation-positive rate for an extended cohort (n=1906), which included first-tier test results, revealed that targeted genetic testing identified 74% of all pathogenic variants. This percentage justified the use of targeted genetic testing as a first-tier assay. Inferred haplotype analysis confirmed the founder status of BRCA2 c.5771_5774del (rs80359535) and c.7934del (rs80359688) and revealed an additional African founder variant (BRCA2 c.582G>A - rs80358810). A risk-benefit analysis using a questionnaire-based survey was performed in parallel to determine genetic professionals' views regarding POC testing. This was done to bridge the clinical implementation gap between haplotype analysis and POC testing as a first-tier screen during risk stratification of breast and ovarian cancer patients. The results reflected high acceptance (94%) of BRCA1/2 POC testing when accompanied by genetic counselling. Establishing the founder status for several recurrent BRCA2 variants across ethnic groups supports unselected use of the BRCA POC assay in all SA breast/ovarian cancer patients by recent local and international public health recommendations. Incorporating POC genotyping into the planned NGS screening algorithm of the Department of Health will ensure optimal use of the country's recourses to adhere to the set standards for optimal care and management for all breast cancer patients.

Past and current biological factors affecting malaria in the low transmission setting of Botswana: A review

Malaria continues to be one of the top infectious agents contributing to morbidity and mortality in sub-Saharan Africa. Annually, Botswana accounts only for a small proportion of cases (<<1%). Despite significantly reduced incidence rate, the country still experiences sporadic outbreaks that hamper the goal of malaria elimination. This review evaluated previous and current biological factors that impact malaria in Botswana, specifically focussing on the vectors, the parasite and the host. This was accomplished via a literature review evaluating these variables in Botswana. Current literature suggests that Anopheles arabiensis is the main malaria vector in the country. Several other potential vectors have been found widely distributed throughout Botswana in high numbers, yet remain largely unstudied with regards to their contribution to the country's malaria burden. We also report the most up to date list of all Anopheles species that have been found in Botswana. Plasmodium falciparum is responsible for the vast majority of symptomatic malaria in the country and some drug resistance markers have been documented for this species. Plasmodium vivax has been reported in asymptomatic subjects, even though a large proportion of the Botswana population appears to be Duffy antigen negative. Very little is known about the true distribution of P. vivax and no point of care testing infrastructure for this species exists in Botswana, making it difficult to tailor treatment to address possible recrudescence or relapse. Due to a genetically diverse population with a substantial Khoisan contribution into the Bantu genetic background, several phenotypes that potentially impact prevalence and severity of malaria exist within the country. These include sickle cell trait, Glucose-6-Phosphate Dehydrogenase deficiency, and Duffy negativity. This review highlights the information that currently exists on malaria in Botswana. It also postulates that a comprehensive understanding of these aforementioned biological factors may help to explain malaria persistence in Botswana.

The genomic prehistory of peoples speaking Khoisan languages

Peoples speaking so-called Khoisan languages-that is, indigenous languages of southern Africa that do not belong to the Bantu family-are culturally and linguistically diverse. They comprise herders, hunter-gatherers as well as groups of mixed modes of subsistence, and their languages are classified into three distinct language families. This cultural and linguistic variation is mirrored by extensive genetic diversity. We here review the recent genomics literature and discuss the genetic evidence for a formerly wider geographic spread of peoples with Khoisan-related ancestry, for the deep divergence among populations speaking Khoisan languages overlaid by more recent gene flow among these groups and for the impact of admixture with immigrant food-producers in their prehistory.

WITHDRAWN: Can Traditional Chinese Medicine provide insights into controlling the COVID-19 pandemic: Serpentinization-induced lithospheric long-wavelength magnetic anomalies in Proterozoic bedrocks in a weakened geomagnetic field mediate the aberrant transformation of biogenic molecules in COVID-19 via magnetic catalysis

This article has been withdrawn at the request of the authors and the editors. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/our-business/policies/article-withdrawal.

Ethnicity, age and disease-associated variation in body fluid-specific CpG sites in a diverse South African cohort

Tissue-specific differential DNA methylation has been an attractive target for the development of markers for discrimination of body fluids found at crime scenes. Though mostly stable, DNA methylation patterns have been shown to vary between different ethnic groups, in different age groups as well as between healthy and diseased individuals. To the best of our knowledge, none of the markers for body fluid identification have been applied to different ethnic groups to ascertain if variability exists. In the present study, saliva and blood were collected to determine the effects of ethnicity (Blacks, Whites, Coloureds and Indians), age (20-30 years, 40-50years and above 60 years) and diabetes on methylation profiles of potential saliva- and blood-specific DMSs. Both DMSs were previously shown to exhibit hypermethylation in their target body fluids at single CpG sites, however in the present study, additional CpG sites flanking the reported sites were also screened. Bisulfite sequencing revealed that Coloureds showed highest methylation levels for both body fluids, and blacks displayed significant differences between other ethnic groups in the blood-specific CpG sites. A decline in methylation for both potential DMRs was observed with increasing age. Heavily methylated CpG sites in different ethnic groups and previously reported DMSs displayed hypomethylation with increasing age and disease status. Diabetic status did not show any significant difference in methylation when compared to healthy counterparts. Thus, the use of methylation markers for forensics needs thorough investigation of influence of external factors and ideally, several CpG sites should be co-analysed instead of a single DMS.

Human leukocyte antigen class I (A, B, C) and class II (DPB1, DQB1, DRB1) allele and haplotype variation in Black South African individuals

South Africa has a population of 58.78 million, of which 80.7% are Black African individuals, representing 9 predominant ethnic/linguistic groups (Zulu, Xhosa, Pedi, Tswana, South Sotho, Tsonga, Swati, Venda and Ndebele). HIV-1 and Mycobacterium tuberculosis infection are the leading causes of death (7.8% and 5.9%, respectively) in this population group. To provide reference HLA allele and haplotype data for studies of gene-associations with infectious/non-infectious diseases or vaccine development, we have updated previously published HLA class I (A, B, C) and class II DRB1 genotypes and determined high-resolution class II (DPB1, DQB1) genotypes for n = 142 healthy, unrelated Black South African individuals.