Search

IntroductionThe 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.MethodsWe 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.ResultsAll 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.DiscussionOverall, 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.

Type 2 diabetes is associated with cognitive impairments, but it is unclear whether common genetic factors influence both type 2 diabetes risk and cognition. 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. 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. 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. 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.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare variants’ (RVs) associations with complex human traits. Variant set analysis is a powerful approach to study RV association, and a key component of it is constructing RV sets for analysis. However, existing methods have limited ability to define analysis units in the noncoding genome. Furthermore, there is a lack of robust pipelines for comprehensive and scalable noncoding RV association analysis. Here we propose a computationally-efficient noncoding RV association-detection framework that uses STAAR (variant-set test for association using annotation information) to group noncoding variants in gene-centric analysis based on functional categories. We also propose SCANG (scan the genome)-STAAR, which uses dynamic window sizes and incorporates multiple functional annotations, in a non-gene-centric analysis. We furthermore develop STAARpipeline to perform flexible noncoding RV association analysis, including gene-centric analysis as well as fixed-window-based and dynamic-window-based non-gene-centric analysis. We apply STAARpipeline to identify noncoding RV sets associated with four quantitative lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several noncoding RV associations in an additional 9,123 TOPMed samples.

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.

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 2 =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.

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.

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.

BACKGROUND: The identification and understanding of therapeutic targets for atherosclerotic cardiovascular disease (ASCVD) is of fundamental importance given its global health and economic burden. Inhibition of angiopoietin-like 3 (ANGPTL3) has demonstrated a cardioprotective effect, showing promise for ASCVD 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 ~1,000 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(2)=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 (LOD = 6.18) to a locus on chromosome 1 at 90 cM, corresponding to the ANGPTL3 gene location. In the genomes of 23 individuals from a single pedigree, we identified a loss of function (LoF) variant, rs398122988 (N121Kfs*2), in ANGPTL3, that was significantly associated with lower ANGPTL3 levels (β=−1.69 SDU, 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 LoF variants in ANGPTL3 have significant associations with the plasma lipidome. These findings further our understanding of ANGPTL3 as a therapeutic target for ASCVD.

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 −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.

High-throughput platforms allow the characterization of thousands of previously known methylation sites. These platforms have great potential for investigating the epigenetic effects that are partially responsible for gene expression control. Methylation sites provide a bridge for the investigation of real-time environmental contributions on genomic events by the alteration of methylation status of those sites. Using the data provided by GAW20’s organization committee, we calculated the heritability estimates of each cytosine-phosphate-guanine (CpG) island before and after the use of fenofibrate, a lipid-control drug. Surprisingly, we detected substantially high heritability estimates before drug usage. This somewhat unexpected high sample correlation was corrected by the use of principal components and the distributions of heritability estimates before and after fenofibrate treatment, which made the distributions comparable. The methylation sites located near a gene were collected and a genetic relationship matrix estimated to represent the overall correlation between samples. We implemented a random-effect association test to screen genes whose methylation patterns partially explain the observable high-density lipoprotein (HDL) heritability. Our leading association was observed for the TMEM52 gene that encodes a transmembrane protein, and is largely expressed in the liver, had not been previously associated with HDL until this manuscript. Using a variance component decomposition framework with the linear mixed model allows the integration of data from different sources, such as methylation, gene expression, metabolomics, and proteomics. The decomposition of the genetic variance component decomposition provides a flexible analytical approach for the challenges of this new omics era.

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(2) = 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

Lipoprotein(a), Lp(a), is a modified low-density lipoprotein particle that contains apolipoprotein(a), encoded by LPA, and is a highly heritable, causal risk factor for cardiovascular diseases that varies in concentrations across ancestries. Here, we use deep-coverage whole genome sequencing in 8392 individuals of European and African ancestry to discover and interpret both single-nucleotide variants and copy number (CN) variation associated with Lp(a). We observe that genetic determinants between Europeans and Africans have several unique determinants. The common variant rs12740374 associated with Lp(a) cholesterol is an eQTL for SORT1 and independent of LDL cholesterol. Observed associations of aggregates of rare non-coding variants are largely explained by LPA structural variation, namely the LPA kringle IV 2 (KIV2)-CN. Finally, we find that LPA risk genotypes confer greater relative risk for incident atherosclerotic cardiovascular diseases compared to directly measured Lp(a), and are significantly associated with measures of subclinical atherosclerosis in African Americans., Circulating lipoprotein(a) is an important risk factor for cardiovascular disease and shows variability between different ethnic groups. Here, Zekavat et al. perform whole-genome sequencing in individuals of European and African ancestries and find ancestry-specific genetic determinants for lipoprotein(a) levels.

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.

Statistical analysis of multinomial data in complex datasets often requires estimation of the multivariate normal (mvn) distribution for models in which the dimensionality can easily reach 10–1000 and higher. Few algorithms for estimating the mvn distribution can offer robust and efficient performance over such a range of dimensions. We report a simulation-based comparison of two algorithms for the mvn that are widely used in statistical genetic applications. The venerable Mendell-Elston approximation is fast but execution time increases rapidly with the number of dimensions, estimates are generally biased, and an error bound is lacking. The correlation between variables significantly affects absolute error but not overall execution time. The Monte Carlo-based approach described by Genz returns unbiased and error-bounded estimates, but execution time is more sensitive to the correlation between variables. For ultra-high-dimensional problems, however, the Genz algorithm exhibits better scale characteristics and greater time-weighted efficiency of estimation.

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

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.

The accurate estimation of age using methylation data has proved a useful and heritable biomarker, with acceleration in epigenetic age predicting a number of age-related phenotypes. Measures of white matter integrity in the brain are also heritable and highly sensitive to both normal and pathological aging processes across adulthood. We consider the phenotypic and genetic interrelationships between epigenetic age acceleration and white matter integrity in humans. Our goal was to investigate processes that underlie interindividual variability in age-related changes in the brain. Using blood taken from a Mexican-American extended pedigree sample (n= 628; age = 23.28–93.11 years), epigenetic age was estimated using the method developed by Horvath (2013). Forn= 376 individuals, diffusion tensor imaging scans were also available. The interrelationship between epigenetic age acceleration and global white matter integrity was investigated with variance decomposition methods. To test for neuroanatomical specificity, 16 specific tracts were additionally considered. We observed negative phenotypic correlations between epigenetic age acceleration and global white matter tract integrity (ρpheno= −0.119,p= 0.028), with evidence of shared genetic (ρgene= −0.463,p= 0.013) but not environmental influences. Negative phenotypic and genetic correlations with age acceleration were also seen for a number of specific white matter tracts, along with additional negative phenotypic correlations between granulocyte abundance and white matter integrity. These findings (i.e., increased acceleration in epigenetic age in peripheral blood correlates with reduced white matter integrity in the brain and shares common genetic influences) provide a window into the neurobiology of aging processes within the brain and a potential biomarker of normal and pathological brain aging.SIGNIFICANCE STATEMENTEpigenetic measures can be used to predict age with a high degree of accuracy and so capture acceleration in biological age, relative to chronological age. The white matter tracts within the brain are also highly sensitive to aging processes. We show that increased biological aging (measured using epigenetic data from blood samples) is correlated with reduced integrity of white matter tracts within the human brain (measured using diffusion tensor imaging) with data from a large sample of Mexican-American families. Given the family design of the sample, we are also able to demonstrate that epigenetic aging and white matter tract integrity also share common genetic influences. Therefore, epigenetic age may be a potential, and accessible, biomarker of brain aging.

Obesity is one of the most prevalent metabolic diseases in the Western world and correlates directly with insulin resistance, which may ultimately culminate in type 2 diabetes (T2D). We sought to ascertain whether the human metalloproteinase A Disintegrin and Metalloproteinase 19 (ADAM19) correlates with parameters of the metabolic syndrome in humans and mice. To determine the potential novel role of ADAM19 in the metabolic syndrome, we first conducted microarray studies on peripheral blood mononuclear cells from a well-characterised human cohort. Secondly, we examined the expression of ADAM19 in liver and gonadal white adipose tissue using an in vivo diet induced obesity mouse model. Finally, we investigated the effect of neutralising ADAM19 on diet induced weight gain, insulin resistance in vivo, and liver TNF-α levels. Significantly, we show that, in humans, ADAM19 strongly correlates with parameters of the metabolic syndrome, particularly BMI, relative fat, HOMA-IR, and triglycerides. Furthermore, we identified that ADAM19 expression was markedly increased in the liver and gonadal white adipose tissue of obese and T2D mice. Excitingly, we demonstrate in our diet induced obesity mouse model that neutralising ADAM19 therapy results in weight loss, improves insulin sensitivity, and reduces liver TNF-α levels. Our novel data suggest that ADAM19 is pro-obesogenic and enhances insulin resistance. Therefore, neutralisation of ADAM19 may be a potential therapeutic approach to treat obesity and T2D.

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.

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.

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 [h2 = 0.614, P = 1.00 × 10−6, standard error (SE) = 0.151] and major depression (h2 = 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.

Type 2 diabetes (T2D) is a common, multifactorial disease that is influenced by genetic and environmental factors and their interactions. However, common variants identified by genome wide association studies (GWAS) explain only about 10% of the total trait variance for T2D and less than 5% of the variance for obesity, indicating that a large proportion of heritability is still unexplained. The transcriptomic approach described here uses quantitative gene expression and disease-related physiological data (deep phenotyping) to measure the direct correlation between the expression of specific genes and physiological traits. Transcriptomic analysis bridges the gulf between GWAS and physiological studies. Recent GWAS studies have utilized very large population samples, numbering in the tens of thousands (or even hundreds of thousands) of individuals, yet establishing causal functional relationships between strongly associated genetic variants and disease remains elusive. In light of the findings described below, it is appropriate to consider how and why transcriptomic approaches in small samples might be capable of identifying complex disease-related genes which are not apparent using GWAS in large samples.

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.

Processing speed is a psychological construct that refers to the speed with which an individual can perform any cognitive operation. Processing speed correlates strongly with general cognitive ability, declines sharply with age and is impaired across a number of neurological and psychiatric disorders. Thus, identifying genes that influence processing speed will likely improve understanding of the genetics of intelligence, biological aging and the etiologies of numerous disorders. Previous genetics studies of processing speed have relied on simple phenotypes (eg, mean reaction time) derived from single tasks. This strategy assumes, erroneously, that processing speed is a unitary construct. In the present study, we aimed to characterize the genetic architecture of processing speed by using a multidimensional model applied to a battery of cognitive tasks. Linkage and QTL-specific association analyses were performed on the factors from this model. The randomly ascertained sample comprised 1291 Mexican-American individuals from extended pedigrees. We found that performance on all three distinct processing-speed factors (Psychomotor Speed; Sequencing and Shifting and Verbal Fluency) were moderately and significantly heritable. We identified a genome-wide significant quantitative trait locus (QTL) on chromosome 3q23 for Psychomotor Speed (LOD = 4.83). Within this locus, we identified a plausible and interesting candidate gene for Psychomotor Speed (Z = 2.90, P = 1.86 × 10-03 ).

Background 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. Methods 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. Results 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. Discussion 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.

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.

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.

The insula and medial prefrontal cortex (mPFC) share functional, histological, transcriptional, and developmental characteristics, and they serve higher cognitive functions of theoretical relevance to schizophrenia and related disorders. Meta-analyses and multivariate analysis of structural magnetic resonance imaging (MRI) scans indicate that gray matter density and volume reductions in schizophrenia are the most consistent and pronounced in a network primarily composed of the insula and mPFC. We used source-based morphometry, a multivariate technique optimized for structural MRI, in a large sample of randomly ascertained pedigrees (N = 887) to derive an insula-mPFC component and to investigate its genetic determinants. Firstly, we replicated the insula-mPFC gray matter component as an independent source of gray matter variation in the general population, and verified its relevance to schizophrenia in an independent case-control sample. Secondly, we showed that the neuroanatomical variation defined by this component is largely determined by additive genetic variation (h(2) = 0.59), and genome-wide linkage analysis resulted in a significant linkage peak at 12q24 (LOD = 3.76). This region has been of significant interest to psychiatric genetics as it contains the Darier's disease locus and other proposed susceptibility genes (e.g., DAO, NOS1), and it has been linked to affective disorders and schizophrenia in multiple populations. Thus, in conjunction with previous clinical studies, our data imply that one or more psychiatric risk variants at 12q24 are co-inherited with reductions in mPFC and insula gray matter concentration. © 2015 Wiley Periodicals, Inc.

Background: Metabolomics may unravel important biological pathways involved in the pathophysiology of childhood obesity. Objectives: We aimed to 1) identify metabolites that differ significantly between nonobese and obese Hispanic children; 2) collapse metabolites into principal components (PCs) associated with obesity and metabolic risk, specifically hyperinsulinemia, hypertriglyceridemia, hyperleptinemia, and hyperuricemia; and 3) identify metabolites associated with energy expenditure and fat oxidation. Design: This trial was a cross-sectional observational study of metabolomics by using gas chromatography–mass spectrometry and ultrahigh-performance liquid chromatography–tandem mass spectrometry analyses performed on fasting plasma samples from 353 nonobese and 450 obese Hispanic children. Results: Branched-chained amino acids (BCAAs) (Leu, Ile, and Val) and their catabolites, propionylcarnitine and butyrylcarnitine, were significantly elevated in obese children. Strikingly lower lysolipids and dicarboxylated fatty acids were seen in obese children. Steroid derivatives were markedly higher in obese children as were markers of inflammation and oxidative stress. PC6 (BCAAs and aromatic AAs) and PC10 (asparagine, glycine, and serine) made the largest contributions to body mass index, and PC10 and PC12 (acylcarnitines) made the largest contributions to adiposity. Metabolic risk factors and total energy expenditure were associated with PC6, PC9 (AA and tricarboxylic acid cycle metabolites), and PC10. Fat oxidation was inversely related to PC8 (lysolipids) and positively related to PC16 (acylcarnitines). Conclusions: Global metabolomic profiling in nonobese and obese children replicates the increased BCAA and acylcarnitine catabolism and changes in nucleotides, lysolipids, and inflammation markers seen in obese adults; however, a strong signature of reduced fatty acid catabolism and increased steroid derivatives may be unique to obese children. Metabolic flexibility in fuel use observed in obese children may occur through the activation of alternative intermediary pathways. Insulin resistance, hyperleptinemia, hypertriglyceridemia, hyperuricemia, and oxidative stress and inflammation evident in obese children are associated with distinct metabolomic profiles.

