Detecting gene-gene interactions from GWAS using diffusion kernel principal components.

Genes and gene products do not function in isolation but as components of complex networks of macromolecules through physical or biochemical interactions. Dependencies of gene mutations on genetic background (i.e., epistasis) are believed to play a role in understanding molecular underpinnings of complex diseases such as inflammatory bowel disease (IBD). However, the process of identifying such interactions is complex due to for instance the curse of high dimensionality, dependencies in the data and non-linearity. Here, we propose a novel approach for robust and computationally efficient epistasis detection. We do so by first reducing dimensionality, per gene via diffusion kernel principal components (kpc). Subsequently, kpc gene summaries are used for downstream analysis including the construction of a gene-based epistasis network. We show that our approach is not only able to recover known IBD associated genes but also additional genes of interest linked to this difficult gastrointestinal disease.

Deletions of pfhrp2 and pfhrp3 genes were uncommon in rapid diagnostic test-negative Plasmodium falciparum isolates from Uganda.

BACKGROUND: Rapid diagnostic tests (RDTs) play a key role in malaria case management. The most widely used RDT identifies Plasmodium falciparum based on immunochromatographic recognition of P. falciparum histidine-rich protein 2 (PfHRP2). Deletion of the paralogous pfhrp2 and pfhrp3 genes leads to false-negative PfHRP2-based RDTs, and has been reported in P. falciparum infections from South America and Africa. However, identification of pfhrp2/pfhrp3 deletions has usually been based only on failure to amplify these genes using PCR, without confirmation based on PfHRP2 protein expression, and understanding of the true prevalence of deletions is incomplete. METHODS: Deletions of pfhrp2/pfhrp3 in blood samples were investigated from cross-sectional surveys in 2012-13 in three regions of varied malaria transmission intensity in Uganda. Samples with positive Giemsa-stained thick blood smears, but negative PfHRP2-based RDTs were evaluated by PCR amplification of conserved subunit ribosomal DNA for Plasmodium species, PCR amplification of pfhrp2 and pfhrp3 genes to identify deletions, and bead-based immunoassays for expression of PfHRP2. RESULTS: Of 3516 samples collected in cross-sectional surveys, 1493 (42.5%) had positive blood smears, of which 96 (6.4%) were RDT-negative. Of these 96 RDT-negative samples, P. falciparum DNA was identified by PCR in 56 (58%) and only non-falciparum plasmodial DNA in 40 (42%). In all 56 P. falciparum-positive samples there was a failure to amplify pfhrp2 or pfhrp3: in 25 (45%) pfhrp2 was not amplified, in 39 (70%) pfhrp3 was not amplified, and in 19 (34%) neither gene was amplified. For the 39 P. falciparum-positive, RDT-negative samples available for analysis of protein expression, PfHRP2 was not identified by immunoassay in only four samples (10.3%); these four samples all had failure to amplify both pfhrp2 and pfhrp3 by PCR. Thus, only four of 96 (4.2%) smear-positive, RDT-negative samples had P. falciparum infections with deletion of pfhrp2 and pfhrp3 confirmed by failure to amplify the genes by PCR and lack of expression of PfHRP2 demonstrated by immunoassay. CONCLUSION: False negative RDTs were uncommon. Deletions in pfhrp2 and pfhrp3 explained some of these false negatives, but most false negatives were not due to deletion of the pfhrp2 and pfhrp3 genes.

Associations between red blood cell variants and malaria among children and adults from three areas of Uganda: a prospective cohort study.

BACKGROUND: Multiple red blood cell (RBC) variants appear to offer protection against the most severe forms of Plasmodium falciparum malaria. Associations between these variants and uncomplicated malaria are less clear. METHODS: Data from a longitudinal cohort study conducted in 3 sub-counties in Uganda was used to quantify associations between three red blood cell variants Hb [AA, AS, S (rs334)], alpha thalassaemia 3.7 kb deletion, and glucose-6-phosphate dehydrogenase deficiency A-(G6PD 202A genotype) and malaria incidence, parasite prevalence, parasite density (a measure of anti-parasite immunity) and body temperature adjusted for parasite density (a measure of anti-disease immunity). All analyses were adjusted for age, average household entomological inoculation rate, and study site. Results for all variants were compared to those for wild type genotypes. RESULTS: In children, HbAS was associated, compared to wild type, with a lower incidence of malaria (IRR = 0.78, 95% CI 0.66-0.92, p = 0.003), lower parasite density upon infection (PR = 0.66, 95% CI 0.51-0.85, p = 0.001), and lower body temperature for any given parasite density (- 0.13 â„ƒ, 95% CI - 0.21, - 0.05, p = 0.002). In children, HbSS was associated with a lower incidence of malaria (IRR = 0.17, 95% CI 0.04-0.71, p = 0.02) and lower parasite density upon infection (PR = 0.31, 95% CI 0.18-0.54, p < 0.001). α-/αα thalassaemia, was associated with higher parasite prevalence in both children and adults (RR = 1.23, 95% CI 1.06-1.43, p = 0.008 and RR = 1.52, 95% CI 1.04-2.23, p = 0.03, respectively). G6PD deficiency was associated with lower body temperature for any given parasite density only among male hemizygote children (- 0.19 â„ƒ, 95% CI - 0.31, - 0.06, p = 0.003). CONCLUSION: RBC variants were associated with non-severe malaria outcomes. Elucidation of the mechanisms by which they confer protection will improve understanding of genetic protection against malaria.

