Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa.

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

APOL1 Bi- and Monoallelic Variants and Chronic Kidney Disease in West Africans.

BACKGROUND: Apolipoprotein L1 gene (APOL1) variants are risk factors for chronic kidney disease (CKD) among Black Americans. Data are sparse on the genetic epidemiology of CKD and the clinical association of APOL1 variants with CKD in West Africans, a major group in the Black population. METHODS: We conducted a case-control study involving participants from Ghana and Nigeria who had CKD stages 2 through 5, biopsy-proven glomerular disease, or no kidney disease. We analyzed the association of CKD with APOL1 variants among participants with high-risk genotypes (two APOL1 risk alleles) and those with low-risk genotypes (fewer than two APOL1 risk alleles) by fitting logistic-regression models that controlled for covariates, including clinical site, age, and sex. RESULTS: Among 8355 participants (4712 with CKD stages 2 through 5, 866 with glomerular diseases, and 2777 with no kidney disease), the prevalence of monoallelic APOL1 variants was 43.0% and that of biallelic APOL1 variants was 29.7%. Participants with two APOL1 risk alleles had higher odds of having CKD than those with one risk allele or no risk alleles (adjusted odds ratio, 1.25; 95% confidence interval [CI], 1.11 to 1.40), as well as higher odds of focal segmental glomerulosclerosis (adjusted odds ratio, 1.84; 95% CI, 1.30 to 2.61). Participants with one APOL1 risk allele had higher odds of having CKD than those with no risk alleles (adjusted odds ratio, 1.18; 95% CI, 1.04 to 1.33), as well as higher odds of focal segmental glomerulosclerosis (adjusted odds ratio, 1.61; 95% CI, 1.04 to 2.48). The inclusion of covariates did not modify the association of monoallelic and biallelic APOL1 variants with CKD or focal segmental glomerulosclerosis. CONCLUSIONS: In this study, monoallelic APOL1 variants were associated with 18% higher odds of CKD and 61% higher odds of focal segmental glomerulosclerosis; biallelic APOL1 variants were associated with 25% higher odds of CKD and 84% higher odds of focal segmental glomerulosclerosis. (Funded by the National Human Genome Research Institute and others.).

High-depth African genomes inform human migration and health.

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed1. Here we performed whole-genome sequencing analyses of 426 individuals-comprising 50 ethnolinguistic groups, including previously unsampled populations-to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon-but in other genes, variants denoted as 'likely pathogenic' in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health.

Shared genetic risk between major orofacial cleft phenotypes in an African population.

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

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.

Perceptions and beliefs of community gatekeepers about genomic risk information in African cleft research.

BACKGROUND: A fundamental ethical issue in African genomics research is how socio-cultural factors impact perspectives, acceptance, and utility of genomic information, especially in stigmatizing conditions like orofacial clefts (OFCs). Previous research has shown that gatekeepers (e.g., religious, political, family or community leaders) wield considerable influence on the decision-making capabilities of their members, including health issues. Thus, their perspectives can inform the design of engagement strategies and increase exposure to the benefits of genomics testing/research. This is especially important for Africans underrepresented in genomic research. Our study aims to investigate the perspectives of gatekeepers concerning genomic risk information (GRI) in the presence of OFCs in a sub-Saharan African cohort. METHODS: Twenty-five focus group discussions (FGDs) consisting of 214 gatekeepers (religious, community, ethnic leaders, and traditional birth attendants) in Lagos, Nigeria, explored the opinions of participants on genomic risk information (GRI), OFC experience, and the possibility of involvement in collaborative decision-making in Lagos, Nigeria. Transcripts generated from audio recordings were coded and analyzed in NVivo using thematic analysis. RESULTS: Three main themes-knowledge, beliefs, and willingness to act-emerged from exploring the perspective of gatekeepers about GRI in this group. We observed mixed opinions regarding the acceptance of GRI. Many participants believed their role is to guide and support members when they receive results; this is based on the level of trust their members have in them. However, participants felt they would need to be trained by medical experts to do this. Also, religious and cultural beliefs were crucial to determining participants' understanding of OFCs and the acceptance and utilization of GRI. CONCLUSIONS: Incorporating cultural sensitivity into public engagement could help develop appropriate strategies to manage conflicting ideologies surrounding genomic information in African communities. This will allow for more widespread access to the advances in genomics research in underrepresented populations. We also recommend a synergistic relationship between community health specialists/scientists, and community leaders, including spiritual providers to better understand and utilize GRI.