Migraine has been defined as a common disabling primary headache disorder. Epidemiology studies have provided with the undeniable evidence of genetic components as active players in the development of the disease under a polygenic model in which multiple risk alleles exert modest individual effects. Our objective was to test the contribution of a polygenic effect to migraine risk in the Norfolk Island population using a panel of SNPs reported to be disease associated in published migraine GWAS. We also investigated whether individual SNPs were associated with gene expression levels measured in whole blood. Polygenic scores were calculated in a total of 285 related individuals (74 cases, 211 controls) from the Norfolk Island using 51 SNPs previously reported to be associated with migraine in published GWAS. The association between polygenic score and migraine case-control status was tested using logistic regression. Results indicate that a migraine polygenic risk score was associated with migraine case-control status in this population (P = 0.016). This supports the hypothesis that multiple SNPs with weak effects collectively contribute to migraine risk in this population. Amongst the SNPs included in the polygenic model, four were associated with the expression of the USMG5 gene, including rs171251 (P = 0.012). Results from this study provide evidence for a polygenic contribution to migraine risk in an isolated population and highlight specific SNPs that regulate the expression of USMG5, a gene critical for mitochondrial function.

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.

Whipworms (Trichuris spp.) infect a variety of hosts, including domestic animals and humans. Of considerable interest is the porcine whipworm, T. suis, which is particularly prevalent in outdoor production systems. High infection levels may cause growth retardation, anaemia and haemorrhagic diarrhoea. A significant proportion of the variation in Trichuris faecal egg count (FEC) has been attributed to the host's genetic make-up. The aim of the present study was to identify genetic loci associated with resistance to T. suis in pigs. We used single nucleotide polymorphism (SNP) markers to perform a whole-genome scan of an F1 resource population (n = 195) trickle-infected with T. suis. A measured genotype analysis revealed a putative quantitative trait locus (QTL) for T. suis FEC on chromosome 13 covering ∼ 4.5 Mbp, although none of the SNPs reached genome-wide significance. We tested the hypothesis that this region of SSC13 harboured genes with effects on T. suis burden by genotyping three SNPs within the putative QTL in unrelated pigs exposed to either experimental or natural T. suis infections and from which we had FEC (n = 113) or worm counts (n = 178). In these studies, two of the SNPs (rs55618716, ST) were associated with FEC (P < 0.01), thus confirming our initial findings. However, we did not find any of the SNPs to be associated with T. suis worm burden. In conclusion, our study demonstrates that genetic markers for resistance to T. suis as indicated by low FEC can be identified in pigs.

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.

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.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.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.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.

