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.

Africa-specific human genetic variation near CHD1L associates with HIV-1 load.

HIV-1 remains a global health crisis1, highlighting the need to identify new targets for therapies. Here, given the disproportionate HIV-1 burden and marked human genome diversity in Africa2, we assessed the genetic determinants of control of set-point viral load in 3,879 people of African ancestries living with HIV-1 participating in the international collaboration for the genomics of HIV3. We identify a previously undescribed association signal on chromosome 1 where the peak variant associates with an approximately 0.3 log10-transformed copies per ml lower set-point viral load per minor allele copy and is specific to populations of African descent. The top associated variant is intergenic and lies between a long intergenic non-coding RNA (LINC00624) and the coding gene CHD1L, which encodes a helicase that is involved in DNA repair4. Infection assays in iPS cell-derived macrophages and other immortalized cell lines showed increased HIV-1 replication in CHD1L-knockdown and CHD1L-knockout cells. We provide evidence from population genetic studies that Africa-specific genetic variation near CHD1L associates with HIV replication in vivo. Although experimental studies suggest that CHD1L is able to limit HIV infection in some cell types in vitro, further investigation is required to understand the mechanisms underlying our observations, including any potential indirect effects of CHD1L on HIV spread in vivo that our cell-based assays cannot recapitulate.

The African Genome Variation Project shapes medical genetics in Africa.

Given the importance of Africa to studies of human origins and disease susceptibility, detailed characterization of African genetic diversity is needed. The African Genome Variation Project provides a resource with which to design, implement and interpret genomic studies in sub-Saharan Africa and worldwide. The African Genome Variation Project represents dense genotypes from 1,481 individuals and whole-genome sequences from 320 individuals across sub-Saharan Africa. Using this resource, we find novel evidence of complex, regionally distinct hunter-gatherer and Eurasian admixture across sub-Saharan Africa. We identify new loci under selection, including loci related to malaria susceptibility and hypertension. We show that modern imputation panels (sets of reference genotypes from which unobserved or missing genotypes in study sets can be inferred) can identify association signals at highly differentiated loci across populations in sub-Saharan Africa. Using whole-genome sequencing, we demonstrate further improvements in imputation accuracy, strengthening the case for large-scale sequencing efforts of diverse African haplotypes. Finally, we present an efficient genotype array design capturing common genetic variation in Africa.

Genome-wide association study of type 2 diabetes in Africa.

AIMS/HYPOTHESIS: Genome-wide association studies (GWAS) for type 2 diabetes have uncovered >400 risk loci, primarily in populations of European and Asian ancestry. Here, we aimed to discover additional type 2 diabetes risk loci (including African-specific variants) and fine-map association signals by performing genetic analysis in African populations. METHODS: We conducted two type 2 diabetes genome-wide association studies in 4347 Africans from South Africa, Nigeria, Ghana and Kenya and meta-analysed both studies together. Likely causal variants were identified using fine-mapping approaches. RESULTS: The most significantly associated variants mapped to the widely replicated type 2 diabetes risk locus near TCF7L2 (p = 5.3 × 10-13). Fine-mapping of the TCF7L2 locus suggested one type 2 diabetes association signal shared between Europeans and Africans (indexed by rs7903146) and a distinct African-specific signal (indexed by rs17746147). We also detected one novel signal, rs73284431, near AGMO (p = 5.2 × 10-9, minor allele frequency [MAF] = 0.095; monomorphic in most non-African populations), distinct from previously reported signals in the region. In analyses focused on 100 published type 2 diabetes risk loci, we identified 21 with shared causal variants in African and non-African populations. CONCLUSIONS/INTERPRETATION: These results demonstrate the value of performing GWAS in Africans, provide a resource to larger consortia for further discovery and fine-mapping and indicate that additional large-scale efforts in Africa are warranted to gain further insight in to the genetic architecture of type 2 diabetes.

Linear mixed model for heritability estimation that explicitly addresses environmental variation.

