Methylome-wide Analysis Reveals Epigenetic Marks Associated With Resistance to Tuberculosis in Human Immunodeficiency Virus-Infected Individuals From East Africa.

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.

A chromosome 5q31.1 locus associates with tuberculin skin test reactivity in HIV-positive individuals from tuberculosis hyper-endemic regions in east Africa.

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.

A Locus at 5q33.3 Confers Resistance to Tuberculosis in Highly Susceptible Individuals.

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

Human and Helicobacter pylori coevolution shapes the risk of gastric disease.

Helicobacter pylori is the principal cause of gastric cancer, the second leading cause of cancer mortality worldwide. However, H. pylori prevalence generally does not predict cancer incidence. To determine whether coevolution between host and pathogen influences disease risk, we examined the association between the severity of gastric lesions and patterns of genomic variation in matched human and H. pylori samples. Patients were recruited from two geographically distinct Colombian populations with significantly different incidences of gastric cancer, but virtually identical prevalence of H. pylori infection. All H. pylori isolates contained the genetic signatures of multiple ancestries, with an ancestral African cluster predominating in a low-risk, coastal population and a European cluster in a high-risk, mountain population. The human ancestry of the biopsied individuals also varied with geography, with mostly African ancestry in the coastal region (58%), and mostly Amerindian ancestry in the mountain region (67%). The interaction between the host and pathogen ancestries completely accounted for the difference in the severity of gastric lesions in the two regions of Colombia. In particular, African H. pylori ancestry was relatively benign in humans of African ancestry but was deleterious in individuals with substantial Amerindian ancestry. Thus, coevolution likely modulated disease risk, and the disruption of coevolved human and H. pylori genomes can explain the high incidence of gastric disease in the mountain population.

Epigenetic and genetic variation in GATA5 is associated with gastric disease risk.

Gastric cancer incidence varies considerably among populations, even those with comparable rates of Helicobacter pylori infection. To test the hypothesis that genetic variation plays a role in gastric disease, we assessed the relationship between genotypes and gastric histopathology in a Colombian study population, using a genotyping array of immune-related single nucleotide polymorphisms (SNPs). Two synonymous SNPs (rs6061243 and rs6587239) were associated with progression of premalignant gastric lesions in a dominant-effects model after correction for multiple comparisons (p = 2.63E-07 and p = 7.97E-07, respectively); effect sizes were ß = -0.863 and ß = -0.815, respectively, where ß is an estimate of effect on histopathology scores, which ranged from 1 (normal) to 5 (dysplasia). In our replication cohort, a second Colombian population, both SNPs were associated with histopathology when additively modeled (ß = -0.256, 95 % CI = -0.47, -0.039; and ß = -0.239, 95 % CI = -0.45, -0.024), and rs6587239 was significantly associated in a dominant-effects model (ß = -0.330, 95 % CI = -0.66, 0.00). Because promoter methylation of GATA5 has previously been associated with gastric cancer, we also tested for the association of methylation status with more advanced histopathology scores in our samples and found a significant relationship (p = 0.001). A multivariate regression model revealed that the effects of both the promoter methylation and the exonic SNPs in GATA5 were independent. A SNP-by-methylation interaction term was also significant. This interaction between GATA5 variants and GATA5 promoter methylation indicates that the association of either factor with gastric disease progression is modified by the other.

Expression of complement and toll-like receptor pathway genes is associated with malaria severity in Mali: a pilot case control study.

BACKGROUND: The host response to infection by Plasmodium falciparum, the parasite most often responsible for severe malaria, ranges from asymptomatic parasitaemia to death. The clinical trajectory of malaria is influenced by host genetics and parasite load, but the factors determining why some infections produce uncomplicated malaria and some proceed to severe disease remain incompletely understood. METHODS: To identify molecular markers of severe falciparum malaria, human gene expression patterns were compared between children aged 6 months to 5 years with severe and uncomplicated malaria who were enrolled in a case-control study in Bandiagara, Mali. Microarrays were used to obtain expression data on severe cases and uncomplicated controls at the time of acute disease presentation (five uncomplicated and five severe), 1 week after presentation (three uncomplicated and three severe) and treatment initiation, and in the subsequent dry season (late convalescence, four uncomplicated and four severe). This is a pilot study for the first use of microarray technology in Mali. RESULTS: Complement and toll-like receptor (TLR) pathways were differentially expressed, with severe cases showing higher expression of the C1q, TLR2, TLR4, TLR8, and CR1 genes. Other genes previously associated with malaria pathogenesis, GZMB, FOS and HSPA6, were also higher among severe cases. TLR2, TLR4, TLR8, CR1, GZMB, FOS, and HSPA6 genes were expressed at lower levels in severe cases at late convalescence. CONCLUSIONS: Overexpression of genes previously associated with uncomplicated malaria was associated with severe disease. Low baseline expression of these genes may represent candidate markers for severe malaria. Despite the small sample size, results of this pilot study offer promising targets for follow-up analyses.

