Pooled Sequencing of Candidate Genes Implicates Rare Variants in the Development of Asthma Following Severe RSV Bronchiolitis in Infancy

Severe infection with respiratory syncytial virus (RSV) during infancy is strongly associated with the development of asthma. To identify genetic variation that contributes to asthma following severe RSV bronchiolitis during infancy, we sequenced the coding exons of 131 asthma candidate genes in 182 European and African American children with severe RSV bronchiolitis in infancy using anonymous pools for variant discovery, and then directly genotyped a set of 190 nonsynonymous variants. Association testing was performed for physician-diagnosed asthma before the 7^th birthday (asthma) using genotypes from 6,500 individuals from the Exome Sequencing Project (ESP) as controls to gain statistical power. In addition, among patients with severe RSV bronchiolitis during infancy, we examined genetic associations with asthma, active asthma, persistent wheeze, and bronchial hyperreactivity (methacholine PC₂₀) at age 6 years. We identified four rare nonsynonymous variants that were significantly associated with asthma following severe RSV bronchiolitis, including single variants in ADRB2, FLG and NCAM1 in European Americans (p = 4.6x10^-4, 1.9x10^-13 and 5.0x10^-5, respectively), and NOS1 in African Americans (p = 2.3x10^-11). One of the variants was a highly functional nonsynonymous variant in ADRB2 (rs1800888), which was also nominally associated with asthma (p = 0.027) and active asthma (p = 0.013) among European Americans with severe RSV bronchiolitis without including the ESP. Our results suggest that rare nonsynonymous variants contribute to the development of asthma following severe RSV bronchiolitis in infancy, notably in ADRB2. Additional studies are required to explore the role of rare variants in the etiology of asthma and asthma-related traits following severe RSV bronchiolitis.     

Nucleotide variation,haplotype structure,and association with end-stage renal disease of the human interleukin-1 gene cluster

A dense gene-based SNP map was constructed across a 360-kb region containing the interleukin-1 gene cluster (IL1A, IL1B, and IL1RN), focusing on IL1RN. In total, 95 polymorphisms were confirmed or identified primarily by direct sequencing. Polymorphisms were precisely mapped to completed BAC and genomic sequences spanning this region. The polymorphisms were typed in 443 case-control subjects from Caucasian and African American groups. Consecutive pair-wise marker linkage disequilibrium was not strictly correlated with distance and ranged from D'=0.0079 to 1.000 and D'=0.0521 to 1.0000 in Caucasians and African Americans, respectively. Single markers and haplotypes in IL1 cluster genes were evaluated for association with end-stage renal disease (ESRD). Eleven SNPs show some evidence of association with ESRD, with the strongest associations in two IL1A variants, one SNP, rs1516792-3, in intron 5 (p=0.0015) and a 4-bp insertion/deletion within the 3'UTR, rs16347-2 (p=0.0024), among African Americans with non-T2DM-associated ESRD.     

Influence of interleukin-12 receptor β1 polymorphisms on tuberculosis

