Cocktail-party listening and cognitive abilities show strong pleiotropy

Introduction

The cocktail-party problem refers to the difficulty listeners face when trying to attend to relevant sounds that are mixed with irrelevant ones. Previous studies have shown that solving these problems relies on perceptual as well as cognitive processes. Previously, we showed that speech-reception thresholds (SRTs) on a cocktail-party listening task were influenced by genetic factors. Here, we estimated the degree to which these genetic factors overlapped with those influencing cognitive abilities.

Methods

We measured SRTs and hearing thresholds (HTs) in 493 listeners, who ranged in age from 18 to 91 years old. The same individuals completed a cognitive test battery comprising 18 measures of various cognitive domains. Individuals belonged to large extended pedigrees, which allowed us to use variance component models to estimate the narrow-sense heritability of each trait, followed by phenotypic and genetic correlations between pairs of traits.

Results

All traits were heritable. The phenotypic and genetic correlations between SRTs and HTs were modest, and only the phenotypic correlation was significant. By contrast, all genetic SRT-cognition correlations were strong and significantly different from 0. For some of these genetic correlations, the hypothesis of complete pleiotropy could not be rejected.

Discussion

Overall, the results suggest that there was substantial genetic overlap between SRTs and a wide range of cognitive abilities, including abilities without a major auditory or verbal component. The findings highlight the important, yet sometimes overlooked, contribution of higher-order processes to solving the cocktail-party problem, raising an important caveat for future studies aiming to identify specific genetic factors that influence cocktail-party listening.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

The Genetic contribution to solving the cocktail-party problem

Communicating in everyday situations requires solving the cocktail-party problem, or segregating the acoustic mixture into its constituent sounds and attending to those of most interest. Humans show dramatic variation in this ability, leading some to experience real-world problems irrespective of whether they meet criteria for clinical hearing loss. Here, we estimated the genetic contribution to cocktail-party listening by measuring speech-reception thresholds (SRTs) in 425 people from large families and ranging in age from 18 to 91 years. Roughly half the variance of SRTs was explained by genes (h ² = 0.567). The genetic correlation between SRTs and hearing thresholds (HTs) was medium (ρ _G = 0.392), suggesting that the genetic factors influencing cocktail-party listening were partially distinct from those influencing sound sensitivity. Aging and socioeconomic status also strongly influenced SRTs. These findings may represent a first step toward identifying genes for "hidden hearing loss," or hearing problems in people with normal HTs.

Whole Genome Sequence Data From Captive Baboons Implicate RBFOX1 in Epileptic Seizure Risk

In this study, we investigate the genetic determinants that underlie epilepsy in a captive baboon pedigree and evaluate the potential suitability of this non-human primate model for understanding the genetic etiology of human epilepsy. Archived whole-genome sequence data were analyzed using both a candidate gene approach that targeted variants in baboon homologs of 19 genes (n = 20,881 SNPs) previously implicated in genetic generalized epilepsy (GGE) and a more agnostic approach that examined protein-altering mutations genome-wide as assessed by snpEff (n = 36,169). Measured genotype association tests for baboon cases of epileptic seizure were performed using SOLAR, as well as gene set enrichment analyses (GSEA) and protein–protein interaction (PPI) network construction of top association hits genome-wide (p < 0.01; n = 441 genes). The maximum likelihood estimate of heritability for epileptic seizure in the pedigreed baboon sample is 0.76 (SE = 0.77; p = 0.07). Among candidate genes for GGE, a significant association was detected for an intronic SNP in RBFOX1 (p = 5.92 × 10^–6; adjusted p = 0.016). For protein-altering variants, no genome-wide significant results were observed for epilepsy status. However, GSEA revealed significant positive enrichment for genes involved in the extracellular matrix structure (ECM; FDR = 0.0072) and collagen formation (FDR = 0.017), which was reflected in a major PPI network cluster. This preliminary study highlights the potential role of RBFOX1 in the epileptic baboon, a protein involved in transcriptomic regulation of multiple epilepsy candidate genes in humans and itself previously implicated in human epilepsy, both focal and generalized. Moreover, protein-damaging variants from across the genome exhibit a pattern of association that links collagen-containing ECM to epilepsy risk. These findings suggest a shared genetic etiology between baboon and human forms of GGE and lay the foundation for follow-up research.

Identifying the Lipidomic Effects of a Rare Loss-of-Function Deletion in ANGPTL3

BACKGROUND

The identification and understanding of therapeutic targets for atherosclerotic cardiovascular disease is of fundamental importance given its global health and economic burden. Inhibition of ANGPTL3 (angiopoietin-like 3) has demonstrated a cardioprotective effect, showing promise for atherosclerotic cardiovascular disease treatment, and is currently the focus of ongoing clinical trials. Here, we assessed the genetic basis of variation in ANGPTL3 levels in the San Antonio Family Heart Study.

METHODS

We assayed ANGPTL3 protein levels in ≈1000 Mexican Americans from extended pedigrees. By drawing upon existing plasma lipidome profiles and genomic data we conducted analyses to understand the genetic basis to variation in ANGPTL3 protein levels, and accordingly the correlation with the plasma lipidome.

RESULTS

In a variance components framework, we identified that variation in ANGPTL3 was significantly heritable (h²=0.33, P=1.31×10^-16). To explore the genetic basis of this heritability, we conducted a genome-wide linkage scan and identified significant linkage (logarithm of odds =6.18) to a locus on chromosome 1 at 90 centimorgans, corresponding to the ANGPTL3 gene location. In the genomes of 23 individuals from a single pedigree, we identified a loss-of-function variant, rs398122988 (N121Kfs*2), in ANGPTL3, that was significantly associated with lower ANGPTL3 levels (β=-1.69 SD units, P=3.367×10^-13), and accounted for the linkage signal at this locus. Given the known role of ANGPTL3 as an inhibitor of endothelial and lipoprotein lipase, we explored the association of ANGPTL3 protein levels and rs398122988 with the plasma lipidome and related phenotypes, identifying novel associations with phosphatidylinositols.

CONCLUSIONS

Variation in ANGPTL3 protein levels is heritable and under significant genetic control. Both ANGPTL3 levels and loss-of-function variants in ANGPTL3 have significant associations with the plasma lipidome. These findings further our understanding of ANGPTL3 as a therapeutic target for atherosclerotic cardiovascular disease.

Further evidence supporting a potential role for ADH1B in obesity

Insulin is an essential hormone that regulates glucose homeostasis and metabolism. Insulin resistance (IR) arises when tissues fail to respond to insulin, and it leads to serious health problems including Type 2 Diabetes (T2D). Obesity is a major contributor to the development of IR and T2D. We previously showed that gene expression of alcohol dehydrogenase 1B (ADH1B) was inversely correlated with obesity and IR in subcutaneous adipose tissue of Mexican Americans. In the current study, a meta-analysis of the relationship between ADH1B expression and BMI in Mexican Americans, African Americans, Europeans, and Pima Indians verified that BMI was increased with decreased ADH1B expression. Using established human subcutaneous pre-adipocyte cell lines derived from lean (BMI < 30 kg m-2) or obese (BMI ≥ 30 kg m-2) donors, we found that ADH1B protein expression increased substantially during differentiation, and overexpression of ADH1B inhibited fatty acid binding protein expression. Mature adipocytes from lean donors expressed ADH1B at higher levels than obese donors. Insulin further induced ADH1B protein expression as well as enzyme activity. Knockdown of ADH1B expression decreased insulin-stimulated glucose uptake. Our findings suggest that ADH1B is involved in the proper development and metabolic activity of adipose tissues and this function is suppressed by obesity.

