Imputation-based meta-analysis of severe malaria in three African populations

Non-O ABO blood group genotypes differ in their associations with Plasmodium falciparum rosetting and severe malaria

Blood group O is associated with protection against severe malaria and reduced size and stability of P. falciparum-host red blood cell (RBC) rosettes compared to non-O blood groups. Whether the non-O blood groups encoded by the specific ABO genotypes AO, BO, AA, BB and AB differ in their associations with severe malaria and rosetting is unknown. The A and B antigens are host RBC receptors for rosetting, hence we hypothesized that the higher levels of A and/or B antigen on RBCs from AA, BB and AB genotypes compared to AO/BO genotypes could lead to larger rosettes, increased microvascular obstruction and higher risk of malaria pathology. We used a case-control study of Kenyan children and in vitro adhesion assays to test the hypothesis that "double dose" non-O genotypes (AA, BB, AB) are associated with increased risk of severe malaria and larger rosettes than "single dose" heterozygotes (AO, BO). In the case-control study, compared to OO, the double dose genotypes consistently had higher odds ratios (OR) for severe malaria than single dose genotypes, with AB (OR 1.93) and AO (OR 1.27) showing most marked difference (p = 0.02, Wald test). In vitro experiments with blood group A-preferring P. falciparum parasites showed that significantly larger rosettes were formed with AA and AB host RBCs compared to OO, whereas AO and BO genotypes rosettes were indistinguishable from OO. Overall, the data show that ABO genotype influences P. falciparum rosetting and support the hypothesis that double dose non-O genotypes confer a greater risk of severe malaria than AO/BO heterozygosity.

Genome-wide association testing in malaria studies in the presence of overdominance

BACKGROUND

In human genetics, heterozygote advantage (heterosis) has been detected in studies that focused on specific genes but not in genome-wide association studies (GWAS). For example, heterosis is believed to confer resistance to certain strains of malaria in patients heterozygous for the sickle-cell gene, haemoglobin S (HbS). Yet the power of allelic tests can be substantially diminished by heterosis. Since GWAS (and haplotype-associations) also utilize allelic tests, it is unclear to what degree GWAS could underachieve because heterosis is ignored.

METHODS

In this study, a two-step approach to genetic association testing in malaria studies in a GWAS setting that may enhance the power of the tests was proposed, by identifying the underlying genetic model first before applying the association tests. Generalized linear models for dominant, recessive, additive, and heterotic effects were fitted and model selection was performed. This was achieved via tests of significance using the MAX and allelic tests, noting the minimum p-values across all the models and the proportion of tests that a given genetic model was deemed the best. An example dataset, based on 17 SNPs, from a robust genetic association study and simulated genotype datasets, were used to illustrate the method. Case-control genotype data on malaria from Kenya and Gambia were used for validation.

RESULTS AND CONCLUSION

Results showed that the allelic test returned some false negatives under the heterosis model, suggesting reduced power in testing genetic association. Disparities were observed for some chromosomes in the Kenyan and Gambian datasets, including the sex chromosomes. Thus, GWAS and haplotype associations should be treated with caution, unless the underlying genetic model had been determined.

ATP2B4 regulatory genetic variants are associated with mild malaria

BACKGROUND

Genome-wide association studies have identified ATP2B4 as a severe malaria resistance gene. Recently, 8 potential causal regulatory variants have been shown to be associated with severe malaria.

METHODS

Genotyping of rs10900585, rs11240734, rs1541252, rs1541253, rs1541254, rs1541255, rs10751450, rs10751451 and rs10751452 was performed in 154 unrelated individuals (79 controls and 75 mild malaria patients). rs10751450, rs10751451 and rs10751452 were genotyped by Taqman assays, whereas the fragment of the ATP2B4 gene containing the remaining SNPs was sequenced. Logistic regression analysis was used to assess the association between the SNPs and mild malaria.

RESULTS

The results showed that mild malaria was associated with rs10900585, rs11240734, rs1541252, rs1541253, rs1541254, rs1541255, rs10751450, rs10751451 and rs10751452. The homozygous genotypes for the major alleles were associated with an increased risk of mild malaria. Furthermore, the haplotype containing the major alleles and that containing the minor alleles were the most frequent haplotypes. Individuals with the major haplotypes had a significantly higher risk of mild malaria compared to the carriers of the minor allele haplotype.

CONCLUSIONS

ATP2B4 polymorphisms that have been associated with severe malaria are also associated with mild malaria.

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.

Risk score prediction model based on single nucleotide polymorphism for predicting malaria: a machine learning approach

Background

The malaria risk prediction is currently limited to using advanced statistical methods, such as time series and cluster analysis on epidemiological data. Nevertheless, machine learning models have been explored to study the complexity of malaria through blood smear images and environmental data. However, to the best of our knowledge, no study analyses the contribution of Single Nucleotide Polymorphisms (SNPs) to malaria using a machine learning model. More specifically, this study aims to quantify an individual's susceptibility to the development of malaria by using risk scores obtained from the cumulative effects of SNPs, known as weighted genetic risk scores (wGRS).

Results

We proposed an SNP-based feature extraction algorithm that incorporates the susceptibility information of an individual to malaria to generate the feature set. However, it can become computationally expensive for a machine learning model to learn from many SNPs. Therefore, we reduced the feature set by employing the Logistic Regression and Recursive Feature Elimination (LR-RFE) method to select SNPs that improve the efficacy of our model. Next, we calculated the wGRS of the selected feature set, which is used as the model's target variables. Moreover, to compare the performance of the wGRS-only model, we calculated and evaluated the combination of wGRS with genotype frequency (wGRS + GF). Finally, Light Gradient Boosting Machine (LightGBM), eXtreme Gradient Boosting (XGBoost), and Ridge regression algorithms are utilized to establish the machine learning models for malaria risk prediction.

Conclusions

Our proposed approach identified SNP rs334 as the most contributing feature with an importance score of 6.224 compared to the baseline, with an importance score of 1.1314. This is an important result as prior studies have proven that rs334 is a major genetic risk factor for malaria. The analysis and comparison of the three machine learning models demonstrated that LightGBM achieves the highest model performance with a Mean Absolute Error (MAE) score of 0.0373. Furthermore, based on wGRS + GF, all models performed significantly better than wGRS alone, in which LightGBM obtained the best performance (0.0033 MAE score).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04870-0.

Delta-Catenin as a Modulator of Rho GTPases in Neurons

