On the cross-population generalizability of gene expression prediction models.

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function.

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity.

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

NLRP1 variant M1184V decreases inflammasome activation in the context of DPP9 inhibition and asthma severity.

BACKGROUND: NLRP1 is an innate immune sensor that can form cytoplasmic inflammasome complexes. Polymorphisms in NLRP1 are linked to asthma; however, there is currently no functional or mechanistic explanation for this. OBJECTIVE: We sought to clarify the role of NLRP1 in asthma pathogenesis. METHODS: Results from the GALA II cohort study were used to identify a link between NLRP1 and asthma in Mexican Americans. In vitro and in vivo models for NLRP1 activation were applied to investigate the role of this inflammasome in asthma at the molecular level. RESULTS: We document the association of an NLRP1 haplotype with asthma for which the single nucleotide polymorphism rs11651270 (M1184V) individually is the most significant. Surprisingly, M1184V increases NLRP1 activation in the context of N-terminal destabilization, but decreases NLRP1 activation on dipeptidyl peptidase 9 inhibition. In vitro studies demonstrate that M1184V increases binding to dipeptidyl peptidase 9, which can account for its inhibitory role in this context. In addition, in vivo data from a mouse model of airway inflammation reveal a protective role for NLRP1 inflammasome activation reducing eosinophilia in this setting. CONCLUSIONS: Linking our in vitro and in vivo results, we found that the NLRP1 variant M1184V reduces inflammasome activation in the context of dipeptidyl peptidase 9 inhibition and could thereby increase asthma severity. Our studies may have implications for the treatment of asthma in patients carrying this variant of NLRP1.

Sequencing and imputation in GWAS: Cost-effective strategies to increase power and genomic coverage across diverse populations.

A key aim for current genome-wide association studies (GWAS) is to interrogate the full spectrum of genetic variation underlying human traits, including rare variants, across populations. Deep whole-genome sequencing is the gold standard to fully capture genetic variation, but remains prohibitively expensive for large sample sizes. Array genotyping interrogates a sparser set of variants, which can be used as a scaffold for genotype imputation to capture a wider set of variants. However, imputation quality depends crucially on reference panel size and genetic distance from the target population. Here, we consider sequencing a subset of GWAS participants and imputing the rest using a reference panel that includes both sequenced GWAS participants and an external reference panel. We investigate how imputation quality and GWAS power are affected by the number of participants sequenced for admixed populations (African and Latino Americans) and European population isolates (Sardinians and Finns), and identify powerful, cost-effective GWAS designs given current sequencing and array costs. For populations that are well-represented in existing reference panels, we find that array genotyping alone is cost-effective and well-powered to detect common- and rare-variant associations. For poorly represented populations, sequencing a subset of participants is often most cost-effective, and can substantially increase imputation quality and GWAS power.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Ambient air pollution, asthma drug response, and telomere length in African American youth.

BACKGROUND: Telomere length (TL) can serve as a potential biomarker for conditions associated with chronic oxidative stress and inflammation, such as asthma. Air pollution can induce oxidative stress. Understanding the relationship between TL, asthma, and air pollution is important for identifying risk factors contributing to unhealthy aging in children. OBJECTIVES: We sought to investigate associations between exposures to ambient air pollutants and TL in African American children and adolescents and to examine whether African ancestry, asthma status, and steroid medication use alter the association. METHODS: Linear regression was used to examine associations between absolute telomere length (aTL) and estimated annual average residential ozone (O3) and fine particulate matter with a diameter of 2.5 µm or less (PM2.5) exposures in a cross-sectional analysis of 1072 children in an existing asthma case-control study. African ancestry, asthma status, and use of steroid medications were examined as effect modifiers. RESULTS: Participants' aTLs were measured by using quantitative PCR. A 1-ppb and 1 µg/m3 increase in annual average exposure to O3 and PM2.5 were associated with a decrease in aTL of 37.1 kilo-base pair (kb; 95% CI, -66.7 to -7.4 kb) and 57.1 kb (95% CI, -118.1 to 3.9 kb), respectively. African ancestry and asthma were not effect modifiers; however, exposure to steroid medications modified the relationships between TL and pollutants. Past-year exposure to O3 and PM2.5 was associated with shorter TLs in patients without steroid use. CONCLUSION: Exposure to air pollution was associated with shorter TLs in nonasthmatic children and adolescents. This was not the case for asthmatic children as a group, but those receiving steroid medication had less shortening than those not using steroids. Reduced exposure to air pollution in childhood might help to preserve TL.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma.

