Rapid, Phase-free Detection of Long Identity-by-Descent Segments Enables Effective Relationship Classification

Identity-by-descent (IBD) segments are a useful tool for applications ranging from demographic inference to relationship classification, but most detection methods rely on phasing information and therefore require substantial computation time. As genetic datasets grow, methods for inferring IBD segments that scale well will be critical. We developed IBIS, an IBD detector that locates long regions of allele sharing between unphased individuals, and benchmarked it with Refined IBD, GERMLINE, and TRUFFLE on 3,000 simulated individuals. Phasing these with Beagle 5 takes 4.3 CPU days, followed by either Refined IBD or GERMLINE segment detection in 2.9 or 1.1 h, respectively. By comparison, IBIS finishes in 6.8 min or 7.8 min with IBD2 functionality enabled: speedups of 805-946× including phasing time. TRUFFLE takes 2.6 h, corresponding to IBIS speedups of 20.2-23.3×. IBIS is also accurate, inferring ≥7 cM IBD segments at quality comparable to Refined IBD and GERMLINE. With these segments, IBIS classifies first through third degree relatives in real Mexican American samples at rates meeting or exceeding other methods tested and identifies fourth through sixth degree pairs at rates within 0.0%-2.0% of the top method. While allele frequency-based approaches that do not detect segments can infer relationship degrees faster than IBIS, the fastest are biased in admixed samples, with KING inferring 30.8% fewer fifth degree Mexican American relatives correctly compared with IBIS. Finally, we ran IBIS on chromosome 2 of the UK Biobank dataset and estimate its runtime on the autosomes to be 3.3 days parallelized across 128 cores.

Crossover interference and sex-specific genetic maps shape identical by descent sharing in close relatives

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of pairwise IBD proportions by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD sharing ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available from https://github.com/williamslab/ped-sim.

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.

Rare DEGS1 variant significantly alters de novo ceramide synthesis pathway

The de novo ceramide synthesis pathway is essential to human biology and health, but genetic influences remain unexplored. The core function of this pathway is the generation of biologically active ceramide from its precursor, dihydroceramide. Dihydroceramides have diverse, often protective, biological roles; conversely, increased ceramide levels are biomarkers of complex disease. To explore the genetics of the ceramide synthesis pathway, we searched for deleterious nonsynonymous variants in the genomes of 1,020 Mexican Americans from extended pedigrees. We identified a Hispanic ancestry-specific rare functional variant, L175Q, in delta 4-desaturase, sphingolipid 1 (DEGS1), a key enzyme in the pathway that converts dihydroceramide to ceramide. This amino acid change was significantly associated with large increases in plasma dihydroceramides. Indexes of DEGS1 enzymatic activity were dramatically reduced in heterozygotes. CRISPR/Cas9 genome editing of HepG2 cells confirmed that the L175Q variant results in a partial loss of function for the DEGS1 enzyme. Understanding the biological role of DEGS1 variants, such as L175Q, in ceramide synthesis may improve the understanding of metabolic-related disorders and spur ongoing research of drug targets along this pathway.

A QTL on chromosome 3q23 influences processing speed in humans

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

Inferring Identical-by-Descent Sharing of Sample Ancestors Promotes High-Resolution Relative Detection

As genetic datasets increase in size, the fraction of samples with one or more close relatives grows rapidly, resulting in sets of mutually related individuals. We present DRUID-deep relatedness utilizing identity by descent-a method that works by inferring the identical-by-descent (IBD) sharing profile of an ungenotyped ancestor of a set of close relatives. Using this IBD profile, DRUID infers relatedness between unobserved ancestors and more distant relatives, thereby combining information from multiple samples to remove one or more generations between the deep relationships to be identified. DRUID constructs sets of close relatives by detecting full siblings and also uses an approach to identify the aunts/uncles of two or more siblings, recovering 92.2% of real aunts/uncles with zero false positives. In real and simulated data, DRUID correctly infers up to 10.5% more relatives than PADRE when using data from two sets of distantly related siblings, and 10.7%-31.3% more relatives given two sets of siblings and their aunts/uncles. DRUID frequently infers relationships either correctly or within one degree of the truth, with PADRE classifying 43.3%-58.3% of tenth degree relatives in this way compared to 79.6%-96.7% using DRUID.