Small Rho GTPases are molecular switches that are involved in multiple processes including regulation of the actin cytoskeleton. These GTPases are activated (turned on) and inactivated (turned off) through various upstream effector molecules to carry out many cellular functions. One such upstream modulator of small Rho GTPase activity is delta-catenin, which is a protein in the p120-catenin subfamily that is enriched in the central nervous system. Delta-catenin affects small GTPase activity to assist in the developmental formation of dendrites and dendritic spines and to maintain them once they mature. As the dendritic arbor and spine density are crucial for synapse formation and plasticity, delta-catenin's ability to modulate small Rho GTPases is necessary for proper learning and memory. Accordingly, the misregulation of delta-catenin and small Rho GTPases has been implicated in several neurological and non-neurological pathologies. While links between delta-catenin and small Rho GTPases have yet to be studied in many contexts, known associations include some cancers, Alzheimer's disease (AD), Cri-du-chat syndrome, and autism spectrum disorder (ASD). Drawing from established studies and recent discoveries, this review explores how delta-catenin modulates small Rho GTPase activity. Future studies will likely elucidate how PDZ proteins that bind delta-catenin further influence small Rho GTPases, how delta-catenin may affect small GTPase activity at adherens junctions when bound to N-cadherin, mechanisms behind delta-catenin's ability to modulate Rac1 and Cdc42, and delta-catenin's ability to modulate small Rho GTPases in the context of diseases, such as cancer and AD.

Identification of ATP2B4 Regulatory Element Containing Functional Genetic Variants Associated with Severe Malaria

Genome-wide association studies for severe malaria (SM) have identified 30 genetic variants mostly located in non-coding regions. Here, we aimed to identify potential causal genetic variants located in these loci and demonstrate their functional activity. We systematically investigated the regulatory effect of the SNPs in linkage disequilibrium (LD) with the malaria-associated genetic variants. Annotating and prioritizing genetic variants led to the identification of a regulatory region containing five ATP2B4 SNPs in LD with rs10900585. We found significant associations between SM and rs10900585 and our candidate SNPs (rs11240734, rs1541252, rs1541253, rs1541254, and rs1541255) in a Senegalese population. Then, we demonstrated that both individual SNPs and the combination of SNPs had regulatory effects. Moreover, CRISPR/Cas9-mediated deletion of this region decreased ATP2B4 transcript and protein levels and increased Ca²⁺ intracellular concentration in the K562 cell line. Our data demonstrate that severe malaria-associated genetic variants alter the expression of ATP2B4 encoding a plasma membrane calcium-transporting ATPase 4 (PMCA4) expressed on red blood cells. Altering the activity of this regulatory element affects the risk of SM, likely through calcium concentration effect on parasitaemia.

Leveraging Mann-Whitney U test on large-scale genetic variation data for analysing malaria genetic markers

Background

The malaria risk analysis of multiple populations is crucial and of great importance whilst compressing limitations. However, the exponential growth in diversity and accumulation of genetic variation data obtained from malaria-infected patients through Genome-Wide Association Studies opens up unprecedented opportunities to explore the significant differences between genetic markers (risk factors), particularly in the resistance or susceptibility of populations to malaria risk. Thus, this study proposes using statistical tests to analyse large-scale genetic variation data, comprising 20,854 samples from 11 populations within three continents: Africa, Oceania, and Asia.

Methods

Even though statistical tests have been utilized to conduct case–control studies since the 1950s to link risk factors to a particular disease, several challenges faced, including the choice of data (ordinal vs. non-ordinal) and test (parametric vs. non-parametric). This study overcomes these challenges by adopting the Mann–Whitney U test to analyse large-scale genetic variation data; to explore the statistical significance of markers between populations; and to further identify the highly differentiated markers.

Results

The findings of this study revealed a significant difference in the genetic markers between populations (p < 0.01) in all the case groups and most control groups. However, for the highly differentiated genetic markers, a significant difference (p < 0.01) was present for most genetic markers with varying p-values between the populations in the case and control groups. Moreover, several genetic markers were observed to have very significant differences (p < 0.001) across all populations, while others exist between certain specific populations. Also, several genetic markers have no significant differences between populations.

Conclusions

These findings further support that the genetic markers contribute differently between populations towards malaria resistance or susceptibility, thus showing differences in the likelihood of malaria infection. In addition, this study demonstrated the robustness of the Mann–Whitney U test in analysing genetic markers in large-scale genetic variation data, thereby indicating an alternative method to explore genetic markers in other complex diseases. The findings hold great promise for genetic markers analysis, and the pipeline emphasized in this study can fully be reproduced to analyse new data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04104-x.

Bayesian model comparison for rare-variant association studies

Whole-genome sequencing studies applied to large populations or biobanks with extensive phenotyping raise new analytic challenges. The need to consider many variants at a locus or group of genes simultaneously and the potential to study many correlated phenotypes with shared genetic architecture provide opportunities for discovery not addressed by the traditional one variant, one phenotype association study. Here, we introduce a Bayesian model comparison approach called MRP (multiple rare variants and phenotypes) for rare-variant association studies that considers correlation, scale, and direction of genetic effects across a group of genetic variants, phenotypes, and studies, requiring only summary statistic data. We apply our method to exome sequencing data (n = 184,698) across 2,019 traits from the UK Biobank, aggregating signals in genes. MRP demonstrates an ability to recover signals such as associations between PCSK9 and LDL cholesterol levels. We additionally find MRP effective in conducting meta-analyses in exome data. Non-biomarker findings include associations between MC1R and red hair color and skin color, IL17RA and monocyte count, and IQGAP2 and mean platelet volume. Finally, we apply MRP in a multi-phenotype setting; after clustering the 35 biomarker phenotypes based on genetic correlation estimates, we find that joint analysis of these phenotypes results in substantial power gains for gene-trait associations, such as in TNFRSF13B in one of the clusters containing diabetes- and lipid-related traits. Overall, we show that the MRP model comparison approach improves upon useful features from widely used meta-analysis approaches for rare-variant association analyses and prioritizes protective modifiers of disease risk.

Haplotype heterogeneity and low linkage disequilibrium reduce reliable prediction of genotypes for the ‑α 3.7I form of α-thalassaemia using genome-wide microarray data

Background: The -α ^3.7I-thalassaemia deletion is very common throughout Africa because it protects against malaria. When undertaking studies to investigate human genetic adaptations to malaria or other diseases, it is important to account for any confounding effects of α-thalassaemia to rule out spurious associations.

Methods: In this study, we have used direct α-thalassaemia genotyping to understand why GWAS data from a large malaria association study in Kilifi Kenya did not identify the α-thalassaemia signal. We then explored the potential use of a number of new approaches to using GWAS data for imputing α-thalassaemia as an alternative to direct genotyping by PCR.

