Search

The persecutions of the Jews that began with legislation introduced by Italy's fascist government in the year 1938 ("Leggi Razziali" i.e. "Racial Laws") also affected the sphere of anatomic pathology, coming to bear on Italian physicians belonging to the Jewish communities of several cities and universities. The damage caused by the discrimination against them and their removal from their jobs penetrated a public health world that had hitherto been based on a climate of tolerance and integration. Here we recall some emblematic figures involved in those troubled times in Italy's history.

The Penn Medicine BioBank (PMBB) is an electronic health record (EHR)-linked biobank at the University of Pennsylvania (Penn Medicine). A large variety of health-related information, ranging from diagnosis codes to laboratory measurements, imaging data and lifestyle information, is integrated with genomic and biomarker data in the PMBB to facilitate discoveries and translational science. To date, 174,712 participants have been enrolled into the PMBB, including approximately 30% of participants of non-European ancestry, making it one of the most diverse medical biobanks. There is a median of seven years of longitudinal data in the EHR available on participants, who also consent to permission to recontact. Herein, we describe the operations and infrastructure of the PMBB, summarize the phenotypic architecture of the enrolled participants, and use body mass index (BMI) as a proof-of-concept quantitative phenotype for PheWAS, LabWAS, and GWAS. The major representation of African-American participants in the PMBB addresses the essential need to expand the diversity in genetic and translational research. There is a critical need for a “medical biobank consortium” to facilitate replication, increase power for rare phenotypes and variants, and promote harmonized collaboration to optimize the potential for biological discovery and precision medicine.

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Late-Onset Neonatal Sepsis (LOS) is a rare, but life-threatening condition, involving widespread infection, immune disruption, organ dysfunction, and often death. Because exposure to pathogens is not 100% preventable, identifying susceptibility factors is critical to defining neonates at greater risk. Prior studies demonstrated that both genetics and infant sex influence susceptibility. We, therefore, performed a genome wide association study (GWAS) with 224 cases and 273 controls from six European countries to identify genetic risk factors. We tested for association with both autosomal and X-chromosome variants in the total sample and in the samples stratified by sex. In total 71 SNPs associated with neonatal sepsis at p−4in at least one analysis. Most importantly, the sex stratified analyses revealed associations with multiple SNPs (28 SNPs in males and 16 in females), but none of 44 SNPs from single-sex analyses associated with sepsis in the other sex at p

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4,andLY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. AtACE2,we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencingACE2expression in multiple African populations. AtTMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Host genetic factors can confer resistance against malaria, raising the question of whether this has led to evolutionary adaptation of parasite populations. In this study we investigated the correlation between host and parasite genetic variation in 4,171 Gambian and Kenya children ascertained with severe malaria due to Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and variation in three regions of the parasite genome, including nonsynonymous variants in the acyl-CoA synthetase family member PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The HbS-associated parasite alleles are in strong linkage disequilibrium and have frequencies which covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome, and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

SARS-CoV-2 enters host cells by binding angiotensin-converting enzyme 2 (ACE2). Through a genome-wide association study, we show that a rare variant (MAF = 0.3%, odds ratio 0.60, P=4.5x10-13) that down-regulates ACE2 expression reduces risk of COVID-19 disease, providing human genetics support for the hypothesis that ACE2 levels influence COVID-19 risk. Further, we show that common genetic variants define a risk score that predicts severe disease among COVID-19 cases.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease-19 (COVID-19), a respiratory illness that can result in hospitalization or death. We investigated associations between rare genetic variants and seven COVID-19 outcomes in 543,213 individuals, including 8,248 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome-wide or when specifically focusing on (i) 14 interferon pathway genes in which rare deleterious variants have been reported in severe COVID-19 patients; (ii) 167 genes located in COVID-19 GWAS risk loci; or (iii) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, with results publicly browsable athttps://rgc-covid19.regeneron.com.