Neurodegenerative disease biomarkers Aβ1-40, Aβ1-42, tau, and p-tau181 in the vervet monkey cerebrospinal fluid: Relation to normal aging, genetic influences, and cerebral amyloid angiopathy

Background

The Caribbean vervet monkey (Chlorocebus aethiops sabaeus) is a potentially valuable animal model of neurodegenerative disease. However, the trajectory of aging in vervets and its relationship to human disease is incompletely understood.

Methods

To characterize biomarkers associated with neurodegeneration, we measured cerebrospinal fluid (CSF) concentrations of Aβ1-40, Aβ1-42, total tau, and p-tau181 in 329 members of a multigenerational pedigree. Linkage and genome-wide association were used to elucidate a genetic contribution to these traits.

Results

Aβ1-40 concentrations were significantly correlated with age, brain total surface area, and gray matter thickness. Levels of p-tau181 were associated with cerebral volume and brain total surface area. Among the measured analytes, only CSF Aβ1-40 was heritable. No significant linkage (LOD > 3.3) was found, though suggestive linkage was highlighted on chromosomes 4 and 12. Genome-wide association identified a suggestive locus near the chromosome 4 linkage peak.

Conclusions

Overall, these results support the vervet as a non-human primate model of amyloid-related neurodegeneration, such as Alzheimer's disease and cerebral amyloid angiopathy, and highlight Aβ1-40 and p-tau181 as potentially valuable biomarkers of these processes.

Genetic variation and gene expression across multiple tissues and developmental stages in a nonhuman primate

By analyzing multitissue gene expression and genome-wide genetic variation data in samples from a vervet monkey pedigree, we generated a transcriptome resource and produced the first catalog of expression quantitative trait loci (eQTLs) in a nonhuman primate model. This catalog contains more genome-wide significant eQTLs per sample than comparable human resources and identifies sex- and age-related expression patterns. Findings include a master regulatory locus that likely has a role in immune function and a locus regulating hippocampal long noncoding RNAs (lncRNAs), whose expression correlates with hippocampal volume. This resource will facilitate genetic investigation of quantitative traits, including brain and behavioral phenotypes relevant to neuropsychiatric disorders.

Benchmarking Relatedness Inference Methods with Genome-Wide Data from Thousands of Relatives

Inferring relatedness from genomic data is an essential component of genetic association studies, population genetics, forensics, and genealogy. While numerous methods exist for inferring relatedness, thorough evaluation of these approaches in real data has been lacking. Here, we report an assessment of 12 state-of-the-art pairwise relatedness inference methods using a data set with 2485 individuals contained in several large pedigrees that span up to six generations. We find that all methods have high accuracy (92-99%) when detecting first- and second-degree relationships, but their accuracy dwindles to <43% for seventh-degree relationships. However, most identical by descent (IBD) segment-based methods inferred seventh-degree relatives correct to within one relatedness degree for >76% of relative pairs. Overall, the most accurate methods are Estimation of Recent Shared Ancestry (ERSA) and approaches that compute total IBD sharing using the output from GERMLINE and Refined IBD to infer relatedness. Combining information from the most accurate methods provides little accuracy improvement, indicating that novel approaches, such as new methods that leverage relatedness signals from multiple samples, are needed to achieve a sizeable jump in performance.

Bioprotection of the built environment and cultural heritage

The growth of microbial biofilms and various biomineralization phenomena can lead to the formation of stable layers and veneers on rocks known as ‘rock varnishes’ that can stabilize surfaces and protect from further weathering. This article describes the potential application of fungal systems for bioprotection of rock and mineral‐based substrates and the evidence to support this concept of utilizing natural or engineered colonization and metabolic properties of fungi, including lichens.

Genetic correlation of the plasma lipidome with type 2 diabetes, prediabetes and insulin resistance in Mexican American families

Background

Differential plasma concentrations of circulating lipid species are associated with pathogenesis of type 2 diabetes (T2D). Whether the wide inter-individual variability in the plasma lipidome contributes to the genetic basis of T2D is unknown. Here, we investigated the potential overlap in the genetic basis of the plasma lipidome and T2D-related traits.