Addressing population-specific multiple testing burdens in genetic association studies.

The number of effectively independent tests performed in genome-wide association studies (GWAS) varies by population, making a universal P-value threshold inappropriate. We estimated the number of independent SNPs in Phase 3 HapMap samples by: (1) the LD-pruning function in PLINK, and (2) an autocorrelation-based approach. Autocorrelation was also used to estimate the number of independent SNPs in whole genome sequences from 1000 Genomes. Both approaches yielded consistent estimates of numbers of independent SNPs, which were used to calculate new population-specific thresholds for genome-wide significance. African populations had the most stringent thresholds (1.49 × 10(-7) for YRI at r(2) = 0.3), East Asian populations the least (3.75 × 10(-7) for JPT at r(2) = 0.3). We also assessed how using population-specific significance thresholds compared to using a single multiple testing threshold at the conventional 5 × 10(-8) cutoff. Applied to a previously published GWAS of melanoma in Caucasians, our approach identified two additional genes, both previously associated with the phenotype. In a Chinese breast cancer GWAS, our approach identified 48 additional genes, 19 of which were in or near genes previously associated with the phenotype. We conclude that the conventional genome-wide significance threshold generates an excess of Type 2 errors, particularly in GWAS performed on more recently founded populations.

Successful MDR-TB treatment regimens including amikacin are associated with high rates of hearing loss.

BACKGROUND: Aminoglycosides are a critical component of multidrug-resistant tuberculosis (MDR-TB) treatment but data on their efficacy and adverse effects in Botswana is scarce. We determined the effect of amikacin on treatment outcomes and development of hearing loss in MDR-TB patients. METHODS: Patients started on MDR-TB treatment between 2006 and 2012 were included. Multivariate analysis was used to determine the effect of amikacin on treatment outcomes and development of hearing loss. RESULTS: 437 MDR-TB patients were included, 288 (66%) of whom were HIV co-infected. 270 (62%) developed hearing loss, of whom 147 (54%) had audiometry. Of the 313 (72%) patients who completed treatment, 228 (73%) had a good outcome (cure or treatment completion). Good outcome was associated with longer amikacin treatment (adjusted OR [aOR] 1.13, 95% CI 1.06 - 1.21) and higher dosage (aOR 1.90, 95% CI 1.12 - 2.99). Longer amikacin duration (aOR 1.98, 95% CI 1.86 - 2.12) and higher dosage per weight per month (aOR 1.15, 95% CI 1.04 - 1.28) were associated with development of hearing loss. Amikacin treatment duration modified the effect of the dosage on the risk of hearing loss, increasing this risk as the duration increased. CONCLUSIONS: Amikacin was effective for MDR-TB treatment, but was associated with a high incidence of hearing loss especially in our study population. Total treatment duration and average monthly amikacin dose were associated with improved outcomes; however these were also associated with development of hearing loss.

Serologic and Cytokine Profiles of Children with Concurrent Cerebral Malaria and Severe Malarial Anemia Are Distinct from Other Subtypes of Severe Malaria.

We used a protein microarray featuring Plasmodium falciparum field variants of a merozoite surface antigen to examine malaria exposure in Malian children with different severe malaria syndromes. Unlike children with cerebral malaria alone or severe malarial anemia alone, those with concurrent cerebral malaria and severe malarial anemia had serologic responses demonstrating a broader prior parasite exposure pattern than matched controls with uncomplicated disease. Comparison of levels of malaria-related cytokines revealed that children with the concurrent phenotype had elevated levels of interleukin (IL)-6, IL-8, and IL-10. Our results suggest that the pathophysiology of this severe subtype is unique and merits further investigation.

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.