Background Tuberculosis (TB) is the most deadly infectious disease globally and is highly prevalent in the developing world. For individuals infected with both Mycobacterium tuberculosis (Mtb) and human immunodeficiency virus (HIV), the risk of active TB is 10% or more annually. Previously, we identified in a genome-wide association study (GWAS) a region on chromosome 5 associated with resistance to TB, which included epigenetic marks that could influence gene regulation. We hypothesized that HIV-infected individuals exposed to Mtb who remain disease free carry epigenetic changes that strongly protect them from active TB. Methods We conducted a methylome-wide study in HIV-infected, TB-exposed cohorts from Uganda and Tanzania and integrated data from our GWAS. Results We identified 3 regions of interest that included markers that were differentially methylated between TB cases and controls with latent TB infection: chromosome 1 (RNF220, P = 4 × 10–5), chromosome 2 (between COPS8 and COL6A3, P = 2.7 × 10–5), and chromosome 5 (CEP72, P = 1.3 × 10–5). These methylation results co-localized with associated single-nucleotide polymorphisms (SNPs), methylation QTLs, and methylation × SNP interaction effects. These markers were in regions with regulatory markers for cells involved in TB immunity and/or lung. Conclusions Epigenetic regulation is a potential biologic factor underlying resistance to TB in immunocompromised individuals that can act in conjunction with genetic variants.

BackgroundTuberculosis (TB) is the most deadly infectious disease globally and highly prevalent in the developing world, especially sub-Saharan Africa. Even though a third of humans are exposed to Myocbacterium tuberculosis (Mtb), most infected immunocompetent individuals do not develop active TB. In contrast, for individuals infected with both TB and the human immunodeficiency virus (HIV), the risk of active disease is 10% or more per year. Previously, we identified in a genome-wide association study a region on chromosome 5 that was associated with resistance to TB. This region included epigenetic marks that could influence gene regulation so we hypothesized that HIV-infected individuals exposed to Mtb, who remain disease free, carry epigenetic changes that strongly protect them from active TB. To test this hypothesis, we conducted a methylome-wide study in HIV-infected, TB-exposed cohorts from Uganda and Tanzania.ResultsIn 221 HIV-infected adults from Uganda and Tanzania, we identified 3 regions of interest that included markers that were differentially methylated between TB cases and LTBI controls, that also included methylation QTLs and associated SNPs: chromosome 1 (RNF220, p=4×10−5), chromosome 2 (between COPS8 and COL6A3 genes, p=2.7×10−5), and chromosome 5 (CEP72, p=1.3×10−5). These methylation results colocalized with associated SNPs, methylation QTLs, and methylation × SNP interaction effects. These markers were in regions with regulatory markers for cells involved in TB immunity and/or lung.ConclusionEpigenetic regulation is a potential biologic factor underlying resistance to TB in immunocompromised individuals that can act in conjunction with genetic variants.

Tuberculosis (TB) is the leading cause of death from a single infectious agent. According to the WHO, 85% of cases in 2018 were pulmonary tuberculosis (PTB), making it the most prevalent form of the disease. Although the bacillus responsible for disease, Mycobacterium tuberculosis (MTB), is estimated to infect 1.7 billion people worldwide, only a small portion of those infected (5-10%) will transition into active TB. Because such a small fraction of infected people develop active disease, we hypothesized that underlying host genetic variation associates with developing active pulmonary disease. Variation in CLEC4E has been of interest in previous association studies showing either no effect or protection from PTB. For our study we assessed 60 SNPs in 11 immune genes, including CLEC4E, using a case-control study from Guinea-Bissau. The 289 cases and 322 controls differed in age, sex, and ethnicity all of which were included in adjusted models. Initial association analysis with unadjusted logistic regression revealed putative association with nine SNPs (p < 0.1). Any SNP that was significant below a p-value of 0.05 was then assessed in an adjusted model. Of the six SNPs that remained significant, three of them were assigned to the CLEC4E gene (rs12302046, rs10841847, and rs11046143). Of these, only rs10841847 passed FDR adjustment for multiple testing. Adjusted regression analyses showed that the minor allele at rs10841847 associated with higher risk of developing PTB (OR = 1.55, CI = 1.22-1.96, p-value = 0.00036). Based on these initial association tests, CLEC4E seemed to be the predictor of interest for PTB risk in this population. Haplotype analysis (2-SNP and 3-SNP windows) showed that minor alleles in segments including rs10841847 were the only ones to pass the threshold of global significance, compared to other haplotypes (p-value < 0.05). Linkage disequilibrium patterns showed that rs12302046 is in high LD with rs10841847 (r2 = 0.67), and all other SNPs lost significance when adjusted for rs10841847 effects. These findings indicate that rs10841847 in CLEC4E is the single best predictor of pulmonary tuberculosis risk in our study population. These results provide evidence for the hypothesis that genetic variation of CLEC4E influences risk to TB in Guinea-Bissau.

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), remains a major public health threat globally. Several lines of evidence support a role for host genetic factors in resistance/susceptibility to TB disease and MTB infection. However, results across candidate gene and genome-wide association studies (GWAS) are largely inconsistent, so a cohesive genetic model underlying TB risk has not emerged. Despite the difficulties in identifying consistent genetic associations, genetic studies of TB and MTB infection have revealed a few well-documented loci. These well-validated genes are presented in this review, but there remains a large gap in how these genes translate into better understanding of TB. To address this, we present a pathway-based extension of standard association analyses, seeding the results with the best validated genes from candidate gene and GWAS studies. Several pathways were significantly enriched using pathway analyses that may help to explain population patterns of TB risk. In conclusion, we advocate for novel approaches to the study of host genetic analysis of TB that extend traditional association approaches.

