GWAS on predominantly European Ancestry (EA) samples have implicated over 100 loci as being associated at genome-wide significant levels with risk for Schizophrenia (SZ). Some of these risk loci are shared across ancestral groups, including within variably admixed African American (AA) samples. Regulation of mRNA expression may be involved in the etiology for some of these loci, thus suggesting utility of a transcriptomic profiling approach. In both of two previous large transcriptomic profiling studies on lymphoblastoid cell lines (LCLs) from non-overlapping EA SZ case-control samples, using arrays (N=714; Sanders et al., 2013) and RNAseq (N=1,189; Sanders et al., 2017), we found genes differentially expressed by affection status to be enriched for immune-related genes, consistent with hypothesized immune contributions to SZ risk. The single largest reported SZ GWAS on an AA sample is on the Molecular Genetics of SZ (MGS) dataset (Shi et al., 2009), from which we derived the current transcriptomic (RNAseq) sample. After various sample quality control steps including identity verification ensuring concordance between known sex (i.e., dosages of X and Y chromosomes) and RNAseq expression levels of sex-dimorphic expressed genes on chromosomes X and Y, and RNAseq-derived genotype concordance versus previous GWAS genotypes, we moved samples with at least 8 M read depth from 378 SZ cases and 386 controls to analyses. To adjust for known or potentially confounding effects, we collected data on measured (EBV load, growth rate, and energy status), RNAseq batch, ancestry (top two genotypic principal components from GWAS), and other demographic (sex, age) covariates to include in the regression analysis to identify genes differentially expressed by affection status. We will present and discuss results in several main areas: (1) degree of replication of AA transcriptomic results for previous EA findings, (2) pathway and network analyses of differentially expressed genes, (3) gene set enrichment analyses, especially for immune related genes, (4) differences compared to findings for previous EA transcriptomic work, and (5) combined analyses of the current AA and previous EA transcriptomic datasets. This work was supported NIH grant R01MH098059 using samples from the MGS collaboration.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that positive symptoms of SZ, in particular psychosis, are due to disturbed neurotransmission via the dopamine (DA) receptor D2 (DRD2). However, DA is a reactive molecule that yields various oxidative species, and thus has important non-receptor-mediated effects, with empirical evidence of cellular toxicity and neurodegeneration. Here we examine non-receptor-mediated effects of DA on gene co-expression networks and its potential role in SZ pathology. Transcriptomic profiles were measured by RNA-seq in B-cell transformed lymphoblastoid cell lines from 514 SZ cases and 690 controls, both before and after exposure to DA ex vivo (100 μM). Gene co-expression modules were identified using Weighted Gene Co-expression Network Analysis for both baseline and DA-stimulated conditions, with each module characterized for biological function and tested for association with SZ status and SNPs from a genome-wide panel. We identified seven co-expression modules under baseline, of which six were preserved in DA-stimulated data. One module shows significantly increased association with SZ after DA perturbation (baseline: P = 0.023; DA-stimulated: P = 7.8 × 10 -5 ; ΔAIC = −10.5) and is highly enriched for genes related to ribosomal proteins and translation (FDR = 4 × 10 −141 ), mitochondrial oxidative phosphorylation, and neurodegeneration. SNP association testing revealed tentative QTLs underlying module co-expression, notably at FASTKD2 (top P = 2.8 × 10 −6 ), a gene involved in mitochondrial translation. These results substantiate the role of translational machinery in SZ pathogenesis, providing insights into a possible dopaminergic mechanism disrupting mitochondrial function, and demonstrates the utility of disease-relevant functional perturbation in the study of complex genetic etiologies.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that dopaminergic over activity causes psychosis, a central feature of SZ, based on the observation that blocking dopamine (DA) improves psychotic symptoms. DA is known to have both receptor- and non-receptor-mediated effects, including oxidative mechanisms that lead to apoptosis. The role of DA-mediated oxidative processes in SZ has been little studied. Here, we have used a cell perturbation approach and measured transcriptomic profiles by RNAseq to study the effect of DA exposure on transcription in B-cell transformed lymphoblastoid cell lines (LCLs) from 514 SZ cases and 690 controls. We found that DA had widespread effects on both cell growth and gene expression in LCLs. Overall, 1455 genes showed statistically significant differential DA response in SZ cases and controls. This set of differentially expressed genes is enriched for brain expression and for functions related to immune processes and apoptosis, suggesting that DA may play a role in SZ pathogenesis through modulating those systems. Moreover, we observed a non-significant enrichment of genes near genome-wide significant SZ loci and with genes spanned by SZ-associated copy number variants (CNVs), which suggests convergent pathogenic mechanisms detected by both genetic association and gene expression. The study suggests a novel role of DA in the biological processes of immune and apoptosis that may be relevant to SZ pathogenesis. Furthermore, our results show the utility of pathophysiologically relevant perturbation experiments to investigate the biology of complex mental disorders.