Results

We used plasma lipidomic data (1202 pedigreed individuals, 319 lipid species representing 23 lipid classes) from San Antonio Family Heart Study in Mexican Americans. Bivariate trait analyses were used to estimate the genetic and environmental correlation of all lipid species with three T2D-related traits: risk of T2D, presence of prediabetes and homeostatic model of assessment – insulin resistance. We found that 44 lipid species were significantly genetically correlated with one or more of the three T2D-related traits. Majority of these lipid species belonged to the diacylglycerol (DAG, 17 species) and triacylglycerol (TAG, 17 species) classes. Six lipid species (all belonging to the triacylglycerol class and containing palmitate at the first position) were significantly genetically correlated with all the T2D-related traits.

Conclusions

Our results imply that: a) not all plasma lipid species are genetically informative for T2D pathogenesis; b) the DAG and TAG lipid classes partially share genetic basis of T2D; and c) 1-palmitate containing TAGs may provide additional insights into the genetic basis of T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0515-5) contains supplementary material, which is available to authorized users.

Fast Genome-Wide QTL Association Mapping on Pedigree and Population Data

Since most analysis software for genome-wide association studies (GWAS) currently exploit only unrelated individuals, there is a need for efficient applications that can handle general pedigree data or mixtures of both population and pedigree data. Even datasets thought to consist of only unrelated individuals may include cryptic relationships that can lead to false positives if not discovered and controlled for. In addition, family designs possess compelling advantages. They are better equipped to detect rare variants, control for population stratification, and facilitate the study of parent-of-origin effects. Pedigrees selected for extreme trait values often segregate a single gene with strong effect. Finally, many pedigrees are available as an important legacy from the era of linkage analysis. Unfortunately, pedigree likelihoods are notoriously hard to compute. In this paper, we reexamine the computational bottlenecks and implement ultra-fast pedigree-based GWAS analysis. Kinship coefficients can either be based on explicitly provided pedigrees or automatically estimated from dense markers. Our strategy (a) works for random sample data, pedigree data, or a mix of both; (b) entails no loss of power; (c) allows for any number of covariate adjustments, including correction for population stratification; (d) allows for testing SNPs under additive, dominant, and recessive models; and (e) accommodates both univariate and multivariate quantitative traits. On a typical personal computer (six CPU cores at 2.67 GHz), analyzing a univariate HDL (high-density lipoprotein) trait from the San Antonio Family Heart Study (935,392 SNPs on 1,388 individuals in 124 pedigrees) takes less than 2 min and 1.5 GB of memory. Complete multivariate QTL analysis of the three time-points of the longitudinal HDL multivariate trait takes less than 5 min and 1.5 GB of memory. The algorithm is implemented as the Ped-GWAS Analysis (Option 29) in the Mendel statistical genetics package, which is freely available for Macintosh, Linux, and Windows platforms from http://genetics.ucla.edu/software/mendel.

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.

Independent test assessment using the extreme value distribution theory

The new generation of whole genome sequencing platforms offers great possibilities and challenges for dissecting the genetic basis of complex traits. With a very high number of sequence variants, a naïve multiple hypothesis threshold correction hinders the identification of reliable associations by the overreduction of statistical power. In this report, we examine 2 alternative approaches to improve the statistical power of a whole genome association study to detect reliable genetic associations. The approaches were tested using the Genetic Analysis Workshop 19 (GAW19) whole genome sequencing data. The first tested method estimates the real number of effective independent tests actually being performed in whole genome association project by the use of an extreme value distribution and a set of phenotype simulations. Given the familiar nature of the GAW19 data and the finite number of pedigree founders in the sample, the number of correlations between genotypes is greater than in a set of unrelated samples. Using our procedure, we estimate that the effective number represents only 15 % of the total number of independent tests performed. However, even using this corrected significance threshold, no genome-wide significant association could be detected for systolic and diastolic blood pressure traits.