The linear mixed model (LMM) is now routinely used to estimate heritability. Unfortunately, as we demonstrate, LMM estimates of heritability can be inflated when using a standard model. To help reduce this inflation, we used a more general LMM with two random effects-one based on genomic variants and one based on easily measured spatial location as a proxy for environmental effects. We investigated this approach with simulated data and with data from a Uganda cohort of 4,778 individuals for 34 phenotypes including anthropometric indices, blood factors, glycemic control, blood pressure, lipid tests, and liver function tests. For the genomic random effect, we used identity-by-descent estimates from accurately phased genome-wide data. For the environmental random effect, we constructed a covariance matrix based on a Gaussian radial basis function. Across the simulated and Ugandan data, narrow-sense heritability estimates were lower using the more general model. Thus, our approach addresses, in part, the issue of "missing heritability" in the sense that much of the heritability previously thought to be missing was fictional. Software is available at https://github.com/MicrosoftGenomics/FaST-LMM.

Evaluating the Impact of Functional Genetic Variation on HIV-1 Control.

Background: Previous genetic association studies of human immunodeficiency virus-1 (HIV-1) progression have focused on common human genetic variation ascertained through genome-wide genotyping. Methods: We sought to systematically assess the full spectrum of functional variation in protein coding gene regions on HIV-1 progression through exome sequencing of 1327 individuals. Genetic variants were tested individually and in aggregate across genes and gene sets for an influence on HIV-1 viral load. Results: Multiple single variants within the major histocompatibility complex (MHC) region were observed to be strongly associated with HIV-1 outcome, consistent with the known impact of classical HLA alleles. However, no single variant or gene located outside of the MHC region was significantly associated with HIV progression. Set-based association testing focusing on genes identified as being essential for HIV replication in genome-wide small interfering RNA (siRNA) and clustered regularly interspaced short palindromic repeats (CRISPR) studies did not reveal any novel associations. Conclusions: These results suggest that exonic variants with large effect sizes are unlikely to have a major contribution to host control of HIV infection.

Multiple processes drive genetic structure of humpback whale (Megaptera novaeangliae) populations across spatial scales.

Elucidating patterns of population structure for species with complex life histories, and disentangling the processes driving such patterns, remains a significant analytical challenge. Humpback whale (Megaptera novaeangliae) populations display complex genetic structures that have not been fully resolved at all spatial scales. We generated a data set of nuclear markers for 3575 samples spanning the seven breeding stocks and substocks found in the South Atlantic and western and northern Indian Oceans. For the total sample, and males and females separately, we assessed genetic diversity, tested for genetic differentiation between putative populations and isolation by distance, estimated the number of genetic clusters without a priori population information and estimated rates of gene flow using maximum-likelihood and Bayesian approaches. At the ocean basin scale, structure is governed by geographical distance (IBD P < 0.05) and female fidelity to breeding areas, in line with current understanding of the drivers of broadscale population structure. Consistent with previous studies, the Arabian Sea breeding stock was highly genetically differentiated (FST 0.034-0.161; P < 0.01 for all comparisons). However, the breeding stock boundary between west South Africa and east Africa was more porous than expected based on genetic differentiation, cluster and geneflow analyses. Instances of male fidelity to breeding areas and relatively high rates of dispersal for females were also observed between the three substocks in the western Indian Ocean. The relationships between demographic units and current management boundaries may have ramifications for assessments of the status and continued protections of populations still in recovery from commercial whaling.

Global diversity and oceanic divergence of humpback whales (Megaptera novaeangliae).

Humpback whales (Megaptera novaeangliae) annually undertake the longest migrations between seasonal feeding and breeding grounds of any mammal. Despite this dispersal potential, discontinuous seasonal distributions and migratory patterns suggest that humpbacks form discrete regional populations within each ocean. To better understand the worldwide population history of humpbacks, and the interplay of this species with the oceanic environment through geological time, we assembled mitochondrial DNA control region sequences representing approximately 2700 individuals (465 bp, 219 haplotypes) and eight nuclear intronic sequences representing approximately 70 individuals (3700 bp, 140 alleles) from the North Pacific, North Atlantic and Southern Hemisphere. Bayesian divergence time reconstructions date the origin of humpback mtDNA lineages to the Pleistocene (880 ka, 95% posterior intervals 550-1320 ka) and estimate radiation of current Northern Hemisphere lineages between 50 and 200 ka, indicating colonization of the northern oceans prior to the Last Glacial Maximum. Coalescent analyses reveal restricted gene flow between ocean basins, with long-term migration rates (individual migrants per generation) of less than 3.3 for mtDNA and less than 2 for nuclear genomic DNA. Genetic evidence suggests that humpbacks in the North Pacific, North Atlantic and Southern Hemisphere are on independent evolutionary trajectories, supporting taxonomic revision of M. novaeangliae to three subspecies.