Minimal Relationship between Local Gyrification and General Cognitive Ability in Humans

Previous studies suggest that gyrification is associated with superior cognitive abilities in humans, but the strength of this relationship remains unclear. Here, in two samples of related individuals (total N = 2882), we calculated an index of local gyrification (LGI) at thousands of cortical surface points using structural brain images and an index of general cognitive ability (g) using performance on cognitive tests. Replicating previous studies, we found that phenotypic and genetic LGI-g correlations were positive and statistically significant in many cortical regions. However, all LGI-g correlations in both samples were extremely weak, regardless of whether they were significant or nonsignificant. For example, the median phenotypic LGI-g correlation was 0.05 in one sample and 0.10 in the other. These correlations were even weaker after adjusting for confounding neuroanatomical variables (intracranial volume and local cortical surface area). Furthermore, when all LGIs were considered together, at least 89% of the phenotypic variance of g remained unaccounted for. We conclude that the association between LGI and g is too weak to have profound implications for our understanding of the neurobiology of intelligence. This study highlights potential issues when focusing heavily on statistical significance rather than effect sizes in large-scale observational neuroimaging studies.

Neurocognitive impairment in type 2 diabetes: evidence for shared genetic aetiology

AIMS/HYPOTHESIS

Type 2 diabetes is associated with cognitive impairments, but it is unclear whether common genetic factors influence both type 2 diabetes risk and cognition.

METHODS

Using data from 1892 Mexican-American individuals from extended pedigrees, including 402 with type 2 diabetes, we examined possible pleiotropy between type 2 diabetes and cognitive functioning, as measured by a comprehensive neuropsychological test battery.

RESULTS