The second approach implements a biological relevance-driven hypothesis tested by exploiting prior computational predictions on the effect of nonsynonymous genetic variants detected in a whole genome sequencing association study. This guided testing approach was able to identify 2 promising single-nucleotide polymorphisms (SNPs), 1 for each trait, targeting biologically relevant genes that could help shed light on the genesis of the human hypertension. The first gene, PFH14, associated with systolic blood pressure, interacts directly with genes involved in calcium-channel formation and the second gene, MAP4, encodes a microtubule-associated protein and had already been detected by previous genome-wide association study experiments conducted in an Asian population. Our results highlight the necessity of the development of alternative approached to improve the efficiency on the detection of reasonable candidate associations in whole genome sequencing studies.

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.

Omics-squared: human genomic, transcriptomic and phenotypic data for genetic analysis workshop 19

Background

The Genetic Analysis Workshops (GAW) are a forum for development, testing, and comparison of statistical genetic methods and software. Each contribution to the workshop includes an application to a specified data set. Here we describe the data distributed for GAW19, which focused on analysis of human genomic and transcriptomic data.

Methods

GAW19 data were donated by the T2D-GENES Consortium and the San Antonio Family Heart Study and included whole genome and exome sequences for odd-numbered autosomes, measures of gene expression, systolic and diastolic blood pressures, and related covariates in two Mexican American samples. These two samples were a collection of 20 large families with whole genome sequence and transcriptomic data and a set of 1943 unrelated individuals with exome sequence. For each sample, simulated phenotypes were constructed based on the real sequence data. ‘Functional’ genes and variants for the simulations were chosen based on observed correlations between gene expression and blood pressure. The simulations focused primarily on additive genetic models but also included a genotype-by-medication interaction. A total of 245 genes were designated as ‘functional’ in the simulations with a few genes of large effect and most genes explaining < 1 % of the trait variation. An additional phenotype, Q1, was simulated to be correlated among related individuals, based on theoretical or empirical kinship matrices, but was not associated with any sequence variants. Two hundred replicates of the phenotypes were simulated. The GAW19 data are an expansion of the data used at GAW18, which included the family-based whole genome sequence, blood pressure, and simulated phenotypes, but not the gene expression data or the set of 1943 unrelated individuals with exome sequence.

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.

Lipidomic risk score independently and cost-effectively predicts risk of future type 2 diabetes: results from diverse cohorts

Background

Detection of type 2 diabetes (T2D) is routinely based on the presence of dysglycemia. Although disturbed lipid metabolism is a hallmark of T2D, the potential of plasma lipidomics as a biomarker of future T2D is unknown. Our objective was to develop and validate a plasma lipidomic risk score (LRS) as a biomarker of future type 2 diabetes and to evaluate its cost-effectiveness for T2D screening.

Methods

Plasma LRS, based on significantly associated lipid species from an array of 319 lipid species, was developed in a cohort of initially T2D-free individuals from the San Antonio Family Heart Study (SAFHS). The LRS derived from SAFHS as well as its recalibrated version were validated in an independent cohort from Australia – the AusDiab cohort. The participants were T2D-free at baseline and followed for 9197 person-years in the SAFHS cohort (n = 771) and 5930 person-years in the AusDiab cohort (n = 644). Statistically and clinically improved T2D prediction was evaluated with established statistical parameters in both cohorts. Modeling studies were conducted to determine whether the use of LRS would be cost-effective for T2D screening. The main outcome measures included accuracy and incremental value of the LRS over routinely used clinical predictors of T2D risk; validation of these results in an independent cohort and cost-effectiveness of including LRS in screening/intervention programs for T2D.

Results

The LRS was based on plasma concentration of dihydroceramide 18:0, lysoalkylphosphatidylcholine 22:1 and triacyglycerol 16:0/18:0/18:1. The score predicted future T2D independently of prediabetes with an accuracy of 76 %. Even in the subset of initially euglycemic individuals, the LRS improved T2D prediction. In the AusDiab cohort, the LRS continued to predict T2D significantly and independently. When combined with risk-stratification methods currently used in clinical practice, the LRS significantly improved the model fit (p < 0.001), information content (p < 0.001), discrimination (p < 0.001) and reclassification (p < 0.001) in both cohorts. Modeling studies demonstrated that LRS-based risk-stratification combined with metformin supplementation for high-risk individuals was the most cost-effective strategy for T2D prevention.

Conclusions

