Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genome-Wide Association Study of Staphylococcus aureus Carriage in a Community-Based Sample of Mexican-Americans in Starr County, Texas

  • Eric L. Brown,

    Affiliation Center for Infectious Disease, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center, Houston, TX, United States of America

    ⨯
  • Jennifer E. Below,

    Affiliation Human Genetics Center, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center at Houston, Houston, TX, United States of America

    ⨯
  • Rebecca S. B. Fischer,

    Current address: Baylor College of Medicine, National School of Tropical Medicine, Section of Tropical Medicine, Houston, TX, United States of America

    Affiliation Center for Infectious Disease, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center, Houston, TX, United States of America

    ⨯
  • Heather T. Essigmann,

    Affiliation Center for Infectious Disease, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center, Houston, TX, United States of America

    ⨯
  • Hao Hu,

    Affiliation Department of Epidemiology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States of America

    ⨯
  • Chad Huff,

    Affiliation Department of Epidemiology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States of America

    ⨯
  • D. Ashley Robinson,

    Affiliation Department of Microbiology, University of Mississippi Medical Center, Jackson, MS, United States of America

    ⨯
  • Lauren E. Petty,

    Affiliation Human Genetics Center, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center at Houston, Houston, TX, United States of America

    ⨯
  • David Aguilar,

    Affiliation Winters Center for Heart Failure Research and Section of Cardiology, Baylor College of Medicine, Houston, TX, United States of America

    ⨯
  • Graeme I. Bell,

    Affiliation Departments of Medicine and Human Genetics, The University of Chicago, Chicago, United States of America

    ⨯
  • Craig L. Hanis

    craig.l.hanis@uth.tmc.edu

    Affiliation Human Genetics Center, Division of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center at Houston, Houston, TX, United States of America

    ⨯

Abstract

Staphylococcus aureus is the number one cause of hospital-acquired infections. Understanding host pathogen interactions is paramount to the development of more effective treatment and prevention strategies. Therefore, whole exome sequence and chip-based genotype data were used to conduct rare variant and genome-wide association analyses in a Mexican-American cohort from Starr County, Texas to identify genes and variants associated with S. aureus nasal carriage. Unlike most studies of S. aureus that are based on hospitalized populations, this study used a representative community sample. Two nasal swabs were collected from participants (n = 858) 11–17 days apart between October 2009 and December 2013, screened for the presence of S. aureus, and then classified as either persistent, intermittent, or non-carriers. The chip-based and exome sequence-based single variant association analyses identified 1 genome-wide significant region (KAT2B) for intermittent and 11 regions suggestively associated with persistent or intermittent S. aureus carriage. We also report top findings from gene-based burden analyses of rare functional variation. Notably, we observed marked differences between signals associated with persistent and intermittent carriage. In single variant analyses of persistent carriage, 7 of 9 genes in suggestively associated regions and all 5 top gene-based findings are associated with cell growth or tight junction integrity or are structural constituents of the cytoskeleton, suggesting that variation in genes associated with persistent carriage impact cellular integrity and morphology.

Citation: Brown EL, Below JE, Fischer RSB, Essigmann HT, Hu H, Huff C, et al. (2015) Genome-Wide Association Study of Staphylococcus aureus Carriage in a Community-Based Sample of Mexican-Americans in Starr County, Texas. PLoS ONE 10(11): e0142130. https://doi.org/10.1371/journal.pone.0142130

Editor: Patrick M. Schlievert, University of Iowa Carver College of Medicine, UNITED STATES

Received: August 4, 2015; Accepted: October 16, 2015; Published: November 16, 2015

Copyright: © 2015 Brown et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was supported by NIH grants R01 AI085014-01A1 (E.L. Brown and C.L. Hanis), R01 GM104390 (C.D. Huff and H. Hao), HL102830 (C.L. Hanis), DK085501 (C.L. Hanis) and P30DK020595 (G.I. Bell). H. Hao is also supported by the M.D. Anderson Cancer Center Odyssey Program. This work was also partially supported from a grant from the Kleberg Foundation to E.L. Brown. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Infectious diseases result from complex interactions between the microorganism, the host, and the environment. Host genetic factors play a major role in determining differential susceptibility to major infectious diseases of humans, including malaria [1], HIV/AIDS [2], tuberculosis [3], hepatitis B [4], Norovirus diarrhea [5], prion disease [6], Cholera [7], and Helicobacter pylori infections [8]. The first evidence that genetic factors could impact infectious disease outcomes was derived from epidemiological studies that identified differences between human populations exposed to the same infectious organism [9]. This is equally true for S. aureus [10–12], but this pathogen represents a special case because it is an opportunistic pathogen that can colonize humans without causing overt disease [13]. It is therefore an ideal system for examining host pathogen interactions.

Even though humans are exposed to S. aureus at birth, not all are equally susceptible to colonization [9]. Many body sites can be colonized by S. aureus, but nasal decolonization has been shown to be effective in reducing colonization at other body sites, suggesting that the anterior nares is one of the primary S. aureus reservoirs [14, 15]. Human carriage has been classified as either persistent, intermittent, or non-carriage with rates of carriage ranging from 10–35%, 20–75%, and 5–70%, respectively, depending on race, age, gender, and whether the population examined was hospital- or community-based [9, 16–18]. Carriage is not representative of infection, per se. Rather, carriage impacts the risk of acquiring infection, disease presentation, and disease severity [13]. Furthermore, the genotype of the colonizing S. aureus strain, the nature of the immune response elicited following exposure, and underlying host genetic factors may all play a role in susceptibility to colonization and/or infection [9, 19–24]. Like other complex conditions, susceptibility to infectious agents does not typically follow a simple Mendelian pattern of inheritance, largely due to the fact that human immune responses are controlled by complex genetic mechanisms and modifying environmental influences [25, 26].

Candidate gene studies have uncovered associations between specific genes and carriage status [20–23, 27–29]. For example, IL4 and C-reactive protein have been shown to be associated with carriage in the Rotterdam Study [20, 22]. In the same study, a 68% reduction in risk of persistent carriage was observed related to the glucocorticoid receptor gene [30] (S1 Table). Polymorphisms in genes encoding different defensins and MBL (manose binding lectin) have also been associated with S. aureus persistent carriage [20, 31, 32] (S1 Table). The toll-like receptors have also been associated with increased risk of streptococci and enterococci skin and soft tissue infections [21, 33] suggesting that there may be some commonalities in the genetics of susceptibility to infection with different pathogens. No community-based genome-wide association or whole exome sequencing studies have previously been performed in the context of S. aureus carriage, but recently, 2 hospital-based genome-wide association studies of S. aureus infections were conducted [34, 35]. That these studies failed to identify targets with genome wide significance is not necessarily surprising since hospital environments themselves are a significant risk factor for acquiring S. aureus infections and these effects may overwhelm modest genetic influences on risk [36].

The present study was designed to identify genes/markers associated with persistent and intermittent carriage of S. aureus in a community-based sample of 858 Mexican-Americans from Starr County, Texas. Single nucleotide polymorphism (SNP) data from the Affymetrix Genome-Wide SNPArray 6.0 assay imputed out to the complete SNP set in the 1000 Genomes Project [37] and whole exome sequence data were used to conduct single variant and gene-based burden tests. The single variant analyses identified the KAT2B (lysine acetyltransferase 2B) region as significantly associated with intermittent S. aureus carriage. All 5 top genes identified in the gene-based burden test and at least 1 gene in each region suggestively associated with persistent carriage in the single variant analysis are associated in some fashion with maintenance of cellular integrity, the cytoskeleton, or the cell cycle. On the other hand, genes associated with intermittent carriage were largely associated with immune function, adipogeneisis, or inflammation. These analyses identified little evidence of overlap between genes or regions corresponding to different carriage phenotypes suggesting that each carrier state may be distinct.

Materials and Methods

Human subjects

This study and the consenting procedures were approved by the University of Texas Health Science Center Institutional Review Board (HSC-SPH-06-0225). Written informed consent was obtained from all participants before they were enrolled in the study.

Microbiologic testing

Specimens were collected from the nares using dry, unmoistened sterile BBL™ CultureSwabs™ Liquid Stuart swabs. Swabs were inserted into the patient’s nostril approximately 1 inch from the edge from the anterior nares placing the swab in proximity with the inferior and middle concha and rolled several times. Bar-coded specimen tubes were stored and shipped at 4°C to the University of Texas Health Science Center at Houston School of Public Health for processing.

To identify and characterize S. aureus from specimens containing mixed flora, nasal swabs were inoculated on manitol salt agar (MSA) plates (Remel Inc., Lenexa, KS) as described [38]. Following inoculation of primary plates, swabs were broken off into tryptic soy broth for enrichment (TSB) (Remel Inc.). The enrichment broths were vortexed for 10 seconds to ensure that any bacteria still attached to the swab were released into the media and the samples subsequently incubated at 37°C for 48 hours and re-plated on secondary MSA plates. Gram staining of respective colonies that turned MSA plates yellow were used to ensure that selected colonies possessed S. aureus morphology. Presumptive S. aureus colonies were streaked on blood agar (BA) (Quad Five, Ryegate, MT) and TSB agar and incubated at 37°C for 24 h.

Following incubations on BA and TSB agar from the primary and secondary MSA plates, colonies were subjected to catalase (Sigma, St. Louis, MO) and coagulase testing (BactiStaph® Latex 450, Remel Inc.). Positive tests were considered diagnostic for S. aureus. The identification of S. aureus was also confirmed genetically by PCR amplification and sequencing of a fragment of the spa gene for 1598/1662 (96%) of the isolates as done previously [38]. The second MSA plates streaked from the overnight liquid broth cultures were examined for additional growth, and colonies with S. aureus morphology were isolated and tested as above. Once isolates were defined as S. aureus, their respective susceptibilities to methicillin were determined using the E-test® (AB Biodisk, Biomerieux, I’Etoile, France). Methicillin resistance was defined by growth at antibiotic concentrations ≥4 μg/ml. All confirmed S. aureus isolates were stored at -80°C [38].

Definition of the S. aureus carriage phenotypes

Carriage status was determined for individuals from whom nasal swabs were collected at two time points, 2 weeks apart (14±3 days) as described previously [39]. Carriers were defined by S. aureus positive cultures at either visit, and intermittent carriers were S. aureus positive at either the first or second visit but not both. Non carriers were negative for S. aureus at both visits [39].

Genome-wide association studies and generation of whole genome imputation data

Subjects (n = 858) were eligible for this study because of prior participation in genome-wide association studies for diabetes [40]. Genotyping was performed at the Center for Inherited Disease Research using the Affymetrix Genome-wide SNPArray 6.0 assay with sample- and SNP-level genotyping quality control performed as described in Below et al. [40]. Imputations were carried out in the full Starr County sample, cleaned of ethnic outliers and including 1,616 unrelated (pairwise identity by descent ≤ 0.3) [41] individuals of which 858 met inclusion criteria in the present study. A set of autosomal scaffold SNPs were selected to drive imputation by excluding those with: 1) minor allele frequency <1%, 2) Hardy-Weinberg p-values < 10−4 in the full sample 3) missingness >10% in the full sample and 4) all ambiguous strand (AT/CG) SNPs. Individual-level missingness is <5% in all samples. 603,042 scaffold SNPs were carried forward into a two-step imputation strategy: i) pre-phasing using the program SHAPEIT [42] and ii) Imputation from the reference panel into the estimated haplotypes with IMPUTE v2 [43–45]. Imputations were done in conjunction with the T2D-GENES consortium as part of a larger set of some 13,000 multiethnic samples. SNPs with imputation quality ≥ 0.8 and minor allele frequency > 0.05 were carried forward for single variant analyses. Population stratification was evaluated using EIGENSOFT on a subset of directly genotyped SNPs pruned for local and long distance linkage disequilibrium as described in Patterson et al. [46].