Evolutionary forces in diabetes and hypertension pathogenesis in Africans.

Rates of type 2 diabetes (T2D) and hypertension are increasing rapidly in urbanizing sub-Saharan Africa (SSA). While lifestyle factors drive the increases in T2D and hypertension prevalence, evidence across populations shows that genetic variation, which is driven by evolutionary forces including a natural selection that shaped the human genome, also plays a role. Here we report the evidence for the effect of selection in African genomes on mechanisms underlying T2D and hypertension, including energy metabolism, adipose tissue biology, insulin action and salt retention. Selection effects found for variants in genes PPARA and TCF7L2 may have enabled Africans to respond to nutritional challenges by altering carbohydrate and lipid metabolism. Likewise, African-ancestry-specific characteristics of adipose tissue biology (low visceral adipose tissue [VAT], high intermuscular adipose tissue and a strong association between VAT and adiponectin) may have been selected for in response to nutritional and infectious disease challenges in the African environment. Evidence for selection effects on insulin action, including insulin resistance and secretion, has been found for several genes including MPHOSPH9, TMEM127, ZRANB3 and MC3R. These effects may have been historically adaptive in critical conditions, such as famine and inflammation. A strong correlation between hypertension susceptibility variants and latitude supports the hypothesis of selection for salt retention mechanisms in warm, humid climates. Nevertheless, adaptive genomics studies in African populations are scarce. More work is needed, particularly genomics studies covering the wide diversity of African populations in SSA and Africans in diaspora, as well as further functional assessment of established risk loci.

GWAS in Africans identifies novel lipids loci and demonstrates heterogenous association within Africa.

Serum lipids are biomarkers of cardiometabolic disease risk, and understanding genomic factors contributing to their distribution is of interest. Studies of lipids in Africans are rare, though it is expected that such studies could identify novel loci. We conducted a GWAS of 4317 Africans enrolled from Nigeria, Ghana and Kenya. We evaluated linear mixed models of high-density lipoprotein cholesterol (HDLC), low-density lipoprotein cholesterol (LDLC), total cholesterol (CHOL), triglycerides (TG) and TG/HDLC. Replication was attempted in 9542 African Americans (AA). In our main analysis, we identified 28 novel associations in Africans. Of the 18 of these that could be tested in AA, three associations replicated (GPNMB-TG, ENPP1-TG and SMARCA4-LDLC). Five additional novel loci were discovered upon meta-analysis with AA (rs138282551-TG, PGBD5-HDLC, CD80-TG/HDLC, SLC44A1-CHOL and TLL2-CHOL). Analyses considering only those with predominantly West African ancestry (Nigeria, Ghana and AA) yielded new insights: ORC5-LDLC and chr20:60973327-CHOL. Among our novel findings are some loci with known connections to lipids pathways. For instance, rs147706369 (TLL2) alters a regulatory motif for sterol regulatory element-binding proteins, a family of transcription factors that control the expression of a range of enzymes involved in cholesterol, fatty acid and TG synthesis, and rs115749422 (SMARCA4), an independent association near the known LDLR locus that is rare or absent in populations without African ancestry. These findings demonstrate the utility of conducting genomic analyses in Africans for discovering novel loci and provide some preliminary evidence for caution against treating 'African ancestry' as a monolithic category.

Hajdu-Cheney syndrome with atypical cardiovascular abnormalities.

Hajdu-Cheney syndrome is an ultra-rare autosomal dominant disorder caused by a heterozygous variant in NOTCH2 gene. Characteristic features include osteolysis, distinct facial appearance, skull deformity, joint laxity, osteoporosis, and short stature. Associated abnormalities are congenital heart disease, congenital defects of the kidney, and neurological problems. Here, we present the first reported case of an African child with a variant in NOTCH2 gene and features of Hajdu-Cheney syndrome in whom we detected a congenital heart defect that has not been previously reported in association with the syndrome. To appropriately characterize this disease and document correct proportion of cardiovascular malformation associations, echocardiography is recommended for all cases of Hajdu Cheney syndrome.