Genetic Effects on the Correlation Structure of CVD Risk Factors: Exome-Wide Data From a Ghanaian Population.

Plasma concentration of plasminogen activator inhibitor-1 (PAI-1) is highly correlated with several cardiovascular disease (CVD) risk factors. It also plays a direct role in CVD, including myocardial infarction and stroke, by impeding the dissolution of thrombi in the blood. Insofar as PAI-1 links CVD's risk factors to its endpoints, genetic variants modulating the relationship between PAI-1 and risk factors may be of particular clinical and biological interest. The high heritability of PAI-1, which has not been explained by genetic association studies, may also, in large part, be due to this relationship with CVD risk factors. Using exome-wide data from 1,032 Ghanaian study participants, we tested for heterogeneity of correlation by genotype between PAI-1 and 4 CVD risk factors (body mass index, triglycerides, mean arterial pressure, and fasting glucose) under the hypothesis that loci involved in the relationship between PAI-1 and other risk factors will also modify their correlational structure. We found more significant heterogeneities of correlation by genotype than we found marginal effects, with no evidence of type I inflation. The most significant result among all univariate and multivariate tests performed in this study was the heterogeneity of correlation between PAI-1 and mean arterial pressure at rs10738554, near SLC24A2, a gene previously associated with high blood pressure in African Americans.

Genomics of human pulmonary tuberculosis: from genes to pathways.

PURPOSE OF REVIEW: 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. RECENT FINDINGS: 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. SUMMARY: 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.

Draft Genome Sequences of 13 Colombian Helicobacter pylori Strains Isolated from Pacific Coast and Andean Residents.

We present here the draft genomes of 13 Helicobacter pylori strains isolated from Colombian residents on the Pacific coast (n = 6) and in the Andes mountains (n = 7), locations that differ in gastric cancer risk. These 13 strains were obtained from individuals with diagnosed gastric lesions.

Sex-Specific Parental Effects on Offspring Lipid Levels.

BACKGROUND: Plasma lipid levels are highly heritable traits, but known genetic loci can only explain a small portion of their heritability. METHODS AND RESULTS: In this study, we analyzed the role of parental levels of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglycerides (TGs) in explaining the values of the corresponding traits in adult offspring. We also evaluated the contribution of nongenetic factors that influence lipid traits (age, body mass index, smoking, medications, and menopause) alone and in combination with variability at the genetic loci known to associate with TC, LDL-C, HDL-C, and TG levels. We performed comparisons among different sex-specific regression models in 416 families from the Framingham Heart Study and 304 from the SardiNIA cohort. Models including parental lipid levels explain significantly more of the trait variation than models without these measures, explaining up to ≈39% of the total trait variation. Of this variation, the parent-of-origin effect explains as much as ≈15% and it does so in a sex-specific way. This observation is not owing to shared environment, given that spouse-pair correlations were negligible (<1.5% explained variation in all cases) and is distinct from previous genetic and acquired factors that are known to influence serum lipid levels. CONCLUSIONS: These findings support the concept that unknown genetic and epigenetic contributors are responsible for most of the heritable component of the plasma lipid phenotype, and that, at present, the clinical utility of knowing age-matched parental lipid levels in assessing risk of dyslipidemia supersedes individual locus effects. Our results support the clinical utility of knowing parental lipid levels in assessing future risk of dyslipidemia.

Co-infection with HPV types from the same species provides natural cross-protection from progression to cervical cancer.

BACKGROUND: The worldwide administration of bivalent and quadrivalent HPV vaccines has resulted in cross-protection against non-vaccine HPV types. Infection with multiple HPV types may offer similar cross-protection in the natural setting. We hypothesized that infections with two or more HPV types from the same species, and independently, infections with two or more HPV types from different species, associate with protection from high-grade lesions. FINDINGS: We recruited a cohort of 94 HIV, HPV-positive women from Botswana, with Grade 2 or higher cervical intraepithelial neoplasia. Infections with 2 or more HPV types from a single species associated with reduced lesion severity in univariate analysis (OR = 0.41, 95% CI 0.18-0.97, p = 0.042), when adjusted for the presence of HPV 16 or 18 types (OR = 0.41, 95% CI 0.17-1.00, p = 0.049), or all high-risk HPV type infections (OR = 0.38, 95% CI 0.16-0.90, p = 0.028). Infections with 2 or more HPV types from different species did not associate (OR = 0.68, 95% CI 0.25-1.81, p = 0.435). CONCLUSIONS: Our findings show that co-infections with genetically similar HPV types reduce the likelihood of progression to high-grade lesions in HIV positive women, an effect not observed in co-infections with taxonomically different HPV types. This observation is possibly caused by an immune cross-protection through a similar mechanism to that observed after HPV vaccination.