Short-acting ß2-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the United States. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African-American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single-nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1830 Latinos (total = 2779). Finally, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p = 7.69 × 10-9). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958 and rs7081864, p ≤ 5 × 10-8). Our results failed to replicate in three additional populations of 416 Latinos and 1615 African Americans. Our findings indicate that both population-specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

Secondhand smoke exposure and asthma outcomes among African-American and Latino children with asthma.

BACKGROUND: Secondhand smoke (SHS) exposures have been linked to asthma-related outcomes but quantitative dose-responses using biomarkers of exposure have not been widely reported. OBJECTIVES: Assess dose-response relationships between plasma cotinine-determined SHS exposure and asthma outcomes in minority children, a vulnerable population exposed to higher levels of SHS and under-represented in the literature. METHODS: We performed analyses in 1172 Latino and African-American children with asthma from the mainland USA and Puerto Rico. We used logistic regression to assess relationships of cotinine levels ≥0.05 ng/mL with asthma exacerbations (defined as asthma-related hospitalisations, emergency room visits or oral steroid prescription) in the previous year and asthma control. The shape of dose-response relationships was assessed using a continuous exposure variable in generalised additive logistic models with penalised splines. RESULTS: The OR for experiencing asthma exacerbations in the previous year for cotinine levels ≥0.05 ng/mL, compared with <0.05 ng/mL, was 1.40 (95% CI 1.03 to 1.89), while the OR for poor asthma control was 1.53 (95% CI 1.12 to 2.13). Analyses for dose-response relationships indicated increasing odds of asthma outcomes related with increasing exposure, even at cotinine levels associated with light SHS exposures. CONCLUSIONS: Exposure to SHS was associated with higher odds of asthma exacerbations and having poorly controlled asthma with an increasing dose-response even at low levels of exposure. Our results support the conclusion that there are no safe levels of SHS exposures.

Air Pollution and Lung Function in Minority Youth with Asthma in the GALA II (Genes-Environments and Admixture in Latino Americans) and SAGE II (Study of African Americans, Asthma, Genes, and Environments) Studies.

RATIONALE: Adverse effects of exposures to ambient air pollution on lung function are well documented, but evidence in racial/ethnic minority children is lacking. OBJECTIVES: To assess the relationship between air pollution and lung function in minority children with asthma and possible modification by global genetic ancestry. METHODS: The study population consisted of 1,449 Latino and 519 African American children with asthma from five different geographical regions in the mainland United States and Puerto Rico. We examined five pollutants (particulate matter ≤10 µm and ≤2.5 µm in diameter, ozone, nitrogen dioxide, and sulfur dioxide), derived from participant residential history and ambient air monitoring data, and assessed over several time windows. We fit generalized additive models for associations between pollutant exposures and lung function parameters and tested for interaction terms between exposures and genetic ancestry. MEASUREMENTS AND MAIN RESULTS: A 5 µg/m(3) increase in average lifetime particulate matter less than or equal to 2.5 µm in diameter exposure was associated with a 7.7% decrease in FEV1 (95% confidence interval = -11.8 to -3.5%) in the overall study population. Global genetic ancestry did not appear to significantly modify these associations, but percent African ancestry was a significant predictor of lung function. CONCLUSIONS: Early-life particulate exposures were associated with reduced lung function in Latino and African American children with asthma. This is the first study to report an association between exposure to particulates and reduced lung function in minority children in which racial/ethnic status was measured by ancestry-informative markers.

Novel genetic risk factors for asthma in African American children: Precision Medicine and the SAGE II Study.

