Search

Populations and individuals differ in susceptibility to infections because of a number of factors, including host genetic variation. We previously demonstrated that differences in antibody titer, which reflect infection history, are significantly heritable. Here we attempt to identify the genetic factors influencing variation in these serological phenotypes. Blood samples from >1300 Mexican Americans were quantified for IgG antibody level against 12 common infections, selected on the basis of their reported role in cardiovascular disease risk: Chlamydia pneumoniae; Helicobacter pylori; Toxoplasma gondii; cytomegalovirus; herpes simplex I virus; herpes simplex II virus; human herpesvirus 6 (HHV6); human herpesvirus 8 (HHV8); varicella zoster virus; hepatitis A virus (HAV); influenza A virus; and influenza B virus. Pathogen-specific quantitative antibody levels were analyzed, as were three measures of pathogen burden. Genome-wide linkage and joint linkage and association analyses were performed using ~1 million SNPs. Significant linkage (lod scores >3.0) was obtained for HHV6 (on chromosome 7), HHV8 (on chromosome 6), and HAV (on chromosome 13). SNP rs4812712 on chromosome 20 was significantly associated with C. pneumoniae (P=5.3 × 10 -8). However, no genome-wide significant loci were obtained for the other investigated antibodies. We conclude that it is possible to localize host genetic factors influencing some of these antibody traits, but that further larger-scale investigations will be required to elucidate the genetic mechanisms contributing to variation in antibody levels.

Increased immunoglobulin G (IgG) response to dietary antigens can be associated with gastrointestinal dysfunction and autoimmunity. The underlying processes contributing to these adverse reactions remain largely unknown, and it is likely that genetic factors play a role. Here, we estimate heritability and attempt to localize genetic factors influencing IgG antibody levels against food-derived antigens using an integrative genomics approach. IgG antibody levels were determined by ELISA in >1,300 Mexican Americans for the following food antigens: wheat gliadin; bovine casein; and two forms of bovine serum albumin (BSA-a and BSA-b). Pedigree-based variance components methods were used to estimate additive genetic heritability (h2), perform genome-wide association analyses, and identify transcriptional signatures (based on 19,858 transcripts from peripheral blood lymphocytes). Heritability estimates were significant for all traits (0.15–0.53), and shared environment (based on shared residency among study participants) was significant for casein (0.09) and BSA-a (0.33). Genome-wide significant evidence of association was obtained only for antibody to gliadin (P = 8.57 × 10−8), mapping to the human leukocyte antigen II region, with HLA-DRA and BTNL2 as the best candidate genes. Lack of association of known celiac disease risk alleles HLA-DQ2.5 and -DQ8 with antigliadin antibodies in the studied population suggests a separate genetic etiology. Significant transcriptional signatures were found for all IgG levels except BSA-b. These results demonstrate that individual genetic differences contribute to food antigen antibody measures in this population. Further investigations may elucidate the underlying immunological processes involved.

Infection with Epstein-Barr virus (EBV) is highly prevalent worldwide, and it has been associated with infectious mononucleosis and severe diseases including Burkitt lymphoma, Hodgkin lymphoma, nasopharyngeal lymphoma, and lymphoproliferative disorders. Although EBV has been the focus of extensive research, much still remains unknown concerning what makes some individuals more sensitive to infection and to adverse outcomes as a result of infection. Here we use an integrative genomics approach in order to localize genetic factors influencing levels of Epstein Barr virus (EBV) nuclear antigen-1 (EBNA-1) IgG antibodies, as a measure of history of infection with this pathogen, in large Mexican American families. Genome-wide evidence of both significant linkage and association was obtained on chromosome 6 in the human leukocyte antigen (HLA) region and replicated in an independent Mexican American sample of large families (minimum p-value in combined analysis of both datasets is 1.4×10−15 for SNPs rs477515 and rs2516049). Conditional association analyses indicate the presence of at least two separate loci within MHC class II, and along with lymphocyte expression data suggest genes HLA-DRB1 and HLA-DQB1 as the best candidates. The association signals are specific to EBV and are not found with IgG antibodies to 12 other pathogens examined, and therefore do not simply reveal a general HLA effect.We investigated whether SNPs significantly associated with diseases in which EBV is known or suspected to play a role (namely nasopharyngeal lymphoma, Hodgkin lymphoma, systemic lupus erythematosus, and multiple sclerosis) also show evidence of associated with EBNA-1 antibody levels, finding an overlap only for the HLA locus, but none elsewhere in the genome. The significance of this work is that a major locus related to EBV infection has been identified, which may ultimately reveal the underlying mechanisms by which the immune system regulates infection with this pathogen.