Host genetic factors may be important determinants of susceptibility to tuberculosis, and several candidate gene polymorphisms have been shown to date. A series of recent reports concerning rare human deficiencies in the type-1 cytokine pathway suggest that more subtle variants of relevant genes may also contribute to susceptibility to tuberculosis at the general population level. To investigate whether polymorphisms in the interleukin-12 receptor (IL-12R) gene predispose individuals to tuberculosis, we studied these genes by single-strand conformational polymorphism analysis and direct sequencing. Although no common polymorphisms could be identified in the IL-12R beta 2 gene ( IL-12RB2), we confirmed four single nucleotide polymorphisms (SNPs; 641A-->G, 684C-->T, 1094T-->C, and 1132G-->C) causing three missense variants (Q214R, M365T, G378R) and one synonymous substitution in the extracellular domain of the IL-12R beta 1 gene ( IL12RB1). All SNPs were in almost perfect linkage disequilibrium (D'=0.98), and two common haplotypes of IL12RB1(allele 1: Q214-M365-G378; allele 2: R214-T365-R378) were revealed. Polymerase chain reaction/restriction fragment length polymorphism and sequence analyses were used to type IL12RB1polymorphisms in 98 patients with tuberculosis and 197 healthy controls in Japanese populations. In our case-control association study of tuberculosis, the R214-T365-R378 allele (allele 2) was over-represented in patients with tuberculosis, and homozygosity for R214-T365-R378 (the 2/2 genotype) was significantly associated with tuberculosis (odds ratio: 2.45; 95% CI: 1.20-4.99; P=0.013). In healthy subjects, homozygotes for R214-T365-R378 had lower levels of IL-12-induced signaling, according to differences in cellular responses to IL-12 between two haplotypes. These data suggest that the R214-T365-R378 allele, i.e., variation in IL12RB1, contribute to tuberculosis susceptibility in the Japanese population. This genetic variation may predispose individuals to tuberculosis infection by diminishing receptor responsiveness to IL-12 and to IL-23, leading to partial dysfunction of interferon-gamma-mediated immunity.     

Differences in candidate gene association between European ancestry and African American asthmatic children

Background

Candidate gene case-control studies have identified several single nucleotide polymorphisms (SNPs) that are associated with asthma susceptibility. Most of these studies have been restricted to evaluations of specific SNPs within a single gene and within populations from European ancestry. Recently, there is increasing interest in understanding racial differences in genetic risk associated with childhood asthma. Our aim was to compare association patterns of asthma candidate genes between children of European and African ancestry.

Methodology/Principal Findings

Using a custom-designed Illumina SNP array, we genotyped 1,485 children within the Greater Cincinnati Pediatric Clinic Repository and Cincinnati Genomic Control Cohort for 259 SNPs in 28 genes and evaluated their associations with asthma. We identified 14 SNPs located in 6 genes that were significantly associated (p-values <0.05) with childhood asthma in African Americans. Among Caucasians, 13 SNPs in 5 genes were associated with childhood asthma. Two SNPs in IL4 were associated with asthma in both races (p-values <0.05). Gene-gene interaction studies identified race specific sets of genes that best discriminate between asthmatic children and non-allergic controls.

Conclusions/Significance

We identified IL4 as having a role in asthma susceptibility in both African American and Caucasian children. However, while IL4 SNPs were associated with asthma in asthmatic children with European and African ancestry, the relative contributions of the most replicated asthma-associated SNPs varied by ancestry. These data provides valuable insights into the pathways that may predispose to asthma in individuals with European vs. African ancestry.     

Missense mutations of the interleukin-12 receptor beta 1(IL12RB1) and interferon-gamma receptor 1 (IFNGR1) genes are not associated with susceptibility to lepromatous leprosy in Korea

Interleukin-12 receptor beta 1 ( IL12RB1), interleukin-12 receptor beta 2 ( IL12RB2), and interferon gamma receptor 1 ( IFNGR1) perform important roles in the host defense against intracellular pathogens such as Mycobacteria. Several mutations within their genes have been confirmed as associated with increased susceptibility to mycobacterial infection. However, the association between mutations of the IL12RB1, IL12RB2, and IFNGR1 encoding genes and lepromatous leprosy has not been studied. This study screened for polymorphisms within IL12RB1, IL12RB2, and IFNGR1 encoding genes in the Korean populations using polymerase chain reaction (PCR)/single-strand conformation polymorphism (SSCP) DNA sequencing assay, and an association study was performed using the missense mutations of 705 A/G (Q214R), 1196 G/C (G378R), 1637 G/A (A525T), and 1664 C/T (P534S) of the IL12RB1, 83 G/A (V14M), and 1443 T/C (L467P) for the IFNGR1 encoding genes. There were no differences in the genotype and allele frequencies of IL12RB1 and IFNGR1 genes between 93 lepromatous leprosy patients and 94 control subjects. In conclusion, missense mutations of 705 A/G (Q214R), 1196 G/C (G378R), 1637 G/A (A525T), 1664 C/T (P534S) of the IL12RB1, 83 G/A (V14 M), and 1443 T/C (L467P) of the IFNGR1 encoding genes have no association with the susceptibility to lepromatous leprosy in the Korean population.     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

Genome-wide local ancestry approach identifies genes and variants associated with chemotherapeutic susceptibility in African Americans

Chemotherapeutic agents are used in the treatment of many cancers, yet variable resistance and toxicities among individuals limit successful outcomes. Several studies have indicated outcome differences associated with ancestry among patients with various cancer types. Using both traditional SNP-based and newly developed gene-based genome-wide approaches, we investigated the genetics of chemotherapeutic susceptibility in lymphoblastoid cell lines derived from 83 African Americans, a population for which there is a disparity in the number of genome-wide studies performed. To account for population structure in this admixed population, we incorporated local ancestry information into our association model. We tested over 2 million SNPs and identified 325, 176, 240, and 190 SNPs that were suggestively associated with cytarabine-, 5'-deoxyfluorouridine (5'-DFUR)-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-4)). Importantly, some of these variants are found only in populations of African descent. We also show that cisplatin-susceptibility SNPs are enriched for carboplatin-susceptibility SNPs. Using a gene-based genome-wide association approach, we identified 26, 11, 20, and 41 suggestive candidate genes for association with cytarabine-, 5'-DFUR-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-3)). Fourteen of these genes showed evidence of association with their respective chemotherapeutic phenotypes in the Yoruba from Ibadan, Nigeria (p<0.05), including TP53I11, COPS5 and GAS8, which are known to be involved in tumorigenesis. Although our results require further study, we have identified variants and genes associated with chemotherapeutic susceptibility in African Americans by using an approach that incorporates local ancestry information.     