Only few systematic studies on the contribution of copy number variation to gene expression variation have been published to date. Here we identify effects of copy number variable regions (CNVRs) on nearby gene expression by investigating 909 CNVRs and expression levels of 12059 nearby genes in white blood cells from Mexican-American participants of the San Antonio Family Heart Study. We empirically evaluate our ability to detect the contribution of CNVs to proximal gene expression (presumably in cis) at various window sizes (up to a 10 Mb distance) between the gene and CNV. We found a ~1-Mb window size to be optimal for capturing cis effects of CNVs. Up to 10% of the CNVs in this study were found to be significantly associated with the expression of at least one gene within their vicinity. As expected, we find that CNVs that directly overlap gene sequences have the largest effects on gene expression (compared with non-overlapping CNVRs located nearby), with positive correlation (except for a few exceptions) between estimated genomic dosage and expression level. We find that genes whose expression level is significantly influenced by nearby CNVRs are enriched for immunity and autoimmunity related genes. These findings add to the currently limited catalog of CNVRs that are recognized as expression quantitative trait loci, and have implications for future study designs as well as for prioritizing candidate causal variants in genomic regions associated with disease.

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.

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.

Phasing, the process of predicting haplotypes from genotype data, is an important undertaking in genetics and an ongoing area of research. Phasing methods, and associated software, designed specifically for pedigrees are urgently needed. Here we present a new method for phasing genotypes from whole genome sequencing data in pedigrees: PULSAR (Phasing Using Lineage Specific Alleles / Rare variants). The method is built upon the idea that alleles that are specific to a single founding chromosome within a pedigree, which we refer to as lineage-specific alleles, are highly informative for identifying haplotypes that are identical-by-decent between individuals within a pedigree. Through extensive simulation we assess the performance of PULSAR in a variety of pedigree sizes and structures, and we explore the effects of genotyping errors and presence of non-sequenced individuals on its performance. If the genotyping error rate is sufficiently low PULSAR can phase > 99.9% of heterozygous genotypes with a switch error rate below 1 x 10-4 in pedigrees where all individuals are sequenced. We demonstrate that the method is highly accurate and consistently outperforms the long-range phasing approach used for comparison in our benchmarking. The method also holds promise for fixing genotype errors or imputing missing genotypes. The software implementation of this method is freely available.

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.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

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.

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.

