Genetic signatures of gene flow and malaria-driven natural selection in sub-Saharan populations of the "endemic Burkitt Lymphoma belt"

Human leukocyte antigen-DQA1*04:01 and rs2040406 variants are associated with elevated risk of childhood Burkitt lymphoma

Burkitt lymphoma (BL) is responsible for many childhood cancers in sub-Saharan Africa, where it is linked to recurrent or chronic infection by Epstein-Barr virus or Plasmodium falciparum. However, whether human leukocyte antigen (HLA) polymorphisms, which regulate immune response, are associated with BL has not been well investigated, which limits our understanding of BL etiology. Here we investigate this association among 4,645 children aged 0-15 years, 800 with BL, enrolled in Uganda, Tanzania, Kenya, and Malawi. HLA alleles are imputed with accuracy >90% for HLA class I and 85-89% for class II alleles. BL risk is elevated with HLA-DQA1*04:01 (adjusted odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.32-1.97, P = 3.71 × 10-6), with rs2040406(G) in HLA-DQA1 region (OR = 1.43, 95% CI = 1.26-1.63, P = 4.62 × 10-8), and with amino acid Gln at position 53 versus other variants in HLA-DQA1 (OR = 1.36, P = 2.06 × 10-6). The associations with HLA-DQA1*04:01 (OR = 1.29, P = 0.03) and rs2040406(G) (OR = 1.68, P = 0.019) persist in mutually adjusted models. The higher risk rs2040406(G) variant for BL is associated with decreased HLA-DQB1 expression in eQTLs in EBV transformed lymphocytes. Our results support the role of HLA variation in the etiology of BL and suggest that a promising area of research might be understanding the link between HLA variation and EBV control.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Mosaic chromosomal alterations in peripheral blood leukocytes of children in sub-Saharan Africa

In high-income countries, mosaic chromosomal alterations in peripheral blood leukocytes are associated with an elevated risk of adverse health outcomes, including hematologic malignancies. We investigate mosaic chromosomal alterations in sub-Saharan Africa among 931 children with Burkitt lymphoma, an aggressive lymphoma commonly characterized by immunoglobulin-MYC chromosomal rearrangements, 3822 Burkitt lymphoma-free children, and 674 cancer-free men from Ghana. We find autosomal and X chromosome mosaic chromosomal alterations in 3.4% and 1.7% of Burkitt lymphoma-free children, and 8.4% and 3.7% of children with Burkitt lymphoma (P-values = 5.7×10-11 and 3.74×10-2, respectively). Autosomal mosaic chromosomal alterations are detected in 14.0% of Ghanaian men and increase with age. Mosaic chromosomal alterations in Burkitt lymphoma cases include gains on chromosomes 1q and 8, the latter spanning MYC, while mosaic chromosomal alterations in Burkitt lymphoma-free children include copy-neutral loss of heterozygosity on chromosomes 10, 14, and 16. Our results highlight mosaic chromosomal alterations in sub-Saharan African populations as a promising area of research.

Unappreciated subcontinental admixture in Europeans and European Americans and implications for genetic epidemiology studies

European-ancestry populations are recognized as stratified but not as admixed, implying that residual confounding by locus-specific ancestry can affect studies of association, polygenic adaptation, and polygenic risk scores. We integrate individual-level genome-wide data from ~19,000 European-ancestry individuals across 79 European populations and five European American cohorts. We generate a new reference panel that captures ancestral diversity missed by both the 1000 Genomes and Human Genome Diversity Projects. Both Europeans and European Americans are admixed at the subcontinental level, with admixture dates differing among subgroups of European Americans. After adjustment for both genome-wide and locus-specific ancestry, associations between a highly differentiated variant in LCT (rs4988235) and height or LDL-cholesterol were confirmed to be false positives whereas the association between LCT and body mass index was genuine. We provide formal evidence of subcontinental admixture in individuals with European ancestry, which, if not properly accounted for, can produce spurious results in genetic epidemiology studies.