Analyses were conducted by comparing persistent S. aureus carriers to noncarriers or intermittent carriers to noncarriers. Persistent carriers were defined as unrelated [41] individuals passing genotyping quality control and testing positive for colonization of S. aureus at both of two time points, 11 to 17 days apart (n = 141). Genes located within 50 kilobases of signals comprised of at least 4 SNPs and study-based minor allele frequency > 0.05 with a p value <10−5 were considered suggestively significant. For each region showing association, we identified a sentinel marker, defined as the most significant SNP meeting all quality control thresholds (locus zoom plots, Figs A-L, in S1 File).

Associations of the imputed genetic markers with S. aureus carrier status were tested with the program SNPTEST v2 [44] using frequentist association tests, based on an additive model. To control for genotype uncertainty, we used the missing data likelihood score test (the score method). All association analyses corrected for ancestry using the first and second principal components from EIGENSOFT as covariates, and all analyses were run once including diabetes status as a covariate and once excluding diabetes status in the model.

Generation and analysis of whole exome single variants

Whole exome sequence data were available for a subset of 792 participants (131 persistent carriers, 88 intermittent carriers, and 573 non-carriers, as defined above). These were part of a larger group sequenced as part of the T2D-GENES Consortium at the Broad Institute using Agilent Truseq capture reagents on Illumina HiSeq2000 instruments.

Association tests of the 1,011 common (minor allele frequency > 0.05) single variants present in the exome sequence data were performed using logistic regression in the program PLINK v2 [47]. As above, association analyses were corrected for ancestry using the first and second principal components, and all analyses were run including and excluding diabetes status as a covariate in the model. These results were combined with the imputed data results in common Manhattan plots (Figs 1 and 2 and Figs M-N in S1 File).

thumbnail
Download:
Fig 1. Manhattan (a) and QQ plots (b) of results of single variant logistic regression of persistent S. aureus carriage versus non-carrier, including PC1 and PC2 as covariates.

The x-axis represents the chromosome number and each dot represents a single polymorphic variant with minor allele frequency greater than 0.05. QQ plot shows the observed versus expected p-values for the same variants shown in (a). Grey shading indicates the 95% confidence interval, the solid line indicates the expected null distribution, and the dotted line indicates the slope after lambda correction for genomic control. The 1,011 common variants identified by whole exome sequencing are shown as x’s in the Manhattan plots.

https://doi.org/10.1371/journal.pone.0142130.g001

thumbnail
Download:
Fig 2. Manhattan (a) and QQ plots (b) of results of single variant logistic regression of intermittent S. aureus carriage versus non-carrier, including PC1 and PC2 as covariates.

The x-axis represents the chromosome number and each dot represents a single polymorphic variant with minor allele frequency greater than 0.05. QQ plot shows the observed versus expected p-values for the same variants shown in (a). Grey shading indicates the 95% confidence interval, the solid line indicates the expected null distribution, and the dotted line indicates the slope after lambda correction for genomic control. The 1,011 common variants identified by whole exome sequencing are shown as x’s in the Manhattan plots.

https://doi.org/10.1371/journal.pone.0142130.g002

Gene-based analysis of whole exome sequence data

We used the Variant Annotation Analysis and Search Tool (VAAST) to identify genes associated with increased risk of S. aureus colonization [48, 49]. For quality-control purposes, we removed sites with missingness >10% in the full sample. We also used the rate option to set the maximum expected disease allele frequency to 0.05, as we expect to be powered to detect effects of common variants in single variant tests. The top two principle components from EIGENSOFT were used as covariates in all VAAST analyses, and analyses were performed with and without diabetes status as a covariate, as above. Statistical significance was assessed using a covariate-adjusted randomization test as previously described [50, 51]; p-value confidence intervals were calculated using a Poisson approximation based on the number of successes in the randomization test. Genome-wide significance thresholds for the gene-based tests were calculated from the number of genes tested (0.05/18665 = 2.68×10−6).

Results

Population demographics and S. aureus carriage determination

Carriage status was established by collecting and analyzing swabs for the presence of S. aureus on 2 occasions from a single nostril 11–17 days apart on 858 Mexican Americans from Starr County, TX, USA [39]. A summary of demographic information for these individuals, who were eligible due to prior participation in a genome-wide association study for type 2 diabetes, are presented in Table 1 [40]. Participants testing positive for S. aureus on 2 separate occasions were defined as persistent carriers (n = 141), participants testing positive once were defined as intermittent carriers (n = 97), and participants testing negative on both occasions were defined as non-carriers (n = 620) as previously described [38, 39].

thumbnail
Download:
Table 1. Demographic information for study participants by S. aureus carriage phenotype.

https://doi.org/10.1371/journal.pone.0142130.t001

Single variant association tests

Single variant association tests of persistent S. aureus carriage identified 5 loci as suggestively significant (p value ≤ 10−5, as defined in the Methods) are summarized in Fig 1 and Table 2, namely MKLN1 (muskelin 1), SORBS1 (sorbin and SH3 domain containing 1), SLC1A2 (solute carrier family 1) SORBS1, a region intergenic between EPB41L4B (erythrocyte membrane protein band 4.1 like 4B) and PTPN3 (cytoskeletal-associated protein tyrosine phosphatase), and a region downstream of FGF3 (fibroblast growth factor 3). MKLN1 encodes an intracellular mediator of cell morphology and cytoskeletal responses [52, 53]. SORBS1 is involved in insulin signaling and SLC1A2 is a member of the solute transporter family. EPB41L4B and PTPN3 are involved in membrane-cytoskeletal interactions while FGF3 is a member of the fibroblast growth factor family of genes. MKLN1 has been previously associated with childhood asthma [54], SORBS1 with suicide risk (46) and childhood obesity in Hispanics [55], SLC1A2 with fatty acid levels [56], essential tremor [57–59], and other traits [58, 59], EPB41L4B with wound healing [60], PTPN3 with cancer [61], and FGF3 with breast cancer [62] and deafness [63, 64]. Whole exome sequencing identified 1,011 common variants (minor allele frequency > 0.05). These are shown as x’s in the Manhattan plots. In no case did any of these variants reach a suggestive level indicating that it is unlikely that there are common protein-coding variants of substantial effect. LocusZoom [65] plots for each top locus highlight LD (linkage disequilibrium) patterns among the top SNPs and show multiple SNPs in LD blocks being associated (Figs A-E in S1 File).

thumbnail
Download:
Table 2. SNPs reaching suggestive significance in single variant logistic regression for persistent S. aureus carrier vs. non-carrier (top) and intermittent S. aureus carrier vs. non-carrier (bottom), including PC1 and PC2 as covariates.

https://doi.org/10.1371/journal.pone.0142130.t002

In addition, we carried out single variant association tests of intermittent carriage of S. aureus, defined as individuals testing positive for S. aureus colonization at either visit compared to non S. aureus carriers (Fig 2, Table 2). The 7 regions suggestively associated (as defined above) with intermittent carriage include a genome-wide significant finding on chromosome 3 at rs61440199 (p value 8.68 x 10−9) that is intronic to KAT2B (lysine acetyltransferase 2B) (also known as PCAF; p300/CBP-associated factor), a gene associated with post-traumatic stress disorder [66], mean arterial blood pressure [67], adipogenesis [68], development of T regulatory cells [69], and recently shown to be a potential regulator of inflammatory responses following infection with S. aureus in a mouse model of disease (Table 2) [70]. Other signals were at or near UBE2E2 (ubiquitin-conjugating enzyme E2E 2), a gene that has been associated with risk to gestational and type 2 diabetes [71–73], ICK (intestinal cell [MAK-like] kinase), and ROBO1 (roundabout, axon guidance receptor, homolog 1), which encodes a member of the immunoglobulin gene superfamily and plays a role in axon guidance and neuronal precursor cell migration (Table 2). A SNP highly correlated with ROBO1 expression in the brain has been reproducibly associated with reading disabilities [74, 75], and SNPs mapped to ROBO1 have been associated with levels of liver enzymes [76] and other pQTLs [77]. Three sentinel SNPs were intergenic between (RELT-like 1) and PGM2 (phosphoglucomutase 2), between genes LOC283585 and GALC, and between ZNF532 (zinc finger protein 532) and MALT1 (mucosa associated lymphoid tissue lymphoma translocation gene 1) (Table 2). GALC encodes the enzyme β-galactocerebrosidase, mutations in which are responsible for Krabbe disease [78, 79]. Homozygous mutations in MALT1 have been associated with immunodeficiency [80–82] (Table 2). MALT1 has also been associated with multiple-sclerosis [83]. No common (minor allele frequency >0.05) variants in the whole exome sequencing data reached p value < 10−5 (shown as x’s in the Manhattan plot, Fig 2). LocusZoom plots for each top locus highlight LD patterns among top SNPs (Figs F-L in S1 File). It is notable that the signals for persistent carriage of S. aureus appear to be largely independent of signals for intermittent carriage of S. aureus. Of all top findings, only rs61440199 (KAT2B) and rs16993852 (RELL1) show nominal evidence of association in both persistent and intermittent carriage of S. aureus. Diabetes stratified and non-stratified analyses of both persistent and intermittent carriage gave highly concordant results across all analyses (Figs M-N in S1 File).

Gene-based tests of functional variants