Ambiguous genitalia, giant congenital melanocytic nevus and subpulmonic outlet ventricular septal defect in an African child with Neurofibromatosis 1.

Neurofibromatosis type 1 is an autosomal dominant multisystemic disease caused by mutation of the neurofibromin (NF1) gene located on chromosome 17q11. We report a case of Neurofibromatosis 1 with ambiguous genitalia, giant congenital melanocytic nevus, and associated subpulmonic outlet ventricular septal defect, hitherto unreported in sub-Saharan Africa. In addition, a literature review of congenital heart diseases associated with Neurofibromatosis 1 is presented.

HLA and autoantibodies define scleroderma subtypes and risk in African and European Americans and suggest a role for molecular mimicry.

Systemic sclerosis (SSc) is a clinically heterogeneous autoimmune disease characterized by mutually exclusive autoantibodies directed against distinct nuclear antigens. We examined HLA associations in SSc and its autoantibody subsets in a large, newly recruited African American (AA) cohort and among European Americans (EA). In the AA population, the African ancestry-predominant HLA-DRB1*08:04 and HLA-DRB1*11:02 alleles were associated with overall SSc risk, and the HLA-DRB1*08:04 allele was strongly associated with the severe antifibrillarin (AFA) antibody subset of SSc (odds ratio = 7.4). These African ancestry-predominant alleles may help explain the increased frequency and severity of SSc among the AA population. In the EA population, the HLA-DPB1*13:01 and HLA-DRB1*07:01 alleles were more strongly associated with antitopoisomerase (ATA) and anticentromere antibody-positive subsets of SSc, respectively, than with overall SSc risk, emphasizing the importance of HLA in defining autoantibody subtypes. The association of the HLA-DPB1*13:01 allele with the ATA+ subset of SSc in both AA and EA patients demonstrated a transancestry effect. A direct correlation between SSc prevalence and HLA-DPB1*13:01 allele frequency in multiple populations was observed (r = 0.98, P = 3 × 10-6). Conditional analysis in the autoantibody subsets of SSc revealed several associated amino acid residues, mostly in the peptide-binding groove of the class II HLA molecules. Using HLA α/ß allelic heterodimers, we bioinformatically predicted immunodominant peptides of topoisomerase 1, fibrillarin, and centromere protein A and discovered that they are homologous to viral protein sequences from the Mimiviridae and Phycodnaviridae families. Taken together, these data suggest a possible link between HLA alleles, autoantibodies, and environmental triggers in the pathogenesis of SSc.

Refining genome-wide associated loci for serum uric acid in individuals with African ancestry.

OBJECTIVE: Serum uric acid is the end-product of purine metabolism and at high levels is a risk factor for several human diseases including gout and cardiovascular disease. Heritability estimates range from 0.32 to 0.63. Genome-wide association studies (GWAS) provide an unbiased approach to identify loci influencing serum uric acid. Here, we performed the first GWAS for serum uric acid in continental Africans, with replication in African Americans. METHODS: Africans (n = 4126) and African Americans (n = 5007) were genotyped on high-density GWAS arrays. Efficient mixed model association, a variance component approach, was used to perform association testing for a total of ~ 18 million autosomal genotyped and imputed variants. CAVIARBF was used to fine map significant regions. RESULTS: We identified two genome-wide significant loci: 4p16.1 (SLC2A9) and 11q13.1 (SLC22A12). At SLC2A9, the most strongly associated SNP was rs7683856 (P = 1.60 × 10-44). Conditional analysis revealed a second signal indexed by rs6838021 (P = 5.75 × 10-17). Gene expression and regulatory motif data prioritized a single-candidate causal variant for each signal. At SLC22A12, the most strongly associated SNP was rs147647315 (P = 6.65 × 10-25). Conditional analysis and functional annotation prioritized the missense variant rs147647315 (R (Arg) > H (His)) as the sole causal variant. Functional annotation of these three signals implicated processes in skeletal muscle, subcutaneous adipose tissue and the kidneys, respectively. CONCLUSIONS: This first GWAS of serum uric acid in continental Africans identified three associations at two loci, SLC2A9 and SLC22A12. The combination of weak linkage disequilibrium in Africans and functional annotation led to the identification of candidate causal SNPs for all three signals. Each candidate causal variant implicated a different cell type. Collectively, the three associations accounted for 4.3% of the variance of serum uric acid.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations.