Linkage studies of alcoholism have implicated several chromosome regions, leading to the successful identification of susceptibility genes, including ADH4 and GABRA2 on chromosome 4. Quantitative endophenotypes that are potentially closer to gene action than clinical endpoints offer a means of obtaining more refined linkage signals of genes that predispose alcohol use disorders (AUD). In this study we examine a self-reported measure of the maximum number of drinks consumed in a 24-hr period (abbreviated Max Drinks), a significantly heritable phenotype (h2 = 0.32 ± 0.05; P = 4.61 × 10−14) with a strong genetic correlation with AUD (ρg = 0.99 ± 0.13) for the San Antonio Family Study (n = 1,203). Genome-wide SNPs were analyzed using variance components linkage methods in the program SOLAR, revealing a novel, genome-wide significant QTL (LOD = 4.17; P = 5.85 × 10−6) for Max Drinks at chromosome 6p22.3, a region with a number of compelling candidate genes implicated in neuronal function and psychiatric illness. Joint analysis of Max Drinks and AUD status shows that the QTL has a significant non-zero effect on diagnosis (P = 4.04 × 10−3), accounting for 8.6% of the total variation. Significant SNP associations for Max Drinks were also identified at the linkage region, including one, rs7761213 (P = 2.14 × 10−4), obtained for an independent sample of Chinese families. Thus, our study identifies a potential risk locus for AUD at 6p22.3, with significant pleiotropic effects on the heaviness of alcohol consumption that may not be population specific. © 2014 Wiley Periodicals, Inc.