The program VAAST [49, 50] was used to identify genes enriched for functional rare variation in cases based on next generation whole exome sequence data. In the analysis of persistent carriers (131 cases, Table 1) versus non-carriers (573 controls, Table 1) of S. aureus, one gene, FAM123C (APC membrane recruitment protein 3), approached genome-wide significance (p value 6.50 x 10−6) (Table 3). Other top gene-based findings include NGEF (neuronal guanine nucleotide exchange factor, p value 1.22 x 10−5), CCDC69 (coiled-coil domain containing 69, p value 1.40 x 10−5), ERP29 (endoplasmic reticulum protein 29, p value 3.72 x 10−5), and TSGA10IP (testis-specific protein 10-interacting protein, p value 7.45 x 10−5 (Table 3 and Fig O in S1 File). In the analysis of intermittent carriers (88 cases, Table 1) versus non-carriers (573 controls, Table 1) top gene-based findings included SLC4A4 (bicarbonate cotransporter, member 4, p value 2.27 x 10−4), TSPAN11 (tetraspanin 11, p value 1.98 x 10−4), TPO (thyroid peroxidase, p value 4.05 x 10−4), ZNF280D (zinc finger protein 280D, p value 3.76 x 10−4), and CSF2RB (colony stimulating factor 2 receptor, beta, low-affinity, p value 4.15 x 10−4) (Table 3 and Fig P in S1 File). Specific variant enrichment and predicted function of variants driving top gene-based findings are shown in Table 3 and gene functions are discussed below.

thumbnail
Download:
Table 3. Top findings from gene-based burden tests of rare functional variation in VAAST for persistent S. aureus carriers versus non-carriers (top) and intermittent S. aureus carriers versus non-carriers (bottom); including PC1, and PC2 as covariates.

https://doi.org/10.1371/journal.pone.0142130.t003

As in the single variant analysis of 1000 Genomes imputed data and common variants from the exome sequence data, the gene-based findings in analyses of intermittent carriers of S. aureus appear to be largely independent between analyses of persistent versus intermittent carriage groups (Table 3). Only genes CCDC69 and ZNF280D reach nominal levels of significance (p value < 0.05) in both tests. CSF2RB shows suggestive enrichment of missense variation in both analyses, and may constitute a gene involved in general S. aureus carriage susceptibility (Table 3). Shared signals should be interpreted with caution given that the non-carriage control group is the same in both tests, and thus the two tests are not strictly independent. As observed for the single variant association tests, diabetes stratified and non-stratified analyses gave highly concordant results across all analyses (Figs Q-R in S1 File). Manhattan and QQ plots suggest the type 1 error for both single variant and gene-based tests are well controlled (Figs M-R in S1 File).

We used the Disease Association Protein-Protein Link Evaluator (DAPPLE) [84] to identify interactions between proteins encoded by the top 5 candidate genes in the persistent versus non-carrier and intermittent versus non-carrier VAAST runs. DAPPLE searches for protein-protein interactions among a candidate gene list; a significant number of protein-protein interaction may indicate a shared molecular pathway relevant to S. aureus susceptibility. In the analysis of persistent carriers versus non-carriers we did not detect any direct protein-protein interactions. However, among the top 5 genes identified from the intermittent carriers versus non-carriers run, we found that TPO is directly interacting with CSF2RB (Fig S in S1 File). The p-value for observing at least one interaction among the top 5 genes is 0.008; the p-values for observing at least one interacting protein for TPO and CSF2RB are 0.015 and 0.014, respectively.

Replication of previously identified loci

We compared our persistent carriage single variant and gene-based test results to all sites previously reported in genetic analyses of S. aureus (S1 Table) [20–23, 27, 30–32, 34, 35, 85]. When the variant in question was not present in our post-quality control imputed or exome sequenced variant lists, and therefore not analyzed in this study, we identified the best proxy variant by assessing linkage disequilibrium patterns in the Mexican-American (MEX) reference population within the 1000 Genomes Project data (release 27). In these cases, statistics for the variant with the highest linkage disequilibrium r2 are provided in S1 Table.

With the exception of CDK7 (discussed below) our findings do not replicate the genes and variants described in 2 previously conducted genome wide association studies, possibly because of several differences between these prior studies and the current study (described in the Discussion)[34, 35]. We found suggestive evidence of association at rs4918120 (p value 0.034) a SNP previously identified by Nelson et al. [34] in Caucasian inpatients; however, we observed the opposite direction of effect of the T allele (odds ratio 0.70 versus 1.68, see S1 Table). Interrogation of our single variant test results for intermittent carriage at previously reported loci yielded replication at three loci identified by Ye et al. [35]: rs12696090 (p value 0.0214), rs7643377 (p value 0.0081), and rs9867210 (p value 0.0079), however as before; we find opposite direction of effect at each of these loci (S1 Table).

We also examined our gene-based test results for replication of previous findings at genes near previously associated SNPs and genes. CDK7 (cyclin-dependent kinase) (gene-based p value 0.040) replicates findings by Ye et al. [35] who studied genetic risk of hospital-based S. aureus infection in Caucasians and identified CDK7 using gene-based tests in the program VEGAS (S1 Table).

Discussion

This was the first genome-wide association study of S. aureus carriage states in a community-based representative population. This approach is significantly different from previously described genome-wide association studies that were carried out in the context of S. aureus infections [34, 35]. We found genome-wide significance at 1 gene region and 11 other regions meeting suggestive levels of significance for association with persistent and intermittent carriage states by single variant analysis. We also reported the 5 top findings from gene-based tests of persistent and intermittent carriage. The lack of overlap in signals between gene-based tests of rare functional variation and single-variant tests suggested that genome-wide association signals were not driven by coding sequence variation. Non-genic regulatory factors affecting gene expression levels or post-translational modifications may also affect carriage phenotypes.

We found that top signals associated with persistent and intermittent carriage captured genes of different cellular functions. Genome-wide single variant analysis identified 5 gene regions suggestively associated with persistent carriage. Gene-based rare variant analysis identified 5 genes in association with persistent carriage. Near genome-wide significance was observed only for FAM123C (p value < 6.50 x 10−6). Each of these genes (except for TSGA10IP, which has not been previously described to our knowledge) was involved with cellular growth, tissue homeostasis, and/or cancer [86–92]. It should be noted however, that TSGA10IP (TSGA10 interacting protein) interacts with TSGA10, a protein also associated with cancer and that binds cytoskeletal proteins (e.g., vimentin and actin-γ1) [93, 94].

In analyses of persistent S. aureus carriage, all of the top 5 findings from gene-based tests and all regions identified in the single variant analysis harbored at least 1 gene associated with either regulation of cell growth or maintenance of cellular integrity (e.g., tight junctions) [95, 96]. Conversely, a minority of genes identified in previous genome-wide association studies of S. aureus infection were involved in cell cycle, cellular growth, or cellular integrity (S1 Table) [34, 35]. These differences are important for 2 reasons: i) carriage and infection are not mutually exclusive i.e., the S. aureus carriage status of individuals was not established in relation to the infections described in the previous genome-wide association studies, and ii) susceptibility to infections in hospital environments may not accurately reflect an individual's susceptibility to an infectious agent. Hospital environments in and of themselves place patients at increased risk for infections with numerous pathogens including S. aureus, an agent responsible for more healthcare-associated infections and surgical site infections than any other pathogen [97].

Genome-wide association analysis of intermittent carriage identified a different set of genes from those identified in association with persistent carriage. This analysis identified 7 gene regions. The top signal (rs61440199) was genome-wide significant (p value 8.68 x 10−9) and intronic to KAT2B. This gene was of particular interest since its expression in mice was affected by the nature of the infecting S. aureus strain [70]. In addition, KAT2B has been linked to immune function, cancer progression, and adipogenesis [68, 98–100]. The association of KAT2B with cancer progression/cell cycle was also shared by SGOL, ROBO1, and ICK, and represents the only functional overlap with genes identified with persistent and intermittent carriage of S. aureus [101–107]. The other themes observed in the context of genes associated with intermittent carriage were genes associated with both adipogenesis and inflammation/immunity (KAT2B, ZNF532, RELL1, FOXO9, MALT1) [68, 80, 98, 100, 101, 103, 108–112]. In light of sample ascertainment for diabetes in this cohort [40], a gene in 1 region, UBE2E2, was of interest because of prior associations with diabetes risk [71–73], however, stratification for diabetes provided highly concordant results with the unstratified analysis (data not shown). Our gene-based analyses did not model complications that present in diabetic patients (e.g., obesity, immune function, elevated blood glucose levels) that may alter susceptibility to intermittent carriage. The number of adipogenesis genes linked to intermittent carriage may be of significance in light of recent studies that identified a protective role for adipose tissue in a murine model of S. aureus skin infection, suggesting that immune factors produced by adipose tissues (e.g., antimicrobial peptides) may play a role in intermittent carriage [112].

Although gene-based analyses of rare functional variants failed to identify any genome-wide significant differences in association with intermittent carriage, a top signal, CSF2RB, demonstrated concordance of burden in both persistent and intermittent carriers (p value 7.04 x 10−4 and p value < 4.15 x 10−4, respectively). CSF2RB codes for CD131, the common β receptor subunit for IL-3, IL-5, and GM-CSF (granulocyte/monocyte colony stimulating factor) that in mice was shown to play a role in regulating Th2 type immune responses [113]. In addition, CD131 stimulated the recruitment of neutrophils (which are a key innate immune component) and controlled the homeostasis of tissue dendritic cells [114, 115]. In addition, DAPPLE analysis identified a significant protein-protein interaction between the CSF2RB and TPO gene products. TPO is critical to the production of thyroid hormones that can impact immune function and is also associated with mucinosis (myxedema), a disease characterized by increased glycosaminoglycan deposition in the skin [116, 117]. Other than CSF2RB, no other top finding in the gene-based tests were even modestly associated with both persistent and intermittent carriage.

Results from the 2 previously described genome-wide association studies identified a number of loci with statistical significance. However, those associations were for the most part not replicated in our studies or previous work [9, 18–23, 27, 30, 32, 34, 35, 85, 118]. Lack of replication between studies may be due to population differences, the impact of the respective colonizing/infecting S. aureus strains (and their relationship with distinct human genetic determinants), study design (i.e., S. aureus infection versus carriage), and the size of respective populations examined [9, 22]. Replication of 1 gene identified by gene-based tests was observed in the context of persistent carriage that identified CDK7 (p value 0.041) from the VEGAS gene test conducted by Ye et al. (S1 Table). We also assessed gene-based evidence of replication in our analyses of intermittent carriage versus non-carriage and found no support for previously identified genes (data not shown).

Previous colonization studies have suggested that the 3 described staphylococcal carriage phenotypes (persistent, intermittent, and non-carriers) be modified to include only 2 carriage phenotypes: persistent carriers and intermittent/non-carriers [17]. However gene targets identified in the present S. aureus carriage genome-wide association study suggested that each phenotype is distinct. That the genome-wide association and rare variant analyses identified relatively little functional similarity between persistent and intermittent carriers may suggest underlying differences between these 2 carriage states. An alternate explanation is that these studies lacked sufficient power to identify common factors across the carriage states. Despite the recommendation of previous studies to consider intermittent and non-carriers as a single group, this reclassification would require ignoring the differences that exist between these 2 carriage states. It is clear, however, that persistent carriers represent the most distinct carriage state. This is supported by colonization studies that demonstrated that non-carriers (and decolonized intermittent carriers) artificially inoculated with S. aureus in the nares cleared the bacteria over a similar time period (4 days for non-carriers and 14 days for intermittent carriers) compared to persistent carriers (decolonized and then re-inoculated) that still harbored the S. aureus inoculum >154 days later [17]. Persistent carriers also had a different antibody profile against some staphylococcal virulence factors compared to the indistinguishable profile described for non-carriers and intermittent carriers [17]. In addition, persistent carriers that were decolonized and re-inoculated with a heterogeneous mix of S. aureus isolates were more likely to be re-colonized with their original colonizing isolate suggestive of an intimate association between the colonizing strain and the host [17].

This difference between persistent carriers and intermittent carriers (and intermittent carriers and non-carriers) is further accentuated by the function of the genes associated with the respective carriage states. Almost all determinants associated with persistent carriage were associated with cellular integrity, morphology, and growth, functions that directly hold the potential of impacting the host/pathogen interface that establishes environments permissive to persistent carriage.

Attachment to host surfaces is requisite for colonization and infection of host tissues by pathogens. S. aureus possesses an arsenal of adhesins capable of binding an array of host extracellular matrix (ECM) components. These components include fibrinogen, fibronectin, collagen, cytokeratin 10, elastin, heparan sulfate proteoglycans, von Willebrand factor, bone sialoprotein, vitronectin, and prothrombin that all facilitate the colonization of diverse tissues and accounts in part for the myriad of diseases than can result following infection with this pathogen [119–121]. It is not surprising therefore that host polymorphisms potentially affecting cellular integrity, morphology and growth could also impact colonization with different pathogens or strains of the same pathogen.

That various potential genes identified by the genome-wide association study (e.g., ALDH18A1, EPB41L4B, FGF3, and FGF4) and all but 1 gene identified in the rare variant analysis have been shown to possess tumorogenic potential should not be surprising since various genes shown to play roles in the progression of various cancers also play critical roles in wound healing, cellular migration, cellular integrity, and angiogenesis [60, 122, 123]. Polymorphisms in these gene products or any gene products with the potential of altering the structural integrity of the host cell could potentially impact staphylococcal colonization.