Results: We found very low linkage-disequilibrium of the directly typed data with the GWAS SNP markers around α-thalassaemia and across the haemoglobin-alpha ( HBA) gene region, which along with a complex haplotype structure, could explain the lack of an association signal from the GWAS SNP data. Some indirect typing methods gave results that were in broad agreement with those derived from direct genotyping and could identify an association signal, but none were sufficiently accurate to allow correct interpretation compared with direct typing, leading to confusing or erroneous results.

Conclusions: We conclude that going forwards, direct typing methods such as PCR will still be required to account for α-thalassaemia in GWAS studies.

Improving statistical power in severe malaria genetic association studies by augmenting phenotypic precision

Severe falciparum malaria has substantially affected human evolution. Genetic association studies of patients with clinically defined severe malaria and matched population controls have helped characterise human genetic susceptibility to severe malaria, but phenotypic imprecision compromises discovered associations. In areas of high malaria transmission, the diagnosis of severe malaria in young children and, in particular, the distinction from bacterial sepsis are imprecise. We developed a probabilistic diagnostic model of severe malaria using platelet and white count data. Under this model, we re-analysed clinical and genetic data from 2220 Kenyan children with clinically defined severe malaria and 3940 population controls, adjusting for phenotype mis-labelling. Our model, validated by the distribution of sickle trait, estimated that approximately one-third of cases did not have severe malaria. We propose a data-tilting approach for case-control studies with phenotype mis-labelling and show that this reduces false discovery rates and improves statistical power in genome-wide association studies.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

BACKGROUND

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility.

RESULTS

Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb)) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs associated with both severe COVID-19 and other human traits demonstrated colocalization of the GWAS signal at the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN). This finding points to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity.

CONCLUSIONS

Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches. The iCPAGdb web portal is accessible at http://cpag.oit.duke.edu and the software code at https://github.com/tbalmat/iCPAGdb .

Fine scale human genetic structure in three regions of Cameroon reveals episodic diversifying selection

Inferences from genetic association studies rely largely on the definition and description of the underlying populations that highlight their genetic similarities and differences. The clustering of human populations into subgroups (population structure) can significantly confound disease associations. This study investigated the fine-scale genetic structure within Cameroon that may underlie disparities observed with Cameroonian ethnicities in malaria genome-wide association studies in sub-Saharan Africa. Genotype data of 1073 individuals from three regions and three ethnic groups in Cameroon were analyzed using measures of genetic proximity to ascertain fine-scale genetic structure. Model-based clustering revealed distinct ancestral proportions among the Bantu, Semi-Bantu and Foulbe ethnic groups, while haplotype-based coancestry estimation revealed possible longstanding and ongoing sympatric differentiation among individuals of the Foulbe ethnic group, and their Bantu and Semi-Bantu counterparts. A genome scan found strong selection signatures in the HLA gene region, confirming longstanding knowledge of natural selection on this genomic region in African populations following immense disease pressure. Signatures of selection were also observed in the HBB gene cluster, a genomic region known to be under strong balancing selection in sub-Saharan Africa due to its co-evolution with malaria. This study further supports the role of evolution in shaping genomes of Cameroonian populations and reveals fine-scale hierarchical structure among and within Cameroonian ethnicities that may impact genetic association studies in the country.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility. Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb); http://cpag.oit.duke.edu) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs with severe COVID-19 demonstrated colocalization of the GWAS signal of the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN), pointing to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches.

Haplotype heterogeneity and low linkage disequilibrium reduce reliable prediction of genotypes for the ‑α3.7I form of α-thalassaemia using genome-wide microarray data

Background: The -α3.7I-thalassaemia deletion is very common throughout Africa because it protects against malaria. When undertaking studies to investigate human genetic adaptations to malaria or other diseases, it is important to account for any confounding effects of α-thalassaemia to rule out spurious associations. Methods: In this study we have used direct α-thalassaemia genotyping to understand why GWAS data from a large malaria association study in Kilifi Kenya did not identify the α-thalassaemia signal. We then explored the potential use of a number of new approaches to using GWAS data for imputing α-thalassaemia as an alternative to direct genotyping by PCR. Results: We found very low linkage-disequilibrium of the directly typed data with the GWAS SNP markers around α-thalassaemia and across the haemoglobin-alpha (HBA) gene region, which along with a complex haplotype structure, could explain the lack of an association signal from the GWAS SNP data. Some indirect typing methods gave results that were in broad agreement with those derived from direct genotyping and could identify an association signal, but none were sufficiently accurate to allow correct interpretation compared with direct typing, leading to confusing or erroneous results. Conclusions: We conclude that going forwards, direct typing methods such as PCR will still be required to account for α-thalassaemia in GWAS studies.

Human genetics and malaria resistance

Malaria has been the pre-eminent cause of early mortality in many parts of the world throughout much of the last five thousand years and, as a result, it is the strongest force for selective pressure on the human genome yet described. Around one third of the variability in the risk of severe and complicated malaria is now explained by additive host genetic effects. Many individual variants have been identified that are associated with malaria protection, but the most important all relate to the structure or function of red blood cells. They include the classical polymorphisms that cause sickle cell trait, α-thalassaemia, G6PD deficiency, and the major red cell blood group variants. More recently however, with improving technology and experimental design, others have been identified that include the Dantu blood group variant, polymorphisms in the red cell membrane protein ATP2B4, and several variants related to the immune response. Characterising how these genes confer their effects could eventually inform novel therapeutic approaches to combat malaria. Nevertheless, all together, only a small proportion of the heritable component of malaria resistance can be explained by the variants described so far, underscoring its complex genetic architecture and the need for continued research.

Rare protein-altering variants in ANGPTL7 lower intraocular pressure and protect against glaucoma