High-resolution African HLA resource uncovers HLA-DRB1 expression effects underlying vaccine response.

How human genetic variation contributes to vaccine effectiveness in infants is unclear, and data are limited on these relationships in populations with African ancestries. We undertook genetic analyses of vaccine antibody responses in infants from Uganda (n = 1391), Burkina Faso (n = 353) and South Africa (n = 755), identifying associations between human leukocyte antigen (HLA) and antibody response for five of eight tested antigens spanning pertussis, diphtheria and hepatitis B vaccines. In addition, through HLA typing 1,702 individuals from 11 populations of African ancestry derived predominantly from the 1000 Genomes Project, we constructed an imputation resource, fine-mapping class II HLA-DR and DQ associations explaining up to 10% of antibody response variance in our infant cohorts. We observed differences in the genetic architecture of pertussis antibody response between the cohorts with African ancestries and an independent cohort with European ancestry, but found no in silico evidence of differences in HLA peptide binding affinity or breadth. Using immune cell expression quantitative trait loci datasets derived from African-ancestry samples from the 1000 Genomes Project, we found evidence of differential HLA-DRB1 expression correlating with inferred protection from pertussis following vaccination. This work suggests that HLA-DRB1 expression may play a role in vaccine response and should be considered alongside peptide selection to improve vaccine design.

Lipoprotein(a) and risk of coronary, cerebrovascular, and peripheral artery disease: the EPIC-Norfolk prospective population study.

OBJECTIVE: Although the association between circulating levels of lipoprotein(a) [Lp(a)] and risk of coronary artery disease (CAD) and stroke is well established, its role in risk of peripheral arterial disease (PAD) remains unclear. Here, we examine the association between Lp(a) levels and PAD in a large prospective cohort. To contextualize these findings, we also examined the association between Lp(a) levels and risk of stroke and CAD and studied the role of low-density lipoprotein as an effect modifier of Lp(a)-associated cardiovascular risk. METHODS AND RESULTS: Lp(a) levels were measured in apparently healthy participants in the European Prospective Investigation of Cancer (EPIC)-Norfolk cohort. Cox regression was used to quantify the association between Lp(a) levels and risk of PAD, stroke, and CAD outcomes. During 212 981 person-years at risk, a total of 2365 CAD, 284 ischemic stroke, and 596 PAD events occurred in 18 720 participants. Lp(a) was associated with PAD and CAD outcomes but not with ischemic stroke (hazard ratio per 2.7-fold increase in Lp(a) of 1.37, 95% CI 1.25-1.50, 1.13, 95% CI 1.04-1.22 and 0.91, 95% CI 0.79-1.03, respectively). Low-density lipoprotein cholesterol levels did not modify these associations. CONCLUSIONS: Lp(a) levels were associated with future PAD and CAD events. The association between Lp(a) and cardiovascular disease was not modified by low-density lipoprotein cholesterol levels.

Targeting hepatitis B vaccine escape using immunogenetics in Bangladeshi infants.

Hepatitis B virus (HBV) vaccine escape mutants (VEM) are increasingly described, threatening progress in control of this virus worldwide. Here we studied the relationship between host genetic variation, vaccine immunogenicity and viral sequences implicating VEM emergence. In a cohort of 1,096 Bangladeshi children, we identified human leukocyte antigen (HLA) variants associated with response vaccine antigens. Using an HLA imputation panel with 9,448 south Asian individuals DPB1*04:01 was associated with higher HBV antibody responses (p=4.5×10-30). The underlying mechanism is a result of higher affinity binding of HBV surface antigen epitopes to DPB1*04:01 dimers. This is likely a result of evolutionary pressure at the HBV surface antigen 'a-determinant' segment incurring VEM specific to HBV. Prioritizing pre-S isoform HBV vaccines may tackle the rise of HBV vaccine evasion.

Transcriptomics and chromatin accessibility in multiple African population samples.

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

Effectiveness of scat-detection dogs in determining species presence in a tropical savanna landscape.