Focal adhesions represent large, multi-protein complexes that are closely associated with cell surface integrins that span the eukaryotic plasma membrane linking the cellular cytoskeleton to the ECM (surrounding the cell) [124]. Most integrins and their respective focal adhesions are expressed in the epidermis and regulate epithelial cell homeostasis by mediating cell adhesion processes (and signaling) critical to tissue repair following injury [124]. Of the gene targets identified in association with S. aureus persistent carriage, EHM2, PTPN3, SORS1, and MKLN1 can impact the integrity of focal adhesions that in turn alters the cytoskeleton [53, 120, 124–133].

SORBS1 encodes CAP (Cbl-Associated Protein) [129, 132, 134] that affects insulin receptor signaling and also functions as a cytoskeletal regulatory protein [129]. In fibroblasts, when CAP associates with actin stress fibers, focal adhesion kinase binds CAP, and CAP over-expression induces the development of actin stress fibers and focal adhesions that physically link intracellular actin bundles to the extracellular substrates of many cell types [127, 135]. Various pathogens like S. aureus usurp focal adhesions as a means of triggering their uptake by various non-professional antigen presenting cells, including epithelia/endothelial cells, osteoclasts, kidney cells, fibroblasts and keratinocytes [120]. S. aureus possess various fibronectin binding proteins (e.g., FnbpA, FnbpB, ClfA, ClfB) that facilitate coating the bacterial surface with this matrix molecule that in turn binds to α5β1 integrins resulting in the formation of a molecular bridge linking S. aureus to the host cell [125]. This interaction triggers the recruitment of focal adhesion proteins that further alter the cytoskeleton facilitating attachment, invasion, and the ability to persist in their hosts [125]. The importance of this interaction for the successful attachment/invasion of human cells by staphylococci was demonstrated by generating fnbpA/fnbpB-deficient S. aureus that less effectively infected epithelial cells and in a mastitis model caused less severe disease [136, 137]. Furthermore, cells unable to form focal adhesions were resistant to integrin α5β1-mediated cellular invasion by S. aureus [120, 127].