We conducted a genome-wide association study of host resistance to severePlasmodium falciparummalaria in over 17,000 individuals from 11 malaria-endemic countries, undertaking a wide ranging analysis which identifies five replicable associations with genome-wide levels of evidence. Our findings include a newly implicated variant on chromosome 6 associated with risk of cerebral malaria, and the discovery of an erythroid-specific transcription start site underlying the association inATP2B4. Previously reported HLA associations cannot be replicated in this dataset. We estimate substantial heritability of severe malaria (h2~ 23%), of which around 10% is explained by the currently identified associations. Our dataset will provide a major building block for future research on the genetic determinants of disease in these diverse human populations.

Aminoglycosides such as amikacin continue to be part of the backbone of treatment of multidrug-resistant tuberculosis (MDR-TB). We measured amikacin concentrations in 28 MDR-TB patients in Botswana receiving amikacin therapy together with oral levofloxacin, ethionamide, cycloserine, and pyrazinamide and calculated areas under the concentration-time curves from 0 to 24 h (AUC 0–24 ). The patients were followed monthly for sputum culture conversion based on liquid cultures. The median duration of amikacin therapy was 184 (range, 28 to 866) days, at a median dose of 17.30 (range 11.11 to 19.23) mg/kg. Only 11 (39%) patients had sputum culture conversion during treatment; the rest failed. We utilized classification and regression tree analyses (CART) to examine all potential predictors of failure, including clinical and demographic features, comorbidities, and amikacin peak concentrations ( C max ), AUC 0–24 , and trough concentrations. The primary node for failure had two competing variables, C max of 0–24 of 41 kg was a secondary node with a score of 35% relative to the primary node. The area under the receiver operating characteristic curve for the CART model was an R 2 = 0.90 on posttest. In patients weighing >41 kg, sputum conversion was 3/3 (100%) in those with an amikacin C max of ≥67 mg/liter versus 3/15 (20%) in those with a C max of C max and AUC 0–24 below the threshold, 7/7 (100%) failed, compared to 7/15 (47%) of those who had these parameters above threshold (RR = 2.14; 95% CI, 1.25 to 43.68). These amikacin dose-schedule patterns and exposures are virtually the same as those identified in the hollow-fiber system model.

Genetic studies of both the human host and Mycobacterium tuberculosis (MTB) demonstrate independent association with tuberculosis (TB) risk. However, neither explains a large portion of disease risk or severity. Based on studies in other infectious diseases and animal models of TB, we hypothesized that the genomes of the two interact to modulate risk of developing active TB or increasing the severity of disease, when present. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which there were 3 MTB lineages of which L4-Ugandan (L4.6) is the most recent. TB severity, measured using the Bandim TBscore, was modeled as a function of host SNP genotype, MTB lineage, and their interaction, within two independent cohorts of TB cases, N = 113 and 121. No association was found between lineage and severity, but association between multiple polymorphisms in IL12B and TBscore was replicated in two independent cohorts (most significant rs3212227, combined p = 0.0006), supporting previous associations of IL12B with TB susceptibility. We also observed significant interaction between a single nucleotide polymorphism (SNP) in SLC11A1 and the L4-Ugandan lineage in both cohorts (rs17235409, meta p = 0.0002). Interestingly, the presence of the L4-Uganda lineage in the presence of the ancestral human allele associated with more severe disease. These findings demonstrate that IL12B is associated with severity of TB in addition to susceptibility, and that the association between TB severity and human genetics can be due to an interaction between genes in the two species, consistent with host-pathogen coevolution in TB., Author summary Susceptibility to tuberculosis (TB) is affected by genetic variation in both the human host and the causative bacterium, Mycobacterium tuberculosis. However, prior studies of the genetics of each species have not explained a large part of TB risk. The possibility exists that risk can be better estimated from patterns of variation in the two species as a unit, such that some combinations provide increased risk, or in the presence of TB, increased disease severity. We hypothesized that alleles in the two species that have co-existed for long periods are more likely to reduce disease severity so as to promote prolonged co-occurrence. We tested this by studying TB severity in two patient cohorts from Uganda for which paired MTB-human DNA were available. We examined severity, as measured by the Bandim TBscore, and assessed whether there was an interaction between MTB lineage and SNPs in the host with this metric. Our results indicate that the most recent TB lineage (L4.6/Uganda) when found together with an ancestral allele in SLC11A1 resulted in more severe disease. This finding is consistent with the conclusion that MTB and human have coevolved to modulate TB severity.