Dense sampling of ethnic groups within African countries reveals fine-scale genetic structure and extensive historical admixture

Previous studies have highlighted how African genomes have been shaped by a complex series of historical events. Despite this, genome-wide data have only been obtained from a small proportion of present-day ethnolinguistic groups. By analyzing new autosomal genetic variation data of 1333 individuals from over 150 ethnic groups from Cameroon, Republic of the Congo, Ghana, Nigeria, and Sudan, we demonstrate a previously underappreciated fine-scale level of genetic structure within these countries, for example, correlating with historical polities in western Cameroon. By comparing genetic variation patterns among populations, we infer that many northern Cameroonian and Sudanese groups share genetic links with multiple geographically disparate populations, likely resulting from long-distance migrations. In Ghana and Nigeria, we infer signatures of intermixing dated to over 2000 years ago, corresponding to reports of environmental transformations possibly related to climate change. We also infer recent intermixing signals in multiple African populations, including Congolese, that likely relate to the expansions of Bantu language-speaking peoples.

Genetic variations in human ATP2B4 gene alter Plasmodium falciparum in vitro growth in RBCs from Gambian adults

BACKGROUND

Polymorphisms in ATP2B4 coding for PMCA4b, the primary regulator of erythrocyte calcium concentration, have been shown by GWAS and cross-sectional studies to protect against severe malaria but the mechanism remains unknown.

METHODS

Using a recall-by-genotype design, we investigated the impact of a common haplotype variant in ATP2B4 using in vitro assays that model erythrocyte stage malaria pathogenesis. Ninety-six donors representing homozygote (carriers of the minor allele, C/C), heterozygote (T/C) and wildtype (T/T) carriers of the tagging SNP rs1541252 were selected from a cohort of over 12,000 participants in the Keneba Biobank.

RESULTS

Red blood cells (RBCs) from homozygotes showed reduced PMCA4b protein expression (mean fluorescence intensities (MFI = 2428 ± 124, 3544 ± 159 and 4261 ± 283], for homozygotes, heterozygotes and wildtypes respectively, p < 0.0001) and slower rates of calcium expulsion (calcium t_½ ± SD = 4.7 ± 0.5, 1.8 ± 0.3 and 1.9 ± 0.4 min, p < 0.0001). Growth of a Plasmodium falciparum laboratory strain (FCR3) and two Gambian field isolates was decreased in RBCs from homozygotes compared to heterozygotes and wildtypes (p < 0.01). Genotype group did not affect parasite adhesion in vitro or var-gene expression in malaria-infected RBCs. Parasite growth was inhibited by a known inhibitor of PMCA4b, aurintricarboxylic acid (IC₅₀ = 122uM CI: 110-134) confirming its sensitivity to calcium channel blockade.

CONCLUSION

The data support the hypothesis that this ATP2B4 genotype, common in The Gambia and other malaria-endemic areas, protects against severe malaria through the suppression of parasitaemia during an infection. Reduction in parasite density plays a pivotal role in disease outcome by minimizing all aspects of malaria pathogenesis. Follow up studies are needed to further elucidate the mechanism of protection and to determine if this ATP2B4 genotype carries a fitness cost or increases susceptibility to other human disease.

Putative pathogen-selected polymorphisms in the PKLR gene are associated with mycobacterial susceptibility in Brazilian and African populations