Negative phenotypic correlations (ρp) were observed between type 2 diabetes and measures of attention (Continuous Performance Test [CPT d']: ρp = -0.143, p = 0.001), verbal memory (California Verbal Learning Test [CVLT] recall: ρp = -0.111, p = 0.004) and face memory (Penn Face Memory Test [PFMT]: ρp = -0.127, p = 0.002; PFMT Delayed: ρp = -0.148, p = 2 × 10-4), replicating findings of cognitive impairment in type 2 diabetes. Negative genetic correlations (ρg) were also observed between type 2 diabetes and measures of attention (CPT d': ρg = -0.401, p = 0.001), working memory (digit span backward test: ρg = -0.380, p = 0.005), and face memory (PFMT: ρg = -0.476, p = 2 × 10-4; PFMT Delayed: ρg = -0.376, p = 0.005), suggesting that the same genetic factors underlying risk for type 2 diabetes also influence poor cognitive performance in these domains. Performance in these domains was also associated with type 2 diabetes risk using an endophenotype ranking value approach. Specifically, on measures of attention (CPT d': β = -0.219, p = 0.005), working memory (digit span backward: β = -0.326, p = 0.035), and face memory (PFMT: β = -0.171, p = 0.023; PFMT Delayed: β = -0.215, p = 0.005), individuals with type 2 diabetes showed the lowest performance, while unaffected/unrelated individuals showed the highest performance, and those related to an individual with type 2 diabetes performed at an intermediate level.

CONCLUSIONS/INTERPRETATION

These findings suggest that cognitive impairment may be a useful endophenotype of type 2 diabetes and, therefore, help to elucidate the pathophysiological underpinnings of this chronic disease.

DATA AVAILABILITY

The data analysed in this study is available in dbGaP: www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs001215.v2.p2.

Genotype phasing in pedigrees using whole-genome sequence data

Phasing is the process of inferring haplotypes from genotype data. Efficient algorithms and associated software for accurate phasing in pedigrees are needed, especially for populations lacking reference panels of sequenced individuals. We present a novel method for phasing genotypes from whole-genome sequence data in pedigrees, called PULSAR (Phasing Using Lineage Specific Alleles/Rare variants). The method is based on the property that alleles specific to a single founding chromosome within a pedigree are highly informative for identifying haplotypes that are shared identical by descent. Simulation studies are used to assess the performance of PULSAR with various pedigree sizes and structures, and the effect of genotyping errors and the presence of nonsequenced individuals is investigated. In pedigrees with complete sequencing and realistic genotyping error rates, PULSAR correctly phases >99.9% of heterozygous genotypes, excluding sites at which all individuals are heterozygous, and does so with a switch error rate frequently below 10^-4. PULSAR is highly accurate, capable of genotype error correction and imputation, and computationally competitive with alternative phasing software applicable to pedigrees. Our method has the significant advantage of not requiring reference panels that are essential for other population-based phasing algorithms. A software implementation of PULSAR is freely available.

TRAK2, a novel regulator of ABCA1 expression, cholesterol efflux and HDL biogenesis

Aims

The recent failures of HDL-raising therapies have underscored our incomplete understanding of HDL biology. Therefore there is an urgent need to comprehensively investigate HDL metabolism to enable the development of effective HDL-centric therapies. To identify novel regulators of HDL metabolism, we performed a joint analysis of human genetic, transcriptomic, and plasma HDL-cholesterol (HDL-C) concentration data and identified a novel association between trafficking protein, kinesin binding 2 (TRAK2) and HDL-C concentration. Here we characterize the molecular basis of the novel association between TRAK2 and HDL-cholesterol concentration.

Methods and results

Analysis of lymphocyte transcriptomic data together with plasma HDL from the San Antonio Family Heart Study (n = 1240) revealed a significant negative correlation between TRAK2 mRNA levels and HDL-C concentration, HDL particle diameter and HDL subspecies heterogeneity. TRAK2 siRNA-mediated knockdown significantly increased cholesterol efflux to apolipoprotein A-I and isolated HDL from human macrophage (THP-1) and liver (HepG2) cells by increasing the mRNA and protein expression of the cholesterol transporter ATP-binding cassette, sub-family A member 1 (ABCA1). The effect of TRAK2 knockdown on cholesterol efflux was abolished in the absence of ABCA1, indicating that TRAK2 functions in an ABCA1-dependent efflux pathway. TRAK2 knockdown significantly increased liver X receptor (LXR) binding at the ABCA1 promoter, establishing TRAK2 as a regulator of LXR-mediated transcription of ABCA1.

Conclusion

We show, for the first time, that TRAK2 is a novel regulator of LXR-mediated ABCA1 expression, cholesterol efflux, and HDL biogenesis. TRAK2 may therefore be an important target in the development of anti-atherosclerotic therapies.

Family-based analyses reveal novel genetic overlap between cytokine interleukin-8 and risk for suicide attempt

BACKGROUND

Suicide is major public health concern. It is imperative to find robust biomarkers so that at-risk individuals can be identified in a timely and reliable manner. Previous work suggests mechanistic links between increased cytokines and risk for suicide, but questions remain regarding the etiology of this association, as well as the roles of sex and BMI.

METHODS

Analyses were conducted using a randomly-ascertained extended-pedigree sample of 1882 Mexican-American individuals (60% female, mean age = 42.04, range = 18-97). Genetic correlations were calculated using a variance components approach between the cytokines TNF-α, IL-6 and IL-8, and Lifetime Suicide Attempt and Current Suicidal Ideation. The potentially confounding effects of sex and BMI were considered.

RESULTS

159 individuals endorse a Lifetime Suicide Attempt. IL-8 and IL-6 shared significant genetic overlap with risk for suicide attempt (ρg = 0.49, pFDR = 7.67 × 10-03; ρg = 0.53, pFDR = 0.01), but for IL-6 this was attenuated when BMI was included as a covariate (ρg = 0.37, se = 0.23, pFDR = 0.12). Suicide attempts were significantly more common in females (pFDR = 0.01) and the genetic overlap between IL-8 and risk for suicide attempt was significant in females (ρg = 0.56, pFDR = 0.01), but not in males (ρg = 0.44, pFDR = 0.30).

DISCUSSION

These results demonstrate that: IL-8 shares genetic influences with risk for suicide attempt; females drove this effect; and BMI should be considered when assessing the association between IL-6 and suicide. This finding represents a significant advancement in knowledge by demonstrating that cytokine alterations are not simply a secondary manifestation of suicidal behavior, but rather, the pathophysiology of suicide attempts is, at least partly, underpinned by the same biological mechanisms responsible for regulating inflammatory response.

Disentangling the genetic overlap between cholesterol and suicide risk

Suicide is major public health concern; one million individuals worldwide die by suicide each year of which there are many more attempts. Thus, it is imperative that robust and reliable indicators, or biomarkers, of suicide risk be identified so that individuals at risk can be identified and provided appropriate interventions as quickly as possible. Previous work has revealed a relationship between low levels of circulating cholesterol and suicide risk, implicating cholesterol level as one such potential biomarker, but the factors underlying this relationship remain unknown. In the present study, we applied a combination of bivariate polygenic and coefficient-of-relatedness analysis, followed by mediation analysis, in a large sample of Mexican-American individuals from extended pedigrees [N = 1897; 96 pedigrees (average size = 19.17 individuals, range = 2-189) 60% female; mean age = 42.58 years, range = 18-97 years, sd = 15.75 years] with no exclusion criteria for any given psychiatric disorder. We observed that total esterified cholesterol measured at the time of psychiatric assessment shared a significant genetic overlap with risk for suicide attempt (ρ_g = -0.64, p = 1.24 × 10^-04). We also found that total unesterified cholesterol measured around 20 years prior to assessment varied as a function of genetic proximity to an affected individual (h² = 0.21, se = 0.10, p = 8.73 × 10^-04; β_suicide = -0.70, se = 0.25, p = 8.90 × 10^-03). Finally, we found that the relationship between total unesterified cholesterol and suicide risk was significantly mediated by ABCA-1-specific cholesterol efflux capacity (β_{suicide-efflux} = -0.45, p = 0.039; β_{efflux-cholexterol} = -0.34, p < 0.0001; β_indirect = -0.15, p = 0.044). These findings suggest that the relatively well-delineated process of cholesterol metabolism and associated molecular pathways will be informative for understanding the neurobiological underpinnings of risk for suicide attempt.

Contribution of Inbred Singletons to Variance Component Estimation of Heritability and Linkage

OBJECTIVES

An interesting consequence of consanguinity is that the inbred singleton becomes informative for genetic variance. We determine the contribution of an inbred singleton to variance component analysis of heritability and linkage.

METHODS

Statistical theory for the power of variance component analysis of quantitative traits is used to determine the expected contribution of an inbred singleton to likelihood-ratio tests of heritability and linkage.

RESULTS

In variance component models, an inbred singleton contributes relatively little to a test of heritability but can contribute substantively to a test of linkage. For small-to-moderate quantitative trait locus (QTL) effects and a level of inbreeding comparable to matings between first cousins (the preferred form of union in many human populations), an inbred singleton can carry nearly 25% of the information of a non-inbred sib pair. In more highly inbred contexts available with experimental animal populations, nonhuman primate colonies, and some human subpopulations, the contribution of an inbred singleton relative to a sib pair can exceed 50%.

CONCLUSIONS

Inbred individuals, even in isolation from other members of a sample, can contribute to variance component estimation and tests of heritability and linkage. Under certain conditions, the informativeness of the inbred singleton can approach that of a non-inbred sib pair.

Transcriptomic signatures of schizophrenia revealed by dopamine perturbation in an ex vivo model

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.

Identity-by-Descent Mapping Identifies Major Locus for Serum Triglycerides in Amerindians Largely Explained by an APOC3 Founder Mutation

BACKGROUND

Identity-by-descent mapping using empirical estimates of identity-by-descent allele sharing may be useful for studies of complex traits in founder populations, where hidden relationships may augment the inherent genetic information that can be used for localization.

METHODS AND RESULTS

Through identity-by-descent mapping, using ≈400 000 single-nucleotide polymorphisms (SNPs), of serum lipid profiles, we identified a major linkage signal for triglycerides in 1007 Pima Indians (LOD=9.23; P=3.5×10^-11 on chromosome 11q). In subsequent fine-mapping and replication association studies in ≈7500 Amerindians, we determined that this signal reflects effects of a loss-of-function Ala43Thr substitution in APOC3 (rs147210663) and 3 established functional SNPs in APOA5. The association with rs147210663 was particularly strong; each copy of the Thr allele conferred 42% lower triglycerides (β=-0.92±0.059 SD unit; P=9.6×10^-55 in 4668 Pimas and 2793 Southwest Amerindians combined). The Thr allele is extremely rare in most global populations but has a frequency of 2.5% in Pimas. We further demonstrated that 3 APOA5 SNPs with established functional impact could explain the association with the most well-replicated SNP (rs964184) for triglycerides identified by genome-wide association studies. Collectively, these 4 SNPs account for 6.9% of variation in triglycerides in Pimas (and 4.1% in Southwest Amerindians), and their inclusion in the original linkage model reduced the linkage signal to virtually null.

CONCLUSIONS

APOC3/APOA5 constitutes a major locus for serum triglycerides in Amerindians, especially the Pimas, and these results provide an empirical example for the concept that population-based linkage analysis is a useful strategy to identify complex trait variants.

Shared Genetic Factors Influence Head Motion During MRI and Body Mass Index

Head movements are typically viewed as a nuisance to functional magnetic resonance imaging (fMRI) analysis, and are particularly problematic for resting state fMRI. However, there is growing evidence that head motion is a behavioral trait with neural and genetic underpinnings. Using data from a large randomly ascertained extended pedigree sample of Mexican Americans (n = 689), we modeled the genetic structure of head motion during resting state fMRI and its relation to 48 other demographic and behavioral phenotypes. A replication analysis was performed using data from the Human Connectome Project, which uses an extended twin design (n = 864). In both samples, head motion was significantly heritable (h2 = 0.313 and 0.427, respectively), and phenotypically correlated with numerous traits. The most strongly replicated relationship was between head motion and body mass index, which showed evidence of shared genetic influences in both data sets. These results highlight the need to view head motion in fMRI as a complex neurobehavioral trait correlated with a number of other demographic and behavioral phenotypes. Given this, when examining individual differences in functional connectivity, the confounding of head motion with other traits of interest needs to be taken into consideration alongside the critical important of addressing head motion artifacts.

Exome sequences of multiplex, multigenerational families reveal schizophrenia risk loci with potential implications for neurocognitive performance

Schizophrenia is a serious mental illness, involving disruptions in thought and behavior, with a worldwide prevalence of about one percent. Although highly heritable, much of the genetic liability of schizophrenia is yet to be explained. We searched for susceptibility loci in multiplex, multigenerational families affected by schizophrenia, targeting protein-altering variation with in silico predicted functional effects. Exome sequencing was performed on 136 samples from eight European-American families, including 23 individuals diagnosed with schizophrenia or schizoaffective disorder. In total, 11,878 non-synonymous variants from 6,396 genes were tested for their association with schizophrenia spectrum disorders. Pathway enrichment analyses were conducted on gene-based test results, protein-protein interaction (PPI) networks, and epistatic effects. Using a significance threshold of FDR < 0.1, association was detected for rs10941112 (p = 2.1 × 10^-5 ; q-value = 0.073) in AMACR, a gene involved in fatty acid metabolism and previously implicated in schizophrenia, with significant cis effects on gene expression (p = 5.5 × 10^-4 ), including brain tissue data from the Genotype-Tissue Expression project (minimum p = 6.0 × 10^-5 ). A second SNP, rs10378 located in TMEM176A, also shows risk effects in the exome data (p = 2.8 × 10^-5 ; q-value = 0.073). PPIs among our top gene-based association results (p < 0.05; n = 359 genes) reveal significant enrichment of genes involved in NCAM-mediated neurite outgrowth (p = 3.0 × 10^-5 ), while exome-wide SNP-SNP interaction effects for rs10941112 and rs10378 indicate a potential role for kinase-mediated signaling involved in memory and learning. In conclusion, these association results implicate AMACR and TMEM176A in schizophrenia risk, whose effects may be modulated by genes involved in synaptic plasticity and neurocognitive performance.

Association of TMTC2 With Human Nonsyndromic Sensorineural Hearing Loss

IMPORTANCE

Sensorineural hearing loss (SNHL) is commonly caused by conditions that affect cochlear structures or the auditory nerve, and the genes identified as causing SNHL to date only explain a fraction of the overall genetic risk for this debilitating disorder. It is likely that other genes and mutations also cause SNHL.

OBJECTIVE

To identify a candidate gene that causes bilateral, symmetric, progressive SNHL in a large multigeneration family of Northern European descent.

DESIGN, SETTING, AND PARTICIPANTS

In this prospective genotype and phenotype study performed from January 1, 2006, through April 1, 2016, a 6-generation family of Northern European descent with 19 individuals having reported early-onset hearing loss suggestive of an autosomal dominant inheritance were studied at a tertiary academic medical center. In addition, 179 unrelated adult individuals with SNHL and 186 adult individuals reporting nondeafness were examined.

MAIN OUTCOMES AND MEASURES

Sensorineural hearing loss.

RESULTS

Nine family members (5 women [55.6%]) provided clinical audiometric and medical records that documented hearing loss. The hearing loss is characterized as bilateral, symmetric, progressive SNHL that reached severe to profound loss in childhood. Audiometric configurations demonstrated a characteristic dip at 1000 to 2000 Hz. All affected family members wear hearing aids or have undergone cochlear implantation. Exome sequencing and linkage and association analyses identified a fully penetrant sequence variant (rs35725509) on chromosome 12q21 (logarithm of odds, 3.3) in the TMTC2 gene region that segregates with SNHL in this family. This gene explains the SNHL occurrence in this family. The variant is also associated with SNHL in a cohort of 363 unrelated individuals (179 patients with confirmed SNHL and 184 controls, P = 7 × 10-4).

CONCLUSIONS AND RELEVANCE

A previously uncharacterized gene, TMTC2, has been identified as a candidate for causing progressive SNHL in humans. This finding identifies a novel locus that causes autosomal dominant SNHL and therefore a more detailed understanding of the genetic basis of SNHL. Because TMTC2 has not been previously reported to regulate auditory function, the discovery reveals a potentially new, uncharacterized mechanism of hearing loss.

Transcriptome sequencing study implicates immune-related genes differentially expressed in schizophrenia: new data and a meta-analysis

We undertook an RNA sequencing (RNAseq)-based transcriptomic profiling study on lymphoblastoid cell lines of a European ancestry sample of 529 schizophrenia cases and 660 controls, and found 1058 genes to be differentially expressed by affection status. These differentially expressed genes were enriched for involvement in immunity, especially the 697 genes with higher expression in cases. Comparing the current RNAseq transcriptomic profiling to our previous findings in an array-based study of 268 schizophrenia cases and 446 controls showed a highly significant positive correlation over all genes. Fifteen (18%) of the 84 genes with significant (false discovery rate<0.05) expression differences between cases and controls in the previous study and analyzed here again were differentially expressed by affection status here at a genome-wide significance level (Bonferroni P<0.05 adjusted for 8141 analyzed genes in total, or P<~6.1 × 10^-6), all with the same direction of effect, thus providing corroborative evidence despite each sample of fully independent subjects being studied by different technological approaches. Meta-analysis of the RNAseq and array data sets (797 cases and 1106 controls) showed 169 additional genes (besides those found in the primary RNAseq-based analysis) to be differentially expressed, and provided further evidence of immune gene enrichment. In addition to strengthening our previous array-based gene expression differences in schizophrenia cases versus controls and providing transcriptomic support for some genes implicated by other approaches for schizophrenia, our study detected new genes differentially expressed in schizophrenia. We highlight RNAseq-based differential expression of various genes involved in neurodevelopment and/or neuronal function, and discuss caveats of the approach.

The lipidome in major depressive disorder: Shared genetic influence for ether-phosphatidylcholines, a plasma-based phenotype related to inflammation, and disease risk

BACKGROUND

The lipidome is rapidly garnering interest in the field of psychiatry. Recent studies have implicated lipidomic changes across numerous psychiatric disorders. In particular, there is growing evidence that the concentrations of several classes of lipids are altered in those diagnosed with MDD. However, for lipidomic abnormalities to be considered potential treatment targets for MDD (rather than secondary manifestations of the disease), a shared etiology between lipid concentrations and MDD should be demonstrated.

METHODS

In a sample of 567 individuals from 37 extended pedigrees (average size 13.57 people, range=3-80), we used mass spectrometry lipidomic measures to evaluate the genetic overlap between twenty-three biologically distinct lipid classes and a dimensional scale of MDD.

RESULTS

We found that the lipid class with the largest endophenotype ranking value (ERV, a standardized parametric measure of pleiotropy) were ether-phosphodatidylcholines (alkylphosphatidylcholine, PC(O) and alkenylphosphatidylcholine, PC(P) subclasses). Furthermore, we examined the cluster structure of the twenty-five species within the top-ranked lipid class, and the relationship of those clusters with MDD. This analysis revealed that species containing arachidonic acid generally exhibited the greatest degree of genetic overlap with MDD.

CONCLUSIONS

This study is the first to demonstrate a shared genetic etiology between MDD and ether-phosphatidylcholine species containing arachidonic acid, an omega-6 fatty acid that is a precursor to inflammatory mediators, such as prostaglandins. The study highlights the potential utility of the well-characterized linoleic/arachidonic acid inflammation pathway as a diagnostic marker and/or treatment target for MDD.

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

The genetic basis of the comorbidity between cannabis use and major depression

BACKGROUND AND AIMS

While the prevalence of major depression is elevated among cannabis users, the role of genetics in this pattern of comorbidity is not clear. This study aimed to estimate the heritability of cannabis use and major depression, quantify the genetic overlap between these two traits and localize regions of the genome that segregate in families with cannabis use and major depression.

DESIGN

Family-based univariate and bivariate genetic analysis.

SETTING

San Antonio, Texas, USA.

PARTICIPANTS

Genetics of Brain Structure and Function study (GOBS) participants: 1284 Mexican Americans from 75 large multi-generation families and an additional 57 genetically unrelated spouses.

MEASUREMENTS

Phenotypes of life-time history of cannabis use and major depression, measured using the semistructured MINI-Plus interview. Genotypes measured using ~1 M single nucleotide polymorphisms (SNPs) on Illumina BeadChips. A subselection of these SNPs were used to build multi-point identity-by-descent matrices for linkage analysis.

FINDINGS

Both cannabis use [h²  = 0.614, P = 1.00 × 10^-6 , standard error (SE) = 0.151] and major depression (h²  = 0.349, P = 1.06 × 10^-5 , SE = 0.100) are heritable traits, and there is significant genetic correlation between the two (ρ_g  = 0.424, P = 0.0364, SE = 0.195). Genome-wide linkage scans identify a significant univariate linkage peak for major depression on chromosome 22 [logarithm of the odds (LOD) = 3.144 at 2 centimorgans (cM)], with a suggestive peak for cannabis use on chromosome 21 (LOD = 2.123 at 37 cM). A significant pleiotropic linkage peak influencing both cannabis use and major depression was identified on chromosome 11 using a bivariate model (LOD = 3.229 at 112 cM). Follow-up of this pleiotropic signal identified a SNP 20 kb upstream of NCAM1 (rs7932341) that shows significant bivariate association (P = 3.10 × 10^-5 ). However, this SNP is rare (seven minor allele carriers) and does not drive the linkage signal observed.

CONCLUSIONS

There appears to be a significant genetic overlap between cannabis use and major depression among Mexican Americans, a pleiotropy that appears to be localized to a region on chromosome 11q23 that has been linked previously to these phenotypes.

Novel genetic loci underlying human intracranial volume identified through genome-wide association

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Genome-wide significant loci for addiction and anxiety

BACKGROUND

Psychiatric comorbidity is common among individuals with addictive disorders, with patients frequently suffering from anxiety disorders. While the genetic architecture of comorbid addictive and anxiety disorders remains unclear, elucidating the genes involved could provide important insights into the underlying etiology.

METHODS

Here we examine a sample of 1284 Mexican-Americans from randomly selected extended pedigrees. Variance decomposition methods were used to examine the role of genetics in addiction phenotypes (lifetime history of alcohol dependence, drug dependence or chronic smoking) and various forms of clinically relevant anxiety. Genome-wide univariate and bivariate linkage scans were conducted to localize the chromosomal regions influencing these traits.

RESULTS

Addiction phenotypes and anxiety were shown to be heritable and univariate genome-wide linkage scans revealed significant quantitative trait loci for drug dependence (14q13.2-q21.2, LOD=3.322) and a broad anxiety phenotype (12q24.32-q24.33, LOD=2.918). Significant positive genetic correlations were observed between anxiety and each of the addiction subtypes (ρg=0.550-0.655) and further investigation with bivariate linkage analyses identified significant pleiotropic signals for alcohol dependence-anxiety (9q33.1-q33.2, LOD=3.054) and drug dependence-anxiety (18p11.23-p11.22, LOD=3.425).

CONCLUSIONS

This study confirms the shared genetic underpinnings of addiction and anxiety and identifies genomic loci involved in the etiology of these comorbid disorders. The linkage signal for anxiety on 12q24 spans the location of TMEM132D, an emerging gene of interest from previous GWAS of anxiety traits, whilst the bivariate linkage signal identified for anxiety-alcohol on 9q33 peak coincides with a region where rare CNVs have been associated with psychiatric disorders. Other signals identified implicate novel regions of the genome in addiction genetics.

GWAS and transcriptional analysis prioritize ITPR1 and CNTN4 for a serum uric acid 3p26 QTL in Mexican Americans

Background

The variation in serum uric acid concentrations is under significant genetic influence. Elevated SUA concentrations have been linked to increased risk for gout, kidney stones, chronic kidney disease, and cardiovascular disease whereas reduced serum uric acid concentrations have been linked to multiple sclerosis, Parkinson’s disease and Alzheimer’s disease. Previously, we identified a novel locus on chromosome 3p26 affecting serum uric acid concentrations in Mexican Americans from San Antonio Family Heart Study. As a follow up, we examined genome-wide single nucleotide polymorphism data in an extended cohort of 1281 Mexican Americans from multigenerational families of the San Antonio Family Heart Study and the San Antonio Family Diabetes/Gallbladder Study. We used a linear regression-based joint linkage/association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component.

Results

Univariate genetic analysis indicated serum uric acid concentrations to be significant heritable (h² = 0.50 ± 0.05, p < 4 × 10⁻³⁵), and linkage analysis of serum uric acid concentrations confirmed our previous finding of a novel locus on 3p26 (LOD = 4.9, p < 1 × 10⁻⁵) in the extended sample. Additionally, we observed strong association of serum uric acid concentrations with variants in following candidate genes in the 3p26 region; inositol 1,4,5-trisphosphate receptor, type 1 (ITPR1), contactin 4 (CNTN4), decapping mRNA 1A (DCP1A); transglutaminase 4 (TGM4) and rho guanine nucleotide exchange factor (GEF) 26 (ARHGEF26) [p < 3 × 10⁻⁷; minor allele frequencies ranged between 0.003 and 0.42] and evidence of cis-regulation for ITPR1 transcripts.

Conclusion

Our results confirm the importance of the chromosome 3p26 locus and genetic variants in this region in the regulation of serum uric acid concentrations.

Genome- and epigenome-wide association study of hypertriglyceridemic waist in Mexican American families

BACKGROUND

There is growing interest in the hypertriglyceridemic waist (HTGW) phenotype, defined as high waist circumference (≥95 cm in males and ≥80 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) as a marker of type 2 diabetes (T2D) and cardiovascular disease. However, the prevalence of this phenotype in high-risk populations, its association with T2D, and the genetic or epigenetic influences on HTGW are not well explored. Using data from large, extended families of Mexican Americans (a high-risk minority population in the USA) we aimed to: (1) estimate the prevalence of this phenotype, (2) test its association with T2D and related traits, and (3) dissect out the genetic and epigenetic associations with this phenotype using genome-wide and epigenome-wide studies, respectively.

RESULTS

Data for this study was from 850 Mexican American participants (representing 39 families) recruited under the ongoing San Antonio Family Heart Study, 26 % of these individuals had HTGW. This phenotype was significantly heritable (h (2) r = 0.52, p = 1.1 × 10(-5)) and independently associated with T2D as well as fasting glucose levels and insulin resistance. We conducted genome-wide association analyses using 759,809 single nucleotide polymorphisms (SNPs) and epigenome-wide association analyses using 457,331 CpG sites. There was no evidence of any SNP associated with HTGW at the genome-wide level but two CpG sites (cg00574958 and cg17058475) in CPT1A and one CpG site (cg06500161) in ABCG1 were significantly associated with HTGW and remained significant after adjusting for the closely related components of metabolic syndrome. CPT1A holds a cardinal position in the metabolism of long-chain fatty acids while ABCG1 plays a role in triglyceride metabolism.

CONCLUSIONS

Our results reemphasize the value of HTGW as a marker of T2D. This phenotype shows association with DNA methylation within CPT1A and ABCG1, genes involved in fatty acid and triglyceride metabolism. Our results underscore the importance of epigenetics in a clinically informative phenotype.

The transcriptional landscape of age in human peripheral blood

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the ‘transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Ageing increases the risk of many diseases. Here the authors compare blood cell transcriptomes of over 14,000 individuals and identify a set of about 1,500 genes that are differently expressed with age, shedding light on transcriptional programs linked to the ageing process and age-associated diseases.

Recurrent major depression and right hippocampal volume: A bivariate linkage and association study

Previous work has shown that the hippocampus is smaller in the brains of individuals suffering from major depressive disorder (MDD) than those of healthy controls. Moreover, right hippocampal volume specifically has been found to predict the probability of subsequent depressive episodes. This study explored the utility of right hippocampal volume as an endophenotype of recurrent MDD (rMDD). We observed a significant genetic correlation between the two traits in a large sample of Mexican American individuals from extended pedigrees (ρg = -0.34, p = 0.013). A bivariate linkage scan revealed a significant pleiotropic quantitative trait locus on chromosome 18p11.31-32 (LOD = 3.61). Bivariate association analysis conducted under the linkage peak revealed a variant (rs574972) within an intron of the gene SMCHD1 meeting the corrected significance level (χ(2) = 19.0, p = 7.4 × 10(-5)). Univariate association analyses of each phenotype separately revealed that the same variant was significant for right hippocampal volume alone, and also revealed a suggestively significant variant (rs12455524) within the gene DLGAP1 for rMDD alone. The results implicate right-hemisphere hippocampal volume as a possible endophenotype of rMDD, and in so doing highlight a potential gene of interest for rMDD risk.

Novel epigenetic determinants of type 2 diabetes in Mexican-American families

Although DNA methylation is now recognized as an important mediator of complex diseases, the extent to which the genetic basis of such diseases is accounted for by DNA methylation is unknown. In the setting of large, extended families representing a minority, high-risk population of the USA, we aimed to characterize the role of epigenome-wide DNA methylation in type 2 diabetes (T2D). Using Illumina HumanMethylation450 BeadChip arrays, we tested for association of DNA methylation at 446 356 sites with age, sex and phenotypic traits related to T2D in 850 pedigreed Mexican-American individuals. Robust statistical analyses showed that (i) 15% of the methylome is significantly heritable, with a median heritability of 0.14; (ii) DNA methylation at 14% of CpG sites is associated with nearby sequence variants; (iii) 22% and 3% of the autosomal CpG sites are associated with age and sex, respectively; (iv) 53 CpG sites were significantly associated with liability to T2D, fasting blood glucose and insulin resistance; (v) DNA methylation levels at five CpG sites, mapping to three well-characterized genes (TXNIP, ABCG1 and SAMD12) independently explained 7.8% of the heritability of T2D (vi) methylation at these five sites was unlikely to be influenced by neighboring DNA sequence variation. Our study has identified novel epigenetic indicators of T2D risk in Mexican Americans who have increased risk for this disease. These results provide new insights into potential treatment targets of T2D.

Transcriptome outlier analysis implicates schizophrenia susceptibility genes and enriches putatively functional rare genetic variants

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.

Transcriptomic identification of ADH1B as a novel candidate gene for obesity and insulin resistance in human adipose tissue in Mexican Americans from the Veterans Administration Genetic Epidemiology Study (VAGES)

Type 2 diabetes (T2D) is a complex metabolic disease that is more prevalent in ethnic groups such as Mexican Americans, and is strongly associated with the risk factors obesity and insulin resistance. The goal of this study was to perform whole genome gene expression profiling in adipose tissue to detect common patterns of gene regulation associated with obesity and insulin resistance. We used phenotypic and genotypic data from 308 Mexican American participants from the Veterans Administration Genetic Epidemiology Study (VAGES). Basal fasting RNA was extracted from adipose tissue biopsies from a subset of 75 unrelated individuals, and gene expression data generated on the Illumina BeadArray platform. The number of gene probes with significant expression above baseline was approximately 31,000. We performed multiple regression analysis of all probes with 15 metabolic traits. Adipose tissue had 3,012 genes significantly associated with the traits of interest (false discovery rate, FDR ≤ 0.05). The significance of gene expression changes was used to select 52 genes with significant (FDR ≤ 10(-4)) gene expression changes across multiple traits. Gene sets/Pathways analysis identified one gene, alcohol dehydrogenase 1B (ADH1B) that was significantly enriched (P < 10(-60)) as a prime candidate for involvement in multiple relevant metabolic pathways. Illumina BeadChip derived ADH1B expression data was consistent with quantitative real time PCR data. We observed significant inverse correlations with waist circumference (2.8 x 10(-9)), BMI (5.4 x 10(-6)), and fasting plasma insulin (P < 0.001). These findings are consistent with a central role for ADH1B in obesity and insulin resistance and provide evidence for a novel genetic regulatory mechanism for human metabolic diseases related to these traits.

Genome-wide genetic investigation of serological measures of common infections

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.

Genome-wide genetic and transcriptomic investigation of variation in antibody response to dietary antigens

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 (h(2) ), 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.

Pleiotropic locus for emotion recognition and amygdala volume identified using univariate and bivariate linkage

OBJECTIVE

The role of the amygdala in emotion recognition is well established, and amygdala volume and emotion recognition performance have each been shown separately to be highly heritable traits, but the potential role of common genetic influences on both traits has not been explored. The authors investigated the pleiotropic influences of amygdala volume and emotion recognition performance.

METHOD

In a sample of randomly selected extended pedigrees (N=858), the authors used a combination of univariate and bivariate linkage to investigate pleiotropy between amygdala volume and emotion recognition performance and followed up with association analysis.

RESULTS

The authors found a pleiotropic region for amygdala volume and emotion recognition performance on chromosome 4q26 (LOD score=4.40). Association analysis conducted in the region underlying the bivariate linkage peak revealed a variant meeting the corrected significance level (Bonferroni-corrected p=5.01×10(-5)) within an intron of PDE5A (rs2622497, p=4.4×10(-5)) as being jointly influential on both traits. PDE5A has been implicated previously in recognition-memory deficits and is expressed in subcortical structures that are thought to underlie memory ability, including the amygdala.

CONCLUSIONS

This study extends our understanding of the shared etiology between the amygdala and emotion recognition by showing that the overlap between amygdala volume and emotion recognition performance is due at least in part to common genetic influences. Moreover, this study identifies a pleiotropic locus for the two traits and an associated variant, which localizes the genetic signal even more precisely. These results, when taken in the context of previous research, highlight the potential utility of PDE5 inhibitors for ameliorating emotion recognition deficits in individuals suffering from mental or neurodegenerative illness.

Common genetic variants influence human subcortical brain structures

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

Discovering schizophrenia endophenotypes in randomly ascertained pedigrees

BACKGROUND

Although case-control approaches are beginning to disentangle schizophrenia's complex polygenic burden, other methods will likely be necessary to fully identify and characterize risk genes. Endophenotypes, traits genetically correlated with an illness, can help characterize the impact of risk genes by providing genetically relevant traits that are more tractable than the behavioral symptoms that classify mental illness. Here, we present an analytic approach for discovering and empirically validating endophenotypes in extended pedigrees with very few affected individuals. Our approach indexes each family member's risk as a function of shared genetic kinship with an affected individual, often referred to as the coefficient of relatedness. To demonstrate the utility of this approach, we search for neurocognitive and neuroanatomic endophenotypes for schizophrenia in large unselected multigenerational pedigrees.

METHODS

A fixed-effects test within the variance component framework was performed on neurocognitive and cortical surface area traits in 1606 Mexican-American individuals from large, randomly ascertained extended pedigrees who participated in the Genetics of Brain Structure and Function study. As affecteds were excluded from analyses, results were not influenced by disease state or medication usage.

RESULTS

Despite having sampled just 6 individuals with schizophrenia, our sample provided 233 individuals at various levels of genetic risk for the disorder. We identified three neurocognitive measures (digit-symbol substitution, facial memory, and emotion recognition) and six medial temporal and prefrontal cortical surfaces associated with liability for schizophrenia.

CONCLUSIONS

With our novel analytic approach, one can discover and rank endophenotypes for schizophrenia, or any heritable disease, in randomly ascertained pedigrees.

Influence of age, sex and genetic factors on the human brain

We report effects of age, age(2), sex and additive genetic factors on variability in gray matter thickness, surface area and white matter integrity in 1,010 subjects from the Genetics of Brain Structure and Function Study. Age was more strongly associated with gray matter thickness and fractional anisotropy of water diffusion in white matter tracts, while sex was more strongly associated with gray matter surface area. Widespread heritability of neuroanatomic traits was observed, suggesting that brain structure is under strong genetic control. Furthermore, our findings indicate that neuroimaging-based measurements of cerebral variability are sensitive to genetic mediation. Fundamental studies of genetic influence on the brain will help inform gene discovery initiatives in both clinical and normative samples.

Whole-genome sequencing to understand the genetic architecture of common gene expression and biomarker phenotypes

Initial results from sequencing studies suggest that there are relatively few low-frequency (<5%) variants associated with large effects on common phenotypes. We performed low-pass whole-genome sequencing in 680 individuals from the InCHIANTI study to test two primary hypotheses: (i) that sequencing would detect single low-frequency–large effect variants that explained similar amounts of phenotypic variance as single common variants, and (ii) that some common variant associations could be explained by low-frequency variants. We tested two sets of disease-related common phenotypes for which we had statistical power to detect large numbers of common variant–common phenotype associations—11 132 cis-gene expression traits in 450 individuals and 93 circulating biomarkers in all 680 individuals. From a total of 11 657 229 high-quality variants of which 6 129 221 and 5 528 008 were common and low frequency (<5%), respectively, low frequency–large effect associations comprised 7% of detectable cis-gene expression traits [89 of 1314 cis-eQTLs at P < 1 × 10⁻⁰⁶ (false discovery rate ∼5%)] and one of eight biomarker associations at P < 8 × 10⁻¹⁰. Very few (30 of 1232; 2%) common variant associations were fully explained by low-frequency variants. Our data show that whole-genome sequencing can identify low-frequency variants undetected by genotyping based approaches when sample sizes are sufficiently large to detect substantial numbers of common variant associations, and that common variant associations are rarely explained by single low-frequency variants of large effect.

Human plasma lipidome is pleiotropically associated with cardiovascular risk factors and death

BACKGROUND

Cardiovascular disease (CVD) is the most common cause of death in the United States and is associated with a high economic burden. Prevention of CVD focuses on controlling or improving the lipid profile of patients at risk. The human lipidome is made up of thousands of ubiquitous lipid species. By studying biologically simple canonical lipid species, we investigated whether the lipidome is genetically redundant and whether its genetic influences can provide clinically relevant clues of CVD risk.

METHODS AND RESULTS

We performed a genetic study of the human lipidome in 1212 individuals from 42 extended Mexican American families. High-throughput mass spectrometry enabled rapid capture of precise lipidomic profiles, providing 319 unique species. Using variance component-based heritability analyses and bivariate trait analyses, we detected significant genetic influences on each lipid assayed. Median heritability of the plasma lipid species was 0.37. Hierarchical clustering based on complex genetic correlation patterns identified 12 genetic clusters that characterized the plasma lipidome. These genetic clusters were differentially but consistently associated with risk factors of CVD, including central obesity, obesity, type 2 diabetes mellitus, raised serum triglycerides, and metabolic syndrome. Also, these clusters consistently predicted occurrence of cardiovascular deaths during follow-up.

CONCLUSIONS

The human plasma lipidome is heritable. Shared genetic influences reduce the dimensionality of the human lipidome into clusters that are associated with risk factors of CVD.

Plasma levels of soluble interleukin 1 receptor accessory protein are reduced in obesity

CONTEXT

Adipokines actuate chronic, low-grade inflammation through a complex network of immune markers, but the current understanding of these networks is incomplete. The soluble isoform of the IL-1 receptor accessory protein (sIL1RAP) occupies an important position in the inflammatory pathways involved in obesity. The pathogenetic and clinical influences of sIL1RAP are unknown.

OBJECTIVE

The objective of the study was to elucidate whether plasma levels of sIL1RAP are reduced in obesity, using affluent clinical, biochemical, and genetic data from two diverse cohorts.

DESIGN, SETTING, AND PARTICIPANTS

The study was conducted in two cohorts: the San Antonio Family Heart Study (n = 1397 individuals from 42 families) and South Asians living in Mauritius, n = 230).

MAIN OUTCOME MEASURES

Plasma sIL1RAP levels were measured using an ELISA. The genetic basis of sIL1RAP levels were investigated using both a large-scale gene expression profiling study and a genome-wide association study.

RESULTS

A significant decrease in plasma sIL1RAP levels were observed in obese subjects, even after adjustment for age and sex. The sIL1RAP levels demonstrated a strong inverse association with obesity measures in both populations. All associations were more significant in females. Plasma sIL1RAP levels were significantly heritable, correlated with IL1RAP transcript levels (NM_134470), showed evidence for shared genetic influences with obesity measures and were significantly associated with the rs2885373 single-nucleotide polymorphism (P = 6.7 × 10(-23)) within the IL1RAP gene.

CONCLUSIONS

Plasma sIL1RAP levels are reduced in obesity and can potentially act as biomarkers of obesity. Mechanistic studies are required to understand the exact contribution of sIL1RAP to the pathogenesis of obesity.

The ENIGMA Consortium: large-scale collaborative analyses of neuroimaging and genetic data

The Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium is a collaborative network of researchers working together on a range of large-scale studies that integrate data from 70 institutions worldwide. Organized into Working Groups that tackle questions in neuroscience, genetics, and medicine, ENIGMA studies have analyzed neuroimaging data from over 12,826 subjects. In addition, data from 12,171 individuals were provided by the CHARGE consortium for replication of findings, in a total of 24,997 subjects. By meta-analyzing results from many sites, ENIGMA has detected factors that affect the brain that no individual site could detect on its own, and that require larger numbers of subjects than any individual neuroimaging study has currently collected. ENIGMA's first project was a genome-wide association study identifying common variants in the genome associated with hippocampal volume or intracranial volume. Continuing work is exploring genetic associations with subcortical volumes (ENIGMA2) and white matter microstructure (ENIGMA-DTI). Working groups also focus on understanding how schizophrenia, bipolar illness, major depression and attention deficit/hyperactivity disorder (ADHD) affect the brain. We review the current progress of the ENIGMA Consortium, along with challenges and unexpected discoveries made on the way.

Heritable changes in regional cortical thickness with age

It is now well established that regional indices of brain structure such as cortical thickness, surface area or grey matter volume exhibit spatially variable patterns of heritability. However, a recent study found these patterns to change with age during development, a result supported by gene expression studies. Changes in heritability have not been investigated in adulthood so far and could have important implications in the study of heritability and genetic correlations in the brain as well as in the discovery of specific genes explaining them. Herein, we tested for genotype by age (G ×A) interactions, an extension of genotype by environment interactions, through adulthood and healthy aging in 902 subjects from the Genetics of Brain Structure (GOBS) study. A "jackknife" based method for the analysis of stable cortical thickness clusters (JASC) and scale selection is also introduced. Although additive genetic variance remained constant throughout adulthood, we found evidence for incomplete pleiotropy across age in the cortical thickness of paralimbic and parieto-temporal areas. This suggests that different genetic factors account for cortical thickness heritability at different ages in these regions.

Plasma lipidome is independently associated with variability in metabolic syndrome in Mexican American families

Plasma lipidome is now increasingly recognized as a potentially important marker of chronic diseases, but the exact extent of its contribution to the interindividual phenotypic variability in family studies is unknown. Here, we used the rich data from the ongoing San Antonio Family Heart Study (SAFHS) and developed a novel statistical approach to quantify the independent and additive value of the plasma lipidome in explaining metabolic syndrome (MS) variability in Mexican American families recruited in the SAFHS. Our analytical approach included two preprocessing steps: principal components analysis of the high-resolution plasma lipidomics data and construction of a subject-subject lipidomic similarity matrix. We then used the Sequential Oligogenic Linkage Analysis Routines software to model the complex family relationships, lipidomic similarities, and other important covariates in a variance components framework. Our results suggested that even after accounting for the shared genetic influences, indicators of lipemic status (total serum cholesterol, TGs, and HDL cholesterol), and obesity, the plasma lipidome independently explained 22% of variability in the homeostatic model of assessment-insulin resistance trait and 16% to 22% variability in glucose, insulin, and waist circumference. Our results demonstrate that plasma lipidomic studies can additively contribute to an understanding of the interindividual variability in MS.

Genome-wide association analysis confirms and extends the association of SLC2A9 with serum uric acid levels to Mexican Americans

Increased serum uric acid (SUA) is a risk factor for gout and renal and cardiovascular disease (CVD). The purpose of this study was to identify genetic factors that affect the variation in SUA in 632 Mexican Americans participants of the San Antonio Family Heart Study (SAFHS). A genome-wide association (GWA) analysis was performed using the Illumina Human Hap 550K single nucleotide polymorphism (SNP) microarray. We used a linear regression-based association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component. All analyses were performed in the software package SOLAR. SNPs rs6832439, rs13131257, and rs737267 in solute carrier protein 2 family, member 9 (SLC2A9) were associated with SUA at genome-wide significance (p < 1.3 × 10⁻⁷). The minor alleles of these SNPs had frequencies of 36.2, 36.2, and 38.2%, respectively, and were associated with decreasing SUA levels. All of these SNPs were located in introns 3–7 of SLC2A9, the location of the previously reported associations in European populations. When analyzed for association with cardiovascular-renal disease risk factors, conditional on SLC2A9 SNPs strongly associated with SUA, significant associations were found for SLC2A9 SNPs with BMI, body weight, and waist circumference (p < 1.4 × 10⁻³) and suggestive associations with albumin-creatinine ratio and total antioxidant status (TAS). The SLC2A9 gene encodes an urate transporter that has considerable influence on variation in SUA. In addition to the primary association locus, suggestive evidence (p < 1.9 × 10⁻⁶) for joint linkage/association (JLA) was found at a previously-reported urate quantitative trait locus (Logarithm of odds score = 3.6) on 3p26.3. In summary, our GWAS extends and confirms the association of SLC2A9 with SUA for the first time in a Mexican American cohort and also shows for the first time its association with cardiovascular-renal disease risk factors.

QTL-based association analyses reveal novel genes influencing pleiotropy of metabolic syndrome (MetS)

OBJECTIVE

Metabolic Syndrome (MetS) is a phenotype cluster predisposing to type 2 diabetes and cardiovascular disease. We conducted a study to elucidate the genetic basis underlying linkage signals for multiple representative traits of MetS that we had previously identified at two significant QTLs on chromosomes 3q27 and 17p12.

DESIGN AND METHODS

We performed QTL-specific genomic and transcriptomic analyses in 1,137 individuals from 85 extended families that contributed to the original linkage. We tested in SOLAR association of MetS phenotypes with QTL-specific haplotype-tagging SNPs as well as transcriptional profiles of peripheral blood mononuclear cells (PBMCs).

RESULTS

SNPs significantly associated with MetS phenotypes under the prior hypothesis of linkage mapped to seven genes at 3q27 and seven at 17p12. Prioritization based on biologic relevance, SNP association, and expression analyses identified two genes: insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) at 3q27 and tumor necrosis factor receptor 13B (TNFRSF13B) at 17p12. Prioritized genes could influence cell-cell adhesion and adipocyte differentiation, insulin/glucose responsiveness, cytokine effectiveness, plasma lipid levels, and lipoprotein densities.