IMPORTANCE: Hereditary transthyretin (TTR) amyloid cardiomyopathy (hATTR-CM) due to the TTR V122I variant is an autosomal-dominant disorder that causes heart failure in elderly individuals of African ancestry. The clinical associations of carrying the variant, its effect in other African ancestry populations including Hispanic/Latino individuals, and the rates of achieving a clinical diagnosis in carriers are unknown. OBJECTIVE: To assess the association between the TTR V122I variant and heart failure and identify rates of hATTR-CM diagnosis among carriers with heart failure. DESIGN, SETTING, AND PARTICIPANTS: Cross-sectional analysis of carriers and noncarriers of TTR V122I of African ancestry aged 50 years or older enrolled in the Penn Medicine Biobank between 2008 and 2017 using electronic health record data from 1996 to 2017. Case-control study in participants of African and Hispanic/Latino ancestry with and without heart failure in the Mount Sinai BioMe Biobank enrolled between 2007 and 2015 using electronic health record data from 2007 to 2018. EXPOSURES: TTR V122I carrier status. MAIN OUTCOMES AND MEASURES: The primary outcome was prevalent heart failure. The rate of diagnosis with hATTR-CM among TTR V122I carriers with heart failure was measured. RESULTS: The cross-sectional cohort included 3724 individuals of African ancestry with a median age of 64 years (interquartile range, 57-71); 1755 (47%) were male, 2896 (78%) had a diagnosis of hypertension, and 753 (20%) had a history of myocardial infarction or coronary revascularization. There were 116 TTR V122I carriers (3.1%); 1121 participants (30%) had heart failure. The case-control study consisted of 2307 individuals of African ancestry and 3663 Hispanic/Latino individuals; the median age was 73 years (interquartile range, 68-80), 2271 (38%) were male, 4709 (79%) had a diagnosis of hypertension, and 1008 (17%) had a history of myocardial infarction or coronary revascularization. There were 1376 cases of heart failure. TTR V122I was associated with higher rates of heart failure (cross-sectional cohort: n = 51/116 TTR V122I carriers [44%], n = 1070/3608 noncarriers [30%], adjusted odds ratio, 1.7 [95% CI, 1.2-2.4], P = .006; case-control study: n = 36/1376 heart failure cases [2.6%], n = 82/4594 controls [1.8%], adjusted odds ratio, 1.8 [95% CI, 1.2-2.7], P = .008). Ten of 92 TTR V122I carriers with heart failure (11%) were diagnosed as having hATTR-CM; the median time from onset of symptoms to clinical diagnosis was 3 years. CONCLUSIONS AND RELEVANCE: Among individuals of African or Hispanic/Latino ancestry enrolled in 2 academic medical center–based biobanks, the TTR V122I genetic variant was significantly associated with heart failure.