Considering the novelty, incremental value and cost-effectiveness of LRS it should be used for risk-stratification of future T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12944-016-0234-3) contains supplementary material, which is available to authorized users.

Soluble Forms of Intercellular and Vascular Cell Adhesion Molecules Independently Predict Progression to Type 2 Diabetes in Mexican American Families

Objective

While the role of type 2 diabetes (T2D) in inducing endothelial dysfunction is fairly well-established the etiological role of endothelial dysfunction in the onset of T2D is still a matter of debate. In the light of conflicting evidence in this regard, we conducted a prospective study to determine the association of circulating levels of soluble intercellular adhesion molecule 1 (sICAM-1) and soluble vessel cell adhesion molecule 1 (sVCAM-1) with incident T2D.

Methods

Data from this study came from 1,269 Mexican Americans of whom 821 initially T2D-free individuals were longitudinally followed up in the San Antonio Family Heart Study. These individuals were followed for 9752.95 person-years for development of T2D. Prospective association of sICAM-1 and sVCAM-1 with incident T2D was studied using Kaplan-Meier survival plots and mixed effects Cox proportional hazards modeling to account for relatedness among study participants. Incremental value of adhesion molecule biomarkers was studied using integrated discrimination improvement (IDI) and net reclassification improvement (NRI) indexes.

Results

Decreasing median values for serum concentrations of sICAM-1 and sVCAM-1 were observed in the following groups in this order: individuals with T2D at baseline, individuals who developed T2D during follow-up, individuals with prediabetes at baseline and normal glucose tolerant (NGT) individuals who remained T2D-free during follow-up. Top quartiles for sICAM-1 and sVCAM-1 were strongly and significantly associated with homeostatic model of assessment—insulin resistance (HOMA-IR). Mixed effects Cox proportional hazards modeling revealed that after correcting for important clinical confounders, high sICAM-1 and sVCAM-1 concentrations were associated with 2.52 and 1.99 times faster progression to T2D as compared to low concentrations, respectively. Individuals with high concentrations for both sICAM-1 and sVCAM-1 progressed to T2D 3.42 times faster than those with low values for both sICAM-1 and sVCAM-1. The results were similar in women in reproductive age group and the remainder of the cohort. Inclusion of sICAM-1 and sVCAM-1 in predictive models significantly improved reclassification and discrimination. The majority of these results were seen even when the analyses were restricted to NGT individuals.

Conclusion

Serum concentrations of sICAM-1 and sVCAM-1 independently and additively predict future T2D and represent important candidate biomarkers of T2D.

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.

Lack of Association between SLC30A8 Variants and Type 2 Diabetes in Mexican American Families

SLC30A8 encodes zinc transporter 8 which is involved in packaging and release of insulin. Evidence for the association of SLC30A8 variants with type 2 diabetes (T2D) is inconclusive. We interrogated single nucleotide polymorphisms (SNPs) around SLC30A8 for association with T2D in high-risk, pedigreed individuals from extended Mexican American families. This study of 118 SNPs within 50 kb of the SLC30A8 locus tested the association with eight T2D-related traits at four levels: (i) each SNP using measured genotype approach (MGA); (ii) interaction of SNPs with age and sex; (iii) combinations of SNPs using Bayesian Quantitative Trait Nucleotide (BQTN) analyses; and (iv) entire gene locus using the gene burden test. Only one SNP (rs7817754) was significantly associated with incident T2D but a summary statistic based on all T2D-related traits identified 11 novel SNPs. Three SNPs and one SNP were weakly but interactively associated with age and sex, respectively. BQTN analyses could not demonstrate any informative combination of SNPs over MGA. Lastly, gene burden test results showed that at best the SLC30A8 locus could account for only 1-2% of the variability in T2D-related traits. Our results indicate a lack of association of the SLC30A8 SNPs with T2D in Mexican American families.

Genome-wide significant linkage of schizophrenia-related neuroanatomical trait to 12q24

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.

The Quantitative-MFG Test: A Linear Mixed Effect Model to Detect Maternal-Offspring Gene Interactions