Protein-altering variants that are protective against human disease provide in vivo validation of therapeutic targets. Here we use genotyping data from UK Biobank (n = 337,151 unrelated White British individuals) and FinnGen (n = 176,899) to conduct a search for protein-altering variants conferring lower intraocular pressure (IOP) and protection against glaucoma. Through rare protein-altering variant association analysis, we find a missense variant in ANGPTL7 in UK Biobank (rs28991009, p.Gln175His, MAF = 0.8%, genotyped in 82,253 individuals with measured IOP and an independent set of 4,238 glaucoma patients and 250,660 controls) that significantly lowers IOP (β = -0.53 and -0.67 mmHg for heterozygotes, -3.40 and -2.37 mmHg for homozygotes, P = 5.96 x 10-9 and 1.07 x 10-13 for corneal compensated and Goldman-correlated IOP, respectively) and is associated with 34% reduced risk of glaucoma (P = 0.0062). In FinnGen, we identify an ANGPTL7 missense variant at a greater than 50-fold increased frequency in Finland compared with other populations (rs147660927, p.Arg220Cys, MAF Finland = 4.3%), which was genotyped in 6,537 glaucoma patients and 170,362 controls and is associated with a 29% lower glaucoma risk (P = 1.9 x 10-12 for all glaucoma types and also protection against its subtypes including exfoliation, primary open-angle, and primary angle-closure). We further find three rarer variants in UK Biobank, including a protein-truncating variant, which confer a strong composite lowering of IOP (P = 0.0012 and 0.24 for Goldman-correlated and corneal compensated IOP, respectively), suggesting the protective mechanism likely resides in the loss of interaction or function. Our results support inhibition or down-regulation of ANGPTL7 as a therapeutic strategy for glaucoma.

Genome-wide heritability analysis of severe malaria resistance reveals evidence of polygenic inheritance

BACKGROUND

Estimating single nucleotide polymorphism (SNP)-heritability (h2g) of severe malaria resistance and its distribution across the genome might shed new light in to the underlying biology.

METHOD

We investigated h2g of severe malaria resistance from a genome-wide association study (GWAS) dataset (sample size = 11 657). We estimated the h2g and partitioned in to chromosomes, allele frequencies and annotations using the genetic relationship-matrix restricted maximum likelihood approach. We further examined non-cell type-specific and cell type-specific enrichments from GWAS-summary statistics.

RESULTS

The h2g of severe malaria resistance was estimated at 0.21 (se = 0.05, P = 2.7 × 10-5), 0.20 (se = 0.05, P = 7.5 × 10-5) and 0.17 (se = 0.05, P = 7.2 × 10-4) in Gambian, Kenyan and Malawi populations, respectively. A comparable range of h2g [0.21 (se = 0.02, P < 1 × 10-5)] was estimated from GWAS-summary statistics meta-analysed across the three populations. Partitioning analysis from raw genotype data showed significant enrichment of h2g in genic SNPs while summary statistics analysis suggests evidences of enrichment in multiple categories. Supporting the polygenic inheritance, the h2g of severe malaria resistance is distributed across the chromosomes and allelic frequency spectrum. However, the h2g is disproportionately concentrated on three chromosomes (chr 5, 11 and 20), suggesting cost-effectiveness of targeting these chromosomes in future malaria genomic sequencing studies.

CONCLUSION

We report for the first time that the heritability of malaria resistance is largely ascribed by common SNPs and the causal variants are overrepresented in protein coding regions of the genome. Further studies with larger sample sizes are needed to better understand the underpinning genetics of severe malaria resistance.

Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania

The human genetic factors that affect resistance to infectious disease are poorly understood. Here we report a genome-wide association study in 17,000 severe malaria cases and population controls from 11 countries, informed by sequencing of family trios and by direct typing of candidate loci in an additional 15,000 samples. We identify five replicable associations with genome-wide levels of evidence including a newly implicated variant on chromosome 6. Jointly, these variants account for around one-tenth of the heritability of severe malaria, which we estimate as ~23% using genome-wide genotypes. We interrogate available functional data and discover an erythroid-specific transcription start site underlying the known association in ATP2B4, but are unable to identify a likely causal mechanism at the chromosome 6 locus.  Previously reported HLA associations do not replicate in these samples. This large dataset will provide a foundation for further research on the genetic determinants of malaria resistance in diverse populations.

Four genome-wide associated loci are currently known for malaria susceptibility. Here, the authors expand on earlier work by combining data from 11 malaria-endemic countries and additional population sequencing informing an African-enriched imputation reference panel, with findings including a previously unreported association on chromosome 6.

Analysis of five deep-sequenced trio-genomes of the Peninsular Malaysia Orang Asli and North Borneo populations

Background

Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated.

Results

We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10^− 8 – 1.33 × 10^− 8, 1.0 × 10^− 9 – 2.9 × 10^− 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples.

Conclusion

Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.

Gene expression profiling in blood from cerebral malaria patients and mild malaria patients living in Senegal

Background

Plasmodium falciparum malaria remains a major health problem in Africa. The mechanisms of pathogenesis are not fully understood. Transcriptomic studies may provide new insights into molecular pathways involved in the severe form of the disease.

Methods

Blood transcriptional levels were assessed in patients with cerebral malaria, non-cerebral malaria, or mild malaria by using microarray technology to look for gene expression profiles associated with clinical status. Multi-way ANOVA was used to extract differentially expressed genes. Network and pathways analyses were used to detect enrichment for biological pathways.

Results

We identified a set of 443 genes that were differentially expressed in the three patient groups after applying a false discovery rate of 10%. Since the cerebral patients displayed a particular transcriptional pattern, we focused our analysis on the differences between cerebral malaria patients and mild malaria patients. We further found 842 differentially expressed genes after applying a false discovery rate of 10%. Unsupervised hierarchical clustering of cerebral malaria-informative genes led to clustering of the cerebral malaria patients. The support vector machine method allowed us to correctly classify five out of six cerebral malaria patients and six of six mild malaria patients. Furthermore, the products of the differentially expressed genes were mapped onto a human protein-protein network. This led to the identification of the proteins with the highest number of interactions, including GSK3B, RELA, and APP. The enrichment analysis of the gene functional annotation indicates that genes involved in immune signalling pathways play a role in the occurrence of cerebral malaria. These include BCR-, TCR-, TLR-, cytokine-, FcεRI-, and FCGR- signalling pathways and natural killer cell cytotoxicity pathways, which are involved in the activation of immune cells. In addition, our results revealed an enrichment of genes involved in Alzheimer’s disease.

Conclusions

In the present study, we examine a set of genes whose expression differed in cerebral malaria patients and mild malaria patients. Moreover, our results provide new insights into the potential effect of the dysregulation of gene expression in immune pathways. Host genetic variation may partly explain such alteration of gene expression. Further studies are required to investigate this in African populations.

First genome-wide association study of non-severe malaria in two birth cohorts in Benin