Pyruvate kinase (PK), encoded by the PKLR gene, is a key player in glycolysis controlling the integrity of erythrocytes. Due to Plasmodium selection, mutations for PK deficiency, which leads to hemolytic anemia, are associated with resistance to malaria in sub-Saharan Africa and with susceptibility to intracellular pathogens in experimental models. In this case-control study, we enrolled 4,555 individuals and investigated whether PKLR single nucleotide polymorphisms (SNPs) putatively selected for malaria resistance are associated with susceptibility to leprosy across Brazil (Manaus-North; Salvador-Northeast; Rondonópolis-Midwest and Rio de Janeiro-Southeast) and with tuberculosis in Mozambique. Haplotype T/G/G (rs1052176/rs4971072/rs11264359) was associated with leprosy susceptibility in Rio de Janeiro (OR = 2.46, p = 0.00001) and Salvador (OR = 1.57, p = 0.04), and with tuberculosis in Mozambique (OR = 1.52, p = 0.07). This haplotype downregulates PKLR expression in nerve and skin, accordingly to GTEx, and might subtly modulate ferritin and haptoglobin levels in serum. Furthermore, we observed genetic signatures of positive selection in the HCN3 gene (xpEHH>2 -recent selection) in Europe but not in Africa, involving 6 SNPs which are PKLR/HCN3 eQTLs. However, this evidence was not corroborated by the other tests (FST, Tajima's D and iHS). Altogether, we provide evidence that a common PKLR locus in Africans contribute to mycobacterial susceptibility in African descent populations and also highlight, for first, PKLR as a susceptibility gene for leprosy and TB.

Modulation of the Human Erythroid Plasma Membrane Calcium Pump (PMCA4b) Expression by Polymorphic Genetic Variants

In the human ATP2B4 gene, coding for the plasma membrane calcium pump PMCA4b, a minor haplotype results in the decreased expression of this membrane protein in erythroid cells. The presence of this haplotype and the consequently reduced PMCA4b expression have been suggested to affect red blood cell hydration and malaria susceptibility. By using dual-luciferase reporter assays, we have localized the erythroid-specific regulatory region within the haplotype of the ATP2B4 gene, containing predicted GATA1 binding sites that are affected by SNPs in the minor haplotype. Our results show that, in human erythroid cells, the regulation of ATP2B4 gene expression is significantly affected by GATA1 expression, and we document the role of specific SNPs involved in predicted GATA1 binding. Our findings provide a mechanistic explanation at the molecular level for the reduced erythroid-specific PMCA4b expression in carriers of ATP2B4 gene polymorphic variants.

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host-parasite coevolutionary dynamics.

A Review of Cancer Genetics and Genomics Studies in Africa

Cancer is the second leading cause of death globally and is projected to overtake infectious disease as the leading cause of mortality in Africa within the next two decades. Cancer is a group of genomic diseases that presents with intra- and inter-population unique phenotypes, with Black populations having the burden of morbidity and mortality for most types. At large, the prevention and treatment of cancers have been propelled by the understanding of the genetic make-up of the disease of mostly non-African populations. By the same token, there is a wide knowledge gap in understanding the underlying genetic causes of, and genomic alterations associated with, cancer among black Africans. Accordingly, we performed a review of the literature to survey existing studies on cancer genetics/genomics and curated findings pertaining to publications across multiple cancer types conducted on African populations. We used PubMed MeSH terms to retrieve the relevant publications from 1990 to December 2019. The metadata of these publications were extracted using R text mining packages: RISmed and Pubmed.mineR. The data showed that only 0.329% of cancer publications globally were on Africa, and only 0.016% were on cancer genetics/genomics from Africa. Although the most prevalent cancers in Africa are cancers of the breast, cervix, uterus, and prostate, publications representing breast, colorectal, liver, and blood cancers were the most frequent in our review. The most frequently reported cancer genes were BRCA1, BRCA2, and TP53. Next, the genes reported in the reviewed publications’ abstracts were extracted and annotated into three gene ontology classes. Genes in the cellular component class were mostly associated with cell part and organelle part, while those in biological process and molecular function classes were mainly associated with cell process, biological regulation, and binding, and catalytic activity, respectively. Overall, this review highlights the paucity of research on cancer genomics on African populations, identified gaps, and discussed the need for concerted efforts to encourage more research on cancer genomics in Africa.

Pediatric Aggressive Mature B-Cell Lymphomas, Version 2.2020, NCCN Clinical Practice Guidelines in Oncology

Pediatric aggressive mature B-cell lymphomas are the most common types of non-Hodgkin lymphoma in children, and they include Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL). These diseases are highly aggressive but curable, the treatment is complex, and patients may have many complicated supportive care issues. The NCCN Guidelines for Pediatric Aggressive Mature B-Cell Lymphomas provide guidance regarding pathology and diagnosis, staging, initial treatment, disease reassessment, surveillance, therapy for relapsed/refractory disease, and supportive care for clinicians who treat sporadic pediatric BL and DLBCL.