Most protected areas are too small to sustain populations of wide-ranging mammals; thus, identification and conservation of high-quality habitat for those animals outside parks is often a high priority, particularly for regions where extensive land conversion is occurring. This is the case in the vicinity of Emas National Park, a small protected area in the Brazilian Cerrado. Over the last 40 years the native vegetation surrounding the park has been converted to agriculture, but the region still supports virtually all of the animals native to the area. We determined the effectiveness of scat-detection dogs in detecting presence of five species of mammals threatened with extinction by habitat loss: maned wolf (Chrysocyon brachyurus), puma (Puma concolor), jaguar (Panthera onca), giant anteater (Myrmecophaga tridactyla), and giant armadillo (Priodontes maximus). The probability of scat detection varied among the five species and among survey quadrats of different size, but was consistent across team, season, and year. The probability of occurrence, determined from the presence of scat, in a randomly selected site within the study area ranged from 0.14 for jaguars, which occur primarily in the forested areas of the park, to 0.91 for maned wolves, the most widely distributed species in our study area. Most occurrences of giant armadillos in the park were in open grasslands, but in the agricultural matrix they tended to occur in riparian woodlands. At least one target species occurred in every survey quadrat, and giant armadillos, jaguars, and maned wolves were more likely to be present in quadrats located inside than outside the park. The effort required for detection of scats was highest for the two felids. We were able to detect the presence for each of five wide-ranging species inside and outside the park and to assign occurrence probabilities to specific survey sites. Thus, scat dogs provide an effective survey tool for rare species even when accurate detection likelihoods are required. We believe the way we used scat-detection dogs to determine the presence of species can be applied to the detection of other mammalian species in other ecosystems.

An Ancient Baboon Genome Demonstrates Long-Term Population Continuity in Southern Africa.

Baboons are one of the most abundant large nonhuman primates and are widely studied in biomedical, behavioral, and anthropological research. Despite this, our knowledge of their evolutionary and demographic history remains incomplete. Here, we report a 0.9-fold coverage genome sequence from a 5800-year-old baboon from the site of Ha Makotoko in Lesotho. The ancient baboon is closely related to present-day Papio ursinus individuals from southern Africa-indicating a high degree of continuity in the southern African baboon population. This level of population continuity is rare in recent human populations but may provide a good model for the evolution of Homo and other large primates over similar timespans in structured populations throughout Africa.

Distinct genetic architectures and environmental factors associate with host response to the γ2-herpesvirus infections.

Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr Virus (EBV) establish life-long infections and are associated with malignancies. Striking geographic variation in incidence and the fact that virus alone is insufficient to cause disease, suggests other co-factors are involved. Here we present epidemiological analysis and genome-wide association study (GWAS) in 4365 individuals from an African population cohort, to assess the influence of host genetic and non-genetic factors on virus antibody responses. EBV/KSHV co-infection (OR = 5.71(1.58-7.12)), HIV positivity (OR = 2.22(1.32-3.73)) and living in a more rural area (OR = 1.38(1.01-1.89)) are strongly associated with immunogenicity. GWAS reveals associations with KSHV antibody response in the HLA-B/C region (p = 6.64 × 10-09). For EBV, associations are identified for VCA (rs71542439, p = 1.15 × 10-12). Human leucocyte antigen (HLA) and trans-ancestry fine-mapping substantiate that distinct variants in HLA-DQA1 (p = 5.24 × 10-44) are driving associations for EBNA-1 in Africa. This study highlights complex interactions between KSHV and EBV, in addition to distinct genetic architectures resulting in important differences in pathogenesis and transmission.

The ferroportin Q248H mutation protects from anemia, but not malaria or bacteremia.

Iron acquisition is critical for life. Ferroportin (FPN) exports iron from mature erythrocytes, and deletion of the Fpn gene results in hemolytic anemia and increased fatality in malaria-infected mice. The FPN Q248H mutation (glutamine to histidine at position 248) renders FPN partially resistant to hepcidin-induced degradation and was associated with protection from malaria in human studies of limited size. Using data from cohorts including over 18,000 African children, we show that the Q248H mutation is associated with modest protection against anemia, hemolysis, and iron deficiency, but we found little evidence of protection against severe malaria or bacteremia. We additionally observed no excess Plasmodium growth in Q248H erythrocytes ex vivo, nor evidence of selection driven by malaria exposure, suggesting that the Q248H mutation does not protect from malaria and is unlikely to deprive malaria parasites of iron essential for their growth.