Asthma, an inflammatory disorder of the airways, is the most common chronic disease of children worldwide. There are significant racial/ethnic disparities in asthma prevalence, morbidity, and mortality among US children. This trend is mirrored in obesity, which may share genetic and environmental risk factors with asthma. The majority of asthma biomedical research has been performed in populations of European decent. We sought to identify genetic risk factors for asthma in African American children. We also assessed the generalizability of genetic variants associated with asthma in European and Asian populations to African American children. Our study population consisted of 1227 (812 asthma cases, 415 controls) African American children with genome-wide single nucleotide polymorphism (SNP) data. Logistic regression was used to identify associations between SNP genotype and asthma status. We identified a novel variant in the PTCHD3 gene that is significantly associated with asthma (rs660498, p = 2.2 × 10(-7)) independent of obesity status. Approximately 5 % of previously reported asthma genetic associations identified in European populations replicated in African Americans. Our identification of novel variants associated with asthma in African American children, coupled with our inability to replicate the majority of findings reported in European Americans, underscores the necessity for including diverse populations in biomedical studies of asthma.

Diversity in Clinical and Biomedical Research: A Promise Yet to Be Fulfilled.

Esteban Gonzalez Burchard and colleagues explore how making medical research more diverse would aid not only social justice but scientific quality and clinical effectiveness, too.

A single nucleotide polymorphism in SLC7A5 is associated with gastrointestinal toxicity after high-dose melphalan and autologous stem cell transplantation for multiple myeloma.

Multiple myeloma is the most frequent indication for high-dose melphalan (HDM) chemotherapy with autologous stem cell transplantation (ASCT). Gastrointestinal symptoms represent the most significant nonhematological toxicity of HDM. However, specific, especially genetic, predictors of their incidence or clinical severity are lacking. The amino acid transporters LAT1 and LAT2 encoded by the SLC7A5 and SLC7A8 genes, respectively, are the principal mediators of melphalan uptake into cells. To determine whether genetic variability at these loci contributed to interindividual differences in the development of gastrointestinal complications of HDM, we analyzed single nucleotide polymorphisms (SNPs) in these genes in 135 patients with multiple myeloma treated with HDM and ASCT and correlated these with the need for total parenteral nutrition (TPN). Seven SNPs in SLC7A5 and 20 in SLC7A8 were genotyped. Multiple analyses indicated that 1 SNP in the first intron of SLC7A5, rs4240803, was significantly associated with TPN use (odds ratio = .45, 95% confidence interval, .25 to .79; P = .007). Further, every haplotype that correlated with TPN requirement included this SNP. These results suggest that variability in melphalan transport affects mucosal injury after HDM. This finding could help in individualizing the dose of this effective and widely used chemotherapeutic agent for multiple myeloma.

Increased variance in germline allele-specific expression of APC associates with colorectal cancer.

BACKGROUND & AIMS: Germline variations in allele-specific expression (ASE) are associated with highly penetrant familial cancers, but their role in common sporadic cancers is unclear. ASE of adenomatous polyposis coli (APC) is associated with pathogenesis of familial adenomatous polyposis. We investigated whether moderate variations in ASE of APC contribute to common forms of colorectal cancer (CRC). METHODS: Denaturing high-performance liquid chromatography was used to analyze germline ASE of APC in blood samples from patients with CRC (cases, n = 53) and controls (n = 68). Means, medians, and variances of ASE were compared. Variants in the APC gene region also were analyzed. RESULTS: The distribution of ASE differed significantly between groups; cases had significantly larger amounts of variance than controls (P = .0004). Risk for CRC increased proportionally with the degree of deviation from the mean. The odds ratio for individuals with levels of ASE that deviated more than 1 standard deviation from the mean was 3.97 (95% confidence interval, 1.71-9.24; P = .001); for those with levels greater than 1.645 standard deviations, the odds ratio was 13.46 (95% confidence interval, 1.76-609.40; P = .005). Sequence analysis revealed that a patient with a high level of ASE who did not have a family history of CRC carried a nonsense mutation in APC (p.Arg216X). Genotype analysis of APC associated multiple single-nucleotide polymorphisms with ASE values and/or variance among cases, but not controls. Cis variants, therefore, might account for some of the variance in ASE of APC. CONCLUSIONS: Patients with CRC have a larger variance in germline levels of ASE in APC than controls; large distances from the mean ASE were associated with risk for common forms of CRC.

Diagnostic tools for hypertension and salt sensitivity testing.