Sequence variants in the genes for the interleukin-23 receptor (IL23R) and its ligand (IL12B) confer protection against psoriasis

Psoriasis is an inflammatory skin disorder that is inherited as a multifactorial trait. Genetic analyses have repeatedly identified a primary disease susceptibility locus lying within the major histocompatibility complex (MHC), on chromosome 6p21. A small number of non-MHC susceptibility loci have also been identified. These regions tend to overlap with susceptibility intervals for Crohn’s disease and atopic dermatitis, suggesting the possibility that genetic variants affecting inflammatory pathways may contribute to the pathogenesis of multiple disorders. Here, we report a genetic analysis of the interleukin 23 receptor gene (IL23R), which was recently identified as a susceptibility determinant for Crohn’s disease. We initially examined the results of a whole-genome association scan, carried out on 318 cases and 288 controls. We observed a significant increase of a non-synonymous substitution (p.Arg381Gln) among controls (P = 0.00036). We validated this finding by extending our cohort to include a further 519 cases and 528 controls. In the overall sample, the frequency of the 381Gln allele was 3.6% in cases and 7% in controls, yielding a P value of 0.00014. Next, we examined genetic variation at the IL12RB1, IL23A and IL12B genes, respectively, encoding the second subunit of the IL23R receptor and the two subunits of its ligand. This analysis identified independent associations for IL12B SNPs rs10045431 (P value for the extended dataset = 0.0001) and rs3212227 (P = 0.036). Altogether, these findings indicate that genes participating in IL23 signalling play a significant role in the pathogenesis of chronic epithelial inflammation. Francesca Capon and Paola Di Meglio have contributed equally to this work.     

Racial differences in selected cytokine allelic and genotypic frequencies among healthy,pregnant women in North Carolina