Differences in the frequency of genetic variants associated with iron imbalance among global populations

Iron deficiency anaemia is a major health problem affecting approximately 1.2 billion people worldwide. Young children, women of reproductive age and pregnant women living in sub-Saharan Africa are the most vulnerable. It is estimated that iron deficiency accounts for half of anaemia cases. Apart from nutritional deficiency, infection, inflammation and genetic factors are the major drivers of anaemia. However, the role of genetic risk factors has not been thoroughly investigated. This is particularly relevant in African populations, as they carry high genetic diversity and have a high prevalence of anaemia. Multiple genetic variations in iron regulatory genes have been linked to impaired iron status. Here we conducted a literature review to identify genetic variants associated with iron imbalance among global populations. We compare their allele frequencies and risk scores and we investigated population-specific selection among populations of varying geographic origin using data from the Keneba Biobank representing individuals in rural Gambia and the 1000 Genomes Project. We identified a significant lack of data on the genetic determinants of iron status in sub-Saharan Africa. Most of the studies on genetic determinants of iron status have been conducted in Europeans. Also, we identified population differences in allele frequencies in candidate putative genetic risk factors. Given the disproportionately high genetic diversity in African populations coupled with their high prevalence of iron deficiency, there is need to investigate the genetic influences of low iron status in Sub-Saharan Africa. The resulting insights may inform the future implementation of iron intervention strategies.

Origins, Admixture Dynamics, and Homogenization of the African Gene Pool in the Americas

The Transatlantic Slave Trade transported more than 9 million Africans to the Americas between the early 16th and the mid-19th centuries. We performed a genome-wide analysis using 6,267 individuals from 25 populations to infer how different African groups contributed to North-, South-American, and Caribbean populations, in the context of geographic and geopolitical factors, and compared genetic data with demographic history records of the Transatlantic Slave Trade. We observed that West-Central Africa and Western Africa-associated ancestry clusters are more prevalent in northern latitudes of the Americas, whereas the South/East Africa-associated ancestry cluster is more prevalent in southern latitudes of the Americas. This pattern results from geographic and geopolitical factors leading to population differentiation. However, there is a substantial decrease in the between-population differentiation of the African gene pool within the Americas, when compared with the regions of origin from Africa, underscoring the importance of historical factors favoring admixture between individuals with different African origins in the New World. This between-population homogenization in the Americas is consistent with the excess of West-Central Africa ancestry (the most prevalent in the Americas) in the United States and Southeast-Brazil, with respect to historical-demography expectations. We also inferred that in most of the Americas, intercontinental admixture intensification occurred between 1750 and 1850, which correlates strongly with the peak of arrivals from Africa. This study contributes with a population genetics perspective to the ongoing social, cultural, and political debate regarding ancestry, admixture, and the mestizaje process in the Americas.

Mean platelet counts are relatively decreased with malaria but relatively increased with endemic Burkitt Lymphoma in Uganda, Tanzania, and Kenya