PURPOSE OF REVIEW: One-third of the world's population has hypertension and it is responsible for almost 50% of deaths from stroke or coronary heart disease. These statistics do not distinguish salt-sensitive from salt-resistant hypertension or include normotensives who are salt-sensitive even though salt sensitivity, independent of blood pressure, is a risk factor for cardiovascular and other diseases, including cancer. This review describes new personalized diagnostic tools for salt sensitivity. RECENT FINDINGS: The relationship between salt intake and cardiovascular risk is not linear, but rather fits a J-shaped curve relationship. Thus, a low-salt diet may not be beneficial to everyone and may paradoxically increase blood pressure in some individuals. Current surrogate markers of salt sensitivity are not adequately sensitive or specific. Tests in the urine that could be surrogate markers of salt sensitivity with a quick turn-around time include renal proximal tubule cells, exosomes, and microRNA shed in the urine. SUMMARY: Accurate testing of salt sensitivity is not only laborious but also expensive, and with low patient compliance. Patients who have normal blood pressure but are salt-sensitive cannot be diagnosed in an office setting and there are no laboratory tests for salt sensitivity. Urinary surrogate markers for salt sensitivity are being developed.

PTX3 genetic variation and dizygotic twinning in the Gambia: could pleiotropy with innate immunity explain common dizygotic twinning in Africa?

Dizygotic (DZ) twinning has a genetic component and is common among sub-Saharan Africans; in The Gambia its frequency is up to 3% of live births. Variation in PTX3, encoding Pentraxin 3, a soluble pattern recognition receptor that plays an important role both in innate immunity and in female fertility, has been associated with resistance to Mycobacterium tuberculosis pulmonary disease and to Pseudomonas aeruginosa infection in cystic fibrosis patients. We tested whether PTX3 variants in Gambian women associate with DZ twinning, by genotyping five PTX3 single nucleotide polymorphisms (SNPs) in 130 sister pairs (96 full sibs and 34 half sibs) who had DZ twins. Two, three and five SNP haplotypes differed in frequency between twinning mothers and those without a history of twinning (from P = 0.006 to 3.03e-06 for two SNP and three SNP haplotypes, respectively). Twinning mothers and West African tuberculosis-controls from a previous study shared several frequent haplotypes. Most importantly, our data are consistent with an independently reported association of PTX3 and female fertility in a sample from Ghana. Taken together, these results indicate that selective pressure on PTX3 variants that affect the innate immune response to infectious agents, could also produce the observed high incidence of DZ twinning in Gambians.

Native American Ancestry and Air Pollution Interact to Impact Bronchodilator Response in Puerto Rican Children with Asthma.

Objective: Asthma is the most common chronic disease in children. Short-acting bronchodilator medications are the most commonly prescribed asthma treatment worldwide, regardless of disease severity. Puerto Rican children display the highest asthma morbidity and mortality of any US population. Alarmingly, Puerto Rican children with asthma display poor bronchodilator drug response (BDR). Reduced BDR may explain, in part, the increased asthma morbidity and mortality observed in Puerto Rican children with asthma. Gene-environment interactions may explain a portion of the heritability of BDR. We aimed to identify gene-environment interactions associated with BDR in Puerto Rican children with asthma. Setting: Genetic, environmental, and psycho-social data from the Genes-environments and Admixture in Latino Americans (GALA II) case-control study. Participants: Our discovery dataset consisted of 658 Puerto Rican children with asthma; our replication dataset consisted of 514 Mexican American children with asthma. Main Outcome Measures: We assessed the association of pairwise interaction models with BDR using ViSEN (Visualization of Statistical Epistasis Networks). Results: We identified a non-linear interaction between Native American genetic ancestry and air pollution significantly associated with BDR in Puerto Rican children with asthma. This interaction was robust to adjustment for age and sex but was not significantly associated with BDR in our replication population. Conclusions: Decreased Native American ancestry coupled with increased air pollution exposure was associated with increased BDR in Puerto Rican children with asthma. Our study acknowledges BDR's phenotypic complexity, and emphasizes the importance of integrating social, environmental, and biological data to further our understanding of complex disease.

CLEC4E (Mincle) genetic variation associates with pulmonary tuberculosis in Guinea-Bissau (West Africa).

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 seven SNPs (p < 0.05). All SNPs were 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.