Egg allergy is one of the most common food allergies of childhood. There is a lack of information on the immunologic basis of egg allergy beyond the role of IgE.
To use transcriptional profiling as a novel approach to uncover immunologic processes associated with different phenotypes of egg allergy.
Peripheral blood mononuclear cells (PBMCs) were obtained from egg-allergic children who were defined as reactive (BER) or tolerant (BET) to baked egg, and from food allergic controls (AC) who were egg non-allergic. PBMCs were stimulated with egg white protein. Gene transcription was measured by microarray after 24 h, and cytokine secretion by multiplex assay after 5 days.
The transcriptional response of PBMCs to egg protein differed between BER and BET versus AC subjects. Compared to the AC group, the BER group displayed increased expression of genes associated with allergic inflammation as well as corresponding increased secretion of IL-5, IL-9 and TNF-α. A similar pattern was observed for the BET group. Further similarities in gene expression patterns between BER and BET groups, as well as some important differences, were revealed using a novel Immune Annotation resource developed for this project. This approach identified several novel processes not previously associated with egg allergy, including positive associations with TLR4-stimulated myeloid cells and activated NK cells, and negative associations with an induced Treg signature. Further pathway analysis of differentially expressed genes comparing BER to BET subjects showed significant enrichment of IFN-α and IFN-γ response genes, as well as genes associated with virally-infected DCs.
Transcriptional profiling identified several novel pathways and processes that differed when comparing the response to egg allergen in BET, BER, and AC groups. We conclude that this approach is a useful hypothesis-generating mechanism to identify novel immune processes associated with allergy and tolerance to forms of egg.
Motivation: Underrepresentation of racial groups represents an important challenge and major gap in phenomics research. Most of the current human phenomics research is based primarily on European populations; hence it is an important challenge to expand it to consider other population groups. One approach is to utilize data from EMR databases that contain patient data from diverse demographics and ancestries. The implications of this racial underrepresentation of data can be profound regarding effects on the healthcare delivery and actionability. To the best of our knowledge, our work is the first attempt to perform comparative, population-scale analyses of disease networks across three different populations, namely Caucasian (EA), African American (AA) and Hispanic/Latino (HL).
Results: We compared susceptibility profiles and temporal connectivity patterns for 1988 diseases and 37 282 disease pairs represented in a clinical population of 1 025 573 patients. Accordingly, we revealed appreciable differences in disease susceptibility, temporal patterns, network structure and underlying disease connections between EA, AA and HL populations. We found 2158 significantly comorbid diseases for the EA cohort, 3265 for AA and 672 for HL. We further outlined key disease pair associations unique to each population as well as categorical enrichments of these pairs. Finally, we identified 51 key ‘hub’ diseases that are the focal points in the race-centric networks and of particular clinical importance. Incorporating race-specific disease comorbidity patterns will produce a more accurate and complete picture of the disease landscape overall and could support more precise understanding of disease relationships and patient management towards improved clinical outcomes.
Contacts: firstname.lastname@example.org or email@example.com
Supplementary data are available at Bioinformatics online.
To investigate the association between African admixture and glaucoma prevalence amongst African American women.
Design, Setting, Participants
Participants included 11616 African American women from the Women’s Health Initiative Study (WHI) for whom admixture information was available and included 2548 who self-reported a diagnosis of glaucoma.
Main Outcome Measures
Significant association was observed between self-identified glaucoma status and admixture. However, this association was not significant in a model that included neighborhood socioeconomic status (NSES), hypertension, diabetes and body mass index (BMI). Self-identified glaucoma status was associated with diabetes that persisted after adjustment for admixture, NSES, hypertension, and BMI. Lower NSES was also associated with higher glaucoma risk but this association was marginal in the fully adjusted model and neither hypertension nor BMI showed association. When glaucoma status was limited to those reporting use or no use of appropriate ophthalmologic medication, no associations were observed in any of the models.
This study failed to find an independent association of glaucoma status and African admixture and these findings suggest that the higher frequency glaucoma in African Americans may be largely due to other factors.
Admixture; Glaucoma; African American women
Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10-21), A928V (rs35018800, OR = 0.53, P = 1.2x10-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, Pomnibus = 6x10-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; Pomnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.
Most studies suggest that hysterectomies are more common in African Americans than in other ethnic groups. To assess this ethnic surgical disparity in a novel way, our main goal was to determine whether admixture (the proportion of sub-Saharan African or European origin in individuals) is associated with hysterectomy frequency in African American women in the Women’s Health Initiative (WHI).
In this retrospective study, we used ancestry informative single nucleotide polymorphisms (SNPs) to estimate admixture proportions in >10,000 African American women from the WHI. Logistic regression models were used to assess the association between admixture and self-reported history of hysterectomy with and without controlling for relevant covariates. Multinomial logistic regression models were used to assess the association between admixture and self-reported age of hysterectomy. We also considered other potential risk factors (adiposity, hypertension, and education) for hysterectomy accounting for admixture.
African admixture was a strong risk factor after adjusting for multiple covariates (OR 1.85, P<.0001). The admixture risk for hysterectomy was highest for those performed in the 35–39 age range (OR 3.08, P<.0001) and least evident in oldest ages (45 or older). Our analyses also suggest that adiposity, hypertension and education were independently associated with hysterectomy in this population group.
These results suggest that higher African admixture is associated with higher frequencies of hysterectomy and that genetic studies specifically targeting African American women and diseases associated with hysterectomy may be especially useful in understanding the pathogenesis and underlying cause of this disparity in health outcome.
ancestry; ancestry informative markers; fibroid tumors; Women’s Health Initiative
Previous work has demonstrated that northern and southern European ancestries are associated with specific systemic lupus erythematosus (SLE) manifestations. Here, 1855 SLE cases of European descent were genotyped for 4965 single nucleotide polymorphisms and principal components analysis of genotype information was used to define population substructure. The first principal component (PC1) distinguished northern from southern European ancestry, PC2 differentiated eastern from western European ancestry, and PC3 delineated Ashkenazi Jewish ancestry. Compared to northern European ancestry, southern European ancestry was associated with autoantibody production (OR=1.40, 95% CI 1.07-1.83) and renal involvement (OR 1.41, 95% CI 1.06-1.87), and was protective for discoid rash (OR=0.51, 95% CI 0.32-0.82) and photosensitivity (OR=0.74, 95% CI 0.56-0.97). Both serositis (OR=1.46, 95% CI 1.12-1.89) and autoantibody production (OR=1.38, 95% CI 1.06-1.80) were associated with Western compared to Eastern European ancestry. Ashkenazi Jewish ancestry was protective against neurologic manifestations of SLE (OR=0.62, 95% CI 0.40-0.94). Homogeneous clusters of cases defined by multiple PCs demonstrated stronger phenotypic associations. Genetic ancestry may contribute to the development of SLE endophenotypes and should be accounted for in genetic studies of disease characteristics.
Systemic lupus erythematosus; epidemiology; population substructure; genetics
The objective of this study is to comprehensively define the genetic basis of Early Onset Myasthenia Gravis.
We have carried out a two-stage genome-wide association study on a total of 649 North European EOMG patients. Cases were matched 1:4 with controls of European ancestry. We performed imputation and conditional analyses across the major histocompatibility complex, as well as in the top regions of association outside the HLA region.
We observed the strongest association in the HLA class I region at rs7750641 (p = 1.2 × 10−92, OR = 6.25). By imputation and conditional analyses, HLA-B*08 proves to be the major associated allele (p = 2.87 × 10−113, OR = 6.41). In addition to the expected association with PTPN22 (rs2476601, OR =1.71, p = 8.2 ×10−10), an imputed coding variant (rs2233290) at position 151 (Pro→Ala) in the TNFAIP3-interacting protein 1, TNIP1, confers even stronger risk than PTPN22 (OR = 1.91, p = 3.2 × 10−10).
The association at TNIP1 in EOMG implies disease mechanisms involving ubiquitin-dependent dysregulation of NF-κB signaling. The localization of the major HLA signal to the HLA-B*08 allele suggests that CD8+ T-cells may play a key role in disease initiation or pathogenesis.
To further characterize the genetic basis of primary biliary cirrhosis (PBC), we genotyped 2426 PBC patients and 5731 unaffected controls from three independent cohorts using a single nucleotide polymorphism (SNP) array (Immunochip) enriched for autoimmune disease risk loci. Meta-analysis of the genotype data sets identified a novel disease-associated locus near the TNFSF11 gene at 13q14, provided evidence for association at six additional immune-related loci not previously implicated in PBC and confirmed associations at 19 of 22 established risk loci. Results of conditional analyses also provided evidence for multiple independent association signals at four risk loci, with haplotype analyses suggesting independent SNP effects at the 2q32 and 16p13 loci, but complex haplotype driven effects at the 3q25 and 6p21 loci. By imputing classical HLA alleles from this data set, four class II alleles independently contributing to the association signal from this region were identified. Imputation of genotypes at the non-HLA loci also provided additional associations, but none with stronger effects than the genotyped variants. An epistatic interaction between the IL12RB2 risk locus at 1p31and the IRF5 risk locus at 7q32 was also identified and suggests a complementary effect of these loci in predisposing to disease. These data expand the repertoire of genes with potential roles in PBC pathogenesis that need to be explored by follow-up biological studies.
Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with HLA region polymorphisms. To determine if associations can be explained by classical HLA determinants we studied Italian 676 cases and 1440 controls with genotyped with dense single nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (p = 1.59 × 10−11) was the strongest predisposing allele where as DRB1*11 (p = 1.42 × 10−10) was protective. Additionally DRB1*14 and the DPB1 association (DPB1*03:01) (p = 9.18 × 10−7) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.
genetic risk; risk allele; imputation; antigen binding pocket; autoimmune disease
We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10−10). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.
To investigate whether differences in admixture in African American (AFA) and Hispanic American (HA) adult women are associated with adiposity and adipose distribution.
The proportion of European, sub– Saharan African and Amerindian admixture was estimated for AFA and HA women in the Women's Heath Initiative using 92 ancestry informative markers. Analyses assessed the relationship between admixture and adiposity indices.
11712 AFA and 5088 HA self– identified post– menopausal women.
There was a significant positive association between body mass index (BMI) and African admixture when BMI was considered as a continuous variable, and age, education, physical activity, parity, family income and smoking were included covariates (p < 10− 4). A dichotomous model (upper and lower BMI quartiles) showed that African admixture was associated with a high odds ratio [OR = 3.27 (for 100% admixture compared to 0% admixture), 95% confidence interval (CI) 2.08 – 5.15]. For HA there was no association between BMI and admixture. In contrast, when waist to hip ratio (WHR) was used as a measure of adipose distribution, there was no significant association between WHR and admixture in AFA but there was a strong association in HA (p<10− 4; OR Amerindian admixture = 5.93, CI = 3.52 – 9.97).
These studies show that 1) African admixture is associated with BMI in AFA women; 2) Amerindian admixture is associated with WHR but not BMI in HA women; and 3) it may be important to consider different measurements of adiposity and adipose distribution in different ethnic population groups.
European admixture; African admixture; Amerindian admixture; post–menopausal women; waist to hip ratio; body mass index
A genome-wide association screen for primary biliary cirrhosis risk alleles was performed in an Italian cohort. The results from the Italian cohort replicated IL12A and IL12RB associations, and a combined meta-analysis using a Canadian dataset identified newly associated loci at SPIB (P = 7.9 × 10–11, odds ratio (OR) = 1.46), IRF5-TNPO3 (P = 2.8 × 10–10, OR = 1.63) and 17q12-21 (P = 1.7 × 10–10, OR = 1.38).
To provide a resource for assessing continental ancestry in a wide variety of genetic studies we identified, validated and characterized a set of 128 ancestry informative markers (AIMs). The markers were chosen for informativeness, genome-wide distribution, and genotype reproducibility on two platforms (TaqMan® assays and Illumina arrays). We analyzed genotyping data from 825 subjects with diverse ancestry, including European, East Asian, Amerindian, African, South Asian, Mexican, and Puerto Rican. A comprehensive set of 128 AIMs and subsets as small as 24 AIMs are shown to be useful tools for ascertaining the origin of subjects from particular continents, and to correct for population stratification in admixed population sample sets. Our findings provide general guidelines for the application of specific AIM subsets as a resource for wide application. We conclude that investigators can use TaqMan assays for the selected AIMs as a simple and cost efficient tool to control for differences in continental ancestry when conducting association studies in ethnically diverse populations.
population structure; continental ancestry; population stratification; ancestry informative markers
Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 × 10−8): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P ≤ 1 × 10−5. A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 × 10−3) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.
The definition of European population genetic substructure and its application to understanding complex phenotypes is becoming increasingly important. In the current study using over 4,000 subjects genotyped for 300,000 single-nucleotide polymorphisms (SNPs), we provide further insight into relationships among European population groups and identify sets of SNP ancestry informative markers (AIMs) for application in genetic studies. In general, the graphical description of these principal components analyses (PCA) of diverse European subjects showed a strong correspondence to the geographical relationships of specific countries or regions of origin. Clearer separation of different ethnic and regional populations was observed when northern and southern European groups were considered separately and the PCA results were influenced by the inclusion or exclusion of different self-identified population groups including Ashkenazi Jewish, Sardinian, and Orcadian ethnic groups. SNP AIM sets were identified that could distinguish the regional and ethnic population groups. Moreover, the studies demonstrated that most allele frequency differences between different European groups could be controlled effectively in analyses using these AIM sets. The European substructure AIMs should be widely applicable to ongoing studies to confirm and delineate specific disease susceptibility candidate regions without the necessity of performing additional genome-wide SNP studies in additional subject sets.
Case-control genetic studies of complex human diseases can be confounded by population stratification. This issue can be addressed using panels of ancestry informative markers (AIMs) that can provide substantial population substructure information. Previously, we described a panel of 128 SNP AIMs that were designed as a tool for ascertaining the origins of subjects from Europe, Sub-Saharan Africa, Americas, and East Asia.
In this study, genotypes from Human Genome Diversity Panel populations were used to further evaluate a 93 SNP AIM panel, a subset of the 128 AIMS set, for distinguishing continental origins. Using both model-based and relatively model-independent methods, we here confirm the ability of this AIM set to distinguish diverse population groups that were not previously evaluated. This study included multiple population groups from Oceana, South Asia, East Asia, Sub-Saharan Africa, North and South America, and Europe. In addition, the 93 AIM set provides population substructure information that can, for example, distinguish Arab and Ashkenazi from Northern European population groups and Pygmy from other Sub-Saharan African population groups.
These data provide additional support for using the 93 AIM set to efficiently identify continental subject groups for genetic studies, to identify study population outliers, and to control for admixture in association studies.
Accounting for population genetic substructure is important in reducing type 1 errors in genetic studies of complex disease. As efforts to understand complex genetic disease are expanded to different continental populations the understanding of genetic substructure within these continents will be useful in design and execution of association tests. In this study, population differentiation (Fst) and Principal Components Analyses (PCA) are examined using >200 K genotypes from multiple populations of East Asian ancestry. The population groups included those from the Human Genome Diversity Panel [Cambodian, Yi, Daur, Mongolian, Lahu, Dai, Hezhen, Miaozu, Naxi, Oroqen, She, Tu, Tujia, Naxi, Xibo, and Yakut], HapMap [ Han Chinese (CHB) and Japanese (JPT)], and East Asian or East Asian American subjects of Vietnamese, Korean, Filipino and Chinese ancestry. Paired Fst (Wei and Cockerham) showed close relationships between CHB and several large East Asian population groups (CHB/Korean, 0.0019; CHB/JPT, 00651; CHB/Vietnamese, 0.0065) with larger separation with Filipino (CHB/Filipino, 0.014). Low levels of differentiation were also observed between Dai and Vietnamese (0.0045) and between Vietnamese and Cambodian (0.0062). Similarly, small Fst's were observed among different presumed Han Chinese populations originating in different regions of mainland of China and Taiwan (Fst's <0.0025 with CHB). For PCA, the first two PC's showed a pattern of relationships that closely followed the geographic distribution of the different East Asian populations. PCA showed substructure both between different East Asian groups and within the Han Chinese population. These studies have also identified a subset of East Asian substructure ancestry informative markers (EASTASAIMS) that may be useful for future complex genetic disease association studies in reducing type 1 errors and in identifying homogeneous groups that may increase the power of such studies.
We cloned two genes, KIN1 and KIN2, encoding kinesin-II homologues from the ciliate Tetrahymena thermophila and constructed strains lacking either KIN1 or KIN2 or both genes. Cells with a single disruption of either gene showed partly overlapping sets of defects in cell growth, motility, ciliary assembly, and thermoresistance. Deletion of both genes resulted in loss of cilia and arrests in cytokinesis. Mutant cells were unable to assemble new cilia or to maintain preexisting cilia. Double knockout cells were not viable on a standard medium but could be grown on a modified medium on which growth does not depend on phagocytosis. Double knockout cells could be rescued by transformation with a gene encoding an epitope-tagged Kin1p. In growing cells, epitope-tagged Kin1p preferentially accumulated in cilia undergoing active assembly. Kin1p was also detected in the cell body but did not show any association with the cleavage furrow. The cell division arrests observed in kinesin-II knockout cells appear to be induced by the loss of cilia and resulting cell paralysis.