OBJECTIVES: Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11. MATERIALS AND METHODS: We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools. RESULTS: Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11. CONCLUSION: This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P).

OBJECTIVE: Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset. METHODS: The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P. RESULTS: We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16. CONCLUSION: Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Damaging Mutations in AFDN Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate.

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

Plasma cell immunoglobulin heavy chain repertoire dynamics before and after tetanus booster vaccination.

Antibody sequence repertoire analysis of plasma cells (PC) isolated before and 1 week after a vaccine provides time-specific snapshots of the antibody response. Comparison of the immunoglobulin (Ig) sequences pre- and post-vaccination allows analysis of maturation over time and identification of antigen specific Ig. Here we compare the Ig heavy chain (Ig-H) repertoire of circulating PCs isolated from 109 peripheral blood mononuclear cells (PBMC) collected by apheresis 1 week after a tetanus toxoid vaccine booster with the Ig-H repertoire of PCs collected 2 and 11 weeks prior to the booster. A total of 21,060 unique Ig nucleotide sequences encoding 14,307 unique heavy chain complementarity determining region 3 (CDR-H3) amino acid sequences, also called clonotypes, were identified. Only 466 clonotypes (3.3%) were present at all 3 time points. In contrast, 90% of the 30 highest frequency CDR-H3 regions at +1w were also identified at another time point and 50% were present at all time points, suggesting the rapid expansion of a memory B cell population. The tetanus toxoid specificity of the CDR-H3 region with the 7th highest frequency at +1w was confirmed using immunoprecipitation and mass spectroscopy, and two public tetanus toxoid-specific CDR-H3 regions were also overrepresented at +1w. In summary, we have used the tetanus vaccine model system to demonstrate that bulk PC Ig repertoire analysis can identify PC populations that expand and mature following antigen exposure. The application of this approach before and after clinical infections should advance our understanding of clinical protection and facilitate vaccine design.

Prevalence and predictors of vitamin D deficiency in young African children.

BACKGROUND: Children living in sub-Saharan Africa have a high burden of rickets and infectious diseases, conditions that are linked to vitamin D deficiency. However, data on the vitamin D status of young African children and its environmental and genetic predictors are limited. We aimed to examine the prevalence and predictors of vitamin D deficiency in young African children. METHODS: We measured 25-hydroxyvitamin D (25(OH)D) and typed the single nucleotide polymorphisms, rs4588 and rs7041, in the GC gene encoding the vitamin D binding protein (DBP) in 4509 children aged 0-8 years living in Kenya, Uganda, Burkina Faso, The Gambia and South Africa. We evaluated associations between vitamin D status and country, age, sex, season, anthropometric indices, inflammation, malaria and DBP haplotypes in regression analyses. RESULTS: Median age was 23.9 months (interquartile range [IQR] 12.3, 35.9). Prevalence of vitamin D deficiency using 25(OH)D cut-offs of < 30 nmol/L and < 50 nmol/L was 0.6% (95% CI 0.4, 0.9) and 7.8% (95% CI 7.0, 8.5), respectively. Overall median 25(OH)D level was 77.6 nmol/L (IQR 63.6, 94.2). 25(OH)D levels were lower in South Africa, in older children, during winter or the long rains, and in those with afebrile malaria, and higher in children with inflammation. 25(OH)D levels did not vary by stunting, wasting or underweight in adjusted regression models. The distribution of Gc variants was Gc1f 83.3%, Gc1s 8.5% and Gc2 8.2% overall and varied by country. Individuals carrying the Gc2 variant had lower median 25(OH)D levels (72.4 nmol/L (IQR 59.4, 86.5) than those carrying the Gc1f (77.3 nmol/L (IQR 63.5, 92.8)) or Gc1s (78.9 nmol/L (IQR 63.8, 95.5)) variants. CONCLUSIONS: Approximately 0.6% and 7.8% of young African children were vitamin D deficient as defined by 25(OH)D levels < 30 nmol/L and < 50 nmol/L, respectively. Latitude, age, season, and prevalence of inflammation and malaria should be considered in strategies to assess and manage vitamin D deficiency in young children living in Africa.