BACKGROUND: Genetic susceptibility to diseases is likely influenced by common DNA variants in the form of single nucleotide polymorphisms (SNPs). The value of SNPs for linkage and association mapping studies may depend on the distribution of SNP allele frequencies across populations. OBJECTIVES: To establish the SNP allelic frequencies among Caucasian and African American women for tumor necrosis factor (TNF)alpha, transforming growth factor (TGF)beta1, interleukin-10 (IL10), interleukin-6 (IL6), and interferon (IFN)gamma. MATERIALS AND METHODS: DNA was extracted from whole blood from 123 healthy, pregnant women. PCR-based genotyping was performed for the genes encoding TNFalpha (-308G/A), TGFbeta1 (codon 10C/T, codon 25C/G), IL10 (-1082A/G, -819T/C, -592A/C), IL6 (-174C/G) and IFNgamma (874T/A). Allele frequencies were determined by Hardy-Weinberg Equilibrium and Linkage Disequilibrium tests. Differences in the SNP allelic frequencies between Caucasians and African Americans were assessed by the chi(2) of Amitage trend test. RESULTS: SNP allelic and genotypic frequencies for IL6 and IFNgamma, but not for TNFalpha, TGFbeta1, and IL10, differed significantly between the Caucasian and African American women. CONCLUSIONS: Recognition of racial differences in SNP allelic and genotypic frequencies for selected cytokines is important for designing and powering future linkage and association mapping studies investigating the role of cytokines in human disease.     

IL12RB2 Gene Is Associated with the Age of Type 1 Diabetes Onset in Croatian Family Trios

Background

Common complex diseases are influenced by both genetic and environmental factors. Many genetic factors overlap between various autoimmune diseases. The aim of the present study is to determine whether four genetic variants known to be risk variants for several autoimmune diseases could be associated with an increased susceptibility to type 1 diabetes mellitus.

Methods and Findings

We genotyped four genetic variants (rs2358817, rs1049550, rs6679356, rs9865818) within VTCN1, ANXA11, IL12RB2 and LPP genes respectively, in 265 T1DM family trios in Croatian population. We did not detect association of these polymorphisms with T1DM. However, quantitative transmission disequilibrium test (QTDT, orthogonal model) revealed a significant association between the age of onset of T1DM and IL12RB2 rs6679356 variant. An earlier onset of T1DM was associated with the rs6679356 minor dominant allele C (p = 0.005). The association remained significant even after the Bonferroni correction for multiple testing and permutation.

Conclusions

Variants originally associated with juvenile idiopathic arthritis (VTCN1 gene), sarcoidosis (ANXA11 gene), primary biliary cirrhosis (IL12RB2 gene) and celiac disease (LPP gene) were not associated with type 1 diabetes in our dataset. Nevertheless, association of IL12RB2 rs6679356 polymorphism with the age of T1DM onset suggests that this gene plays a role in defining the time of disease onset.     

Unifying candidate gene and GWAS Approaches in Asthma

The first genome wide association study (GWAS) for childhood asthma identified a novel major susceptibility locus on chromosome 17q21 harboring the ORMDL3 gene, but the role of previous asthma candidate genes was not specifically analyzed in this GWAS. We systematically identified 89 SNPs in 14 candidate genes previously associated with asthma in >3 independent study populations. We re-genotyped 39 SNPs in these genes not covered by GWAS performed in 703 asthmatics and 658 reference children. Genotyping data were compared to imputation data derived from Illumina HumanHap300 chip genotyping. Results were combined to analyze 566 SNPs covering all 14 candidate gene loci. Genotyped polymorphisms in ADAM33, GSTP1 and VDR showed effects with p-values <0.0035 (corrected for multiple testing). Combining genotyping and imputation, polymorphisms in DPP10, EDN1, IL12B, IL13, IL4, IL4R and TNF showed associations at a significance level between p = 0.05 and p = 0.0035. These data indicate that (a) GWAS coverage is insufficient for many asthma candidate genes, (b) imputation based on these data is reliable but incomplete, and (c) SNPs in three previously identified asthma candidate genes replicate in our GWAS population with significance after correction for multiple testing in 14 genes.     

Deep resequencing identifies candidate functional genes in leprosy GWAS loci

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.     

Whole-exome sequencing reveals insights into genetic susceptibility to Congenital Zika Syndrome