EHM2 is a member of the 4.1R, ezrin, radixin, moesin (FERM) protein superfamily consisting of over 40 proteins that contain the characteristic 3-lobed FERM domain on the N-terminus that binds various cell membrane-associated proteins and lipids and the spectrin/actin binding domain (SABD) at the C-terminus [126]. The PTPN3 gene product also belongs to the FERM family and is a protein phospatase that is a structural constituent of the cytoskeletal shown to play a role in T cell activation, maintenance of tight junction integrity (between the cell membrane and the cytoskeleton) and both EHM2 and PTPN3 gene products are associated with focal adhesions [95, 128, 130, 133, 138]. EHM2 expression has been observed on wounds undergoing healing (primarily at the wound's leading edge) functioning as a positive regulator of keratinocyte adhesion and motility in addition to affecting the rates of cellular invasion and adhesion to collagen via regulation of matrix metaloprotease 9 (MMP9) i.e., EHM2 knockdown cells expressed significantly reduced levels of MMP9. This is of interest in the context of S. aureus since up- or down-regulation of MMP9 levels has been shown to affect disease progression resulting from S. aureus infections, that is, MMP9 levels that are either too high or too low can negatively affect wound healing and MMP9-deficient mice poorly controlled S. aureus infections [60, 126, 131, 139–143]. In addition, MMPs play critical roles in tissue remodeling (including the maintenance of the ECM), altering immune cell migration and infiltration patterns, and impacted inflammation by exerting effects on cytokines and chemokines [143, 144]. As it relates to S. aureus colonization, a role for MMP9 has yet to be described; however, staphylococcal lipoteichoic acid has been shown to increase production of MMP9 in middle ear epithelial cells suggesting that increased MMP9 levels could be involved in progression of otitis media [141].

Unlike EHM2, PTPN3, and SORS1; the MKLN1 gene product muskelin mediates ECM binding via complex mechanisms involving interactions between different thrombospondin-1 (TSP-1) domains and various ligands (expressed by different cell types) including integrins, proteoglycans, or integrin-associated proteins. Alterations to muskelin expression levels altered attachment to TSP-1 in association with subtle changes to the organization of focal contacts [53]. Since TSP-1 has also been shown to serve as a ligand for S. aureus, polymorphisms in MKLN1 could alter staphylococcal binding or prevent clearance of S. aureus since TSP-1 breakdown products function as antimicrobial peptides (AMPs) that have broad antibacterial properties affecting both Gram-positive and -negative bacteria [145–148].

Homozygous mutations in ALDH18A1 (or other genes e.g., PYCR1, ATP6V0A2) can result in a heteogenous group of rare diseases characterized by loose or wrinkly skin known as cutis laxa [149–151]. Histologic analysis of skin from cutis laxa patients identified reduced elastin levels with less-well defined collagen fibers lacking the characteristic wavy morphology. In addition, collagen I and III levels were significantly reduced, and fibroblasts harvested from the dermis presented with reduced growth rates [149]. The majority of studies that have examined genes resulting in this rare condition have only described case reports of patients with homozygous mutations, making it difficult to interpret how polymorphisms with a less pronounced phenotype present at the cytoskeletal level.

Although adherence to host surfaces also represents a component of intermittent carriage (i.e., the organism has to attach to host tissues even if this association is transient), the intermittent periods of carriage, carriage of different strains over time, carriage of multiple strains, and the reduced S. aureus inocula recoverable from the nares of intermittent carriers suggests that different determinants are associated with this phenotype [17]. This is emphasized by the observation that the majority of gene targets associated with intermittent carriage were also associated with immune function/inflammation.

Our data suggested that determinants associated with persistent carriage and intermittent carriage differed. A limitation to the present study was the analysis of only two nasal swabs to establish carriage. Even though Nouwen et al. established that the 'two-culture' rule was 93.6% reliable [39] and numerous studies have used this approach to establish S. aureus carriage phenotypes [20, 38, 118, 152–156] there exists room for classification error. Second, because only one nostril was sampled some participants may have been misclassified as intermittent or non carriers based on one study that described differences in S. aureus carriage between colonization [157]; however, two other studies did not identify any differences [158, 159]. It should be noted that samples that were collected and analyzed for the present study were of the ciliated pseudostratified columnar epithelium associate the inferior and middle concha and not the nonkeratinized, squamous epithelium present in the anterior nares and used to establish S. aureus carriage by other studies. Furthermore, due to population differences and power we should be cautious in making assumptions with regard to specific genes associated with respective carriage states. We should therefore further dissect the observation that persistent carriage of S. aureus is affected primarily by polymorphisms at the host/pathogen interface and that intermittent carriage is more likely impacted by environmental factors combined with the heterogeneity of the host immune response.

Supporting Information

S1 File.

Figs A-E. LocusZoom plots of each top finding in the single variant association analyses of persistent S. aureus carriage versus non-carrier. (A) EPB41L4B, (B) LINC-PINT, (C) SORBS1, ALDH18A1, (D) SLC1A2, and (E) FGF4, FGF3. Figs F-L. LocusZoom plots of each top finding in the single variant association analyses of persistent S. aureus carriage versus non-carrier. (F) KAT2B, (G) UBE2E2, MIR548AC, (H) ROBO1, (I) RELL1, (J) GSTA4, ICK, FBXO9, (K) LOC283585, GALC, and (L) ZNF532. Fig M. Manhattan (a) and QQ plots (b) of results of single variant logistic regression of persistent S. aureus carriage versus non-carrier, including diabetes, PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represents a single polymorphic variant with minor allele frequency greater than 0.05. QQ plot shows the observed versus expected p-values for the same variants shown in (a). Grey shading indicates the 95% confidence interval, the solid line indicates the expected null distribution, and the dotted line indicates the slope after lambda correction for genomic control. The 1,011 common variants identified by whole exome sequencing are shown as x’s in the Manhattan plots. Fig N. Manhattan (a) and QQ plots (b) of results of single variant logistic regression of intermittent S. aureus carriage versus non-carrier, including diabetes, PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represents a single polymorphic variant with minor allele frequency greater than 0.05. QQ plot shows the observed versus expected p-values for the same variants shown in (a). Grey shading indicates the 95% confidence interval, the solid line indicates the expected null distribution, and the dotted line indicates the slope after lambda correction for genomic control. The 1,011 common variants identified by whole exome sequencing are shown as x’s in the Manhattan plots. Figs O. Manhattan (a) and QQ plots (b) of results of gene-based burden tests of rare functional variation in VAAST for persistent S. aureus carriage versus non-carrier including PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represent one protein-coding gene. QQ plot shows the observed versus expected p-values for all protein-coding genes, grey shading represents 95% confidence interval, the red line indicates the null distribution of p-values. Fig P. Manhattan (a) and QQ plots (b) of results of gene-based burden tests of rare functional variation in VAAST for intermittent S. aureus carriage versus non-carrier including PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represent one protein-coding gene. QQ plot shows the observed versus expected p-values for all protein-coding genes, grey shading represents 95% confidence interval, the red line indicates the null distribution of p-values. Fig Q. Manhattan (a) and QQ plots (b) of results of gene-based burden tests of rare functional variation in VAAST for persistent S. aureus carriage versus non-carrier including diabetes, PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represent one protein-coding gene. QQ plot shows the observed versus expected p-values for all protein-coding genes, grey shading represents 95% confidence interval, the red line indicates the null distribution of p-values. Fig R. Manhattan (a) and QQ plots (b) of results of gene-based burden tests of rare functional variation in VAAST for intermittent S. aureus carriage versus non-carrier including diabetes, PC1, and PC2 as covariates. The x-axis represents the chromosome number and each dot represent one protein-coding gene. QQ plot shows the observed versus expected p-values for all protein-coding genes, grey shading represents 95% confidence interval, the red line indicates the null distribution of p-values. Fig S. Protein-protein interactions among top-5 candidate genes in the gene-based test of intermittent carriersversus non-carriers analysis. Red: genes that encode proteins with direct interactions to another top-5 candidate; blue: genes that encode proteins with second-degree interactions to another top-5 candidate; grey: genes that are not top-5 candidates, but encode proteins interacting with at least two top-5 candidates. The figure was generated using DAPPLE software.

https://doi.org/10.1371/journal.pone.0142130.s001

(DOCX)

S1 Table. Previous genes and SNPs associated with S. aureus carriage or infection.

https://doi.org/10.1371/journal.pone.0142130.s002

(XLSX)

Acknowledgments

This work was supported by NIH grants R01 AI085014-01A1 (E.L. Brown and C.L. Hanis), R01 GM104390 (C.D. Huff and H. Hao), HL102830 (C.L. Hanis), DK085501 (C.L. Hanis) and P30DK020595 (G.I. Bell). H. Hao is also supported by the M.D. Anderson Cancer Center Odyssey Program. This work was also partially supported from a grant from the Kleberg Foundation to E.L.B. Genotyping imputation and whole exome sequencing were performed as part of our involvement in the T2D-GENES Consortium and we acknowledge those efforts. We also express appreciation to the field staff in Starr County who contacted and collected the necessary participant data and samples. Lastly, we thank the participants for their generous and willing participation.

Author Contributions

Conceived and designed the experiments: ELB CLH. Performed the experiments: ELB CLH HTE RSF. Analyzed the data: ELB JEB HTE RSF CH HH LEP GIB CLH DA DAR. Contributed reagents/materials/analysis tools: ELB CLH JEB CH HH DAR. Wrote the paper: ELB JEB CLH HH CH DAR RSF HTE.

References

  1. 1. Kwiatkowski DP. How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet. 2005;77(2):171–92. pmid:16001361.
  2. 2. Liu R, Paxton WA, Choe S, Ceradini D, Martin SR, Horuk R, et al. Homozygous defect in HIV-1 coreceptor accounts for resistance of some multiply-exposed individuals to HIV-1 infection. Cell. 1996;86(3):367–77. pmid:8756719.
  3. 3. Bellamy R. Genetic susceptibility to tuberculosis. Clin Chest Med. 2005;26(2):233–46, vi. pmid:15837108.
  4. 4. Frodsham AJ, Zhang L, Dumpis U, Taib NA, Best S, Durham A, et al. Class II cytokine receptor gene cluster is a major locus for hepatitis B persistence. Proc Natl Acad Sci U S A. 2006;103(24):9148–53. pmid:16757563.
  5. 5. Lindesmith L, Moe C, Marionneau S, Ruvoen N, Jiang X, Lindblad L, et al. Human susceptibility and resistance to Norwalk virus infection. Nat Med. 2003;9(5):548–53. pmid:12692541.
  6. 6. Zeidler M, Stewart G, Cousens SN, Estibeiro K, Will RG. Codon 129 genotype and new variant CJD. Lancet. 1997;350(9078):668. pmid:9288076.
  7. 7. Levine MM, Nalin DR, Rennels MB, Hornick RB, Sotman S, Van Blerk G, et al. Genetic susceptibility to cholera. Ann Hum Biol. 1979;6(4):369–74. pmid:394667.
  8. 8. Malaty HM, Engstrand L, Pedersen NL, Graham DY. Helicobacter pylori infection: genetic and environmental influences. A study of twins. Ann Intern Med. 1994;120(12):982–6. pmid:8185146.
  9. 9. Sollid JU, Furberg AS, Hanssen AM, Johannessen M. Staphylococcus aureus: determinants of human carriage. Infect Genet Evol. 2014;21:531–41. pmid:23619097.
  10. 10. Embil J, Ramotar K, Romance L, Alfa M, Conly J, Cronk S, et al. Methicillin-resistant Staphylococcus aureus in tertiary care institutions on the Canadian prairies 1990–1992. Infection control and hospital epidemiology: the official journal of the Society of Hospital Epidemiologists of America. 1994;15(10):646–51. pmid:7844335.
  11. 11. Hill PC, Wong CG, Voss LM, Taylor SL, Pottumarthy S, Drinkovic D, et al. Prospective study of 125 cases of Staphylococcus aureus bacteremia in children in New Zealand. Pediatr Infect Dis J. 2001;20(9):868–73. pmid:11734766.
  12. 12. Tong SY, Bishop EJ, Lilliebridge RA, Cheng AC, Spasova-Penkova Z, Holt DC, et al. Community-associated strains of methicillin-resistant Staphylococcus aureus and methicillin-susceptible S. aureus in indigenous Northern Australia: epidemiology and outcomes. J Infect Dis. 2009;199(10):1461–70. pmid:19392622.
  13. 13. Lowy FD. Staphylococcus aureus infections. N Engl J Med. 1998;339:520–32. pmid:9709046
  14. 14. Kalmeijer MD, Coertjens H, van Nieuwland-Bollen PM, Bogaers-Hofman D, de Baere GA, Stuurman A, et al. Surgical site infections in orthopedic surgery: the effect of mupirocin nasal ointment in a double-blind, randomized, placebo-controlled study. Clin Infect Dis. 2002;35(4):353–8. pmid:12145715.
  15. 15. von Eiff C, Becker K, Machka K, Stammer H, Peters G. Nasal carriage as a source of Staphylococcus aureus bacteremia. Study Group. N Engl J Med. 2001;344(1):11–6. pmid:11136954.
  16. 16. Wertheim HF, Vos MC, Ott A, van Belkum A, Voss A, Kluytmans JA, et al. Risk and outcome of nosocomial Staphylococcus aureus bacteraemia in nasal carriers versus non-carriers. Lancet. 2004;364(9435):703–5. pmid:15325835.
  17. 17. van Belkum A, Verkaik NJ, de Vogel CP, Boelens HA, Verveer J, Nouwen JL, et al. Reclassification of Staphylococcus aureus nasal carriage types. J Infect Dis. 2009;199(12):1820–6. Epub 2009/05/08. pmid:19419332.
  18. 18. Vandenbergh MF, Verbrugh HA. Carriage of Staphylococcus aureus: epidemiology and clinical relevance. J Lab Clin Med. 1999;133(6):525–34. pmid:10360626.
  19. 19. von Kockritz-Blickwede M, Rohde M, Oehmcke S, Miller LS, Cheung AL, Herwald H, et al. Immunological mechanisms underlying the genetic predisposition to severe Staphylococcus aureus infection in the mouse model. Am J Pathol. 2008;173(6):1657–68. pmid:18974303.
  20. 20. van den Akker EL, Nouwen JL, Melles DC, van Rossum EF, Koper JW, Uitterlinden AG, et al. Staphylococcus aureus nasal carriage is associated with glucocorticoid receptor gene polymorphisms. J Infect Dis. 2006;194(6):814–8. pmid:16941349.
  21. 21. van Belkum A, Emonts M, Wertheim H, de Jongh C, Nouwen J, Bartels H, et al. The role of human innate immune factors in nasal colonization by Staphylococcus aureus. Microbes Infect. 2007;9(12–13):1471–7. pmid:17913546.
  22. 22. Emonts M, Uitterlinden AG, Nouwen JL, Kardys I, Maat MP, Melles DC, et al. Host polymorphisms in interleukin 4, complement factor H, and C-reactive protein associated with nasal carriage of Staphylococcus aureus and occurrence of boils. J Infect Dis. 2008;197(9):1244–53. pmid:18422436.
  23. 23. Emonts M, de Jongh CE, Houwing-Duistermaat JJ, van Leeuwen WB, de Groot R, Verbrugh HA, et al. Association between nasal carriage of Staphylococcus aureus and the human complement cascade activator serine protease C1 inhibitor (C1INH) valine vs. methionine polymorphism at amino acid position 480. FEMS Immunol Med Microbiol. 2007;50(3):330–2. pmid:17498209.
  24. 24. Cole AM, Tahk S, Oren A, Yoshioka D, Kim YH, Park A, et al. Determinants of Staphylococcus aureus nasal carriage. Clin Diagn Lab Immunol. 2001;8(6):1064–9. pmid:11687441.
  25. 25. Schumann RR. Host cell-pathogen interface: molecular mechanisms and genetics. Vaccine. 2004;22 Suppl 1:S21–4. pmid:15576197.
  26. 26. Clementi M, Forabosco P, Amadori A, Zamarchi R, De Silvestro G, Di Gianantonio E, et al. CD4 and CD8 T lymphocyte inheritance. Evidence for major autosomal recessive genes. Hum Genet. 1999;105(4):337–42. pmid:10543402.
  27. 27. Panierakis C, Goulielmos G, Mamoulakis D, Maraki S, Papavasiliou E, Galanakis E. Staphylococcus aureus nasal carriage might be associated with vitamin D receptor polymorphisms in type 1 diabetes. Int J Infect Dis. 2009. pmid:19411183.
  28. 28. O'Brien LM, Walsh EJ, Massey RC, Peacock SJ, Foster TJ. Staphylococcus aureus clumping factor B (ClfB) promotes adherence to human type I cytokeratin 10: implications for nasal colonization. Cell Microbiol. 2002;4(11):759–70. pmid:12427098.
  29. 29. Moore CE, Segal S, Berendt AR, Hill AV, Day NP. Lack of association between Toll-like receptor 2 polymorphisms and susceptibility to severe disease caused by Staphylococcus aureus. Clin Diagn Lab Immunol. 2004;11(6):1194–7. pmid:15539529.
  30. 30. Ruimy R, Angebault C, Djossou F, Dupont C, Epelboin L, Jarraud S, et al. Are host genetics the predominant determinant of persistent nasal Staphylococcus aureus carriage in humans? J Infect Dis. 2010;202(6):924–34. pmid:20677941.
  31. 31. Fode P, Larsen AR, Feenstra B, Jespersgaard C, Skov RL, Stegger M, et al. Genetic variability in beta-defensins is not associated with susceptibility to Staphylococcus aureus bacteremia. PLoS One. 2012;7(2):e32315. pmid:22384213; PubMed Central PMCID: PMC3285211.
  32. 32. Nurjadi D, Herrmann E, Hinderberger I, Zanger P. Impaired beta-defensin expression in human skin links DEFB1 promoter polymorphisms with persistent Staphylococcus aureus nasal carriage. J Infect Dis. 2013;207(4):666–74. pmid:23204181.
  33. 33. Stappers MH, Thys Y, Oosting M, Plantinga TS, Ioana M, Reimnitz P, et al. TLR1, TLR2, and TLR6 gene polymorphisms are associated with increased susceptibility to complicated skin and skin structure infections. J Infect Dis. 2014;210(2):311–8. pmid:24511099.
  34. 34. Nelson CL, Pelak K, Podgoreanu MV, Ahn SH, Scott WK, Allen AS, et al. A genome-wide association study of variants associated with acquisition of Staphylococcus aureus bacteremia in a healthcare setting. BMC infectious diseases. 2014;14:83. pmid:24524581; PubMed Central PMCID: PMC3928605.
  35. 35. Ye Z, Vasco DA, Carter TC, Brilliant MH, Schrodi SJ, Shukla SK. Genome wide association study of SNP-, gene-, and pathway-based approaches to identify genes influencing susceptibility to Staphylococcus aureus infections. Frontiers in genetics. 2014;5:125. pmid:24847357; PubMed Central PMCID: PMC4023021.
  36. 36. System N. National Nosocomial Infections Surveillance (NNIS) System Report, Data Summary from January 1990-May 1999, issued June 1999. A report from the NNIS System. Am J Infect Control. 1999;27(6):520–32. pmid:10586157.
  37. 37. Hayes MG, Pluzhnikov A, Miyake K, Sun Y, Ng MC, Roe CA, et al. Identification of type 2 diabetes genes in Mexican Americans through genome-wide association studies. Diabetes. 2007;56(12):3033–44. pmid:17846124.
  38. 38. Leung NS, Padgett P, Robinson DA, Brown EL. Prevalence and behavioural risk factors of Staphylococcus aureus nasal colonization in community-based injection drug users. Epidemiology and infection. 2014:1–10. pmid:25434806.
  39. 39. Nouwen JL, Ott A, Kluytmans-Vandenbergh MF, Boelens HA, Hofman A, van Belkum A, et al. Predicting the Staphylococcus aureus nasal carrier state: derivation and validation of a "culture rule". Clin Infect Dis. 2004;39(6):806–11. pmid:15472812.
  40. 40. Below JE, Gamazon ER, Morrison JV, Konkashbaev A, Pluzhnikov A, McKeigue PM, et al. Genome-wide association and meta-analysis in populations from Starr County, Texas, and Mexico City identify type 2 diabetes susceptibility loci and enrichment for expression quantitative trait loci in top signals. Diabetologia. 2011;54(8):2047–55. pmid:21647700; PubMed Central PMCID: PMC3761075.
  41. 41. Staples J, Nickerson DA, Below JE. Utilizing graph theory to select the largest set of unrelated individuals for genetic analysis. Genet Epidemiol. 2013;37(2):136–41. pmid:22996348; PubMed Central PMCID: PMC3770842.
  42. 42. Delaneau O, Marchini J, Zagury JF. A linear complexity phasing method for thousands of genomes. Nature methods. 2012;9(2):179–81. pmid:22138821.
  43. 43. Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nature genetics. 2012;44(8):955–9. pmid:22820512; PubMed Central PMCID: PMC3696580.
  44. 44. Marchini J, Howie B. Genotype imputation for genome-wide association studies. Nature reviews Genetics. 2010;11(7):499–511. pmid:20517342.
  45. 45. Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007;39(7):906–13. pmid:17572673.
  46. 46. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLoS genetics. 2006;2(12):e190. pmid:17194218; PubMed Central PMCID: PMC1713260.
  47. 47. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007;81(3):559–75. pmid:17701901; PubMed Central PMCID: PMC1950838.
  48. 48. Yandell M, Huff C, Hu H, Singleton M, Moore B, Xing J, et al. A probabilistic disease-gene finder for personal genomes. Genome Res. 2011;21(9):1529–42. pmid:21700766; PubMed Central PMCID: PMC3166837.
  49. 49. Hu H, Huff CD, Moore B, Flygare S, Reese MG, Yandell M. VAAST 2.0: improved variant classification and disease-gene identification using a conservation-controlled amino acid substitution matrix. Genet Epidemiol. 2013;37(6):622–34. pmid:23836555; PubMed Central PMCID: PMC3791556.
  50. 50. Hu H, Roach JC, Coon H, Guthery SL, Voelkerding KV, Margraf RL, et al. A unified test of linkage analysis and rare-variant association for analysis of pedigree sequence data. Nat Biotechnol. 2014;32(7):663–9. pmid:24837662; PubMed Central PMCID: PMC4157619.
  51. 51. Epstein MP, Duncan R, Jiang Y, Conneely KN, Allen AS, Satten GA. A permutation procedure to correct for confounders in case-control studies, including tests of rare variation. Am J Hum Genet. 2012;91(2):215–23. pmid:22818855; PubMed Central PMCID: PMC3415546.
  52. 52. Valiyaveettil M, Bentley AA, Gursahaney P, Hussien R, Chakravarti R, Kureishy N, et al. Novel role of the muskelin-RanBP9 complex as a nucleocytoplasmic mediator of cell morphology regulation. The Journal of cell biology. 2008;182(4):727–39. pmid:18710924; PubMed Central PMCID: PMC2518711.
  53. 53. Adams JC, Seed B, Lawler J. Muskelin, a novel intracellular mediator of cell adhesive and cytoskeletal responses to thrombospondin-1. The EMBO journal. 1998;17(17):4964–74. pmid:9724633; PubMed Central PMCID: PMC1170825.
  54. 54. Ding L, Abebe T, Beyene J, Wilke RA, Goldberg A, Woo JG, et al. Rank-based genome-wide analysis reveals the association of ryanodine receptor-2 gene variants with childhood asthma among human populations. Human genomics. 2013;7:16. pmid:23829686; PubMed Central PMCID: PMC3708719.
  55. 55. Comuzzie AG, Cole SA, Laston SL, Voruganti VS, Haack K, Gibbs RA, et al. Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PloS one. 2012;7(12):e51954. pmid:23251661; PubMed Central PMCID: PMC3522587.
  56. 56. Guan W, Steffen BT, Lemaitre RN, Wu JH, Tanaka T, Manichaikul A, et al. Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium. Circulation Cardiovascular genetics. 2014;7(3):321–31. pmid:24823311; PubMed Central PMCID: PMC4123862.
  57. 57. Thier S, Lorenz D, Nothnagel M, Poremba C, Papengut F, Appenzeller S, et al. Polymorphisms in the glial glutamate transporter SLC1A2 are associated with essential tremor. Neurology. 2012;79(3):243–8. pmid:22764253; PubMed Central PMCID: PMC3398434.
  58. 58. Garcia-Martin E, Martinez C, Serrador M, Alonso-Navarro H, Navacerrada F, Agundez JA, et al. SLC1A2 rs3794087 variant and risk for migraine. Journal of the neurological sciences. 2014;338(1–2):92–5. pmid:24412224.
  59. 59. Jimenez-Jimenez FJ, Alonso-Navarro H, Martinez C, Zurdo M, Turpin-Fenoll L, Millan-Pascual J, et al. The solute carrier family 1 (glial high affinity glutamate transporter), member 2 gene, SLC1A2, rs3794087 variant and assessment risk for restless legs syndrome. Sleep medicine. 2014;15(2):266–8. pmid:24424098.
  60. 60. Bosanquet DC, Ye L, Harding KG, Jiang WG. Expressed in high metastatic cells (Ehm2) is a positive regulator of keratinocyte adhesion and motility: The implication for wound healing. Journal of dermatological science. 2013;71(2):115–21. pmid:23664528.
  61. 61. Gao Q, Zhao YJ, Wang XY, Guo WJ, Gao S, Wei L, et al. Activating mutations in PTPN3 promote cholangiocarcinoma cell proliferation and migration and are associated with tumor recurrence in patients. Gastroenterology. 2014;146(5):1397–407. pmid:24503127.
  62. 62. Turnbull C, Ahmed S, Morrison J, Pernet D, Renwick A, Maranian M, et al. Genome-wide association study identifies five new breast cancer susceptibility loci. Nature genetics. 2010;42(6):504–7. pmid:20453838; PubMed Central PMCID: PMC3632836.
  63. 63. Riazuddin S, Ahmed ZM, Hegde RS, Khan SN, Nasir I, Shaukat U, et al. Variable expressivity of FGF3 mutations associated with deafness and LAMM syndrome. BMC medical genetics. 2011;12:21. pmid:21306635; PubMed Central PMCID: PMC3042908.
  64. 64. Sensi A, Ceruti S, Trevisi P, Gualandi F, Busi M, Donati I, et al. LAMM syndrome with middle ear dysplasia associated with compound heterozygosity for FGF3 mutations. American journal of medical genetics Part A. 2011;155A(5):1096–101. pmid:21480479.
  65. 65. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26(18):2336–7. pmid:20634204; PubMed Central PMCID: PMC2935401.
  66. 66. Wolf EJ, Rasmusson AM, Mitchell KS, Logue MW, Baldwin CT, Miller MW. A genome-wide association study of clinical symptoms of dissociation in a trauma-exposed sample. Depression and anxiety. 2014;31(4):352–60. pmid:24677629; PubMed Central PMCID: PMC3984628.
  67. 67. Simino J, Sung YJ, Kume R, Schwander K, Rao DC. Gene-alcohol interactions identify several novel blood pressure loci including a promising locus near SLC16A9. Frontiers in genetics. 2013;4:277. pmid:24376456; PubMed Central PMCID: PMC3860258.
  68. 68. Jin Q, Wang C, Kuang X, Feng X, Sartorelli V, Ying H, et al. Gcn5 and PCAF regulate PPARgamma and Prdm16 expression to facilitate brown adipogenesis. Molecular and cellular biology. 2014;34(19):3746–53. pmid:25071153; PubMed Central PMCID: PMC4187735.
  69. 69. Xiong Y, Svingen PA, Sarmento OO, Smyrk TC, Dave M, Khanna S, et al. Differential coupling of KLF10 to Sin3-HDAC and PCAF regulates the inducibility of the FOXP3 gene. American journal of physiology Regulatory, integrative and comparative physiology. 2014;307(6):R608–20. pmid:24944246; PubMed Central PMCID: PMC4166759.
  70. 70. Modak R, Das Mitra S, Vasudevan M, Krishnamoorthy P, Kumar M, Bhat AV, et al. Epigenetic response in mice mastitis: Role of histone H3 acetylation and microRNA(s) in the regulation of host inflammatory gene expression during Staphylococcus aureus infection. Clinical epigenetics. 2014;6(1):12. pmid:25075227; PubMed Central PMCID: PMC4114167.
  71. 71. Kim JY, Cheong HS, Park BL, Baik SH, Park S, Kim S, et al. Putative association between UBE2E2 polymorphisms and the risk of gestational diabetes mellitus. Gynecological endocrinology: the official journal of the International Society of Gynecological Endocrinology. 2013;29(10):904–8. pmid:23862583.
  72. 72. Yamauchi T, Hara K, Maeda S, Yasuda K, Takahashi A, Horikoshi M, et al. A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B. Nature genetics. 2010;42(10):864–8. pmid:20818381.
  73. 73. DIAbetes Genetics Replication Meta-analysis Consortium, Asian Genetic Epidemiology Network Type 2 Diabetes Consortium, South Asian Type 2 Diabetes Consortium, Mexican American Type 2 Diabetes Consortium, Type 2 Diabetes Genetic Exploration by Nex-generation sequencing in muylti-Ethnic Samples Consortium, Mahajan A, et al. Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility. Nature genetics. 2014;46(3):234–44. pmid:24509480; PubMed Central PMCID: PMC3969612.
  74. 74. Lamminmaki S, Massinen S, Nopola-Hemmi J, Kere J, Hari R. Human ROBO1 regulates interaural interaction in auditory pathways. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2012;32(3):966–71. pmid:22262894.
  75. 75. Tran C, Wigg KG, Zhang K, Cate-Carter TD, Kerr E, Field LL, et al. Association of the ROBO1 gene with reading disabilities in a family-based analysis. Genes, brain, and behavior. 2014;13(4):430–8. pmid:24612512.
  76. 76. Park TJ, Hwang JY, Go MJ, Lee HJ, Jang HB, Choi Y, et al. Genome-wide association study of liver enzymes in korean children. Genomics & Informatics. 2013;11(3):149–54. pmid:24124411; PubMed Central PMCID: PMC3794088.
  77. 77. Melzer D, Perry JR, Hernandez D, Corsi AM, Stevens K, Rafferty I, et al. A genome-wide association study identifies protein quantitative trait loci (pQTLs). PLoS genetics. 2008;4(5):e1000072. pmid:18464913; PubMed Central PMCID: PMC2362067.
  78. 78. Hill CH, Graham SC, Read RJ, Deane JE. Structural snapshots illustrate the catalytic cycle of beta-galactocerebrosidase, the defective enzyme in Krabbe disease. Proceedings of the National Academy of Sciences of the United States of America. 2013;110(51):20479–84. pmid:24297913; PubMed Central PMCID: PMC3870757.
  79. 79. Kardas F, Uzak AS, Hossain MA, Sakai N, Canpolat M, Yikilmaz A. A novel homozygous GALC mutation: very early onset and rapidly progressive Krabbe disease. Gene. 2013;517(1):125–7. pmid:23276707.
  80. 80. Jabara HH, Ohsumi T, Chou J, Massaad MJ, Benson H, Megarbane A, et al. A homozygous mucosa-associated lymphoid tissue 1 (MALT1) mutation in a family with combined immunodeficiency. J Allergy Clin Immunol. 2013;132(1):151–8. pmid:23727036; PubMed Central PMCID: PMC3700575.
  81. 81. McKinnon ML, Rozmus J, Fung SY, Hirschfeld AF, Del Bel KL, Thomas L, et al. Combined immunodeficiency associated with homozygous MALT1 mutations. The Journal of allergy and clinical immunology. 2014;133(5):1458–62, 62 e1–7. pmid:24332264.
  82. 82. Turvey SE, Durandy A, Fischer A, Fung SY, Geha RS, Gewies A, et al. The CARD11-BCL10-MALT1 (CBM) signalosome complex: Stepping into the limelight of human primary immunodeficiency. The Journal of allergy and clinical immunology. 2014;134(2):276–84. pmid:25087226; PubMed Central PMCID: PMC4167767.
  83. 83. International Multiple Sclerosis Genetics C, Wellcome Trust Case Control C, Sawcer S, Hellenthal G, Pirinen M, Spencer CC, et al. Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature. 2011;476(7359):214–9. pmid:21833088; PubMed Central PMCID: PMC3182531.
  84. 84. Rossin EJ, Lage K, Raychaudhuri S, Xavier RJ, Tatar D, Benita Y, et al. Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology. PLoS genetics. 2011;7(1):e1001273. pmid:21249183; PubMed Central PMCID: PMC3020935.
  85. 85. Vuononvirta J, Toivonen L, Grondahl-Yli-Hannuksela K, Barkoff AM, Lindholm L, Mertsola J, et al. Nasopharyngeal bacterial colonization and gene polymorphisms of mannose-binding lectin and toll-like receptors 2 and 4 in infants. PLoS One. 2011;6(10):e26198. pmid:22022564; PubMed Central PMCID: PMC3192769.
  86. 86. Bambang IF, Lee YK, Richardson DR, Zhang D. Endoplasmic reticulum protein 29 regulates epithelial cell integrity during the mesenchymal-epithelial transition in breast cancer cells. Oncogene. 2013;32(10):1240–51. pmid:22543584.
  87. 87. Brauburger K, Akyildiz S, Ruppert JG, Graeb M, Bernkopf DB, Hadjihannas MV, et al. Adenomatous polyposis coli (APC) membrane recruitment 3, a member of the APC membrane recruitment family of APC-binding proteins, is a positive regulator of Wnt-beta-catenin signalling. The FEBS journal. 2014;281(3):787–801. pmid:24251807.
  88. 88. Chen S, Zhang D. Friend or foe: Endoplasmic reticulum protein 29 (ERp29) in epithelial cancer. FEBS open bio. 2015;5:91–8. pmid:25709888; PubMed Central PMCID: PMC4329646.
  89. 89. Nathke IS. The adenomatous polyposis coli protein: the Achilles heel of the gut epithelium. Annual review of cell and developmental biology. 2004;20:337–66. pmid:15473844.
  90. 90. Pal D, Wu D, Haruta A, Matsumura F, Wei Q. Role of a novel coiled-coil domain-containing protein CCDC69 in regulating central spindle assembly. Cell cycle. 2010;9(20):4117–29. pmid:20962590; PubMed Central PMCID: PMC3055196.
  91. 91. Shamah SM, Lin MZ, Goldberg JL, Estrach S, Sahin M, Hu L, et al. EphA receptors regulate growth cone dynamics through the novel guanine nucleotide exchange factor ephexin. Cell. 2001;105(2):233–44. pmid:11336673.
  92. 92. Kawai H, Kobayashi M, Hiramoto-Yamaki N, Harada K, Negishi M, Katoh H. Ephexin4-mediated promotion of cell migration and anoikis resistance is regulated by serine 897 phosphorylation of EphA2. FEBS open bio. 2013;3:78–82. pmid:23772378; PubMed Central PMCID: PMC3668535.
  93. 93. Tanaka R, Ono T, Sato S, Nakada T, Koizumi F, Hasegawa K, et al. Over-expression of the testis-specific gene TSGA10 in cancers and its immunogenicity. Microbiol Immunol. 2004;48(4):339–45. pmid:15107545.
  94. 94. Roghanian A, Jones DC, Pattisapu JV, Wolfe J, Young NT, Behnam B. Filament-associated TSGA10 protein is expressed in professional antigen presenting cells and interacts with vimentin. Cell Immunol. 2010;265(2):120–6. pmid:20797700.
  95. 95. Yang Q, Tonks NK. Isolation of a cDNA clone encoding a human protein-tyrosine phosphatase with homology to the cytoskeletal-associated proteins band 4.1, ezrin, and talin. Proc Natl Acad Sci U S A. 1991;88(14):5949–53. pmid:1648725; PubMed Central PMCID: PMC51999.
  96. 96. Deng X, Shibata H, Ninomiya H, Tashiro N, Iwata N, Ozaki N, et al. Association study of polymorphisms in the excitatory amino acid transporter 2 gene (SLC1A2) with schizophrenia. BMC psychiatry. 2004;4:21. pmid:15296513; PubMed Central PMCID: PMC514708.
  97. 97. Sievert DM, Ricks P, Edwards JR, Schneider A, Patel J, Srinivasan A, et al. Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2009–2010. Infection control and hospital epidemiology: the official journal of the Society of Hospital Epidemiologists of America. 2013;34(1):1–14. pmid:23221186.
  98. 98. Brasacchio D, Noori T, House C, Brennan AJ, Simpson KJ, Susanto O, et al. A functional genomics screen identifies PCAF and ADA3 as regulators of human granzyme B-mediated apoptosis and Bid cleavage. Cell Death Differ. 2014;21(5):748–60. pmid:24464226; PubMed Central PMCID: PMC3978306.
  99. 99. Mizuguchi Y, Specht S, Lunz JG 3rd, Isse K, Corbitt N, Takizawa T, et al. Cooperation of p300 and PCAF in the control of microRNA 200c/141 transcription and epithelial characteristics. PLoS One. 2012;7(2):e32449. pmid:22384255; PubMed Central PMCID: PMC3284570.
  100. 100. van Loosdregt J, Vercoulen Y, Guichelaar T, Gent YY, Beekman JM, van Beekum O, et al. Regulation of Treg functionality by acetylation-mediated Foxp3 protein stabilization. Blood. 2010;115(5):965–74. pmid:19996091.
  101. 101. Chen T, Wu D, Moskaluk CA, Fu Z. Distinct expression patterns of ICK/MAK/MOK protein kinases in the intestine implicate functional diversity. PLoS One. 2013;8(11):e79359. pmid:24244486; PubMed Central PMCID: PMC3820702.
  102. 102. Chetaille P, Preuss C, Burkhard S, Cote JM, Houde C, Castilloux J, et al. Mutations in SGOL1 cause a novel cohesinopathy affecting heart and gut rhythm. Nat Genet. 2014;46(11):1245–9. pmid:25282101.
  103. 103. Fu Z, Kim J, Vidrich A, Sturgill TW, Cohn SM. Intestinal cell kinase, a MAP kinase-related kinase, regulates proliferation and G1 cell cycle progression of intestinal epithelial cells. Am J Physiol Gastrointest Liver Physiol. 2009;297(4):G632–40. pmid:19696144; PubMed Central PMCID: PMC2763805.
  104. 104. Jiang L, Wang Y, Rong Y, Xu L, Chu Y, Zhang Y, et al. miR-1179 promotes cell invasion through SLIT2/ROBO1 axis in esophageal squamous cell carcinoma. International journal of clinical and experimental pathology. 2015;8(1):319–27. pmid:25755718; PubMed Central PMCID: PMC4348900.
  105. 105. Nakayama T, Mizohata E, Yamashita T, Nagatoishi S, Nakakido M, Iwanari H, et al. Structural features of interfacial tyrosine residue in ROBO1 fibronectin domain-antibody complex: Crystallographic, thermodynamic, and molecular dynamic analyses. Protein Sci. 2015;24(3):328–40. pmid:25492858; PubMed Central PMCID: PMC4353359.
  106. 106. Ravnskjaer K, Hogan MF, Lackey D, Tora L, Dent SY, Olefsky J, et al. Glucagon regulates gluconeogenesis through KAT2B- and WDR5-mediated epigenetic effects. J Clin Invest. 2013;123(10):4318–28. pmid:24051374; PubMed Central PMCID: PMC3784539.
  107. 107. Zia A, Bhatti A, John P, Kiani AK. Data interpretation: deciphering the biological function of Type 2 diabetes associated risk loci. Acta diabetologica. 2015. pmid:25585593.
  108. 108. Lee KW, Kwak SH, Ahn BY, Lee HM, Jung HS, Cho YM, et al. F-box only protein 9 is required for adipocyte differentiation. Biochem Biophys Res Commun. 2013;435(2):239–43. pmid:23643813.
  109. 109. Sica GL, Zhu G, Tamada K, Liu D, Ni J, Chen L. RELT, a new member of the tumor necrosis factor receptor superfamily, is selectively expressed in hematopoietic tissues and activates transcription factor NF-kappaB. Blood. 2001;97(9):2702–7. pmid:11313261.
  110. 110. Wei S, Zhang L, Zhou X, Du M, Jiang Z, Hausman GJ, et al. Emerging roles of zinc finger proteins in regulating adipogenesis. Cellular and molecular life sciences: CMLS. 2013;70(23):4569–84. pmid:23760207; PubMed Central PMCID: PMC4100687.
  111. 111. Yang IV, Wade CM, Kang HM, Alper S, Rutledge H, Lackford B, et al. Identification of novel genes that mediate innate immunity using inbred mice. Genetics. 2009;183(4):1535–44. pmid:19805818; PubMed Central PMCID: PMC2787437.
  112. 112. Zhang LJ, Guerrero-Juarez CF, Hata T, Bapat SP, Ramos R, Plikus MV, et al. Innate immunity. Dermal adipocytes protect against invasive Staphylococcus aureus skin infection. Science. 2015;347(6217):67–71. pmid:25554785; PubMed Central PMCID: PMC4318537.
  113. 113. Asquith KL, Ramshaw HS, Hansbro PM, Beagley KW, Lopez AF, Foster PS. The IL-3/IL-5/GM-CSF common receptor plays a pivotal role in the regulation of Th2 immunity and allergic airway inflammation. J Immunol. 2008;180(2):1199–206. pmid:18178860.
  114. 114. Greter M, Helft J, Chow A, Hashimoto D, Mortha A, Agudo-Cantero J, et al. GM-CSF controls nonlymphoid tissue dendritic cell homeostasis but is dispensable for the differentiation of inflammatory dendritic cells. Immunity. 2012;36(6):1031–46. pmid:22749353; PubMed Central PMCID: PMC3498051.
  115. 115. Khajah M, Millen B, Cara DC, Waterhouse C, McCafferty DM. Granulocyte-macrophage colony-stimulating factor (GM-CSF): a chemoattractive agent for murine leukocytes in vivo. J Leukoc Biol. 2011;89(6):945–53. pmid:21393420.
  116. 116. Safer JD. Thyroid hormone action on skin. Dermato-endocrinology. 2011;3(3):211–5. pmid:22110782; PubMed Central PMCID: PMC3219173.
  117. 117. Robinson MV, Obut TA, Melnikova EV, Trufakin VA. Parameters of cellular and humoral immunity in experimental hyperthyroidism and its correction. Bulletin of experimental biology and medicine. 2014;156(4):473–5. pmid:24771430.
  118. 118. Roghmann MC, Johnson JK, Stine OC, Lydecker AD, Ryan KA, Mitchell BD, et al. Persistent Staphylococcus aureus colonization is not a strongly heritable trait in Amish families. PLoS One. 2011;6(2):e17368. pmid:21386985; PubMed Central PMCID: PMC3046241.
  119. 119. Patti JM, Allen BL, McGavin MJ, Höök M. MSCRAMM-mediated adherence of microorganisms to host tissues. Annu Rev Microbiol. 1994;48:585–617. pmid:7826020
  120. 120. Hauck CR, Ohlsen K. Sticky connections: extracellular matrix protein recognition and integrin-mediated cellular invasion by Staphylococcus aureus. Curr Opin Microbiol. 2006;9(1):5–11. pmid:16406780.
  121. 121. Hayashida A, Amano S, Park PW. Syndecan-1 promotes Staphylococcus aureus corneal infection by counteracting neutrophil-mediated host defense. J Biol Chem. 2011;286(5):3288–97. pmid:21127056; PubMed Central PMCID: PMC3030334.
  122. 122. Roy D, Calaf GM. Allelic loss at chromosome 11q13 alters FGF3 gene expression in a human breast cancer progression model. Oncology reports. 2014;32(6):2445–52. pmid:25333703.
  123. 123. Harding JJ, Abou-Alfa GK. Predicting responsiveness to sorafenib: can the determination of FGF3/FGF4 amplifications enrich for clinical benefit? Hepatobiliary surgery and nutrition. 2014;3(4):168–71. pmid:25202691; PubMed Central PMCID: PMC4141298.
  124. 124. Duperret EK, Ridky TW. Focal adhesion complex proteins in epidermis and squamous cell carcinoma. Cell cycle. 2013;12(20):3272–85. pmid:24036537; PubMed Central PMCID: PMC3885638.
  125. 125. Agerer F, Lux S, Michel A, Rohde M, Ohlsen K, Hauck CR. Cellular invasion by Staphylococcus aureus reveals a functional link between focal adhesion kinase and cortactin in integrin-mediated internalisation. Journal of cell science. 2005;118(Pt 10):2189–200. pmid:15855238.
  126. 126. Bosanquet DC, Ye L, Harding KG, Jiang WG. FERM family proteins and their importance in cellular movements and wound healing (review). International journal of molecular medicine. 2014;34(1):3–12. pmid:24820650.
  127. 127. Hauck CR, Agerer F, Muenzner P, Schmitter T. Cellular adhesion molecules as targets for bacterial infection. European journal of cell biology. 2006;85(3–4):235–42. pmid:16546567.
  128. 128. Itoh F, Ikuta S, Hinoda Y, Arimura Y, Ohe M, Adachi M, et al. Expression and chromosomal assignment of PTPH1 gene encoding a cytosolic protein tyrosine phosphatase homologous to cytoskeletal-associated proteins. International journal of cancer Journal international du cancer. 1993;55(6):947–51. pmid:8253532.
  129. 129. Matson SA, Pare GC, Kapiloff MS. A novel isoform of Cbl-associated protein that binds protein kinase A. Biochim Biophys Acta. 2005;1727(2):145–9. pmid:15716063.
  130. 130. Sozio MS, Mathis MA, Young JA, Walchli S, Pitcher LA, Wrage PC, et al. PTPH1 is a predominant protein-tyrosine phosphatase capable of interacting with and dephosphorylating the T cell receptor zeta subunit. J Biol Chem. 2004;279(9):7760–9. pmid:14672952.
  131. 131. Yu H, Ye L, Mansel RE, Zhang Y, Jiang WG. Clinical implications of the influence of Ehm2 on the aggressiveness of breast cancer cells through regulation of matrix metalloproteinase-9 expression. Molecular cancer research: MCR. 2010;8(11):1501–12. pmid:21047774.
  132. 132. Zhang M, Kimura A, Saltiel AR. Cloning and characterization of Cbl-associated protein splicing isoforms. Molecular medicine. 2003;9(1–2):18–25. pmid:12765336; PubMed Central PMCID: PMC1430382.
  133. 133. Zhang SH, Kobayashi R, Graves PR, Piwnica-Worms H, Tonks NK. Serine phosphorylation-dependent association of the band 4.1-related protein-tyrosine phosphatase PTPH1 with 14-3-3beta protein. J Biol Chem. 1997;272(43):27281–7. pmid:9341175.
  134. 134. Mi H, Muruganujan A, Thomas PD. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 2013;41(Database issue):D377–86. pmid:23193289; PubMed Central PMCID: PMC3531194.
  135. 135. Abercrombie M, Dunn GA. Adhesions of fibroblasts to substratum during contact inhibition observed by interference reflection microscopy. Exp Cell Res. 1975;92(1):57–62. pmid:1169157.
  136. 136. Brouillette E, Grondin G, Shkreta L, Lacasse P, Talbot BG. In vivo and in vitro demonstration that Staphylococcus aureus is an intracellular pathogen in the presence or absence of fibronectin-binding proteins. Microb Pathog. 2003;35(4):159–68. pmid:12946328.
  137. 137. Brouillette E, Talbot BG, Malouin F. The fibronectin-binding proteins of Staphylococcus aureus may promote mammary gland colonization in a lactating mouse model of mastitis. Infect Immun. 2003;71(4):2292–5. pmid:12654860; PubMed Central PMCID: PMC152093.
  138. 138. Kurella S, Yaciuk JC, Dozmorov I, Frank MB, Centola M, Farris AD. Transcriptional modulation of TCR, Notch and Wnt signaling pathways in SEB-anergized CD4+ T cells. Genes Immun. 2005;6(7):596–608. pmid:16034473; PubMed Central PMCID: PMC2593626.
  139. 139. Moleirinho S, Tilston-Lunel A, Angus L, Gunn-Moore F, Reynolds PA. The expanding family of FERM proteins. Biochem J. 2013;452(2):183–93. pmid:23662806.
  140. 140. Niebuhr M, Baumert K, Werfel T. TLR-2-mediated cytokine and chemokine secretion in human keratinocytes. Experimental dermatology. 2010;19(10):873–7. pmid:20849532.
  141. 141. Park HY, Song JW, Hong SP, In SM, Kim HJ. Lipoteichoic acid from Staphylococcus aureus induced expression of MMP-9 in human middle ear epithelial cells. International journal of pediatric otorhinolaryngology. 2012;76(4):475–9. pmid:22305687.
  142. 142. Itoh S, Hamada E, Kamoshida G, Takeshita K, Oku T, Tsuji T. Staphylococcal superantigen-like protein 5 inhibits matrix metalloproteinase 9 from human neutrophils. Infect Immun. 2010;78(7):3298–305. pmid:20479083; PubMed Central PMCID: PMC2897379.
  143. 143. Calander AM, Starckx S, Opdenakker G, Bergin P, Quiding-Jarbrink M, Tarkowski A. Matrix metalloproteinase-9 (gelatinase B) deficiency leads to increased severity of Staphylococcus aureus-triggered septic arthritis. Microbes Infect. 2006;8(6):1434–9. pmid:16697230.
  144. 144. Parks WC, Wilson CL, Lopez-Boado YS. Matrix metalloproteinases as modulators of inflammation and innate immunity. Nat Rev Immunol. 2004;4(8):617–29. pmid:15286728.
  145. 145. Hussain M, Haggar A, Peters G, Chhatwal GS, Herrmann M, Flock JI, et al. More than one tandem repeat domain of the extracellular adherence protein of Staphylococcus aureus is required for aggregation, adherence, and host cell invasion but not for leukocyte activation. Infect Immun. 2008;76(12):5615–23. pmid:18794290; PubMed Central PMCID: PMC2583574.
  146. 146. Malmsten M, Davoudi M, Schmidtchen A. Bacterial killing by heparin-binding peptides from PRELP and thrombospondin. Matrix biology: journal of the International Society for Matrix Biology. 2006;25(5):294–300. pmid:16730966.
  147. 147. Herrmann M, Suchard SJ, Boxer LA, Waldvogel FA, Lew PD. Thrombospondin binds to Staphylococcus aureus and promotes staphylococcal adherence to surfaces. Infect Immun. 1991;59(1):279–88. pmid:1987042; PubMed Central PMCID: PMC257738.
  148. 148. Rennemeier C, Hammerschmidt S, Niemann S, Inamura S, Zahringer U, Kehrel BE. Thrombospondin-1 promotes cellular adherence of gram-positive pathogens via recognition of peptidoglycan. FASEB J. 2007;21(12):3118–32. pmid:17507668.
  149. 149. Skidmore DL, Chitayat D, Morgan T, Hinek A, Fischer B, Dimopoulou A, et al. Further expansion of the phenotypic spectrum associated with mutations in ALDH18A1, encoding Delta(1)-pyrroline-5-carboxylate synthase (P5CS). American journal of medical genetics Part A. 2011;155A(8):1848–56. pmid:21739576.
  150. 150. Wolthuis DF, van Asbeck E, Mohamed M, Gardeitchik T, Lim-Melia ER, Wevers RA, et al. Cutis laxa, fat pads and retinopathy due to ALDH18A1 mutation and review of the literature. European journal of paediatric neurology: EJPN: official journal of the European Paediatric Neurology Society. 2014;18(4):511–5. pmid:24767728.
  151. 151. Fischer B, Callewaert B, Schroter P, Coucke PJ, Schlack C, Ott CE, et al. Severe congenital cutis laxa with cardiovascular manifestations due to homozygous deletions in ALDH18A1. Molecular genetics and metabolism. 2014;112(4):310–6. pmid:24913064.
  152. 152. Claassen M, Nouwen J, Fang Y, Ott A, Verbrugh H, Hofman A, et al. Staphylococcus aureus nasal carriage is not associated with known polymorphism in the Vitamin D receptor gene. FEMS Immunol Med Microbiol. 2005;43(2):173–6. pmid:15681147.
  153. 153. Lai CF, Liao CH, Pai MF, Chu FY, Hsu SP, Chen HY, et al. Nasal carriage of methicillin-resistant Staphylococcus aureus is associated with higher all-cause mortality in hemodialysis patients. Clinical journal of the American Society of Nephrology: CJASN. 2011;6(1):167–74. pmid:20947786; PubMed Central PMCID: PMC3022239.
  154. 154. VandenBergh MF, Yzerman EP, van Belkum A, Boelens HA, Sijmons M, Verbrugh HA. Follow-up of Staphylococcus aureus nasal carriage after 8 years: redefining the persistent carrier state. J Clin Microbiol. 1999;37(10):3133–40. pmid:10488166.
  155. 155. Verwer PE, Robinson JO, Coombs GW, Wijesuriya T, Murray RJ, Verbrugh HA, et al. Prevalence of nasal methicillin-resistant Staphylococcus aureus colonization in healthcare workers in a Western Australian acute care hospital. Eur J Clin Microbiol Infect Dis. 2012;31(6):1067–72. pmid:21909648.
  156. 156. Sangvik M, Olsen RS, Olsen K, Simonsen GS, Furberg AS, Sollid JU. Age- and gender-associated Staphylococcus aureus spa types found among nasal carriers in a general population: the Tromso Staph and Skin Study. J Clin Microbiol. 2011;49(12):4213–8. pmid:21998436; PubMed Central PMCID: PMC3232966.
  157. 157. Kildow BJ, Conradie JP, Robson RL. Nostrils of healthy volunteers are independent with regard to Staphylococcus aureus carriage. J Clin Microbiol. 2012;50(11):3744–6. pmid:22915611; PubMed Central PMCID: PMC3486198.
  158. 158. Yan M, Pamp SJ, Fukuyama J, Hwang PH, Cho DY, Holmes S, et al. Nasal microenvironments and interspecific interactions influence nasal microbiota complexity and S. aureus carriage. Cell host & microbe. 2013;14(6):631–40. pmid:24331461; PubMed Central PMCID: PMC3902146.
  159. 159. Gillen AL, Daycock WO, Serafin A. High MRSA carriage rate among nursing microbiolgy students. Advances in Microbiology. 2014;4:871–7.