Recent research efforts to identify genes involved in malaria susceptibility using genome-wide approaches have focused on severe malaria. Here, we present the first GWAS on non-severe malaria designed to identify genetic variants involved in innate immunity or innate resistance mechanisms. Our study was performed on two cohorts of infants from southern Benin (525 and 250 individuals used as discovery and replication cohorts, respectively) closely followed from birth to 18-24 months of age, with an assessment of a space- and time-dependent environmental risk of exposure. Both the recurrence of mild malaria attacks and the recurrence of malaria infections as a whole (symptomatic and asymptomatic) were considered. Post-GWAS functional analyses were performed using positional, eQTL, and chromatin interaction mapping to identify the genes underlying association signals. Our study highlights a role of PTPRT, a tyrosine phosphatase receptor involved in STAT3 pathway, in the protection against both mild malaria attacks and malaria infections (p = 9.70 × 10^-8 and p = 1.78 × 10^-7, respectively, in the discovery cohort). Strong statistical support was also found for a role of MYLK4 (meta-analysis, p = 5.29 × 10^-8 with malaria attacks), and for several other genes, whose biological functions are relevant in malaria infection. Results shows that GWAS on non-severe malaria can successfully identify new candidate genes and inform physiological mechanisms underlying natural protection against malaria.

Interferon lambda 4 impacts the genetic diversity of hepatitis C virus

Hepatitis C virus (HCV) is a highly variable pathogen that frequently establishes chronic infection. This genetic variability is affected by the adaptive immune response but the contribution of other host factors is unclear. Here, we examined the role played by interferon lambda-4 (IFN-λ4) on HCV diversity; IFN-λ4 plays a crucial role in spontaneous clearance or establishment of chronicity following acute infection. We performed viral genome-wide association studies using human and viral data from 485 patients of white ancestry infected with HCV genotype 3a. We demonstrate that combinations of host genetic variants, which determine IFN-λ4 protein production and activity, influence amino acid variation across the viral polyprotein - not restricted to specific viral proteins or HLA restricted epitopes - and modulate viral load. We also observed an association with viral di-nucleotide proportions. These results support a direct role for IFN-λ4 in exerting selective pressure across the viral genome, possibly by a novel mechanism.

Genetics, adaptation to environmental changes and archaic admixture in the pathogenesis of diabetes mellitus in Indigenous Australians

Indigenous Australians are particularly affected by type 2 diabetes mellitus (T2D) due to both their genetic susceptibility and a range of environmental and lifestyle risk factors. Recent genetic studies link predisposition to some diseases, including T2D, to alleles acquired from archaic hominins, such as Neanderthals and Denisovans, which persist in the genomes of modern humans today. Indo-Pacific human populations, including Indigenous Australians, remain extremely underrepresented in genomic research with a paucity of data examining the impact of Denisovan or Neanderthal lineages on human phenotypes in Oceania. The few genetic studies undertaken emphasize the uniqueness and antiquity of Indigenous Australian genomes, with possibly the largest proportion of Denisovan ancestry of any population in the world. In this review, we focus on the potential contributions of ancient genes/pathways to modern human phenotypes, while also highlighting the evolutionary roles of genetic adaptation to dietary and environmental changes associated with an adopted Western lifestyle. We discuss the role of genetic and epigenetic factors in the pathogenesis of T2D in understudied Indigenous Australians, including the potential impact of archaic gene lineages on this disease. Finally, we propose that greater understanding of the underlying genetic predisposition may contribute to the clinical efficacy of diabetes management in Indigenous Australians. We suggest that improved identification of T2D risk variants in Oceania is needed. Such studies promise to clarify how genetic and phenotypic differences vary between populations and, crucially, provide novel targets for personalised medical therapies in currently marginalized groups.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

Genetics of Malaria Inflammatory Responses: A Pathogenesis Perspective

Despite significant progress in combating malaria in recent years the burden of severe disease and death due to Plasmodium infections remains a global public health concern. Only a fraction of infected people develops severe clinical syndromes motivating a longstanding search for genetic determinants of malaria severity. Strong genetic effects have been repeatedly ascribed to mutations and allelic variants of proteins expressed in red blood cells but the role of inflammatory response genes in disease pathogenesis has been difficult to discern. We revisited genetic evidence provided by inflammatory response genes that have been repeatedly associated to malaria, namely TNF, NOS2, IFNAR1, HMOX1, TLRs, CD36, and CD40LG. This highlighted specific genetic variants having opposing roles in the development of distinct malaria clinical outcomes and unveiled diverse levels of genetic heterogeneity that shaped the complex association landscape of inflammatory response genes with malaria. However, scrutinizing genetic effects of individual variants corroborates a pathogenesis model where pro-inflammatory genetic variants acting in early infection stages contribute to resolve infection but at later stages confer increased vulnerability to severe organ dysfunction driven by tissue inflammation. Human genetics studies are an invaluable tool to find genes and molecular pathways involved in the inflammatory response to malaria but their precise roles in disease pathogenesis are still unexploited. Genome editing in malaria experimental models and novel genotyping-by-sequencing techniques are promising approaches to delineate the relevance of inflammatory response gene variants in the natural history of infection thereby will offer new rational angles on adjuvant therapeutics for prevention and clinical management of severe malaria.

Bayesian meta-analysis across genome-wide association studies of diverse phenotypes

Genome-wide association studies (GWAS) are a powerful tool for understanding the genetic basis of diseases and traits, but most studies have been conducted in isolation, with a focus on either a single or a set of closely related phenotypes. We describe MetABF, a simple Bayesian framework for performing integrative meta-analysis across multiple GWAS using summary statistics. The approach is applicable across a wide range of study designs and can increase the power by 50% compared with standard frequentist tests when only a subset of studies have a true effect. We demonstrate its utility in a meta-analysis of 20 diverse GWAS which were part of the Wellcome Trust Case Control Consortium 2. The novelty of the approach is its ability to explore, and assess the evidence for a range of possible true patterns of association across studies in a computationally efficient framework.

Whole genome sequence association with E-selectin levels reveals loss-of-function variant in African Americans

E-selectin mediates the rolling of circulating leukocytes during inflammatory processes. Previous genome-wide association studies in European and Asian individuals have identified the ABO locus associated with E-selectin levels. Using Trans-Omics for Precision Medicine whole genome sequencing data in 2249 African Americans (AAs) from the Jackson Heart Study, we examined genome-wide associations with soluble E-selectin levels. In addition to replicating known signals at ABO, we identified a novel association of a common loss-of-function, missense variant in Fucosyltransferase 6 (FUT6; rs17855739,p.Glu274Lys, P = 9.02 × 10-24) with higher soluble E-selectin levels. This variant is considerably more common in populations of African ancestry compared to non-African ancestry populations. We replicated the association of FUT6 p.Glu274Lys with higher soluble E-selectin in an independent population of 748 AAs from the Women's Health Initiative and identified an additional pleiotropic association with vitamin B12 levels. Despite the broad role of both selectins and fucosyltransferases in various inflammatory, immune and cancer-related processes, we were unable to identify any additional disease associations of the FUT6 p.Glu274Lys variant in an electronic medical record-based phenome-wide association scan of over 9000 AAs.