Changing Antimalarial Drug Resistance Patterns Identified by Surveillance at Three Sites in Uganda.

We assessed Plasmodium falciparum drug resistance markers in parasites collected in 2012, 2013, and 2015 at 3 sites in Uganda. The prevalence and frequency of parasites with mutations in putative transporters previously associated with resistance to aminoquinolines, but increased sensitivity to lumefantrine (pfcrt 76T; pfmdr1 86Y and 1246Y), decreased markedly at all sites. Antifolate resistance mutations were common, with apparent emergence of mutations (pfdhfr 164L; pfdhps 581G) associated with high-level resistance. K13 mutations linked to artemisinin resistance were uncommon and did not increase over time. Changing malaria treatment practices have been accompanied by profound changes in markers of resistance.

Impact of antimalarial treatment and chemoprevention on the drug sensitivity of malaria parasites isolated from ugandan children.

Changing treatment practices may be selecting for changes in the drug sensitivity of malaria parasites. We characterized ex vivo drug sensitivity and parasite polymorphisms associated with sensitivity in 459 Plasmodium falciparum samples obtained from subjects enrolled in two clinical trials in Tororo, Uganda, from 2010 to 2013. Sensitivities to chloroquine and monodesethylamodiaquine varied widely; sensitivities to quinine, dihydroartemisinin, lumefantrine, and piperaquine were generally good. Associations between ex vivo drug sensitivity and parasite polymorphisms included decreased chloroquine and monodesethylamodiaquine sensitivity and increased lumefantrine and piperaquine sensitivity with pfcrt 76T, as well as increased lumefantrine sensitivity with pfmdr1 86Y, Y184, and 1246Y. Over time, ex vivo sensitivity decreased for lumefantrine and piperaquine and increased for chloroquine, the prevalences of pfcrt K76 and pfmdr1 N86 and D1246 increased, and the prevalences of pfdhfr and pfdhps polymorphisms associated with antifolate resistance were unchanged. In recurrent infections, recent prior treatment with artemether-lumefantrine was associated with decreased ex vivo lumefantrine sensitivity and increased prevalence of pfcrt K76 and pfmdr1 N86, 184F, and D1246. In children assigned chemoprevention with monthly dihydroartemisinin-piperaquine with documented circulating piperaquine, breakthrough infections had increased the prevalence of pfmdr1 86Y and 1246Y compared to untreated controls. The noted impacts of therapy and chemoprevention on parasite polymorphisms remained significant in multivariate analysis correcting for calendar time. Overall, changes in parasite sensitivity were consistent with altered selective pressures due to changing treatment practices in Uganda. These changes may threaten the antimalarial treatment and preventive efficacies of artemether-lumefantrine and dihydroartemisinin-piperaquine, respectively.

Integrative computational modeling to unravel novel potential biomarkers in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a major health problem around the world. The management of this disease is complicated by the lack of noninvasive diagnostic tools and the few treatment options available. Better clinical outcomes can be achieved if HCC is detected early, but unfortunately, clinical signs appear when the disease is in its late stages. We aim to identify novel genes that can be targeted for the diagnosis and therapy of HCC. We performed a meta-analysis of transcriptomics data to identify differentially expressed genes and applied network analysis to identify hub genes. Fatty acid metabolism, complement and coagulation cascade, chemical carcinogenesis and retinol metabolism were identified as key pathways in HCC. Furthermore, we integrated transcriptomics data into a reference human genome-scale metabolic model to identify key reactions and subsystems relevant in HCC. We conclude that fatty acid activation, purine metabolism, vitamin D, and E metabolism are key processes in the development of HCC and therefore need to be further explored for the development of new therapies. We provide the first evidence that GABRP, HBG1 and DAK (TKFC) genes are important in HCC in humans and warrant further studies.

Guided extraction of genome-scale metabolic models for the integration and analysis of omics data.