Maternal-offspring gene interactions, aka maternal-fetal genotype (MFG) incompatibilities, are neglected in complex diseases and quantitative trait studies. They are implicated in birth to adult onset diseases but there are limited ways to investigate their influence on quantitative traits. We present the quantitative-MFG (QMFG) test, a linear mixed model where maternal and offspring genotypes are fixed effects and residual correlations between family members are random effects. The QMFG handles families of any size, common or general scenarios of MFG incompatibility, and additional covariates. We develop likelihood ratio tests (LRTs) and rapid score tests and show they provide correct inference. In addition, the LRT's alternative model provides unbiased parameter estimates. We show that testing the association of SNPs by fitting a standard model, which only considers the offspring genotypes, has very low power or can lead to incorrect conclusions. We also show that offspring genetic effects are missed if the MFG modeling assumptions are too restrictive. With genome-wide association study data from the San Antonio Family Heart Study, we demonstrate that the QMFG score test is an effective and rapid screening tool. The QMFG test therefore has important potential to identify pathways of complex diseases for which the genetic etiology remains to be discovered.

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.

Sequencing strategies and characterization of 721 vervet monkey genomes for future genetic analyses of medically relevant traits

Background

We report here the first genome-wide high-resolution polymorphism resource for non-human primate (NHP) association and linkage studies, constructed for the Caribbean-origin vervet monkey, or African green monkey (Chlorocebus aethiops sabaeus), one of the most widely used NHPs in biomedical research. We generated this resource by whole genome sequencing (WGS) of monkeys from the Vervet Research Colony (VRC), an NIH-supported research resource for which extensive phenotypic data are available.

Results

We identified genome-wide single nucleotide polymorphisms (SNPs) by WGS of 721 members of an extended pedigree from the VRC. From high-depth WGS data we identified more than 4 million polymorphic unequivocal segregating sites; by pruning these SNPs based on heterozygosity, quality control filters, and the degree of linkage disequilibrium (LD) between SNPs, we constructed genome-wide panels suitable for genetic association (about 500,000 SNPs) and linkage analysis (about 150,000 SNPs). To further enhance the utility of these resources for linkage analysis, we used a further pruned subset of the linkage panel to generate multipoint identity by descent matrices.

Conclusions

The genetic and phenotypic resources now available for the VRC and other Caribbean-origin vervets enable their use for genetic investigation of traits relevant to human diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0152-2) contains supplementary material, which is available to authorized users.

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.

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.

Data for Genetic Analysis Workshop 18: human whole genome sequence, blood pressure, and simulated phenotypes in extended pedigrees

Genetic Analysis Workshop 18 (GAW18) focused on identification of genes and functional variants that influence complex phenotypes in human sequence data. Data for the workshop were donated by the T2D-GENES Consortium and included whole genome sequences for odd-numbered autosomes in 464 key individuals selected from 20 Mexican American families, a dense set of single-nucleotide polymorphisms in 959 individuals in these families, and longitudinal data on systolic and diastolic blood pressure measured at 1-4 examinations over a period of 20 years. Simulated phenotypes were generated based on the real sequence data and pedigree structures. In the design of the simulation model, gene expression measures from the San Antonio Family Heart Study (not distributed as part of the GAW18 data) were used to identify genes whose mRNA levels were correlated with blood pressure. Observed variants within these genes were designated as functional in the GAW18 simulation if they were nonsynonymous and predicted to have deleterious effects on protein function or if they were noncoding and associated with mRNA levels. Two simulated longitudinal phenotypes were modeled to have the same trait distributions as the real systolic and diastolic blood pressure data, with effects of age, sex, and medication use, including a genotype-medication interaction. For each phenotype, more than 1000 sequence variants in more than 200 genes present on the odd-numbered autosomes individually explained less than 0.01-2.78% of phenotypic variance. Cumulatively, variants in the most influential gene explained 7.79% of trait variance. An additional simulated phenotype, Q1, was designed to be correlated among family members but to not be associated with any sequence variants. Two hundred replicates of the phenotypes were simulated, with each including data for 849 individuals.

Evaluation of estimated genetic values and their application to genome-wide investigation of systolic blood pressure