Disrupted human-pathogen co-evolution: a model for disease.

A major goal in infectious disease research is to identify the human and pathogenic genetic variants that explain differences in microbial pathogenesis. However, neither pathogenic strain nor human genetic variation in isolation has proven adequate to explain the heterogeneity of disease pathology. We suggest that disrupted co-evolution between a pathogen and its human host can explain variation in disease outcomes, and that genome-by-genome interactions should therefore be incorporated into genetic models of disease caused by infectious agents. Genetic epidemiological studies that fail to take both the pathogen and host into account can lead to false and misleading conclusions about disease etiology. We discuss our model in the context of three pathogens, Helicobacter pylori, Mycobacterium tuberculosis and human papillomavirus, and generalize the conditions under which it may be applicable.

Diverse convergent evidence in the genetic analysis of complex disease: coordinating omic, informatic, and experimental evidence to better identify and validate risk factors.

In omic research, such as genome wide association studies, researchers seek to repeat their results in other datasets to reduce false positive findings and thus provide evidence for the existence of true associations. Unfortunately this standard validation approach cannot completely eliminate false positive conclusions, and it can also mask many true associations that might otherwise advance our understanding of pathology. These issues beg the question: How can we increase the amount of knowledge gained from high throughput genetic data? To address this challenge, we present an approach that complements standard statistical validation methods by drawing attention to both potential false negative and false positive conclusions, as well as providing broad information for directing future research. The Diverse Convergent Evidence approach (DiCE) we propose integrates information from multiple sources (omics, informatics, and laboratory experiments) to estimate the strength of the available corroborating evidence supporting a given association. This process is designed to yield an evidence metric that has utility when etiologic heterogeneity, variable risk factor frequencies, and a variety of observational data imperfections might lead to false conclusions. We provide proof of principle examples in which DiCE identified strong evidence for associations that have established biological importance, when standard validation methods alone did not provide support. If used as an adjunct to standard validation methods this approach can leverage multiple distinct data types to improve genetic risk factor discovery/validation, promote effective science communication, and guide future research directions.

Salt sensitivity of blood pressure is associated with polymorphisms in the sodium-bicarbonate cotransporter.

Previous studies have demonstrated that single nucleotide polymorphisms (SNPs) of the sodium-bicarbonate co-transporter gene (SLC4A5) are associated with hypertension. We tested the hypothesis that SNPs in SLC4A5 are associated with salt sensitivity of blood pressure in 185 whites consuming an isocaloric constant diet with a randomized order of 7 days of low Na(+) (10 mmol/d) and 7 days of high Na(+) (300 mmol/d) intake. Salt sensitivity was defined as a ≥ 7-mm Hg increase in mean arterial pressure during a randomized transition between high and low Na(+) diet. A total of 35 polymorphisms in 17 candidate genes were assayed, 25 of which were tested for association. Association analyses with salt sensitivity revealed 3 variants that associated with salt sensitivity, 2 in SLC4A5 (P<0.001) and 1 in GRK4 (P=0.020). Of these, 2 SNPs in SLC4A5 (rs7571842 and rs10177833) demonstrated highly significant results and large effects sizes, using logistic regression. These 2 SNPs had P values of 1.0 × 10(-4) and 3.1 × 10(-4) with odds ratios of 0.221 and 0.221 in unadjusted regression models, respectively, with the G allele at both sites conferring protection. These SNPs remained significant after adjusting for body mass index and age (P=8.9 × 10(-5) and 2.6 × 10(-4) and odds ratios 0.210 and 0.286, respectively). Furthermore, the association of these SNPs with salt sensitivity was replicated in a second hypertensive population. Meta-analysis demonstrated significant associations of both SNPs with salt sensitivity (rs7571842 [P=1.2 × 10(-5)]; rs1017783 [P=1.1 × 10(-4)]). In conclusion, SLC4A5 variants are strongly associated with salt sensitivity of blood pressure in 2 separate white populations.