Omics data can be integrated into a reference model using various model extraction methods (MEMs) to yield context-specific genome-scale metabolic models (GEMs). How to chose the appropriate MEM, thresholding rule and threshold remains a challenge. We integrated mouse transcriptomic data from a Cyp51 knockout mice diet experiment (GSE58271) using five MEMs (GIMME, iMAT, FASTCORE, INIT an tINIT) in a combination with a recently published mouse GEM iMM1865. Except for INIT and tINIT, the size of extracted models varied with the MEM used (t-test: p-value < 0.001). The Jaccard index of iMAT models ranged from 0.27 to 1.0. Out of the three factors under study in the experiment (diet, gender and genotype), gender explained most of the variability ( > 90%) in PC1 for FASTCORE. In iMAT, each of the three factors explained less than 40% of the variability within PC1, PC2 and PC3. Among all the MEMs, FASTCORE captured the most of the true variability in the data by clustering samples by gender. Our results show that for the efficient use of MEMs in the context of omics data integration and analysis, one should apply various MEMs, thresholding rules, and thresholding values to select the MEM and its configuration that best captures the true variability in the data. This selection can be guided by the methodology as proposed and used in this paper. Moreover, we describe certain approaches that can be used to analyse the results obtained with the selected MEM and to put these results in a biological context.

Associations between erythrocyte polymorphisms and risks of uncomplicated and severe malaria in Ugandan children: A case control study.

BACKGROUND: Evidence for association between sickle cell and alpha thalassemia trait and severe malaria is compelling. However, for these polymorphisms associations with uncomplicated malaria, and for G6PD deficiency associations with uncomplicated and severe malaria, findings have been inconsistent. We studied samples from a three-arm case-control study with the objective of determining associations between common host erythrocyte polymorphisms and both uncomplicated and severe malaria, including different severe malaria phenotypes. METHOD: We assessed hemoglobin abnormalities, α-thalassemia, and G6PD deficiency by molecular methods in 325 children with severe malaria age-matched to 325 children with uncomplicated malaria and 325 healthy community controls. Conditional logistic regression was used to measure associations between specified genotypes and malaria outcomes. RESULTS: No tested polymorphisms offered significant protection against uncomplicated malaria. α-thalassemia homozygotes (_α/_α) had increased risk of uncomplicated malaria (OR 2.40; 95%CI 1.15, 5.03, p = 0.020). HbAS and α-thalassemia heterozygous (_α/αα) genotypes protected against severe malaria compared to uncomplicated malaria (HbAS OR 0.46; 0.23, 0.95, p = 0.036; _α/αα OR 0.51; 0.24, 0.77; p = 0.001) or community (HbAS OR 0.23; 0.11, 0.50; p<0.001; _α/αα; OR 0.49; 0.32, 0.76; p = 0.002) controls. The α-thalassemia homozygous (_α/_α) genotype protected against severe malaria when compared to uncomplicated malaria controls (OR 0.34; 95%CI 0.156, 0.73, p = 0.005), but not community controls (OR 1.03; 0.46, 2.27, p = 0.935). Stratifying by the severe malaria phenotype, compared to community controls, the protective effect of HbAS was limited to children with severe anemia (OR 0.17; 95%CI 0.04, 0.65; p = 0.009) and that of _α/αα to those with altered consciousness (OR 0.24; 0.09, 0.59; p = 0.002). A negative epistatic effect was seen between HbAS and _α/αα; protection compared to uncomplicated malaria controls was not seen in individuals with both polymorphisms (OR 0.45; 0.11, 1.84; p = 0.269). G6PD deficiency was not protective against severe malaria. CONCLUSION: Associations were complex, with HbAS principally protective against severe anemia, _α/αα against altered consciousness, and negative epistasis between the two polymorphisms.

Plasmodium Species Infecting Children Presenting with Malaria in Uganda.

Contributions of species other than Plasmodium falciparum to human malaria in sub-Saharan Africa are uncertain. We collected blood from children aged 6 months to 10 years diagnosed with malaria by Giemsa-stained blood smears (176 subjects) or histidine rich protein-2-based rapid diagnostic tests (323 subjects) in 2016; 50 samples from each of 10 sites across Uganda were studied to identify infecting species. Of 499 available samples, 474 demonstrated plasmodial infection by polymerase chain reaction amplification of 18S ribosomal RNA genes, including P. falciparum in 472, Plasmodium malariae in 22, Plasmodium ovale in 15, and Plasmodium vivax in four; 435 were pure P. falciparum, two did not contain P. falciparum, and the remainder were mixed infections including P. falciparum. The prevalence of nonfalciparum species varied geographically. Stratifying based on recent history of indoor residual spraying (IRS) of insecticides, nonfalciparum infections were seen in 27/189 (14.8%) samples from sites that received and 13/285 (4.6%) samples from sites that did not receive IRS since 2010 (P = 0.0013). Overall, 39/474 (8.2%) samples from individuals diagnosed with malaria included nonfalciparum infections. Thus, a substantial proportion of episodes of malaria in Uganda include infections with plasmodial species other than P. falciparum.