Platelet counts are decreased in Plasmodium falciparum malaria, which is aetiologically linked with endemic Burkitt lymphoma (eBL). However, the pattern of platelet counts in eBL cases is unknown. We studied platelet counts in 582 eBL cases and 2 248 controls enrolled in a case-control study in Uganda, Tanzania and Kenya (2010-2016). Mean platelet counts in controls or eBL cases with or without malaria-infection in controls versus eBLcases were compared using Student's t-test. Odds ratios (ORs) and two-sided 95% confidence intervals (95% CIs) were estimated using multiple logistic regression, controlling for age, sex, haemoglobin and white blood cell counts. Platelets were decreased with malaria infection in the controls [263 vs. 339 × 10⁹ platelets/l, P < 0·0001; adjusted OR (aOR) = 3·42, 95% CI: 2·79-4·18] and eBL cases (314 vs. 367 × 10⁹ platelets/l, P-value = 0·002; aOR = 2·36, 95% CI: 1·49-3·73). Unexpectedly, platelets were elevated in eBL cases versus  controls in overall analyses (mean: 353 vs. 307 × 10⁹ platelets/l, P < 0·0001; aOR = 1·41; 95% CI: 1·12-1·77), and when restricted to malaria-positive (mean 314 vs. 263 × 10⁹ platelets/l, P < 0·0001; OR = 2·26; 95% CI: 1·56-3·27) or malaria-negative (mean 367 vs. 339 × 10⁹ platelets/l, P < 0·001; OR = 1·46; 95% CI: 1·17-1·83) subjects. Platelets were decreased with malaria infection in controls and eBL cases but elevated with eBL.

Variation in the Human Leukocyte Antigen system and risk for endemic Burkitt lymphoma in northern Uganda

Endemic Burkitt lymphoma (eBL) is an aggressive childhood B-cell lymphoma associated with Plasmodium falciparum (Pf) malaria and Epstein-Barr virus (EBV) infections. Variation in the Human Leukocyte Antigen (HLA) system is suspected to play a role, but assessments using less accurate serology-based HLA typing techniques in small studies yielded conflicting results. We studied 200 eBL cases and 400 controls aged 0-15 years enrolled in northern Uganda and typed by accurate high-resolution HLA sequencing methods. HLA results were analyzed at one- or two-field resolution. Odds ratios and 95% confidence intervals (aOR, 95% CI) for eBL risk associated with common HLA alleles versus alleles that were rare (<1%) or differed by <2% between the cases and controls as the reference category, were estimated using multiple logistic regression adjusting for age, sex, microgeography, region, malaria positivity and treatment history, and genetic variants associated with eBL. Compared to the controls, eBL cases had a lower frequency of HLA-A*02 (aOR = 0·59, 95% CI 0·38-0·91), HLA-B*41 (aOR = 0·36, 95% CI 0·13-1·00), and HLA-B*58 alleles (aOR = 0·59, 95% CI 0·36-0·97). eBL cases had a lower frequency of HLA-DPB1 homozygosity (aOR = 0·57, 95% CI 0·40-0·82) but a higher frequency of HLA-DQA1 homozygosity (aOR = 2·19, 95% CI 1·42-3·37). Our results suggest that variation in HLA may be associated with eBL risk.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

Genetics of cognitive trajectory in Brazilians: 15 years of follow-up from the Bambuí-Epigen Cohort Study of Aging

Age-related cognitive decline (ACD) is the gradual process of decreasing of cognitive function over age. Most genetic risk factors for ACD have been identified in European populations and there are no reports in admixed Latin American individuals. We performed admixture mapping, genome-wide association analysis (GWAS), and fine-mapping to examine genetic factors associated with 15-year cognitive trajectory in 1,407 Brazilian older adults, comprising 14,956 Mini-Mental State Examination measures. Participants were enrolled as part of the Bambuí-Epigen Cohort Study of Aging. Our admixture mapping analysis identified a genomic region (3p24.2) in which increased Native American ancestry was significantly associated with faster ACD. Fine-mapping of this region identified a single nucleotide polymorphism (SNP) rs142380904 (β = -0.044, SE = 0.01, p = 7.5 × 10-5) associated with ACD. In addition, our GWAS identified 24 associated SNPs, most in genes previously reported to influence cognitive function. The top six associated SNPs accounted for 18.5% of the ACD variance in our data. Furthermore, our longitudinal study replicated previous GWAS hits for cognitive decline and Alzheimer's disease. Our 15-year longitudinal study identified both ancestry-specific and cosmopolitan genetic variants associated with ACD in Brazilians, highlighting the need for more trans-ancestry genomic studies, especially in underrepresented ethnic groups.