CONCLUSIONS

Using an approach combining genomic, transcriptomic, and bioinformatic data we identified novel candidate genes for MetS.

Epidemiology and genetic determinants of progressive deterioration of glycaemia in American Indians: the Strong Heart Family Study

AIMS/HYPOTHESIS

Type 2 diabetes is a chronic, heterogeneous disease and a major risk factor for cardiovascular diseases. The underlying mechanisms leading to progression to type 2 diabetes are not fully understood and genetic tools may help to identify important pathways of glycaemic deterioration.

METHODS

Using prospective data on American Indians from the Strong Heart Family Study, we identified 373 individuals defined as progressors (diabetes incident cases), 566 individuals with transitory impaired fasting glucose (IFG) and 1,011 controls (normal fasting glycaemia at all visits). We estimated the heritability (h(2)) of the traits and the evidence for association with 16 known variants identified in type 2 diabetes genome-wide association studies.

RESULTS

We noted high h(2) for diabetes progression (h(2) = 0.65 ± 0.16, p = 2.7 × 10(-6)) but little contribution of genetic factors to transitory IFG (h(2) = 0.09 ± 0.10, p = 0.19) for models adjusted for multiple risk factors. At least three variants (in WFS1, TSPAN8 and THADA) were nominally associated with diabetes progression in age- and sex-adjusted analyses with estimates showing the same direction of effects as reported in the discovery European ancestry studies.

CONCLUSIONS/INTERPRETATION

Our findings do not exclude these loci for diabetes susceptibility in American Indians and suggest phenotypic heterogeneity of the IFG trait, which may have implications for genetic studies when diagnosis is based on a single time-point measure.