Background: Metabolic diseases such as obesity and diabetes are associated with changes in high-density lipoprotein (HDL) particles, including changes in particle size and protein composition, often resulting in abnormal function. Recent studies suggested that patients with non-alcoholic fatty liver disease (NAFLD), including individuals with nonalcoholic steatohepatitis (NASH), have smaller HDL particles when compared to individuals without liver pathologies. However, no studies have investigated potential changes in HDL particle protein composition in patients with NAFLD, in addition to changes related to obesity, to explore putative functional changes of HDL which may increase the risk of cardiovascular complications. Methods: From a cohort of morbidly obese females who were diagnosed with simple steatosis (SS), NASH, or normal liver histology, we selected five matched individuals from each condition for a preliminary pilot HDL proteome analysis. HDL particles were enriched using size-exclusion chromatography, and the proteome of the resulting fraction was analyzed by liquid chromatography tandem mass spectrometry. Differences in the proteomes between the three conditions (normal, SS, NASH) were assessed using label-free quantitative analysis. Gene ontology term analysis was performed to assess the potential impact of proteomic changes on specific functions of HDL particles. Results: Of the 95 proteins identified, 12 proteins showed nominally significant differences between the three conditions. Gene ontology term analysis revealed that severity of the liver pathology may significantly impact the antithrombotic functions of HDL particles, as suggested by changes in the abundance of HDL-associated proteins such as antithrombin III and plasminogen. Conclusions: The pilot data from this study suggest that changes in the HDL proteome may impact the functionality of HDL particles in NAFLD and NASH patients. These proteome changes may alter cardio-protective properties of HDL, potentially contributing to the increased cardiovascular disease risk in affected individuals. Further validation of these protein changes by orthogonal approaches is key to confirming the role of alterations in the HDL proteome in NAFLD and NASH. This will help elucidate the mechanistic effects of the altered HDL proteome on cardioprotective properties of HDL particles.

The dopaminergic hypothesis of schizophrenia postulates a dopamine D2 receptor (DRD2)-mediated hyperdopaminergic neurotransmission in schizophrenia. However, dopamine is also known to have other, non-receptor-mediated effects, including oxidative mechanisms leading to apoptosis. While oxidative processes have been reported to contribute to dopamine neuron loss in Parkinson's disease, their role in schizophrenia remains little studied.

Currently available methods for interrogating DNA-protein interactions at individual genomic loci have significant limitations, and make it difficult to work with unmodified cells or examine single-copy regions without specific antibodies. In this study, we describe a physiological application of the Hybridization Capture of Chromatin-Associated Proteins for Proteomics (HyCCAPP) methodology we have developed. Both novel and known locus-specific DNA-protein interactions were identified at the ENO2 and GAL1 promoter regions of Saccharomyces cerevisiae, and revealed subgroups of proteins present in significantly different levels at the loci in cells grown on glucose versus galactose as the carbon source. Results were validated using chromatin immunoprecipitation. Overall, our analysis demonstrates that HyCCAPP is an effective and flexible technology that does not require specific antibodies nor prior knowledge of locally occurring DNA-protein interactions and can now be used to identify changes in protein interactions at target regions in the genome in response to physiological challenges.

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.

We searched a gene expression dataset comprised of 634 schizophrenia (SZ) cases and 713 controls for expression outliers (i.e., extreme tails of the distribution of transcript expression values) with SZ cases overrepresented compared with controls. These outlier genes were enriched for brain expression and for genes known to be associated with neurodevelopmental disorders. SZ cases showed higher outlier burden (i.e., total outlier events per subject) than controls for genes within copy number variants (CNVs) associated with SZ or neurodevelopmental disorders. Outlier genes were enriched for CNVs and for rare putative regulatory variants, but this only explained a small proportion of the outlier subjects, highlighting the underlying presence of additional genetic and potentially, epigenetic mechanisms.

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.

Identifying genes that contribute to white matter microstructure should provide insights into the neurobiological processes that regulate white matter development, plasticity and pathology. We detected five significant SNPs using genome-wide association analysis on a global measure of fractional anisotropy in 776 individuals from large extended pedigrees. Genetic correlations and genome-wide association results indicated that the genetic signal was largely homogeneous across white matter regions. Using RNA transcripts derived from lymphocytes in the same individuals, we identified two genes (GNA13 and CCDC91) that are likely to be cis-regulated by top SNPs, and whose expression levels were also genetically correlated with fractional anisotropy. A transcript of HTR7 was phenotypically associated with FA, and was associated with an intronic genome-wide significant SNP. These results encourage further research in the mechanisms by which GNA13, HTR7 and CCDC91 influence brain structure, and emphasize a role for g-protein signaling in the development and maintenance of white matter microstructure in health and disease.