Open-source electronic data capture system offered increased accuracy and cost-effectiveness compared with paper methods in Africa.

OBJECTIVES: Existing electronic data capture options are often financially unfeasible in resource-poor settings or difficult to support technically in the field. To help facilitate large-scale multicenter studies in sub-Saharan Africa, the African Partnership for Chronic Disease Research (APCDR) has developed an open-source electronic questionnaire (EQ). STUDY DESIGN AND SETTING: To assess its relative validity, we compared the EQ against traditional pen-and-paper methods using 200 randomized interviews conducted in an ongoing type 2 diabetes case-control study in South Africa. RESULTS: During its 3-month validation, the EQ had a lower frequency of errors (EQ, 0.17 errors per 100 questions; paper, 0.73 errors per 100 questions; P-value ≤0.001), and a lower monetary cost per correctly entered question, compared with the pen-and-paper method. We found no marked difference in the average duration of the interview between methods (EQ, 5.4 minutes; paper, 5.6 minutes). CONCLUSION: This validation study suggests that the EQ may offer increased accuracy, similar interview duration, and increased cost-effectiveness compared with paper-based data collection methods. The APCDR EQ software is freely available (https://github.com/apcdr/questionnaire).

The world's most isolated and distinct whale population? Humpback whales of the Arabian Sea.

A clear understanding of population structure is essential for assessing conservation status and implementing management strategies. A small, non-migratory population of humpback whales in the Arabian Sea is classified as "Endangered" on the IUCN Red List of Threatened Species, an assessment constrained by a lack of data, including limited understanding of its relationship to other populations. We analysed 11 microsatellite markers and mitochondrial DNA sequences extracted from 67 Arabian Sea humpback whale tissue samples and compared them to equivalent datasets from the Southern Hemisphere and North Pacific. Results show that the Arabian Sea population is highly distinct; estimates of gene flow and divergence times suggest a Southern Indian Ocean origin but indicate that it has been isolated for approximately 70,000 years, remarkable for a species that is typically highly migratory. Genetic diversity values are significantly lower than those obtained for Southern Hemisphere populations and signatures of ancient and recent genetic bottlenecks were identified. Our findings suggest this is the world's most isolated humpback whale population, which, when combined with low population abundance estimates and anthropogenic threats, raises concern for its survival. We recommend an amendment of the status of the population to "Critically Endangered" on the IUCN Red List.

The association between circulating lipoprotein(a) and type 2 diabetes: is it causal?

Epidemiological evidence supports a direct and causal association between lipoprotein(a) [Lp(a)] levels and coronary risk, but the nature of the association between Lp(a) levels and risk of type 2 diabetes (T2D) is unclear. In this study, we assessed the association of Lp(a) levels with risk of incident T2D and tested whether Lp(a) levels are causally linked to T2D. We analyzed data on 18,490 participants from the European Prospective Investigation of Cancer (EPIC)-Norfolk cohort that included adults aged 40-79 years at baseline 1993-1997. During an average 10 years of follow-up, 593 participants developed incident T2D. Cox regression models were used to estimate the association between Lp(a) levels and T2D. In Mendelian randomization analyses, based on EPIC-Norfolk combined with DIAbetes Genetics Replication And Meta-analysis data involving a total of 10,088 diabetes case participants and 68,346 control participants, we used a genetic variant (rs10455872) as an instrument to test whether the association between Lp(a) levels and T2D is causal. In adjusted analyses, there was an inverse association between Lp(a) levels and T2D: hazard ratio was 0.63 (95% CI 0.49-0.81; P trend = 0.003) comparing the top versus bottom quintile of Lp(a). In EPIC-Norfolk, a 1-SD increase in logLp(a) was associated with a lower risk of T2D (odds ratio [OR] 0.88 [95% CI: 0.80-0.95]). However, in Mendelian randomization analyses, a 1-SD increase in logLp(a) due to rs10455872, which explained 26.8% of the variability in Lp(a) levels, was not associated with risk of T2D (OR 1.03 [0.96-1.10]; P = 0.41). These prospective findings demonstrate a strong inverse association of Lp(a) levels with risk of T2D. However, a genetic variant that elevated Lp(a) levels was not associated with risk of T2D, suggesting that elevated Lp(a) levels are not causally associated with a lower risk of T2D.