Aminoglycosides, such as amikacin, are used to treat multidrug-resistant tuberculosis. However, ototoxicity is a common problem and is monitored using peak and trough amikacin concentrations based on World Health Organization recommendations. Our objective was to identify clinical factors predictive of ototoxicity using an agnostic machine learning method. We used classification and regression tree (CART) analyses to identify clinical factors, including amikacin concentration thresholds that predicted audiometry-confirmed ototoxicity among 28 multidrug-resistant pulmonary tuberculosis patients in Botswana. Amikacin concentrations were measured for all patients. The quantitative relationship between predictive factors and the probability of ototoxicity were then identified using probit analyses. The primary predictors of ototoxicity on CART analyses were cumulative days of therapy, followed by cumulative area under the concentration-time curve (AUC), which improved on the primary predictor by 87%. The area under the receiver operating curve was 0.97 on the test set. Peak and trough were not predictors in any tree. When algorithms were forced to pick peak and trough as primary predictors, the area under the receiver operating curve fell to 0.46. Probit analysis revealed that the probability of ototoxicity increased sharply starting after 6 months of therapy to near maximum at 9 months. A 10% probability of ototoxicity occurred with a threshold cumulative AUC of 87,232 days · mg · h/liter, while that of 20% occurred at 120,000 days · mg · h/liter. Thus, cumulative amikacin AUC and duration of therapy, and not peak and trough concentrations, should be used as the primary decision-making parameters to minimize the likelihood of ototoxicity in multidrug-resistant tuberculosis.

Summary Objectives We determined the performance of a sensor array (an electronic nose) made of 8 metalloporphyrins coated quartz microbalances sensors for the diagnosis and prognosis of pulmonary tuberculosis (TB) using exhaled breath samples. Methods TB cases and healthy controls were prospectively enrolled. Signals from volatile organic compounds (VOCs) in breath samples were measured at days 0, 2, 7, 14, and 30 of TB therapy and correlated with clinical and microbiological measurements. Results Fifty one pulmonary TB cases and 20 healthy HIV-uninfected controls were enrolled in the study. 31 (61%) of the 51 pulmonary TB cases were coinfected with HIV. At day 0 (before TB treatment initiation) the sensitivity of our device was estimated at 94.1% (95% confidence interval [CI], 83.8–98.8%) and specificity was 90.0% (95% CI, 68.3–98.8%) for distinguishing TB cases from controls. Time-dependent changes in the breath signals were identified as time on TB treatment progressed. Time-dependent signal changes were more pronounced among HIV-uninfected patients. Conclusion The identification of VOCs' signals in breath samples using a sensor array achieved high sensitivity and specificity for the diagnosis of TB and allowed following signal changes during TB treatment.

The relationship between malaria infection and volatile compounds has been claimed mainly on the basis that they are believed to be an attractant for mosquitoes. However, since the association of emitted molecules with diseases has been observed for many pathologies, malaria-related volatile compounds are a potential diagnostic tool. The recent confirms of this hypothesis prompts the development of sensors for an effective exploitation of these potentialities. On these bases, we investigated the alteration of volatile compounds in a malaria murine model. For the scope, the total “volatilome” of Plasmodium berghei-infected mice was compared with that of non-infected animals. Gas chromatographic analysis of the sampled air reveals the existence of a pattern of compounds that, collectively considered, detects malaria infection. Finally, an array of porphyrins functionalized quartz microbalance gas sensors was applied to sort non-infected from infected mice. The application of a classification model to the sensor data provided more than 80% of correct identification with errors confined to mice with a low parasitemia level. Noteworthy, the sensor array was trained on data collected months before to run the tests. These results provide, although limited to a murine model, a first evidence of the potentialities of gas sensor technology for malaria diagnosis.

Many human genetic associations with resistance to malaria have been reported, but few have been reliably replicated. We collected data on 11,890 cases of severe malaria due to Plasmodium falciparum and 17,441 controls from 12 locations in Africa, Asia and Oceania. We tested 55 SNPs in 27 loci previously reported to associate with severe malaria. There was evidence of association at P < 1 × 10(-4) with the HBB, ABO, ATP2B4, G6PD and CD40LG loci, but previously reported associations at 22 other loci did not replicate in the multicenter analysis. The large sample size made it possible to identify authentic genetic effects that are heterogeneous across populations or phenotypes, with a striking example being the main African form of G6PD deficiency, which reduced the risk of cerebral malaria but increased the risk of severe malarial anemia. The finding that G6PD deficiency has opposing effects on different fatal complications of P. falciparum infection indicates that the evolutionary origins of this common human genetic disorder are more complex than previously supposed.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is frequent in Africa, because it confers resistance to Plasmodium falciparum malaria; however, the nature of the protection and the genotypes associated with it have been controversial. In 1972, Bienzle and others described protection from malaria in West African females heterozygous for G6PD A-. They determined that G6PD A- heterozygotes had lower parasite counts than A- homozygotes, hemizygous males, and normal individuals. However, other studies have reached different conclusions about the protective genotypes. DNA samples from 135 children with severe malaria and 146 children with mild malaria from The Gambia were genotyped for the G6PD A- mutation that is most frequent among Gambians (G6PD 968 T->C); there was a marked deficiency of heterozygotes and an excess of homozygotes with severe malaria, producing a strong deviation from Hardy–Weinberg equilibrium. Our results support the protective effect in G6PD A- heterozygous females and suggest that homozygotes might be more susceptible to severe malaria attacks.