Marked variation in prevalence of malaria-protective human genetic polymorphisms across Uganda.

A number of human genetic polymorphisms are prevalent in tropical populations and appear to offer protection against symptomatic and/or severe malaria. We compared the prevalence of four polymorphisms, the sickle hemoglobin mutation (ß globin E6V), the α-thalassemia 3.7kb deletion, glucose-6-phosphate dehydrogenase deficiency caused by the common African variant (G6PD A-), and the CD36 T188G mutation in 1344 individuals residing in districts in eastern (Tororo), south-central (Jinja), and southwestern (Kanungu) Uganda. Genes of interest were amplified, amplicons subjected to mutation-specific restriction endonuclease digestion (for sickle hemoglobin, G6PD A-, and CD36 T188G), reaction products resolved by electrophoresis, and genotypes determined based on the sizes of reaction products. Mutant genotypes were common, with many more heterozygous than homozygous alleles identified. The prevalences (heterozygotes plus homozygotes) of sickle hemoglobin (28% Tororo, 25% Jinja, 7% Kanungu), α-thalassemia (53% Tororo, 45% Jinja, 18% Kanungu) and G6PD A- (29% Tororo, 18% Jinja, 8% Kanungu) were significantly greater in Tororo and Jinja compared to Kanungu (p<0.0001 for all three alleles); prevalences were also significantly greater in Tororo compared to Jinja for α-thalassemia (p=0.03) and G6PD A- (p<0.0001). For the CD36 T188G mutation, the prevalence was significantly greater in Tororo compared to Jinja or Kanungu (27% Tororo, 17% Jinja, 18% Kanungu; p=0.0004 and 0.0017, respectively). Considering ethnicity of study subjects, based on primary language spoken, the prevalence of mutant genotypes was lower in Bantu compared to non-Bantu language speakers, but in the Jinja cohort, the only study population with a marked diversity of language groups, prevalence did not differ between Bantu and non-Bantu speakers. These results indicate marked differences in human genetic features between populations in different regions of Uganda. These differences might be explained by both ethnic variation and by varied malaria risk in different regions of Uganda.

Intermittent preventive therapy with sulfadoxine-pyrimethamine for malaria in pregnancy: a cross-sectional study from Tororo, Uganda.

BACKGROUND: Intermittent preventive treatment during pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is widely recommended in sub-Saharan Africa to reduce the risk of malaria and improve birth outcomes. However, there are reports that the efficacy of IPTp with SP is waning, especially in parts of Africa where antimalarial resistance to this drug has become widespread. METHODOLOGY/PRINCIPAL FINDINGS: We conducted a cross-sectional study of 565 HIV-uninfected women giving birth at Tororo District Hospital in southeastern Uganda. The primary objective of the study was to measure associations between use of SP during pregnancy from antenatal records and the risk of adverse outcomes including placental malaria, low birth weight, maternal parasitemia and maternal anemia. The proportion of women who reported taking 0, 1, 2, and 3 doses of SP during pregnancy was 5.7%, 35.8%, 56.6% and 2.0% respectively. Overall, the prevalence of placental malaria was 17.5%, 28.1%, and 66.2% by placental smear, PCR, and histopathology, respectively. In multivariate analyses controlling for potential confounders, ≥ 2 doses of SP was associated with non-significant trends towards lower odds of placental malaria by placental smear (OR = 0.75, p = 0.25), placental malaria by PCR (OR = 0.93, p = 0.71), placental malaria by histopathology (OR = 0.75, p = 0.16), low birth weight (OR = 0.63, p = 0.11), maternal parasitemia (OR = 0.88, p = 0.60) and maternal anemia (OR = 0.88, p = 0.48). Using a composite outcome, ≥ 2 doses of SP was associated with a significantly lower odds of placental malaria, low birth weight, maternal parasitemia, or maternal anemia (OR = 0.52, p = 0.01). CONCLUSIONS/SIGNIFICANCE: In this area of Uganda with intense malaria transmission, the prevalence of placental malaria by histopathology was high even among women who reported taking at least 2 doses of SP during pregnancy. The reported use of ≥ 2 doses of SP was not associated with protection against individual birth and maternal outcome measures but did protect against a composite measure of any adverse outcome.