Congenital Zika Syndrome (CZS) is a critical illness with a wide range of severity caused by Zika virus (ZIKV) infection during pregnancy. Life-threatening neurodevelopmental dysfunctions are among the most common phenotypes observed in affected newborns. Risk factors that contribute to susceptibility and response to ZIKV infection may be related to the virus itself, the environment, and maternal genetic background. Nevertheless, the newborn’s genetic contribution to the critical illness is still not elucidated. Here, we aimed to identify possible genetic variants as well as relevant biological pathways that might be associated with CZS phenotypes. For this purpose, we performed a whole-exome sequencing in 40 children born to women with confirmed exposure to ZIKV during pregnancy. We investigated the occurrence of rare harmful single-nucleotide variants (SNVs) possibly associated with inborn errors in genes ontologically related to CZS phenotypes. Moreover, an exome-wide association analysis was also performed using a case-control design (29 CZS cases and 11 controls), for both common and rare variants. Five out of the 29 CZS patients harbored known pathogenic variants likely to contribute to mild to severe manifestations observed. Approximately, 30% of affected individuals carried at least one pathogenic or likely pathogenic SNV in genes candidates to play a role in CZS. Our common variant association analysis detected a suggestive protective effect of the rs2076469 in DISP3 gene (p-value: 1.39 x 10⁻⁵). The IL12RB2 gene (p-value: 2.18x10^-11) also showed an unusual distribution of nonsynonymous rare SNVs in control samples. Finally, genes harboring harmful variants are involved in processes related to CZS phenotypes such as neurological development and immunity. Therefore, both rare and common variations may be likely to contribute as the underlying genetic cause of CZS susceptibility. The variations and pathways identified in this study may also have implications for the development of therapeutic strategies in the future.     

Variation in conserved non-coding sequences on chromosome 5q and susceptibility to asthma and atopy

Background

Evolutionarily conserved sequences likely have biological function.

Methods

To determine whether variation in conserved sequences in non-coding DNA contributes to risk for human disease, we studied six conserved non-coding elements in the Th2 cytokine cluster on human chromosome 5q31 in a large Hutterite pedigree and in samples of outbred European American and African American asthma cases and controls.

Results

Among six conserved non-coding elements (>100 bp, >70% identity; human-mouse comparison), we identified one single nucleotide polymorphism (SNP) in each of two conserved elements and six SNPs in the flanking regions of three conserved elements. We genotyped our samples for four of these SNPs and an additional three SNPs each in the IL13 and IL4 genes. While there was only modest evidence for association with single SNPs in the Hutterite and European American samples (P < 0.05), there were highly significant associations in European Americans between asthma and haplotypes comprised of SNPs in the IL4 gene (P < 0.001), including a SNP in a conserved non-coding element. Furthermore, variation in the IL13 gene was strongly associated with total IgE (P = 0.00022) and allergic sensitization to mold allergens (P = 0.00076) in the Hutterites, and more modestly associated with sensitization to molds in the European Americans and African Americans (P < 0.01).

Conclusion

These results indicate that there is overall little variation in the conserved non-coding elements on 5q31, but variation in IL4 and IL13, including possibly one SNP in a conserved element, influence asthma and atopic phenotypes in diverse populations.     

Identifying host genetic risk factors in the context of public health surveillance for invasive pneumococcal disease