Natural selection and local adaptation of blood pressure regulation and their perspectives on precision medicine in hypertension

Prevalence of hypertension (HTN) varies substantially across different populations. HTN is not only common - affecting at least one third of the world's adult population - but is also the most important driver for cardiovascular diseases. Yet up to a third of hypertensive patients are resistant to therapy, contributed by secondary hypertension but more commonly the hitherto inability to precisely predict response to specific antihypertensive agents. Population and individual genomics information could be useful in guiding the selection and predicting the response to treatment - an approach known as precision medicine. However this cannot be achieved without the knowledge of genetic variations that influence blood pressure (BP). A number of evolutionary factors including population demographics and forces of natural selection may be involved. This article explores some ideas on how natural selection influences BP regulation in ethnically and geographically diverse populations that could lead to them being susceptible to HTN. We explore how such evolutionary factors could impact the implementation of precision medicine in HTN. Finally, in order to ensure the success of precision medicine in HTN, we call for more initiatives to understand the genetic architecture within and between diverse populations with ancestry from different parts of the world, and to precisely classify the intermediate phenotypes of HTN.

NCR3 polymorphism, haematological parameters, and severe malaria in Senegalese patients

Background

Host factors, including host genetic variation, have been shown to influence the outcome of Plasmodium falciparum infection. Genome-wide linkage studies have mapped mild malaria resistance genes on chromosome 6p21, whereas NCR3-412 polymorphism (rs2736191) lying within this region was found to be associated with mild malaria.

Methods

Blood samples were taken from 188 Plasmodium falciparum malaria patients (76 mild malaria patients, 85 cerebral malaria patients, and 27 severe non-cerebral malaria patients). NCR3-412 (rs2736191) was analysed by sequencing, and haematological parameters were measured. Finally, their association with clinical phenotypes was assessed.

Results

We evidenced an association of thrombocytopenia with both cerebral malaria and severe non-cerebral malaria, and of an association of high leukocyte count with cerebral malaria. Additionally, we found no association of NCR3-412 with either cerebral malaria, severe non-cerebral malaria, or severe malaria after grouping cerebral malaria and severe non-cerebral malaria patients.

Conclusions

Our results suggest that NCR3 genetic variation has no effect, or only a small effect on the occurrence of severe malaria, although it has been strongly associated with mild malaria. We discuss the biological meaning of these results. Besides, we confirmed the association of thrombocytopenia and high leukocyte count with severe malaria phenotypes.

Variation and Functional Impact of Neanderthal Ancestry in Western Asia

Neanderthals contributed genetic material to modern humans via multiple admixture events. Initial admixture events presumably occurred in Western Asia shortly after humans migrated out of Africa. Despite being a focal point of admixture, earlier studies indicate lower Neanderthal introgression rates in some Western Asian populations as compared with other Eurasian populations. To better understand the genome-wide and phenotypic impact of Neanderthal introgression in the region, we sequenced whole genomes of nine present-day Europeans, Africans, and the Western Asian Druze at high depth, and analyzed available whole genome data from various other populations, including 16 genomes from present-day Turkey. Our results confirmed previous observations that contemporary Western Asian populations, on an average, have lower levels of Neanderthal-introgressed DNA relative to other Eurasian populations. Modern Western Asians also show comparatively high variability in Neanderthal ancestry, which may be attributed to the complex demographic history of the region. We further replicated the previously described depletion of putatively functional sequences among Neanderthal-introgressed haplotypes. Still, we find dozens of common Neanderthal-introgressed haplotypes in the Turkish sample associated with human phenotypes, including anthropometric and metabolic traits, as well as the immune response. One of these haplotypes is unusually long and harbors variants that affect the expression of members of the CCR gene family and are associated with celiac disease. Overall, our results paint a complex first picture of the genomic impact of Neanderthal introgression in the Western Asian populations.

Adaptive Landscape of Protein Variation in Human Exomes

The human genome contains hundreds of thousands of missense mutations. However, only a handful of these variants are known to be adaptive, which implies that adaptation through protein sequence change is an extremely rare phenomenon in human evolution. Alternatively, existing methods may lack the power to pinpoint adaptive variation. We have developed and applied an Evolutionary Probability Approach (EPA) to discover candidate adaptive polymorphisms (CAPs) through the discordance between allelic evolutionary probabilities and their observed frequencies in human populations. EPA reveals thousands of missense CAPs, which suggest that a large number of previously optimal alleles experienced a reversal of fortune in the human lineage. We explored nonadaptive mechanisms to explain CAPs, including the effects of demography, mutation rate variability, and negative and positive selective pressures in modern humans. Many nonadaptive hypotheses were tested, but failed to explain the data, which suggests that a large proportion of CAP alleles have increased in frequency due to beneficial selection. This suggestion is supported by the fact that a vast majority of adaptive missense variants discovered previously in humans are CAPs, and hundreds of CAP alleles are protective in genotype-phenotype association data. Our integrated phylogenomic and population genetic EPA approach predicts the existence of thousands of nonneutral candidate variants in the human proteome. We expect this collection to be enriched in beneficial variation. The EPA approach can be applied to discover candidate adaptive variation in any protein, population, or species for which allele frequency data and reliable multispecies alignments are available.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Two complement receptor one alleles have opposing associations with cerebral malaria and interact with α+thalassaemia

Malaria has been a major driving force in the evolution of the human genome. In sub-Saharan African populations, two neighbouring polymorphisms in the Complement Receptor One (CR1) gene, named Sl2 and McCb, occur at high frequencies, consistent with selection by malaria. Previous studies have been inconclusive. Using a large case-control study of severe malaria in Kenyan children and statistical models adjusted for confounders, we estimate the relationship between Sl2 and McCb and malaria phenotypes, and find they have opposing associations. The Sl2 polymorphism is associated with markedly reduced odds of cerebral malaria and death, while the McCb polymorphism is associated with increased odds of cerebral malaria. We also identify an apparent interaction between Sl2 and α+thalassaemia, with the protective association of Sl2 greatest in children with normal α-globin. The complex relationship between these three mutations may explain previous conflicting findings, highlighting the importance of considering genetic interactions in disease-association studies.

Genetic susceptibility to infectious diseases: Current status and future perspectives from genome-wide approaches