Schizophrenia genome-wide association studies (GWAS) have identified common SNPs, rare copy number variants (CNVs) and a large polygenic contribution to illness risk, but biological mechanisms remain unclear. Bioinformatic analyses of significantly associated genetic variants point to a large role for regulatory variants. To identify gene expression abnormalities in schizophrenia, we generated whole-genome gene expression profiles using microarrays on lymphoblastoid cell lines (LCLs) from 413 cases and 446 controls. Regression analysis identified 95 transcripts differentially expressed by affection status at a genome-wide false discovery rate (FDR) of 0.05, while simultaneously controlling for confounding effects. These transcripts represented 89 genes with functions such as neurotransmission, gene regulation, cell cycle progression, differentiation, apoptosis, microRNA (miRNA) processing and immunity. This functional diversity is consistent with schizophrenia's likely significant pathophysiological heterogeneity. The overall enrichment of immune-related genes among those differentially expressed by affection status is consistent with hypothesized immune contributions to schizophrenia risk. The observed differential expression of extended major histocompatibility complex (xMHC) region histones (HIST1H2BD, HIST1H2BC, HIST1H2BH, HIST1H2BG and HIST1H4K) converges with the genetic evidence from GWAS, which find the xMHC to be the most significant susceptibility locus. Among the differentially expressed immune-related genes, B3GNT2 is implicated in autoimmune disorders previously tied to schizophrenia risk (rheumatoid arthritis and Graves' disease), and DICER1 is pivotal in miRNA processing potentially linking to miRNA alterations in schizophrenia (e. g. MIR137, the second strongest GWAS finding). Our analysis provides novel candidate genes for further study to assess their potential contribution to schizophrenia.

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.

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.

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 x 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. (C) 2012 Elsevier Ireland Ltd. All rights reserved.

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 1A, 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 that rs3019594 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 10(-75) and 10(-7)) 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.

Abstract Abstract 3366 Background/Aims: The hemophilia A (HA)-causing Factor (F)VIII gene (F8) mutation type is a well-established determinant of risk for the development of alloimmune inhibitors that neutralize replacement FVIII proteins in ∼20% of all HA patients. Studies have investigated variants of immune response genes to determine if they may account for the inter-individual variability in FVIII immunogenicity observed in patients with the same F8 abnormality, e.g., the intron-22 inversion. While some studies have found associations between inhibitor status and promoter polymorphisms in CTLA4, TNFA and/or IL10, others have not. If these promoter polymorphisms are indeed functional and truly influence inhibitor development, their alleles could modulate transcription initiation rates. The goal of this study was to investigate the possibility of cis-acting genotype-specific differences in mean steady-state mRNA levels encoded by CTLA4, TNFA and IL10. Methods: We examined the relationship of lymphocyte CTLA4, TNFA and IL10 mRNA levels with the genotypes of 265 SNPs located across their structural loci in 1189 Mexican American subjects in the San Antonio Family Heart Study. Expression profiles were generated using Illumina's HWG-6 BeadChips and genotypes came from the Illumina OmniExpress-12 BeadChip. Measured genotype association analyses that accounted for non-independence of family members and employed an additive model (in which testing to determine whether gene expression varies by genotype, with the model constrained so that each “dose” of the minor allele raises or lowers gene expression by an equal amount “beta”) were performed using the software package SOLAR. P-values were calculated using a 1 degree-of-freedom chi-square test comparing the likelihood of a model where the change in expression levels by genotype is estimated to the likelihood of a model where beta is constrained to zero. Results: None of the 265 genotyped SNPs within or near these three genes (i.e., 49 SNPs in IL10, 35 SNPs in CTLA4 & 181 SNPs in TNFA) function as cis-acting regulatory variants, as no significant genotype-specific associations with these genes' transcript levels were identified. Conclusions: We observed no evidence for cis-regulation of CTLA4, TNFA or IL10 in Mexican Americans, the largest and most rapidly growing minority population in the United States, despite having genotyped directly the previously implicated promoter polymorphisms in the current analysis (e.g., see Figure). Although Hispanic American HA patients were recently found to have a significantly higher risk for inhibitor development than White HA patients, it is possible that cis-acting functional variants in this minority population are rare and not well-represented by the common GWAS SNPs used for these analyses. Since linkage disequilibrium patterns between markers are population-specific, we are also currently genotyping these SNPs in a large cohort of African and Caucasian American HA patients. Disclosures: No relevant conflicts of interest to declare.