Host genetic factors that modify risk of pneumococcal disease may help target future public health interventions to individuals at highest risk of disease. We linked data from population-based surveillance for invasive pneumococcal disease (IPD) with state-based newborn dried bloodspot repositories to identify biological samples from individuals who developed invasive pneumococcal disease. Genomic DNA was extracted from 366 case and 732 anonymous control samples. TagSNPs were selected in 34 candidate genes thought to be associated with host response to invasive pneumococcal disease, and a total of 326 variants were successfully genotyped. Among 543 European Americans (EA) (182 cases and 361 controls), and 166 African Americans (AA) (53 cases and 113 controls), common variants in surfactant protein D (SFTPD) are consistently underrepresented in IPD. SFTPD variants with the strongest association for IPD are intronic rs17886286 (allelic OR 0.45, 95% confidence interval (CI) [0.25, 0.82], with p = 0.007) in EA and 5′ flanking rs12219080 (allelic OR 0.32, 95%CI [0.13, 0.78], with p = 0.009) in AA. Variants in CD46 and IL1R1 are also associated with IPD in both EA and AA, but with effects in different directions; FAS, IL1B, IL4, IL10, IL12B, SFTPA1, SFTPB, and PTAFR variants are associated (p≤0.05) with IPD in EA or AA. We conclude that variants in SFTPD may protect against IPD in EA and AA and genetic variation in other host response pathways may also contribute to risk of IPD. While our associations are not corrected for multiple comparisons and therefore must be replicated in additional cohorts, this pilot study underscores the feasibility of integrating public health surveillance with existing, prospectively collected, newborn dried blood spot repositories to identify host genetic factors associated with infectious diseases.     

Effector CD4+ T cell expression signatures and immune-mediated disease associated genes

Genome-wide association studies (GWAS) in immune-mediated diseases have identified over 150 associated genomic loci. Many of these loci play a role in T cell responses, and regulation of T cell differentiation plays a critical role in immune-mediated diseases; however, the relationship between implicated disease loci and T cell differentiation is incompletely understood. To further address this relationship, we examined differential gene expression in na?ve human CD4+ T cells, as well as in in vitro differentiated Th1, memory Th17-negative and Th17-enriched CD4+ T cells subsets using microarray and RNASeq. We observed a marked enrichment for increased expression in memory CD4+ compared to na?ve CD4+ T cells of genes contained among immune-mediated disease loci. Within memory T cells, expression of disease-associated genes was typically increased in Th17-enriched compared to Th17-negative cells. Utilizing RNASeq and promoter methylation studies, we identified a differential regulation pattern for genes solely expressed in Th17 cells (IL17A and CCL20) compared to genes expressed in both Th17 and Th1 cells (IL23R and IL12RB2), where high levels of promoter methylation are correlated to near zero RNASeq levels for IL17A and CCL20. These findings have implications for human Th17 celI plasticity and for the regulation of Th17-Th1 expression signatures. Importantly, utilizing RNASeq we found an abundant isoform of IL23R terminating before the transmembrane domain that was enriched in Th17 cells. In addition to molecular resolution, we find that RNASeq provides significantly improved power to define differential gene expression and identify alternative gene variants relative to microarray analysis. The comprehensive integration of differential gene expression between cell subsets with disease-association signals, and functional pathways provides insight into disease pathogenesis.     

Genome-wide Association and Population Genetic Analysis of C-Reactive Protein in African American and Hispanic American Women

C-reactive protein (CRP) is a systemic inflammation marker that predicts future cardiovascular risk. CRP levels are higher in African Americans and Hispanic Americans than in European Americans, but the genetic determinants of CRP in these admixed United States minority populations are largely unknown. We performed genome-wide association studies (GWASs) of 8,280 African American (AA) and 3,548 Hispanic American (HA) postmenopausal women from the Women''s Health Initiative SNP Health Association Resource. We discovered and validated a CRP-associated variant of triggering receptors expressed by myeloid cells 2 (TREM2) in chromosomal region 6p21 (p = 10⁻¹⁰). The TREM2 variant associated with higher CRP is common in Africa but rare in other ancestral populations. In AA women, the CRP region in 1q23 contained a strong admixture association signal (p = 10⁻¹⁷), which appears to be related to several independent CRP-associated alleles; the strongest of these is present only in African ancestral populations and is associated with higher CRP. Of the other genomic loci previously associated with CRP through GWASs of European populations, most loci (LEPR, IL1RN, IL6R, GCKR, NLRP3, HNF1A, HNF4A, and APOC1) showed consistent patterns of association with CRP in AA and HA women. In summary, we have identified a common TREM2 variant associated with CRP in United States minority populations. The genetic architecture underlying the CRP phenotype in AA women is complex and involves genetic variants shared across populations, as well as variants specific to populations of African descent.