The concept of breeding values, an individual's phenotypic deviation from the population mean as a result of the sum of the average effects of the genes they carry, is of great importance in livestock, aquaculture, and cash crop industries where emphasis is placed on an individual's potential to pass desirable phenotypes on to the next generation. As breeding or genetic values (as referred to here) cannot be measured directly, estimated genetic values (EGVs) are based on an individual's own phenotype, phenotype information from relatives, and, increasingly, genetic data. Because EGVs represent additive genetic variation, calculating EGVs in an extended human pedigree is expected to provide a more refined phenotype for genetic analyses. To test the utility of EGVs in genome-wide association, EGVs were calculated for 847 members of 20 extended Mexican American families based on 100 replicates of simulated systolic blood pressure. Calculations were performed in GAUSS to solve a variation on the standard Best Linear Unbiased Predictor (BLUP) mixed model equation with age, sex, and the first 3 principal components of sample-wide genetic variability as fixed effects and the EGV as a random effect distributed around the relationship matrix. Three methods of calculating kinship were considered: expected kinship from pedigree relationships, empirical kinship from common variants, and empirical kinship from both rare and common variants. Genome-wide association analysis was conducted on simulated phenotypes and EGVs using the additive measured genotype approach in the SOLAR software package. The EGV-based approach showed only minimal improvement in power to detect causative loci.

Increased waist circumference is independently associated with hypothyroidism in Mexican Americans: replicative evidence from two large, population-based studies

BACKGROUND

Mexican Americans are at an increased risk of both thyroid dysfunction and metabolic syndrome (MS). Thus it is conceivable that some components of the MS may be associated with the risk of thyroid dysfunction in these individuals. Our objective was to investigate and replicate the potential association of MS traits with thyroid dysfunction in Mexican Americans.

METHODS

We conducted association testing for 18 MS traits in two large studies on Mexican Americans - the San Antonio Family Heart Study (SAFHS) and the National Health and Nutrition Examination Survey (NHANES) 2007-10. A total of 907 participants from 42 families in SAFHS and 1633 unrelated participants from NHANES 2007-10 were included in this study. The outcome measures were prevalence of clinical and subclinical hypothyroidism and thyroid function index (TFI) - a measure of thyroid function. For the SAFHS, we used polygenic regression analyses with multiple covariates to test associations in setting of family studies. For the NHANES 2007-10, we corrected for the survey design variables as needed for association analyses in survey data. In both datasets, we corrected for age, sex and their linear and quadratic interactions.

RESULTS

TFI was an accurate indicator of clinical thyroid status (area under the receiver-operating-characteristic curve to detect clinical hypothyroidism, 0.98) in both SAFHS and NHANES 2007-10. Of the 18 MS traits, waist circumference (WC) showed the most consistent association with TFI in both studies independently of age, sex and body mass index (BMI). In the SAFHS and NHANES 2007-10 datasets, each standard deviation increase in WC was associated with 0.13 (p < 0.001) and 0.11 (p < 0.001) unit increase in the TFI, respectively. In a series of polygenic and linear regression models, central obesity (defined as WC ≥ 102 cm in men and ≥88 cm in women) was associated with clinical and subclinical hypothyroidism independent of age, sex, BMI and type 2 diabetes in both datasets. Estimated prevalence of hypothyroidism was consistently high in those with central obesity, especially below 45y of age.

CONCLUSIONS

WC independently associates with increased risk of thyroid dysfunction. Use of WC to identify Mexican American subjects at high risk of thyroid dysfunction should be investigated in future studies.

Novel QTL at chromosome 6p22 for alcohol consumption: Implications for the genetic liability of alcohol use disorders

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

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.

Common genetic variants and gene expression associated with white matter microstructure in the human brain

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

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.

Gene-by-age effects on BMI from birth to adulthood: the Fels Longitudinal Study

OBJECTIVES

Genome wide association studies have shown 32 loci to influence BMI in European-American adults but replication in other studies is inconsistent and may be attributed to gene-by-age effects. The aims of this study were to determine if the influence of the summed risk score of these 32 loci (GRS) on BMI differed across age from birth to 40 years, and to determine if additive genetic effects other than those in the GRS differed by age.

METHODS

Serial measures of BMI were calculated at 0, 1, 3, 6, 9, 12, 18, and 28 months, and 4, 7, 11, 15, 19, 23, 30, and 40 years for 1,176 (605 females, 571 males) European-American participants in the Fels Longitudinal Study. SOLAR was used for genetic analyses.