Genome-wide association studies (GWASs) have been widely applied to identify genetic factors that affect complex diseases or traits. Presently, the GWAS Catalog includes > 2800 human studies. Of these, only a minority have investigated the susceptibility to infectious diseases or the response to therapies for the treatment or prevention of infections. Despite their limited application in the field, GWASs have provided valuable insights by pinpointing associations to both innate and adaptive immune response loci, as well as novel unexpected risk factors for infection susceptibility. Herein, we discuss some issues and caveats of GWASs for infectious diseases, we review the most recent findings ensuing from these studies, and we provide a brief summary of selected GWASs for infections in non-human mammals. We conclude that, although the general trend in the field of complex traits is to shift from GWAS to next-generation sequencing, important knowledge on infectious disease-related traits can be still gained by GWASs, especially for those conditions that have never been investigated using this approach. We suggest that future studies will benefit from the leveraging of information from the host's and pathogen's genomes, as well as from the exploration of models that incorporate heterogeneity across populations and phenotypes. Interactions within HLA genes or among HLA variants and polymorphisms located outside the major histocompatibility complex may also play an important role in shaping the susceptibility and response to invading pathogens.

•

Relatively few GWASs for infectious diseases were performed.

•

Phenotype heterogeneity and case/control misclassification can affect GWAS power.

•

Adaptive and innate immunity loci were identified in several infectious disease GWASs.

•

Unexpected loci (e.g., lncRNAs) were also associated with infection susceptibility.

•

GWASs should integrate host and pathogen diversity and use complex association models.

Human genetic variation in VAC14 regulates Salmonella invasion and typhoid fever through modulation of cholesterol

Risk, severity, and outcome of infection depend on the interplay of pathogen virulence and host susceptibility. Systematic identification of genetic susceptibility to infection is being undertaken through genome-wide association studies, but how to expeditiously move from genetic differences to functional mechanisms is unclear. Here, we use genetic association of molecular, cellular, and human disease traits and experimental validation to demonstrate that genetic variation affects expression of VAC14, a phosphoinositide-regulating protein, to influence susceptibility to Salmonella enterica serovar Typhi (S Typhi) infection. Decreased VAC14 expression increased plasma membrane cholesterol, facilitating Salmonella docking and invasion. This increased susceptibility at the cellular level manifests as increased susceptibility to typhoid fever in a Vietnamese population. Furthermore, treating zebrafish with a cholesterol-lowering agent, ezetimibe, reduced susceptibility to S Typhi. Thus, coupling multiple genetic association studies with mechanistic dissection revealed how VAC14 regulates Salmonella invasion and typhoid fever susceptibility and may open doors to new prophylactic/therapeutic approaches.

Diversity and inclusion in genomic research: why the uneven progress?

Conducting genomic research in diverse populations has led to numerous advances in our understanding of human history, biology, and health disparities, in addition to discoveries of vital clinical significance. Conducting genomic research in diverse populations is also important in ensuring that the genomic revolution does not exacerbate health disparities by facilitating discoveries that will disproportionately benefit well-represented populations. Despite the general agreement on the need for genomic research in diverse populations in terms of equity and scientific progress, genomic research remains largely focused on populations of European descent. In this article, we describe the rationale for conducting genomic research in diverse populations by reviewing examples of advances facilitated by their inclusion. We also explore some of the factors that perpetuate the disproportionate attention on well-represented populations. Finally, we discuss ongoing efforts to ameliorate this continuing bias. Collaborative and intensive efforts at all levels of research, from the funding of studies to the publication of their findings, will be necessary to ensure that genomic research does not conserve historical inequalities or curtail the contribution that genomics could make to the health of all humanity.

Lack of Association of CD55 Receptor Genetic Variants and Severe Malaria in Ghanaian Children

In a recent report, the cellular receptor CD55 was identified as a molecule essential for the invasion of human erythrocytes by Plasmodium falciparum, the causal agent of the most severe form of malaria. As this invasion process represents a critical step during infection with the parasite, it was hypothesized that genetic variants in the gene could affect severe malaria (SM) susceptibility. We performed high-resolution variant discovery of rare and common genetic variants in the human CD55 gene. Association testing of these variants in over 1700 SM cases and unaffected control individuals from the malaria-endemic Ashanti Region in Ghana, West Africa, were performed on the basis of single variants, combined rare variant analyses, and reconstructed haplotypes. A total of 26 genetic variants were detected in coding and regulatory regions of CD55. Five variants were previously unknown. None of the single variants, rare variants, or haplotypes showed evidence for association with SM or P. falciparum density. Here, we present the first comprehensive analysis of variation in the CD55 gene in the context of SM and show that genetic variants present in a Ghanaian study group appear not to influence susceptibility to the disease.

Building Ancestral Recombination Graphs for Whole Genomes

We propose a heuristic algorithm, called ARG4WG, to build plausible ancestral recombination graphs (ARGs) from thousands of whole genome samples. By using the longest shared end for recombination inference, ARG4WG constructs ARGs with small numbers of recombination events that perform well in association mapping on genome-wide association studies.

Frequency of TNFA, INFG, and IL10 Gene Polymorphisms and Their Association with Malaria Vivax and Genomic Ancestry

Polymorphisms in cytokine genes can alter the production of these proteins and consequently affect the immune response. The trihybrid heterogeneity of the Brazilian population is characterized as a condition for the use of ancestry informative markers. The objective of this study was to evaluate the frequency of -1031T>C, -308G>A and -238G>A TNFA, +874 A>T IFNG and -819C>T, and -592C>A IL10 gene polymorphisms and their association with malaria vivax and genomic ancestry. Samples from 90 vivax malaria-infected individuals and 51 noninfected individuals from northern Brazil were evaluated. Genotyping was carried out by using ASO-PCR or PCR/RFLP. The genomic ancestry of the individuals was classified using 48 insertion/deletion polymorphism biallelic markers. There were no differences in the proportions of African, European, and Native American ancestry between men and women. No significant association was observed for the allele and genotype frequencies of the 6 SNPs between malaria-infected and noninfected individuals. However, there was a trend toward decreasing the frequency of individuals carrying the TNF-308A allele with the increasing proportion of European ancestry. No ethnic-specific SNPs were identified, and there was no allelic or genotype association with susceptibility or resistance to vivax malaria. Understanding the genomic mechanisms by which ancestry influences this association is critical and requires further study.

Common Genetic Polymorphisms Influence Blood Biomarker Measurements in COPD

