Rafal S. Sobota

Human and Helicobacter pylori coevolution shapes the risk of gastric disease

Significance Theory predicts that chronic pathogens with vertical or familial transmission should become less virulent over time because of coevolution. Although transmitted in this way, Helicobacter pylori is the major causative agent of gastric cancer. In two distinct Colombian populations with similar levels of H. pylori infection but different incidences of gastric cancer, we examined human and pathogen ancestry in matched samples to assess whether their genomic variation affects the severity of premalignant lesions. Interaction between human Amerindian ancestry and H. pylori African ancestry accounted for the geographic disparity in clinical presentation. We conclude that coevolutionary relationships are important determinants of gastric disease risk and that the historical colonization of the Americas continues to influence health in modern American populations. 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.

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.

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 r2 = 0.3), East Asian populations the least (3.75 × 10−7 for JPT at r2 = 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.

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.

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

BackgroundAminoglycosides 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.MethodsPatients 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.Results437 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.ConclusionsAmikacin 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.

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.

Cardiovascular Disease Risk Factors in Ghana during the Rural-to-Urban Transition: A Cross-Sectional Study

Populations in sub-Saharan Africa are shifting from rural to increasingly urban. Although the burden of cardiovascular disease is expected to increase with this changing landscape, few large studies have assessed a wide range of risk factors in urban and rural populations, particularly in West Africa. We conducted a cross-sectional, population-based survey of 3317 participants from Ghana (≥18 years old), of whom 2265 (57% female) were from a mid-sized city (Sunyani, population ~250,000) and 1052 (55% female) were from surrounding villages (populations <5000). We measured canonical cardiovascular disease risk factors (BMI, blood pressure, fasting glucose, lipids) and fibrinolytic markers (PAI-1 and t-PA), and assessed how their distributions and related clinical outcomes (including obesity, hypertension and diabetes) varied with urban residence and sex. Urban residence was strongly associated with obesity (OR: 7.8, 95% CI: 5.3–11.3), diabetes (OR 3.6, 95% CI: 2.3–5.7), and hypertension (OR 3.2, 95% CI: 2.6–4.0). Among the quantitative measures, most affected were total cholesterol (+0.81 standard deviations, 95% CI 0.73–0.88), LDL cholesterol (+0.89, 95% CI: 0.79–0.99), and t-PA (+0.56, 95% CI: 0.48–0.63). Triglycerides and HDL cholesterol profiles were similarly poor in both urban and rural environments, but significantly worse among rural participants after BMI-adjustment. For most of the risk factors, the strength of the association with urban residence did not vary with sex. Obesity was a major exception, with urban women at particularly high risk (26% age-standardized prevalence) compared to urban men (7%). Overall, urban residents had substantially worse cardiovascular risk profiles, with some risk factors at levels typically seen in the developed world.

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.

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.

Plasminogen Activator Inhibitor‐1 and Diagnosis of the Metabolic Syndrome in a West African Population

Background Metabolic syndrome (MetS) is diagnosed by the presence of at least 3 of the following: obesity, hypertension, hyperglycemia, hypertriglyceridemia, and low high‐density lipoprotein. Individuals with MetS also typically have elevated plasma levels of the antifibrinolytic factor, plasminogen activator inhibitor‐1 (PAI‐1), but the relationships between PAI‐1 and MetS diagnostic criteria are not clear. Understanding these relationships can elucidate the relevance of MetS to cardiovascular disease risk, because PAI‐1 is associated with ischemic events and directly involved in thrombosis. Methods and Results In a cross‐sectional analysis of 2220 Ghanaian men and women from urban and rural locales, we found the age‐standardized prevalence of MetS to be as high as 21.4% (urban women). PAI‐1 level increased exponentially as the number of diagnostic criteria increased linearly (P<10−13), supporting the conclusion that MetS components have a joint effect that is stronger than their additive contributions. Body mass index, triglycerides, and fasting glucose were more strongly correlated with PAI‐1 than with canonical MetS criteria, and this pattern did not change when pair‐wise correlations were conditioned on all other risk factors, supporting an independent role for PAI‐1 in MetS. Finally, whereas the correlations between conventional risk factors did not vary significantly by sex or across urban and rural environments, correlations with PAI‐1 were generally stronger among urban participants. Conclusions MetS prevalence in the West African population we studied was comparable to that of the industrialized West. PAI‐1 may serve as a key link between MetS, as currently defined, and the endpoints with which it is associated. Whether this association is generalizable will require follow‐up.