RESULTS

GRS was significant (P < 0.05) at ages: 6, 9 months, 4-15 years, and 23-40 years. Remaining additive genetic effects independently influenced BMI (P < 5.3 × 10(-5) , 0.40 < h(2) < 0.76). Some genetic correlations between ages were not significant. Differential GRS effects did not retain significance after multiple comparisons adjustments.

CONCLUSIONS

While well-known BMI variants do not appear to have significant differential effects, other additive genes differ over the lifespan.

Plasma dihydroceramide species associate with waist circumference in Mexican American families

OBJECTIVE

Waist circumference (WC), the clinical marker of central obesity, is gaining popularity as a screening tool for type 2 diabetes (T2D). While there is epidemiologic evidence favoring the WC-T2D association, its biological substantiation is generally weak. Our objective was to determine the independent association of plasma lipid repertoire with WC.

METHODS

Samples and data from the San Antonio Family Heart Study of 1208 Mexican Americans from 42 extended families were used. Association of plasma lipidomic profiles with the cross-sectionally assessed WC was determined. Plasma lipidomic profiling entailed liquid chromatography with mass spectrometry. Statistical analyses included multivariable polygenic regression models and bivariate trait analyses using the SOLAR software.

RESULTS

After adjusting for age and sex interactions, body mass index, homeostasis model of assessment-insulin resistance, total cholesterol, triglycerides, high density lipoproteins and use of lipid lowering drugs, dihydroceramides as a class were associated with WC. Dihydroceramide species 18:0, 20:0, 22:0, and 24:1 were significantly associated and genetically correlated with WC. Two sphingomyelin species (31:1 and 41:1) were also associated with WC.

CONCLUSIONS

Plasma dihydroceramide levels independently associate with WC. Thus, high resolution plasma lipidomic studies can provide further credence to the biological underpinnings of the association of WC with T2D.

Genome-wide significant localization for working and spatial memory: Identifying genes for psychosis using models of cognition

It is well established that risk for developing psychosis is largely mediated by the influence of genes, but identifying precisely which genes underlie that risk has been problematic. Focusing on endophenotypes, rather than illness risk, is one solution to this problem. Impaired cognition is a well-established endophenotype of psychosis. Here we aimed to characterize the genetic architecture of cognition using phenotypically detailed models as opposed to relying on general IQ or individual neuropsychological measures. In so doing we hoped to identify genes that mediate cognitive ability, which might also contribute to psychosis risk. Hierarchical factor models of genetically clustered cognitive traits were subjected to linkage analysis followed by QTL region-specific association analyses in a sample of 1,269 Mexican American individuals from extended pedigrees. We identified four genome wide significant QTLs, two for working and two for spatial memory, and a number of plausible and interesting candidate genes. The creation of detailed models of cognition seemingly enhanced the power to detect genetic effects on cognition and provided a number of possible candidate genes for psychosis.

Waist circumference is genetically correlated with incident Type 2 diabetes in Mexican-American families

AIMS

We aimed to determine the genetic and environmental correlation between various anthropometric indexes and incident Type 2 diabetes with a focus on waist circumference.

METHODS

We used the data on extended Mexican-American families (808 subjects, 7617.92 person-years follow-up) from the San Antonio Family Heart Study and estimated the genetic and environmental correlations of 16 anthropometric indexes with the genetic liability of incident Type 2 diabetes. We performed bivariate trait analyses using the solar software package.

RESULTS

All 16 anthropometric indexes were significantly heritable (range of heritabilities 0.24-0.99). Thirteen indexes were found to have significant environmental correlation with the liability of incident Type 2 diabetes. In contrast, only anthropometric indexes consisting of waist circumference (waist circumference, waist-hip ratio and waist-height ratio) were significantly genetically correlated (genetic correlation coefficients: 0.45, 0.55 and 0.44, respectively) with the liability of incident Type 2 diabetes. We did not observe such a correlation for BMI.

CONCLUSIONS

Waist circumference as a predictor of future Type 2 diabetes is supported by the finding that they share common genetic influences.

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.