Implementing precision medicine for complex diseases such as chronic obstructive lung disease (COPD) will require extensive use of biomarkers and an in-depth understanding of how genetic, epigenetic, and environmental variations contribute to phenotypic diversity and disease progression. A meta-analysis from two large cohorts of current and former smokers with and without COPD [SPIROMICS (N = 750); COPDGene (N = 590)] was used to identify single nucleotide polymorphisms (SNPs) associated with measurement of 88 blood proteins (protein quantitative trait loci; pQTLs). PQTLs consistently replicated between the two cohorts. Features of pQTLs were compared to previously reported expression QTLs (eQTLs). Inference of causal relations of pQTL genotypes, biomarker measurements, and four clinical COPD phenotypes (airflow obstruction, emphysema, exacerbation history, and chronic bronchitis) were explored using conditional independence tests. We identified 527 highly significant (p < 8 X 10-10) pQTLs in 38 (43%) of blood proteins tested. Most pQTL SNPs were novel with low overlap to eQTL SNPs. The pQTL SNPs explained >10% of measured variation in 13 protein biomarkers, with a single SNP (rs7041; p = 10-392) explaining 71%-75% of the measured variation in vitamin D binding protein (gene = GC). Some of these pQTLs [e.g., pQTLs for VDBP, sRAGE (gene = AGER), surfactant protein D (gene = SFTPD), and TNFRSF10C] have been previously associated with COPD phenotypes. Most pQTLs were local (cis), but distant (trans) pQTL SNPs in the ABO blood group locus were the top pQTL SNPs for five proteins. The inclusion of pQTL SNPs improved the clinical predictive value for the established association of sRAGE and emphysema, and the explanation of variance (R2) for emphysema improved from 0.3 to 0.4 when the pQTL SNP was included in the model along with clinical covariates. Causal modeling provided insight into specific pQTL-disease relationships for airflow obstruction and emphysema. In conclusion, given the frequency of highly significant local pQTLs, the large amount of variance potentially explained by pQTL, and the differences observed between pQTLs and eQTLs SNPs, we recommend that protein biomarker-disease association studies take into account the potential effect of common local SNPs and that pQTLs be integrated along with eQTLs to uncover disease mechanisms. Large-scale blood biomarker studies would also benefit from close attention to the ABO blood group.

Genetic alterations of δ-catenin/NPRAP/Neurojungin (CTNND2): functional implications in complex human diseases

Some genes involved in complex human diseases are particularly vulnerable to genetic variations such as single nucleotide polymorphism, copy number variations, and mutations. For example, Ras mutations account for over 30 % of all human cancers. Additionally, there are some genes that can display different variations with functional impact in different diseases that are unrelated. One such gene stands out: δ-catenin/NPRAP/Neurojungin with gene designation as CTNND2 on chromosome 5p15.2. Recent advances in genome wide association as well as molecular biology approaches have uncovered striking involvement of δ-catenin gene variations linked to complex human disorders. These disorders include cancer, bipolar disorder, schizophrenia, autism, Cri-du-chat syndrome, myopia, cortical cataract-linked Alzheimer's disease, and infectious diseases. This list has rapidly grown longer in recent years, underscoring the pivotal roles of δ-catenin in critical human diseases. δ-Catenin is an adhesive junction-associated protein in the delta subfamily of the β-catenin superfamily. δ-Catenin functions in Wnt signaling to regulate gene expression and modulate Rho GTPases of the Ras superfamily in cytoskeletal reorganization. δ-Catenin likely lies where Wnt signaling meets Rho GTPases and is a unique and vulnerable common target for mutagenesis in different human diseases.

Admixture into and within sub-Saharan Africa

Similarity between two individuals in the combination of genetic markers along their chromosomes indicates shared ancestry and can be used to identify historical connections between different population groups due to admixture. We use a genome-wide, haplotype-based, analysis to characterise the structure of genetic diversity and gene-flow in a collection of 48 sub-Saharan African groups. We show that coastal populations experienced an influx of Eurasian haplotypes over the last 7000 years, and that Eastern and Southern Niger-Congo speaking groups share ancestry with Central West Africans as a result of recent population expansions. In fact, most sub-Saharan populations share ancestry with groups from outside of their current geographic region as a result of gene-flow within the last 4000 years. Our in-depth analysis provides insight into haplotype sharing across different ethno-linguistic groups and the recent movement of alleles into new environments, both of which are relevant to studies of genetic epidemiology.

DOI:http://dx.doi.org/10.7554/eLife.15266.001

Our genomes contain a record of historical events. This is because when groups of people are separated for generations, the DNA sequence in the two groups’ genomes will change in different ways. Looking at the differences in the genomes of people from the same population can help researchers to understand and reconstruct the historical interactions that brought their ancestors together. The mixing of two populations that were previously separate is known as admixture.

Africa as a continent has few written records of its history. This means that it is somewhat unknown which important movements of people in the past generated the populations found in modern-day Africa. Busby et al. have now attempted to use DNA to look into this and reconstruct the last 4000 years of genetic history in African populations.

As has been shown in other regions of the world, the new analysis showed that all African populations are the result of historical admixture events. However, Busby et al. could characterize these events to unprecedented level of detail. For example, multiple ethnic groups from The Gambia and Mali all show signs of sharing the same set of ancestors from West Africa, Europe and Asia who mixed around 2000 years ago. Evidence of a migration of people from Central West Africa, known as the Bantu expansion, could also be detected, and was shown to carry genes to the south and east. An important next step will be to now look at the consequences of the observed gene-flow, and ask if it has contributed to spreading beneficial, or detrimental, mutations around Africa.

DOI:http://dx.doi.org/10.7554/eLife.15266.002

Polymorphism in a lincRNA Associates with a Doubled Risk of Pneumococcal Bacteremia in Kenyan Children

Bacteremia (bacterial bloodstream infection) is a major cause of illness and death in sub-Saharan Africa but little is known about the role of human genetics in susceptibility. We conducted a genome-wide association study of bacteremia susceptibility in more than 5,000 Kenyan children as part of the Wellcome Trust Case Control Consortium 2 (WTCCC2). Both the blood-culture-proven bacteremia case subjects and healthy infants as controls were recruited from Kilifi, on the east coast of Kenya. Streptococcus pneumoniae is the most common cause of bacteremia in Kilifi and was thus the focus of this study. We identified an association between polymorphisms in a long intergenic non-coding RNA (lincRNA) gene (AC011288.2) and pneumococcal bacteremia and replicated the results in the same population (p combined = 1.69 × 10(-9); OR = 2.47, 95% CI = 1.84-3.31). The susceptibility allele is African specific, derived rather than ancestral, and occurs at low frequency (2.7% in control subjects and 6.4% in case subjects). Our further studies showed AC011288.2 expression only in neutrophils, a cell type that is known to play a major role in pneumococcal clearance. Identification of this novel association will further focus research on the role of lincRNAs in human infectious disease.