The majority of studies of genetic association with disease have been performed in Europeans. This European bias has important implications for risk prediction of diseases across global populations. In this commentary, we justify the need to study more diverse populations using both empirical examples and theoretical reasoning.

Background Epistasis has been historically used to describe the phenomenon that the effect of a given gene on a phenotype can be dependent on one or more other genes, and is an essential element for understanding the association between genetic and phenotypic variations. Quantifying epistasis of orders higher than two is very challenging due to both the computational complexity of enumerating all possible combinations in genome-wide data and the lack of efficient and effective methodologies. Objectives In this study, we propose a fast, non-parametric, and model-free measure for three-way epistasis. Methods Such a measure is based on information gain, and is able to separate all lower order effects from pure three-way epistasis. Results Our method was verified on synthetic data and applied to real data from a candidate-gene study of tuberculosis in a West African population. In the tuberculosis data, we found a statistically significant pure three-way epistatic interaction effect that was stronger than any lower-order associations. Conclusion Our study provides a methodological basis for detecting and characterizing high-order gene-gene interactions in genetic association studies.

One in three people has been infected with Mycobacterium tuberculosis (MTB), and the risk for MTB infection in HIV-infected individuals is even higher. We hypothesized that HIV-positive individuals living in tuberculosis-endemic regions who do not get infected by Mycobacterium tuberculosis are genetically resistant. Using an “experiment of nature” design that proved successful in our previous work, we performed a genome-wide association study of tuberculin skin test positivity using 469 HIV-positive patients from prospective study cohorts of tuberculosis from Tanzania and Uganda to identify genetic loci associated with MTB infection in the context of HIV-infection. Among these individuals, 244 tested were tuberculin skin test (TST) positive either at enrollment or during the >8 year follow up, while 225 were not. We identified a genome-wide significant association between a dominant model of rs877356 and binary TST status in the combined cohort (Odds ratio = 0.2671, p = 1.22x10-8). Association was replicated with similar significance when examining TST induration as a continuous trait. The variant lies in the 5q31.1 region, 57kb downstream from IL9. Two-locus analyses of association of variants near rs877356 showed a haplotype comprised of rs877356 and an IL9 missense variant, rs2069885, had the most significant association (p = 1.59x10-12). We also replicated previously linked loci on chromosomes 2, 5, and 11. IL9 is a cytokine produced by mast cells and TH2 cells during inflammatory responses, providing a possible link between airway inflammation and protection from MTB infection. Our results indicate that studying uninfected, HIV-positive participants with extensive exposure increases the power to detect associations in complex infectious disease., Author summary Approximately one-third of the world’s population has been exposed to Mycobacterium tuberculosis, the bacterium that causes tuberculosis. A small number of those infected develop active disease; however, there is a substantial portion of exposed people who do not even show evidence of an immunological response. These people who appear to resist infection, as measured by a negative tuberculin skin test, represent a subpopulation from which we can learn about resistance. We used a genome-wide approach to study the genetic basis of this resistance in unique cohorts of hypervulnerable, HIV-positive individuals from Uganda and Tanzania, in which exposure was virtually assured. We identified one locus that was highly significantly associated and conferred more than 70% protection from infection. The most significant variant, rs8773656, was near IL9 and SLC25A48, and a haplotype including this variant and a missense mutation in IL9 was even more significantly associated with negative skin tests. Although it is impossible based solely on our data to determine the causal variant or genes, IL9 is an attractive candidate as its product has previously been associated with bronchial hyperresponsiveness, thereby providing a possible link between inflammation and protection from Mycobacterium tuberculosis infection.