Background: Intima-media thickness (IMT) of the common and internal carotid arteries is an established surrogate for atherosclerosis and predicts risk of stroke and myocardial infarction. Often IMT is measured as the average of these 2 arteries; yet, they are believed to result from separate biological mechanisms. The aim of this study was to conduct a family-based genome-wide association study (GWAS) for IMT to identify polymorphisms influencing IMT and to determine if distinct carotid artery segments are influenced by different genetic components. Methods and Results: IMT for the common and internal carotid arteries was determined through B-mode ultrasound in 772 Mexican Americans from the San Antonio Family Heart Study. A GWAS using 931219 single-nucleotide polymorphisms was undertaken with 6 internal and common carotid artery IMT phenotypes using an additive measured genotype model. The most robust association detected was for 2 single-nucleotide polymorphisms (rs16983261, rs6113474; P=1.60e−7) in complete linkage disequilibrium on chromosome 20p11 for the internal carotid artery near wall, next to the gene PAX1. We also replicated previously reported GWAS regions on chromosomes 19q13 and 7q22. We found no overlapping associations between internal and common carotid artery phenotypes at P<5.0e−6. The genetic correlation between the 2 carotid IMT arterial segments was 0.51. Conclusions: This study represents the first large-scale GWAS of carotid IMT in a non–European population and identified several novel loci. We do not detect any shared GWAS signals between common and internal carotid arterial segments, but the moderate genetic correlation implies both common and unique genetic components.

The delta-5 and delta-6 desaturases (D5D and D6D), encoded by fatty acid desaturase 1 (FADS1) and 2 (FADS2) genes, respectively, are rate-limiting enzymes in the metabolism of ω-3 and ω-6 fatty acids. The objective of this study was to identify genes influencing variation in estimated D5D and D6D activities in plasma and erythrocytes in Alaskan Eskimos (n = 761) participating in the genetics of coronary artery disease in Alaska Natives (GOCADAN) study. Desaturase activity was estimated by product: precursor ratio of polyunsaturated fatty acids. We found evidence of linkage for estimated erythrocyte D5D (eD5D) on chromosome 11q12-q13 (logarithm of odds score = 3.5). The confidence interval contains candidate genes FADS1, FADS2, 7-dehydrocholesterol reductase (DHCR7), and carnitine palmitoyl transferase 1 A, liver (CPT1A). Measured genotype analysis found association between CPT1A, FADS1, and FADS2 single-nucleotide polymorphisms (SNPs) and estimated eD5D activity (p-values between 10-28 and 10-5). A Bayesian quantitative trait nucleotide analysis showed thatrs3019594 in CPT1A, rs174541 in FADS1, and rs174568 in FADS2 had posterior probabilities > 0.8, thereby demonstrating significant statistical support for a functional effect on eD5D activity. Highly significant associations of FADS1, FADS2, and CPT1A transcripts with their respective SNPs (p-values between 1075 and 17) in Mexican Americans of the San Antonio Family Heart Study corroborated our results. These findings strongly suggest a functional role for FADS1, FADS2, and CPT1A SNPs in the variation in eD5D activity.

Individual differences in biological ageing (i.e., the rate of physiological response to the passage of time) may be due in part to genotype-specific variation in gene action. However, the sources of heritable variation in human age-related gene expression profiles are largely unknown. We have profiled genome-wide expression in peripheral blood mononuclear cells from 1240 individuals in large families and found 4472 human autosomal transcripts, representing ~4349 genes, significantly correlated with age. We identified 623 transcripts that show genotype by age interaction in addition to a main effect of age, defining a large set of novel candidates for characterization of the mechanisms of differential biological ageing. We applied a novel SNP genotype × age interaction test to one of these candidates, the ubiquilin-like gene UBQLNL, and found evidence of joint cis-association and genotype by age interaction as well as trans-genotype by age interaction for UBQLNL expression. Both UBQLNL expression levels at recruitment and cis genotype are associated with longitudinal cancer risk in our study cohort.