Laura J. Scott

Christian Fuchsberger, Gonçalo Abecasis and colleagues describe a new web-based imputation service that enables rapid imputation of large numbers of samples and allows convenient access to large reference panels of sequenced individuals. Their state space reduction provides a computationally efficient solution for genotype imputation with no loss in imputation accuracy. Genotype imputation is a key component of genetic association studies, where it increases power, facilitates meta-analysis, and aids interpretation of signals. Genotype imputation is computationally demanding and, with current tools, typically requires access to a high-performance computing cluster and to a reference panel of sequenced genomes. Here we describe improvements to imputation machinery that reduce computational requirements by more than an order of magnitude with no loss of accuracy in comparison to standard imputation tools. We also describe a new web-based service for imputation that facilitates access to new reference panels and greatly improves user experience and productivity.

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

Identifying the genetic variants that increase the risk of type 2 diabetes (T2D) in humans has been a formidable challenge. Adopting a genome-wide association strategy, we genotyped 1161 Finnish T2D cases and 1174 Finnish normal glucose-tolerant (NGT) controls with &gt;315,000 single-nucleotide polymorphisms (SNPs) and imputed genotypes for an additional &gt;2 million autosomal SNPs. We carried out association analysis with these SNPs to identify genetic variants that predispose to T2D, compared our T2D association results with the results of two similar studies, and genotyped 80 SNPs in an additional 1215 Finnish T2D cases and 1258 Finnish NGT controls. We identify T2D-associated variants in an intergenic region of chromosome 11p12, contribute to the identification of T2D-associated variants near the genes IGF2BP2 and CDKAL1 and the region of CDKN2A and CDKN2B , and confirm that variants near TCF7L2 , SLC30A8 , HHEX , FTO , PPARG , and KCNJ11 are associated with T2D risk. This brings the number of T2D loci now confidently identified to at least 10.

We describe a reference panel of 64,976 human haplotypes at 39,235,157 SNPs constructed using whole-genome sequence data from 20 studies of predominantly European ancestry. Using this resource leads to accurate genotype imputation at minor allele frequencies as low as 0.1% and a large increase in the number of SNPs tested in association studies, and it can help to discover and refine causal loci. We describe remote server resources that allow researchers to carry out imputation and phasing consistently and efficiently.

Hutchinson-Gilford progeria syndrome (HGPS) is a rare genetic disorder characterized by features reminiscent of marked premature ageing. Here, we present evidence of mutations in lamin A (LMNA) as the cause of this disorder. The HGPS gene was initially localized to chromosome 1q by observing two cases of uniparental isodisomy of 1q-the inheritance of both copies of this material from one parent-and one case with a 6-megabase paternal interstitial deletion. Sequencing of LMNA, located in this interval and previously implicated in several other heritable disorders, revealed that 18 out of 20 classical cases of HGPS harboured an identical de novo (that is, newly arisen and not inherited) single-base substitution, G608G(GGC > GGT), within exon 11. One additional case was identified with a different substitution within the same codon. Both of these mutations result in activation of a cryptic splice site within exon 11, resulting in production of a protein product that deletes 50 amino acids near the carboxy terminus. Immunofluorescence of HGPS fibroblasts with antibodies directed against lamin A revealed that many cells show visible abnormalities of the nuclear membrane. The discovery of the molecular basis of this disease may shed light on the general phenomenon of human ageing.

Mark McCarthy, Michael Boehnke, Andrew Morris and colleagues perform large-scale association analyses using the Metabochip to gain insights into the genetic architecture of type 2 diabetes. They report several new susceptibility loci, including two that show sex-differentiated effects on disease risk. To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip, including 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two showing sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of additional common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signaling and cell cycle regulation, in diabetes pathogenesis.

Genome-wide association (GWA) studies have identified multiple loci at which common variants modestly but reproducibly influence risk of type 2 diabetes (T2D). Established associations to common and rare variants explain only a small proportion of the heritability of T2D. As previously published analyses had limited power to identify variants with modest effects, we carried out meta-analysis of three T2D GWA scans comprising 10,128 individuals of European descent and approximately 2.2 million SNPs (directly genotyped and imputed), followed by replication testing in an independent sample with an effective sample size of up to 53,975. We detected at least six previously unknown loci with robust evidence for association, including the JAZF1 (P = 5.0 x 10(-14)), CDC123-CAMK1D (P = 1.2 x 10(-10)), TSPAN8-LGR5 (P = 1.1 x 10(-9)), THADA (P = 1.1 x 10(-9)), ADAMTS9 (P = 1.2 x 10(-8)) and NOTCH2 (P = 4.1 x 10(-8)) gene regions. Our results illustrate the value of large discovery and follow-up samples for gaining further insights into the inherited basis of T2D.

By combining genome-wide association data from 8,130 individuals with type 2 diabetes (T2D) and 38,987 controls of European descent and following up previously unidentified meta-analysis signals in a further 34,412 cases and 59,925 controls, we identified 12 new T2D association signals with combined P<5x10(-8). These include a second independent signal at the KCNQ1 locus; the first report, to our knowledge, of an X-chromosomal association (near DUSP9); and a further instance of overlap between loci implicated in monogenic and multifactorial forms of diabetes (at HNF1A). The identified loci affect both beta-cell function and insulin action, and, overall, T2D association signals show evidence of enrichment for genes involved in cell cycle regulation. We also show that a high proportion of T2D susceptibility loci harbor independent association signals influencing apparently unrelated complex traits.

To identify genetic variants influencing plasma lipid concentrations, we first used genotype imputation and meta-analysis to combine three genome-wide scans totaling 8,816 individuals and comprising 6,068 individuals specific to our study (1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables) and 2,758 individuals from the Diabetes Genetics Initiative, reported in a companion study in this issue. We subsequently examined promising signals in 11,569 additional individuals. Overall, we identify strongly associated variants in eleven loci previously implicated in lipid metabolism (ABCA1, the APOA5-APOA4-APOC3-APOA1 and APOE-APOC clusters, APOB, CETP, GCKR, LDLR, LPL, LIPC, LIPG and PCSK9) and also in several newly identified loci (near MVK-MMAB and GALNT2, with variants primarily associated with high-density lipoprotein (HDL) cholesterol; near SORT1, with variants primarily associated with low-density lipoprotein (LDL) cholesterol; near TRIB1, MLXIPL and ANGPTL3, with variants primarily associated with triglycerides; and a locus encompassing several genes near NCAN, with variants strongly associated with both triglycerides and LDL cholesterol). Notably, the 11 independent variants associated with increased LDL cholesterol concentrations in our study also showed increased frequency in a sample of coronary artery disease cases versus controls.

We conducted a combined genome-wide association study (GWAS) of 7,481 individuals with bipolar disorder (cases) and 9,250 controls as part of the Psychiatric GWAS Consortium. Our replication study tested 34 SNPs in 4,496 independent cases with bipolar disorder and 42,422 independent controls and found that 18 of 34 SNPs had P < 0.05, with 31 of 34 SNPs having signals with the same direction of effect (P = 3.8 × 10(-7)). An analysis of all 11,974 bipolar disorder cases and 51,792 controls confirmed genome-wide significant evidence of association for CACNA1C and identified a new intronic variant in ODZ4. We identified a pathway comprised of subunits of calcium channels enriched in bipolar disorder association intervals. Finally, a combined GWAS analysis of schizophrenia and bipolar disorder yielded strong association evidence for SNPs in CACNA1C and in the region of NEK4-ITIH1-ITIH3-ITIH4. Our replication results imply that increasing sample sizes in bipolar disorder will confirm many additional loci.

Blood low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglyceride levels are risk factors for cardiovascular disease. To dissect the polygenic basis of these traits, we conducted genome-wide association screens in 19,840 individuals and replication in up to 20,623 individuals. We identified 30 distinct loci associated with lipoprotein concentrations (each with P < 5 x 10(-8)), including 11 loci that reached genome-wide significance for the first time. The 11 newly defined loci include common variants associated with LDL cholesterol near ABCG8, MAFB, HNF1A and TIMD4; with HDL cholesterol near ANGPTL4, FADS1-FADS2-FADS3, HNF4A, LCAT, PLTP and TTC39B; and with triglycerides near AMAC1L2, FADS1-FADS2-FADS3 and PLTP. The proportion of individuals exceeding clinical cut points for high LDL cholesterol, low HDL cholesterol and high triglycerides varied according to an allelic dosage score (P < 10(-15) for each trend). These results suggest that the cumulative effect of multiple common variants contributes to polygenic dyslipidemia.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10⁻⁹ to P = 1.8 × 10⁻⁴⁰) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10⁻³ to P = 1.2 × 10⁻¹³). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

To identify previously unknown genetic loci associated with fasting glucose concentrations, we examined the leading association signals in ten genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95% CI = 0.06-0.08) mmol/l in fasting glucose levels (P = 3.2 x 10(-50)) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P = 1.1 x 10(-15)). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05-1.12), per G allele P = 3.3 x 10(-7)) in a meta-analysis of 13 case-control studies totaling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P = 1.1 x 10(-57)) and GCK (rs4607517, P = 1.0 x 10(-25)) loci.

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P < 5 × 10-8), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N = 38,580) informative for adult waist circumference (WC) and waist-hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P = 1.9x10(-11)) and MSRA (WC, P = 8.9x10(-9)). A third locus, near LYPLAL1, was associated with WHR in women only (P = 2.6x10(-8)). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.

Identifying genetic variants that influence human height will advance our understanding of skeletal growth and development. Several rare genetic variants have been convincingly and reproducibly associated with height in mendelian syndromes, and common variants in the transcription factor gene HMGA2 are associated with variation in height in the general population. Here we report genome-wide association analyses, using genotyped and imputed markers, of 6,669 individuals from Finland and Sardinia, and follow-up analyses in an additional 28,801 individuals. We show that common variants in the osteoarthritis-associated locus GDF5-UQCC contribute to variation in height with an estimated additive effect of 0.44 cm (overall P < 10(-15)). Our results indicate that there may be a link between the genetic basis of height and osteoarthritis, potentially mediated through alterations in bone growth and development.

We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3–50, P < 2.14 × 10−6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-d-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach. Whole-exome sequencing identifies ten risk genes for schizophrenia implicated by rare protein-coding variants, a subset of which overlap with risk genes in other neurodevelopmental disorders.

In this trans-ethnic multi-omic study, we reinterpret the genetic architecture of blood pressure to identify genes, tissues, phenomes and medication contexts of blood pressure homeostasis. We discovered 208 novel common blood pressure SNPs and 53 rare variants in genome-wide association studies of systolic, diastolic and pulse pressure in up to 776,078 participants from the Million Veteran Program (MVP) and collaborating studies, with analysis of the blood pressure clinical phenome in MVP. Our transcriptome-wide association study detected 4,043 blood pressure associations with genetically predicted gene expression of 840 genes in 45 tissues, and mouse renal single-cell RNA sequencing identified upregulated blood pressure genes in kidney tubule cells.

The extent to which variation in chromatin structure and transcription factor binding may influence gene expression, and thus underlie or contribute to variation in phenotype, is unknown. To address this question, we cataloged both individual-to-individual variation and differences between homologous chromosomes within the same individual (allele-specific variation) in chromatin structure and transcription factor binding in lymphoblastoid cells derived from individuals of geographically diverse ancestry. Ten percent of active chromatin sites were individual-specific; a similar proportion were allele-specific. Both individual-specific and allele-specific sites were commonly transmitted from parent to child, which suggests that they are heritable features of the human genome. Our study shows that heritable chromatin status and transcription factor binding differ as a result of genetic variation and may underlie phenotypic variation in humans.

Bipolar disorder (BP) is a disabling and often life-threatening disorder that affects approximately 1% of the population worldwide. To identify genetic variants that increase the risk of BP, we genotyped on the Illumina HumanHap550 Beadchip 2,076 bipolar cases and 1,676 controls of European ancestry from the National Institute of Mental Health Human Genetics Initiative Repository, and the Prechter Repository and samples collected in London, Toronto, and Dundee. We imputed SNP genotypes and tested for SNP-BP association in each sample and then performed meta-analysis across samples. The strongest association P value for this 2-study meta-analysis was 2.4 x 10(-6). We next imputed SNP genotypes and tested for SNP-BP association based on the publicly available Affymetrix 500K genotype data from the Wellcome Trust Case Control Consortium for 1,868 BP cases and a reference set of 12,831 individuals. A 3-study meta-analysis of 3,683 nonoverlapping cases and 14,507 extended controls on >2.3 M genotyped and imputed SNPs resulted in 3 chromosomal regions with association P approximately 10(-7): 1p31.1 (no known genes), 3p21 (>25 known genes), and 5q15 (MCTP1). The most strongly associated nonsynonymous SNP rs1042779 (OR = 1.19, P = 1.8 x 10(-7)) is in the ITIH1 gene on chromosome 3, with other strongly associated nonsynonymous SNPs in GNL3, NEK4, and ITIH3. Thus, these chromosomal regions harbor genes implicated in cell cycle, neurogenesis, neuroplasticity, and neurosignaling. In addition, we replicated the reported ANK3 association results for SNP rs10994336 in the nonoverlapping GSK sample (OR = 1.37, P = 0.042). Although these results are promising, analysis of additional samples will be required to confirm that variant(s) in these regions influence BP risk.

Insulin secretion has a crucial role in glucose homeostasis, and failure to secrete sufficient insulin is a hallmark of type 2 diabetes. Genome-wide association studies (GWAS) have identified loci contributing to insulin processing and secretion; however, a substantial fraction of the genetic contribution remains undefined. To examine low-frequency (minor allele frequency (MAF) 0.5-5%) and rare (MAF < 0.5%) nonsynonymous variants, we analyzed exome array data in 8,229 nondiabetic Finnish males using the Illumina HumanExome Beadchip. We identified low-frequency coding variants associated with fasting proinsulin concentrations at the SGSM2 and MADD GWAS loci and three new genes with low-frequency variants associated with fasting proinsulin or insulinogenic index: TBC1D30, KANK1 and PAM. We also show that the interpretation of single-variant and gene-based tests needs to consider the effects of noncoding SNPs both nearby and megabases away. This study demonstrates that exome array genotyping is a valuable approach to identify low-frequency variants that contribute to complex traits.

OBJECTIVE Recent genome-wide association studies have revealed loci associated with glucose and insulin-related traits. We aimed to characterize 19 such loci using detailed measures of insulin processing, secretion, and sensitivity to help elucidate their role in regulation of glucose control, insulin secretion and/or action. RESEARCH DESIGN AND METHODS We investigated associations of loci identified by the Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) with circulating proinsulin, measures of insulin secretion and sensitivity from oral glucose tolerance tests (OGTTs), euglycemic clamps, insulin suppression tests, or frequently sampled intravenous glucose tolerance tests in nondiabetic humans (n = 29,084). RESULTS The glucose-raising allele in MADD was associated with abnormal insulin processing (a dramatic effect on higher proinsulin levels, but no association with insulinogenic index) at extremely persuasive levels of statistical significance (P = 2.1 x 10(-71)). Defects in insulin processing and insulin secretion were seen in glucose-raising allele carriers at TCF7L2, SCL30A8, GIPR, and C2CD4B. Abnormalities in early insulin secretion were suggested in glucose-raising allele carriers at MTNR1B, GCK, FADS1, DGKB, and PROX1 (lower insulinogenic index; no association with proinsulin or insulin sensitivity). Two loci previously associated with fasting insulin (GCKR and IGF1) were associated with OGTT-derived insulin sensitivity indices in a consistent direction. CONCLUSIONS Genetic loci identified through their effect on hyperglycemia and/or hyperinsulinemia demonstrate considerable heterogeneity in associations with measures of insulin processing, secretion, and sensitivity. Our findings emphasize the importance of detailed physiological characterization of such loci for improved understanding of pathways associated with alterations in glucose homeostasis and eventually type 2 diabetes.

Common diseases such as type 2 diabetes are phenotypically heterogeneous. Obesity is a major risk factor for type 2 diabetes, but patients vary appreciably in body mass index. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m²) compared to obese cases (BMI≥30 Kg/m²). We performed two case-control genome-wide studies using two accepted cut-offs for defining individuals as overweight or obese. We used 2,112 lean type 2 diabetes cases (BMI<25 kg/m²) or 4,123 obese cases (BMI≥30 kg/m²), and 54,412 un-stratified controls. Replication was performed in 2,881 lean cases or 8,702 obese cases, and 18,957 un-stratified controls. To assess the effects of known signals, we tested the individual and combined effects of SNPs representing 36 type 2 diabetes loci. After combining data from discovery and replication datasets, we identified two signals not previously reported in Europeans. A variant (rs8090011) in the LAMA1 gene was associated with type 2 diabetes in lean cases (P = 8.4×10⁻⁹, OR = 1.13 [95% CI 1.09-1.18]), and this association was stronger than that in obese cases (P = 0.04, OR = 1.03 [95% CI 1.00-1.06]). A variant in HMG20A--previously identified in South Asians but not Europeans--was associated with type 2 diabetes in obese cases (P = 1.3×10⁻⁸, OR = 1.11 [95% CI 1.07-1.15]), although this association was not significantly stronger than that in lean cases (P = 0.02, OR = 1.09 [95% CI 1.02-1.17]). For 36 known type 2 diabetes loci, 29 had a larger odds ratio in the lean compared to obese (binomial P = 0.0002). In the lean analysis, we observed a weighted per-risk allele OR = 1.13 [95% CI 1.10-1.17], P = 3.2×10⁻¹⁴. This was larger than the same model fitted in the obese analysis where the OR = 1.06 [95% CI 1.05-1.08], P = 2.2×10⁻¹⁶. This study provides evidence that stratification of type 2 diabetes cases by BMI may help identify additional risk variants and that lean cases may have a stronger genetic predisposition to type 2 diabetes.

More than 90% of people who attempt suicide have a psychiatric diagnosis; however, twin and family studies suggest that the genetic etiology of suicide attempt is partially distinct from that of the psychiatric disorders themselves. The authors present the largest genome-wide association study (GWAS) on suicide attempt, using cohorts of individuals with major depressive disorder, bipolar disorder, and schizophrenia from the Psychiatric Genomics Consortium.The samples comprised 1,622 suicide attempters and 8,786 nonattempters with major depressive disorder; 3,264 attempters and 5,500 nonattempters with bipolar disorder; and 1,683 attempters and 2,946 nonattempters with schizophrenia. A GWAS on suicide attempt was performed by comparing attempters to nonattempters with each disorder, followed by a meta-analysis across disorders. Polygenic risk scoring was used to investigate the genetic relationship between suicide attempt and the psychiatric disorders.Three genome-wide significant loci for suicide attempt were found: one associated with suicide attempt in major depressive disorder, one associated with suicide attempt in bipolar disorder, and one in the meta-analysis of suicide attempt in mood disorders. These associations were not replicated in independent mood disorder cohorts from the UK Biobank and iPSYCH. No significant associations were found in the meta-analysis of all three disorders. Polygenic risk scores for major depression were significantly associated with suicide attempt in major depressive disorder (R2=0.25%), bipolar disorder (R2=0.24%), and schizophrenia (R2=0.40%).This study provides new information on genetic associations and demonstrates that genetic liability for major depression increases risk for suicide attempt across psychiatric disorders. Further collaborative efforts to increase sample size may help to robustly identify genetic associations and provide biological insights into the etiology of suicide attempt.

Significance The majority of genetic variants associated with type 2 diabetes (T2D) are located outside of genes in noncoding regions that may regulate gene expression in disease-relevant tissues, like pancreatic islets. Here, we present the largest integrated analysis to date of high-resolution, high-throughput human islet molecular profiling data to characterize the genome (DNA), epigenome (DNA packaging), and transcriptome (gene expression). We find that T2D genetic variants are enriched in regions of the genome where transcription Regulatory Factor X (RFX) is predicted to bind in an islet-specific manner. Genetic variants that increase T2D risk are predicted to disrupt RFX binding, providing a molecular mechanism to explain how the genome can influence the epigenome, modulating gene expression and ultimately T2D risk.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power. Exome-wide sequencing studies of populations in Finland identified 26 deleterious alleles associated with 64 quantitative traits that are clinically relevant to cardiovascular and metabolic diseases.

Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10−8) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits. Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10−8) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits.

We integrate comeasured gene expression and DNA methylation (DNAme) in 265 human skeletal muscle biopsies from the FUSION study with >7 million genetic variants and eight physiological traits: height, waist, weight, waist-hip ratio, body mass index, fasting serum insulin, fasting plasma glucose, and type 2 diabetes. We find hundreds of genes and DNAme sites associated with fasting insulin, waist, and body mass index, as well as thousands of DNAme sites associated with gene expression (eQTM). We find that controlling for heterogeneity in tissue/muscle fiber type reduces the number of physiological trait associations, and that long-range eQTMs (>1 Mb) are reduced when controlling for tissue/muscle fiber type or latent factors. We map genetic regulators (quantitative trait loci; QTLs) of expression (eQTLs) and DNAme (mQTLs). Using Mendelian randomization (MR) and mediation techniques, we leverage these genetic maps to predict 213 causal relationships between expression and DNAme, approximately two-thirds of which predict methylation to causally influence expression. We use MR to integrate FUSION mQTLs, FUSION eQTLs, and GTEx eQTLs for 48 tissues with genetic associations for 534 diseases and quantitative traits. We identify hundreds of genes and thousands of DNAme sites that may drive the reported disease/quantitative trait genetic associations. We identify 300 gene expression MR associations that are present in both FUSION and GTEx skeletal muscle and that show stronger evidence of MR association in skeletal muscle than other tissues, which may partially reflect differences in power across tissues. As one example, we find that increased RXRA muscle expression may decrease lean tissue mass.

<h3>Importance</h3> Most previous genome-wide association studies (GWAS) of depression have used data from individuals of European descent. This limits the understanding of the underlying biology of depression and raises questions about the transferability of findings between populations. <h3>Objective</h3> To investigate the genetics of depression among individuals of East Asian and European descent living in different geographic locations, and with different outcome definitions for depression. <h3>Design, Setting, and Participants</h3> Genome-wide association analyses followed by meta-analysis, which included data from 9 cohort and case-control data sets comprising individuals with depression and control individuals of East Asian descent. This study was conducted between January 2019 and May 2021. <h3>Exposures</h3> Associations of genetic variants with depression risk were assessed using generalized linear mixed models and logistic regression. The results were combined across studies using fixed-effects meta-analyses. These were subsequently also meta-analyzed with the largest published GWAS for depression among individuals of European descent. Additional meta-analyses were carried out separately by outcome definition (clinical depression vs symptom-based depression) and region (East Asian countries vs Western countries) for East Asian ancestry cohorts. <h3>Main Outcomes and Measures</h3> Depression status was defined based on health records and self-report questionnaires. <h3>Results</h3> There were a total of 194 548 study participants (approximate mean age, 51.3 years; 62.8% women). Participants included 15 771 individuals with depression and 178 777 control individuals of East Asian descent. Five novel associations were identified, including 1 in the meta-analysis for broad depression among those of East Asian descent: rs4656484 (β = −0.018, SE = 0.003,<i>P</i> = 4.43x10⁻⁸) at 1q24.1. Another locus at 7p21.2 was associated in a meta-analysis restricted to geographically East Asian studies (β = 0.028, SE = 0.005,<i>P</i> = 6.48x10⁻⁹for rs10240457). The lead variants of these 2 novel loci were not associated with depression risk in European ancestry cohorts (β = −0.003, SE = 0.005,<i>P</i> = .53 for rs4656484 and β = −0.005, SE = 0.004,<i>P</i> = .28 for rs10240457). Only 11% of depression loci previously identified in individuals of European descent reached nominal significance levels in the individuals of East Asian descent. The transancestry genetic correlation between cohorts of East Asian and European descent for clinical depression was<i>r</i> = 0.413 (SE = 0.159). Clinical depression risk was negatively genetically correlated with body mass index in individuals of East Asian descent (<i>r</i> = −0.212, SE = 0.084), contrary to findings for individuals of European descent. <h3>Conclusions and Relevance</h3> These results support caution against generalizing findings about depression risk factors across populations and highlight the need to increase the ancestral and geographic diversity of samples with consistent phenotyping.

We report results from the Bipolar Exome (BipEx) collaboration analysis of whole-exome sequencing of 13,933 patients with bipolar disorder (BD) matched with 14,422 controls. We find an excess of ultra-rare protein-truncating variants (PTVs) in patients with BD among genes under strong evolutionary constraint in both major BD subtypes. We find enrichment of ultra-rare PTVs within genes implicated from a recent schizophrenia exome meta-analysis (SCHEMA; 24,248 cases and 97,322 controls) and among binding targets of CHD8. Genes implicated from genome-wide association studies (GWASs) of BD, however, are not significantly enriched for ultra-rare PTVs. Combining gene-level results with SCHEMA, AKAP11 emerges as a definitive risk gene (odds ratio (OR) = 7.06, P = 2.83 × 10-9). At the protein level, AKAP-11 interacts with GSK3B, the hypothesized target of lithium, a primary treatment for BD. Our results lend support to BD's polygenicity, demonstrating a role for rare coding variation as a significant risk factor in BD etiology.

Few studies have explored the impact of rare variants (minor allele frequency < 1%) on highly heritable plasma metabolites identified in metabolomic screens. The Finnish population provides an ideal opportunity for such explorations, given the multiple bottlenecks and expansions that have shaped its history, and the enrichment for many otherwise rare alleles that has resulted. Here, we report genetic associations for 1391 plasma metabolites in 6136 men from the late-settlement region of Finland. We identify 303 novel association signals, more than one third at variants rare or enriched in Finns. Many of these signals identify genes not previously implicated in metabolite genome-wide association studies and suggest mechanisms for diseases and disease-related traits.

Background Sex differences in incidence and/or presentation of schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP) are pervasive. Previous evidence for shared genetic risk and sex differences in brain abnormalities across disorders suggest possible shared sex-dependent genetic risk. Methods We conducted the largest to date genome-wide genotype-by-sex (G×S) interaction of risk for these disorders using 85,735 cases (33,403 SCZ, 19,924 BIP, and 32,408 MDD) and 109,946 controls from the PGC (Psychiatric Genomics Consortium) and iPSYCH. Results Across disorders, genome-wide significant single nucleotide polymorphism–by-sex interaction was detected for a locus encompassing NKAIN2 (rs117780815, p = 3.2 × 10−8), which interacts with sodium/potassium-transporting ATPase (adenosine triphosphatase) enzymes, implicating neuronal excitability. Three additional loci showed evidence (p < 1 × 10−6) for cross-disorder G×S interaction (rs7302529, p = 1.6 × 10−7; rs73033497, p = 8.8 × 10−7; rs7914279, p = 6.4 × 10−7), implicating various functions. Gene-based analyses identified G×S interaction across disorders (p = 8.97 × 10−7) with transcriptional inhibitor SLTM. Most significant in SCZ was a MOCOS gene locus (rs11665282, p = 1.5 × 10−7), implicating vascular endothelial cells. Secondary analysis of the PGC-SCZ dataset detected an interaction (rs13265509, p = 1.1 × 10−7) in a locus containing IDO2, a kynurenine pathway enzyme with immunoregulatory functions implicated in SCZ, BIP, and MDD. Pathway enrichment analysis detected significant G×S interaction of genes regulating vascular endothelial growth factor receptor signaling in MDD (false discovery rate-corrected p < .05). Conclusions In the largest genome-wide G×S analysis of mood and psychotic disorders to date, there was substantial genetic overlap between the sexes. However, significant sex-dependent effects were enriched for genes related to neuronal development and immune and vascular functions across and within SCZ, BIP, and MDD at the variant, gene, and pathway levels.

Abstract Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

The Finland-United States Investigation Of NIDDM Genetics (FUSION) study aims to identify genetic variants that predispose to type 2 diabetes by studying affected sibling pair families from Finland. Chromosome 20 showed our strongest initial evidence for linkage. It currently has a maximum logarithm of odds (LOD) score of 2.48 at 70 cM in a set of 495 families. In this study, we searched for diabetes susceptibility variant(s) at 20q13 by genotyping single nucleotide polymorphism (SNP) markers in case and control DNA pools. Of 291 SNPs successfully typed in a 7.5-Mb interval, the strongest association confirmed by individual genotyping was with SNP rs2144908, located 1.3 kb downstream of the primary β-cell promoter P2 of hepatocyte nuclear factor-4α (HNF4A). This SNP showed association with diabetes disease status (odds ratio [OR] 1.33, 95% CI 1.06–1.65, P = 0.011) and with several diabetes-related traits. Most of the evidence for linkage at 20q13 could be attributed to the families carrying the risk allele. We subsequently found nine additional associated SNPs spanning a 64-kb region, including the P2 and P1 promoters and exons 1–3. Our results and the independent observation of association of SNPs near the P2 promoter with diabetes in a separate study population of Ashkenazi Jewish origin suggests that variant(s) located near or within HNF4A increases susceptibility to type 2 diabetes.

Transcription factor 7-like 2 (TCF7L2) is part of the Wnt signaling pathway. Genetic variants within TCF7L2 on chromosome 10q were recently reported to be associated with type 2 diabetes in Icelandic, Danish, and American (U.S.) samples. We previously observed a modest logarithm of odds score of 0.61 on chromosome 10q, ∼1 Mb from TCF7L2, in the Finland-United States Investigation of NIDDM Genetics study. We tested the five associated TCF7L2 single nucleotide polymorphism (SNP) variants in a Finnish sample of 1,151 type 2 diabetic patients and 953 control subjects. We confirmed the association with the same risk allele (P value &amp;lt;0.05) for all five SNPs. Our strongest results were for rs12255372 (odds ratio [OR] 1.36 [95% CI 1.15–1.61], P = 0.00026) and rs7903146 (1.33 [1.14–1.56], P = 0.00042). Based on the CEU HapMap data, we selected and tested 12 additional SNPs to tag SNPs in linkage disequilibrium with rs12255372. None of these SNPs showed stronger evidence of association than rs12255372 or rs7903146 (OR ≤1.26, P ≥ 0.0054). Our results strengthen the evidence that one or more variants in TCF7L2 are associated with increased risk of type 2 diabetes.

Identifying cis-regulatory elements is important to understanding how human pancreatic islets modulate gene expression in physiologic or pathophysiologic (e.g., diabetic) conditions. We conducted genome-wide analysis of DNase I hypersensitive sites, histone H3 lysine methylation modifications (K4me1, K4me3, K79me2), and CCCTC factor (CTCF) binding in human islets. This identified ∼18,000 putative promoters (several hundred unannotated and islet-active). Surprisingly, active promoter modifications were absent at genes encoding islet-specific hormones, suggesting a distinct regulatory mechanism. Of 34,039 distal (nonpromoter) regulatory elements, 47% are islet unique and 22% are CTCF bound. In the 18 type 2 diabetes (T2D)-associated loci, we identified 118 putative regulatory elements and confirmed enhancer activity for 12 of 33 tested. Among six regulatory elements harboring T2D-associated variants, two exhibit significant allele-specific differences in activity. These findings present a global snapshot of the human islet epigenome and should provide functional context for noncoding variants emerging from genetic studies of T2D and other islet disorders.

Epidemiological studies consistently show that circulating sex hormone binding globulin (SHBG) levels are lower in type 2 diabetes patients than non-diabetic individuals, but the causal nature of this association is controversial. Genetic studies can help dissect causal directions of epidemiological associations because genotypes are much less likely to be confounded, biased or influenced by disease processes. Using this Mendelian randomization principle, we selected a common single nucleotide polymorphism (SNP) near the SHBG gene, rs1799941, that is strongly associated with SHBG levels. We used data from this SNP, or closely correlated SNPs, in 27 657 type 2 diabetes patients and 58 481 controls from 15 studies. We then used data from additional studies to estimate the difference in SHBG levels between type 2 diabetes patients and controls. The SHBG SNP rs1799941 was associated with type 2 diabetes [odds ratio (OR) 0.94, 95% CI: 0.91, 0.97; P = 2 × 10−5], with the SHBG raising allele associated with reduced risk of type 2 diabetes. This effect was very similar to that expected (OR 0.92, 95% CI: 0.88, 0.96), given the SHBG-SNP versus SHBG levels association (SHBG levels are 0.2 standard deviations higher per copy of the A allele) and the SHBG levels versus type 2 diabetes association (SHBG levels are 0.23 standard deviations lower in type 2 diabetic patients compared to controls). Results were very similar in men and women. There was no evidence that this variant is associated with diabetes-related intermediate traits, including several measures of insulin secretion and resistance. Our results, together with those from another recent genetic study, strengthen evidence that SHBG and sex hormones are involved in the aetiology of type 2 diabetes.

Peroxisome proliferator-activated receptor gamma (PPARG) is a master transcriptional regulator of adipocyte differentiation and a canonical target of antidiabetic thiazolidinedione medications. In rare families, loss-of-function (LOF) mutations in PPARG are known to cosegregate with lipodystrophy and insulin resistance; in the general population, the common P12A variant is associated with a decreased risk of type 2 diabetes (T2D). Whether and how rare variants in PPARG and defects in adipocyte differentiation influence risk of T2D in the general population remains undetermined. By sequencing PPARG in 19,752 T2D cases and controls drawn from multiple studies and ethnic groups, we identified 49 previously unidentified, nonsynonymous PPARG variants (MAF < 0.5%). Considered in aggregate (with or without computational prediction of functional consequence), these rare variants showed no association with T2D (OR = 1.35; P = 0.17). The function of the 49 variants was experimentally tested in a novel high-throughput human adipocyte differentiation assay, and nine were found to have reduced activity in the assay. Carrying any of these nine LOF variants was associated with a substantial increase in risk of T2D (OR = 7.22; P = 0.005). The combination of large-scale DNA sequencing and functional testing in the laboratory reveals that approximately 1 in 1,000 individuals carries a variant in PPARG that reduces function in a human adipocyte differentiation assay and is associated with a substantial risk of T2D.

In genome-wide association studies of binary traits, investigators typically use logistic regression to test common variants for disease association within studies, and combine association results across studies using meta-analysis. For common variants, logistic regression tests are well calibrated, and meta-analysis of study-specific association results is only slightly less powerful than joint analysis of the combined individual-level data. In recent sequencing and dense chip based association studies, investigators increasingly test low-frequency variants for disease association. In this paper, we seek to (1) identify the association test with maximal power among tests with well controlled type I error rate and (2) compare the relative power of joint and meta-analysis tests. We use analytic calculation and simulation to compare the empirical type I error rate and power of four logistic regression based tests: Wald, score, likelihood ratio, and Firth bias-corrected. We demonstrate for low-count variants (roughly minor allele count [MAC] < 400) that: (1) for joint analysis, the Firth test has the best combination of type I error and power; (2) for meta-analysis of balanced studies (equal numbers of cases and controls), the score test is best, but is less powerful than Firth test based joint analysis; and (3) for meta-analysis of sufficiently unbalanced studies, all four tests can be anti-conservative, particularly the score test. We also establish MAC as the key parameter determining test calibration for joint and meta-analysis.

A detailed understanding of the genome-wide variability of single-nucleotide germline mutation rates is essential to studying human genome evolution. Here, we use ~36 million singleton variants from 3560 whole-genome sequences to infer fine-scale patterns of mutation rate heterogeneity. Mutability is jointly affected by adjacent nucleotide context and diverse genomic features of the surrounding region, including histone modifications, replication timing, and recombination rate, sometimes suggesting specific mutagenic mechanisms. Remarkably, GC content, DNase hypersensitivity, CpG islands, and H3K36 trimethylation are associated with both increased and decreased mutation rates depending on nucleotide context. We validate these estimated effects in an independent dataset of ~46,000 de novo mutations, and confirm our estimates are more accurate than previously published results based on ancestrally older variants without considering genomic features. Our results thus provide the most refined portrait to date of the factors contributing to genome-wide variability of the human germline mutation rate.

Gene set enrichment testing can enhance the biological interpretation of ChIP-seq data. Here, we develop a method, ChIP-Enrich, for this analysis which empirically adjusts for gene locus length (the length of the gene body and its surrounding non-coding sequence). Adjustment for gene locus length is necessary because it is often positively associated with the presence of one or more peaks and because many biologically defined gene sets have an excess of genes with longer or shorter gene locus lengths. Unlike alternative methods, ChIP-Enrich can account for the wide range of gene locus length-to-peak presence relationships (observed in ENCODE ChIP-seq data sets). We show that ChIP-Enrich has a well-calibrated type I error rate using permuted ENCODE ChIP-seq data sets; in contrast, two commonly used gene set enrichment methods, Fisher's exact test and the binomial test implemented in Genomic Regions Enrichment of Annotations Tool (GREAT), can have highly inflated type I error rates and biases in ranking. We identify DNA-binding proteins, including CTCF, JunD and glucocorticoid receptor α (GRα), that show different enrichment patterns for peaks closer to versus further from transcription start sites. We also identify known and potential new biological functions of GRα. ChIP-Enrich is available as a web interface (http://chip-enrich.med.umich.edu) and Bioconductor package.

Genome-wide association studies have revealed numerous risk loci associated with diverse diseases. However, identification of disease-causing variants within association loci remains a major challenge. Divergence in gene expression due to cis-regulatory variants in noncoding regions is central to disease susceptibility. We show that integrative computational analysis of phylogenetic conservation with a complexity assessment of co-occurring transcription factor binding sites (TFBS) can identify cis-regulatory variants and elucidate their mechanistic role in disease. Analysis of established type 2 diabetes risk loci revealed a striking clustering of distinct homeobox TFBS. We identified the PRRX1 homeobox factor as a repressor of PPARG2 expression in adipose cells and demonstrate its adverse effect on lipid metabolism and systemic insulin sensitivity, dependent on the rs4684847 risk allele that triggers PRRX1 binding. Thus, cross-species conservation analysis at the level of co-occurring TFBS provides a valuable contribution to the translation of genetic association signals to disease-related molecular mechanisms.

Abstract Type 2 diabetes (T2D) results from the combined effects of genetic and environmental factors on multiple tissues over time. Of the &gt;100 variants associated with T2D and related traits in genome-wide association studies (GWAS), &gt;90% occur in non-coding regions, suggesting a strong regulatory component to T2D risk. Here to understand how T2D status, metabolic traits and genetic variation influence gene expression, we analyse skeletal muscle biopsies from 271 well-phenotyped Finnish participants with glucose tolerance ranging from normal to newly diagnosed T2D. We perform high-depth strand-specific mRNA-sequencing and dense genotyping. Computational integration of these data with epigenome data, including ATAC-seq on skeletal muscle, and transcriptome data across diverse tissues reveals that the tissue-specific genetic regulatory architecture of skeletal muscle is highly enriched in muscle stretch/super enhancers, including some that overlap T2D GWAS variants. In one such example, T2D risk alleles residing in a muscle stretch/super enhancer are linked to increased expression and alternative splicing of muscle-specific isoforms of ANK1 .

<h3>Importance</h3> Complex disorders, such as bipolar disorder (BD), likely result from the influence of both common and rare susceptibility alleles. While common variation has been widely studied, rare variant discovery has only recently become feasible with next-generation sequencing. <h3>Objective</h3> To utilize a combined family-based and case-control approach to exome sequencing in BD using multiplex families as an initial discovery strategy, followed by association testing in a large case-control meta-analysis. <h3>Design, Setting, and Participants</h3> We performed exome sequencing of 36 affected members with BD from 8 multiplex families and tested rare, segregating variants in 3 independent case-control samples consisting of 3541 BD cases and 4774 controls. <h3>Main Outcomes and Measures</h3> We used penalized logistic regression and 1-sided gene-burden analyses to test for association of rare, segregating damaging variants with BD. Permutation-based analyses were performed to test for overall enrichment with previously identified gene sets. <h3>Results</h3> We found 84 rare (frequency &lt;1%), segregating variants that were bioinformatically predicted to be damaging. These variants were found in 82 genes that were enriched for gene sets previously identified in de novo studies of autism (19 observed vs. 10.9 expected,<i>P</i> = .0066) and schizophrenia (11 observed vs. 5.1 expected,<i>P</i> = .0062) and for targets of the fragile X mental retardation protein (FMRP) pathway (10 observed vs. 4.4 expected,<i>P</i> = .0076). The case-control meta-analyses yielded 19 genes that were nominally associated with BD based either on individual variants or a gene-burden approach. Although no gene was individually significant after correction for multiple testing, this group of genes continued to show evidence for significant enrichment of de novo autism genes (6 observed vs 2.6 expected,<i>P</i> = .028). <h3>Conclusions and Relevance</h3> Our results are consistent with the presence of prominent locus and allelic heterogeneity in BD and suggest that very large samples will be required to definitively identify individual rare variants or genes conferring risk for this disorder. However, we also identify significant associations with gene sets composed of previously discovered de novo variants in autism and schizophrenia, as well as targets of the FRMP pathway, providing preliminary support for the overlap of potential autism and schizophrenia risk genes with rare, segregating variants in families with BD.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Excessive alcohol consumption is one of the main causes of death and disability worldwide. Alcohol consumption is a heritable complex trait. Here we conducted a meta-analysis of genome-wide association studies of alcohol consumption (g d-1) from the UK Biobank, the Alcohol Genome-Wide Consortium and the Cohorts for Heart and Aging Research in Genomic Epidemiology Plus consortia, collecting data from 480,842 people of European descent to decipher the genetic architecture of alcohol intake. We identified 46 new common loci and investigated their potential functional importance using magnetic resonance imaging data and gene expression studies. We identify genetic pathways associated with alcohol consumption and suggest genetic mechanisms that are shared with neuropsychiatric disorders such as schizophrenia.

Doxorubicin (also called Adriamycin) is effective in treating a wide range of human cancers and currently considered as one of the most important drugs in cancer chemotherapeutics. The clinical use of doxorubicin is, however, associated with dosage-dependent cardiotoxicity and development of heart failure, which diminish the therapeutic index of this widely used anticancer drug. This article first surveys key research findings on doxorubicin redox biology that may impact its cardiotoxicity as well as anticancer activity. It then discusses emerging concepts, especially the topoisomerase IIb?p53?mitochondrion axis that may lead to the development of mechanistically based novel strategies to protect against cardiotoxicity and enhance the effectiveness of doxorubicin therapy.

COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain.In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry.Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression.Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Cigarette smoking and poor glycemic control are risk factors for diabetic nephropathy in type 1 diabetes. However, the specifics of the relation of these risk factors to the onset of this complication have not been elucidated. To investigate these issues, we followed for 4 years 943 Joslin Clinic patients aged 15-44 years with type 1 diabetes who had normoalbuminuria during the 2-year baseline period. Microalbuminuria developed in 109 of the 943 individuals, giving an incidence rate of 3.3/100 person-years. The risk of onset of microalbuminuria was predicted somewhat more precisely by the measurements during the 1st and 2nd years preceding onset than by all the measurements during the longer (4-year) interval, suggesting attenuation of the impact of past hyperglycemia over time. Point estimates of the incidence rate (per 100 person-years) according to quartiles of HbA(1c) during the 1st and 2nd years preceding the outcome were 1.3 in the 1st, 1.5 in the 2nd, 3.1 in the 3rd, and 6.9 in the 4th (P = 1.3 x 10(-9)). Point estimates of the incidence rate (per 100 person-years) according to smoking status were 7.9 for current smokers, 1.8 for past smokers, and 2.2 for those who had never smoked (P = 2.0 x 10(-7)). In a multiple logistic model, the independent effects of HbA(1c) level and cigarette smoking remained highly significant, but their magnitudes were reduced. Using the same covariates in a generalized additive model, we examined the shape of the relationship between HbA(1c) and onset of microalbuminuria and found significant nonlinearity in the logarithm of odds scale (P = 0.04). The slope was steeper with HbA(1c) >8% than <8%. Furthermore, the change in the slope was magnified among current smokers. In conclusion, patients with type 1 diabetes who smoke and have an HbA(1c) >8% have the highest risk of onset of microalbuminuria.

Abstract Genome‐wide association (GWA) studies require genotyping hundreds of thousands of markers on thousands of subjects, and are expensive at current genotyping costs. To conserve resources, many GWA studies are adopting a staged design in which a proportion of the available samples are genotyped on all markers in stage 1, and a proportion of these markers are genotyped on the remaining samples in stage 2. We describe a strategy for designing cost‐effective two‐stage GWA studies. Our strategy preserves much of the power of the corresponding one‐stage design and minimizes the genotyping cost of the study while allowing for differences in per genotyping cost between stages 1 and 2. We show that the ratio of stage 2 to stage 1 per genotype cost can strongly influence both the optimal design and the genotyping cost of the study. Increasing the stage 2 per genotype cost shifts more of the genotyping and study cost to stage 1, and increases the cost of the study. This higher cost can be partially mitigated by adopting a design with reduced power while preserving the false positive rate or by increasing the false positive rate while preserving power. For example, reducing the power preserved in the two‐stage design from 99 to 95% that of the one‐stage design decreases the two‐stage study cost by ∼15%. Alternatively, the same cost savings can be had by relaxing the false positive rate by 2.5‐fold, for example from 1/300,000 to 2.5/300,000, while retaining the same power. Genet. Epidemiol . 2007. © 2007 Wiley‐Liss, Inc.

More than 120 published reports have described associations between single nucleotide polymorphisms (SNPs) and type 2 diabetes. However, multiple studies of the same variant have often been discordant. From a literature search, we identified previously reported type 2 diabetes–associated SNPs. We initially genotyped 134 SNPs on 786 index case subjects from type 2 diabetes families and 617 control subjects with normal glucose tolerance from Finland and excluded from analysis 20 SNPs in strong linkage disequilibrium (r2 &amp;gt; 0.8) with another typed SNP. Of the 114 SNPs examined, we followed up the 20 most significant SNPs (P &amp;lt; 0.10) on an additional 384 case subjects and 366 control subjects from a population-based study in Finland. In the combined data, we replicated association (P &amp;lt; 0.05) for 12 SNPs: PPARG Pro12Ala and His447, KCNJ11 Glu23Lys and rs5210, TNF −857, SLC2A2 Ile110Thr, HNF1A/TCF1 rs2701175 and GE117881_360, PCK1 −232, NEUROD1 Thr45Ala, IL6 −598, and ENPP1 Lys121Gln. The replication of 12 SNPs of 114 tested was significantly greater than expected by chance under the null hypothesis of no association (P = 0.012). We observed that SNPs from genes that had three or more previous reports of association were significantly more likely to be replicated in our sample (P = 0.03), although we also replicated 4 of 58 SNPs from genes that had only one previous report of association.

Type 2 diabetes is a common complex disorder with environmental and genetic components. We used a candidate gene-based approach to identify single nucleotide polymorphism (SNP) variants in 222 candidate genes that influence susceptibility to type 2 diabetes.In a case-control study of 1,161 type 2 diabetic subjects and 1,174 control Finns who are normal glucose tolerant, we genotyped 3,531 tagSNPs and annotation-based SNPs and imputed an additional 7,498 SNPs, providing 99.9% coverage of common HapMap variants in the 222 candidate genes. Selected SNPs were genotyped in an additional 1,211 type 2 diabetic case subjects and 1,259 control subjects who are normal glucose tolerant, also from Finland.Using SNP- and gene-based analysis methods, we replicated previously reported SNP-type 2 diabetes associations in PPARG, KCNJ11, and SLC2A2; identified significant SNPs in genes with previously reported associations (ENPP1 [rs2021966, P = 0.00026] and NRF1 [rs1882095, P = 0.00096]); and implicated novel genes, including RAPGEF1 (rs4740283, P = 0.00013) and TP53 (rs1042522, Arg72Pro, P = 0.00086), in type 2 diabetes susceptibility.Our study provides an effective gene-based approach to association study design and analysis. One or more of the newly implicated genes may contribute to type 2 diabetes pathogenesis. Analysis of additional samples will be necessary to determine their effect on susceptibility.

Identifying the genetic variants that regulate fasting glucose concentrations may further our understanding of the pathogenesis of diabetes. We therefore investigated the association of fasting glucose levels with SNPs in 2 genome-wide scans including a total of 5,088 nondiabetic individuals from Finland and Sardinia. We found a significant association between the SNP rs563694 and fasting glucose concentrations (P = 3.5 x 10(-7)). This association was further investigated in an additional 18,436 nondiabetic individuals of mixed European descent from 7 different studies. The combined P value for association in these follow-up samples was 6.9 x 10(-26), and combining results from all studies resulted in an overall P value for association of 6.4 x 10(-33). Across these studies, fasting glucose concentrations increased 0.01-0.16 mM with each copy of the major allele, accounting for approximately 1% of the total variation in fasting glucose. The rs563694 SNP is located between the genes glucose-6-phosphatase catalytic subunit 2 (G6PC2) and ATP-binding cassette, subfamily B (MDR/TAP), member 11 (ABCB11). Our results in combination with data reported in the literature suggest that G6PC2, a glucose-6-phosphatase almost exclusively expressed in pancreatic islet cells, may underlie variation in fasting glucose, though it is possible that ABCB11, which is expressed primarily in liver, may also contribute to such variation.

Common variants in the transcription factor 7-like 2 (TCF7L2) gene have been identified as the strongest genetic risk factors for type 2 diabetes (T2D). However, the mechanisms by which these non-coding variants increase risk for T2D are not well-established. We used 13 expression assays to survey mRNA expression of multiple TCF7L2 splicing forms in up to 380 samples from eight types of human tissue (pancreas, pancreatic islets, colon, liver, monocytes, skeletal muscle, subcutaneous adipose tissue and lymphoblastoid cell lines) and observed a tissue-specific pattern of alternative splicing. We tested whether the expression of TCF7L2 splicing forms was associated with single nucleotide polymorphisms (SNPs), rs7903146 and rs12255372, located within introns 3 and 4 of the gene and most strongly associated with T2D. Expression of two splicing forms was lower in pancreatic islets with increasing counts of T2D-associated alleles of the SNPs: a ubiquitous splicing form (P = 0.018 for rs7903146 and P = 0.020 for rs12255372) and a splicing form found in pancreatic islets, pancreas and colon but not in other tissues tested here (P = 0.009 for rs12255372 and P = 0.053 for rs7903146). Expression of this form in glucose-stimulated pancreatic islets correlated with expression of proinsulin (r(2) = 0.84-0.90, P < 0.00063). In summary, we identified a tissue-specific pattern of alternative splicing of TCF7L2. After adjustment for multiple tests, no association between expression of TCF7L2 in eight types of human tissue samples and T2D-associated genetic variants remained significant. Alternative splicing of TCF7L2 in pancreatic islets warrants future studies. GenBank Accession Numbers: FJ010164-FJ010174.

Identifying differentially expressed (DE) genes from RNA sequencing (RNAseq) studies is among the most common analyses in genomics. However, RNAseq DE analysis presents several statistical and computational challenges, including over-dispersed read counts and, in some settings, sample non-independence. Previous count-based methods rely on simple hierarchical Poisson models (e.g. negative binomial) to model independent over-dispersion, but do not account for sample non-independence due to relatedness, population structure and/or hidden confounders. Here, we present a Poisson mixed model with two random effects terms that account for both independent over-dispersion and sample non-independence. We also develop a scalable sampling-based inference algorithm using a latent variable representation of the Poisson distribution. With simulations, we show that our method properly controls for type I error and is generally more powerful than other widely used approaches, except in small samples (n <15) with other unfavorable properties (e.g. small effect sizes). We also apply our method to three real datasets that contain related individuals, population stratification or hidden confounders. Our results show that our method increases power in all three data compared to other approaches, though the power gain is smallest in the smallest sample (n = 6). Our method is implemented in MACAU, freely available at www.xzlab.org/software.html.

Molecular mechanisms remain unknown for most type 2 diabetes genome-wide association study identified loci. Variants associated with type 2 diabetes and fasting glucose levels reside in introns of ADCY5, a gene that encodes adenylate cyclase 5. Adenylate cyclase 5 catalyzes the production of cyclic AMP, which is a second messenger molecule involved in cell signaling and pancreatic β-cell insulin secretion. We demonstrated that type 2 diabetes risk alleles are associated with decreased ADCY5 expression in human islets and examined candidate variants for regulatory function. rs11708067 overlaps a predicted enhancer region in pancreatic islets. The type 2 diabetes risk rs11708067-A allele showed fewer H3K27ac ChIP-seq reads in human islets, lower transcriptional activity in reporter assays in rodent β-cells (rat 832/13 and mouse MIN6), and increased nuclear protein binding compared with the rs11708067-G allele. Homozygous deletion of the orthologous enhancer region in 832/13 cells resulted in a 64% reduction in expression level of Adcy5, but not adjacent gene Sec22a, and a 39% reduction in insulin secretion. Together, these data suggest that rs11708067-A risk allele contributes to type 2 diabetes by disrupting an islet enhancer, which results in reduced ADCY5 expression and impaired insulin secretion.

Detecting and estimating DNA sample contamination are important steps to ensure high-quality genotype calls and reliable downstream analysis. Existing methods rely on population allele frequency information for accurate estimation of contamination rates. Correctly specifying population allele frequencies for each individual in early stage of sequence analysis is impractical or even impossible for large-scale sequencing centers that simultaneously process samples from multiple studies across diverse populations. On the other hand, incorrectly specified allele frequencies may result in substantial bias in estimated contamination rates. For example, we observed that existing methods often fail to identify 10% contaminated samples at a typical 3% contamination exclusion threshold when genetic ancestry is misspecified. Such an incomplete screening of contaminated samples substantially inflates the estimated rate of genotyping errors even in deeply sequenced genomes and exomes. We propose a robust statistical method that accurately estimates DNA contamination and is agnostic to genetic ancestry of the intended or contaminating sample. Our method integrates the estimation of genetic ancestry and DNA contamination in a unified likelihood framework by leveraging individual-specific allele frequencies projected from reference genotypes onto principal component coordinates. Our method can also be used for estimating genetic ancestries, similar to LASER or TRACE, but simultaneously accounting for potential contamination. We demonstrate that our method robustly estimates contamination rates and genetic ancestries across populations and contamination scenarios. We further demonstrate that, in the presence of contamination, genetic ancestry inference can be substantially biased with existing methods that ignore contamination, while our method corrects for such biases.

Abstract In genetic association analysis, a joint test of multiple distinct phenotypes can increase power to identify sets of trait‐associated variants within genes or regions of interest. Existing multiphenotype tests for rare variants make specific assumptions about the patterns of association with underlying causal variants, and the violation of these assumptions can reduce power to detect association. Here, we develop a general framework for testing pleiotropic effects of rare variants on multiple continuous phenotypes using multivariate kernel regression (Multi‐SKAT). Multi‐SKAT models affect sizes of variants on the phenotypes through a kernel matrix and perform a variance component test of association. We show that many existing tests are equivalent to specific choices of kernel matrices with the Multi‐SKAT framework. To increase power of detecting association across tests with different kernel matrices, we developed a fast and accurate approximation of the significance of the minimum observed P value across tests. To account for related individuals, our framework uses random effects for the kinship matrix. Using simulated data and amino acid and exome‐array data from the METabolic Syndrome In Men (METSIM) study, we show that Multi‐SKAT can improve power over single‐phenotype SKAT‐O test and existing multiple‐phenotype tests, while maintaining Type I error rate.

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits. Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits.

To facilitate positional cloning of complex trait susceptibility loci, we are investigating methods to reduce the effort required to identify trait-associated alleles. We examined primer extension analysis by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry to screen single-nucleotide polymorphisms (SNPs) for association by using DNA pools. We tested whether this method can accurately estimate allele frequency differences between pools while maintaining the high-throughput nature of assay design, sample handling, and scoring. We follow up interesting allele frequency differences in pools by genotyping individuals. We tested DNA pools of 182, 228, and 499 individuals using 16 SNPs with minor allele frequencies 0.026–0.486 and allele frequency differences 0.001–0.108 that we had genotyped previously on individuals and 381 SNPs that we had not. Precision, as measured by the average standard deviation among 16 semidependent replicates, was 0.021 ± 0.011 for the 16 SNPs and 0.018 ± 0.008 for the 291/381 SNPs used in further analysis. For the 16 SNPs, the average absolute error in predicting allele frequency differences between pools was 0.009; the largest errors were 0.031, 0.028, and 0.027. We determined that compensating for unequal peak heights in heterozygotes improved precision of allele frequency estimates but had only a very minor effect on accuracy of allele frequency differences between pools. Based on these data and assuming pools of 500 individuals, we conclude that at significance level 0.05 we would have 95% (82%) power to detect population allele frequency differences of 0.07 for control allele frequencies of 0.10 (0.50).

While small clinical trials have shown that improved glycemic control reduces the risk of progression of microalbuminuria to proteinuria, two recent clinical trials did not confirm this finding. We sought to reconcile the contradictory evidence by examining the dose-response relationship between hyperglycemia and progression of microalbuminuria to proteinuria in individuals with type 1 diabetes and microalbuminuria (n = 312) who were followed for 4 years with repeated assessments of urinary albumin excretion. Since 33 patients did not participate in follow-up (10.6%), data for 279 patients were analyzed. Urinary albumin excretion level worsened to proteinuria in 40 (4.1 per 100 person-years). To examine the dose-response relationship, baseline HbA1c was divided into four roughly equal groups using the cut points 8, 9, and 10%. The incidence rate varied significantly among the four groups (P = 0.008). Among those with HbA1c <8.0%, the incidence rate of progression was only 1.3 per 100 person-years, while it was 5.1, 4.2, and 6.7 per 100 person-years in the three other groups. We used generalized additive models to examine the dose-response curve using HbA1c as a continuous variable and found that the risk of progression rises steeply between an HbA1c of 7.5-8.5% and then remains approximately constant across higher levels. In conclusion, the results of this study suggest that, in patients with microalbuminuria, the risk of progression to overt proteinuria can be reduced by improved glycemic control only if the HbA1c is maintained below 8.5%. Moreover, below that value, the risk declines as the level of HbA1c decreases.

The goal of this study was to examine the association between known polymorphisms in the apolipoprotein E gene (APOE) and diabetic nephropathy (DN) in type 1 diabetes. We used both a case-control comparison and a family-based study design known as the transmission/disequilibrium test (TDT). For the case-control comparison, we collected DNA from 223 subjects with clinically diagnosed DN and 196 control subjects with normoalbuminuria and long-duration type 1 diabetes (> or = 15 years). For the family-based study, we obtained DNA from both parents of 154 DN subjects and 81 control subjects. The frequency of the epsilon2 allele of exon 4 of APOE was significantly higher in DN subjects than in control subjects. The risk of DN was 3.1 times higher (95% CI 1.6-5.9) in carriers of this allele than in noncarriers. In the family study, heterozygous parents for the E2 allele preferentially transmitted epsilon2 to DN offspring (64 vs. 36%, P < 0.03). Four additional polymorphisms (i.e., -491 A/T, -219 G/T, IE1 G/C, and APOCI insertion/deletion [I/D]) that flank the APOE locus were not associated with DN in either the case-control comparison or in the family-based study. In conclusion, the results of the case-control as well as the family-based study provide evidence that the epsilon2 allele of APOE increases the risk of DN in type 1 diabetes. The molecular mechanisms underlying this risk are unclear at present.

Resistin is a peptide hormone produced by adipocytes that is present at high levels in sera of obese mice and may be involved in glucose homeostasis through regulation of insulin sensitivity. Several studies in humans have found associations between polymorphisms in the resistin gene and obesity, insulin sensitivity and blood pressure. An association between variation in the resistin gene and Type 2 diabetes has been reported in some, but not all studies. The aim of this study was to analyse variants of the resistin gene for association with Type 2 diabetes and related traits in a Finnish sample. In 781 cases with Type 2 diabetes, 187 spouse controls and 222 elderly controls of Finnish origin, we genotyped four previously identified non-coding single-nucleotide polymorphisms (SNPs): -420C>G from the promoter region, +156C>T and +298G>A from intron 2, and +1084G>A from the 3′ untranslated region. We then tested whether these SNPs were associated with Type 2 diabetes and related traits. The SNPs were not significantly associated with Type 2 diabetes. However, SNPs −420C>G, +156C>T and +298G>A and the common haplotype for these three markers were associated with increased values of weight-related traits and diastolic blood pressure in cases, lower weight in elderly control subjects, and lower insulin sensitivity and greater acute insulin response in spouses. Furthermore, the +1084G allele was associated with lower HDL cholesterol in both cases and controls, higher systolic blood pressure and waist circumference in cases, and greater acute insulin response in spouse controls. Our results add to growing evidence that resistin is associated with variation in weight, fat distribution and insulin resistance.

Prior reports have suggested that variants in the genes for maturity-onset diabetes of the young (MODY) may confer susceptibility to type 2 diabetes, but results have been conflicting and coverage of the MODY genes has been incomplete. To complement our previous studies of HNF4A, we examined the other five known MODY genes for association with type 2 diabetes in Finnish individuals. For each of the five genes, we selected 1) nonredundant single nucleotide polymorphisms (SNPs) (r(2)< 0.8 with other SNPs) from the HapMap database or another linkage disequilibrium map, 2) SNPs with previously reported type 2 diabetes association, and 3) nonsynonymous coding SNPs. We tested 128 SNPs for association with type 2 diabetes in 786 index cases from type 2 diabetic families and 619 normal glucose-tolerant control subjects. We followed up 35 of the most significant SNPs by genotyping them on another 384 case subjects and 366 control subjects from Finland. We also supplemented our previous HNF4A results by genotyping 12 SNPs on additional Finnish samples. After correcting for testing multiple correlated SNPs within a gene, we find evidence of type 2 diabetes association with SNPs in five of the six known MODY genes: GCK, HNF1A, HNF1B, NEUROD1, and HNF4A. Our data suggest that common variants in several MODY genes play a modest role in type 2 diabetes susceptibility.

For most associations of common single nucleotide polymorphisms (SNPs) with common diseases, the genetic model of inheritance is unknown.The authors extended and applied a Bayesian meta-analysis approach to data from 19 studies on 17 replicated associations with type 2 diabetes.For 13 SNPs, the data fitted very well to an additive model of inheritance for the diabetes risk allele; for 4 SNPs, the data were consistent with either an additive model or a dominant model; and for 2 SNPs, the data were consistent with an additive or recessive model.Results were robust to the use of different priors and after exclusion of data for which index SNPs had been examined indirectly through proxy markers.The Bayesian meta-analysis model yielded point estimates for the genetic effects that were very similar to those previously reported based on fixed-or random-effects models, but uncertainty about several of the effects was substantially larger.The authors also examined the extent of between-study heterogeneity in the genetic model and found generally small between-study deviation values for the genetic model parameter.Heterosis could not be excluded for 4 SNPs.Information on the genetic model of robustly replicated association signals derived from genome-wide association studies may be useful for predictive modeling and for designing biologic and functional experiments.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

For aggregation tests of genes or regions, the set of included variants often have small total minor allele counts (MACs), and this is particularly true when the most deleterious sets of variants are considered. When MAC is low, commonly used asymptotic tests are not well calibrated for binary phenotypes and can have conservative or anti-conservative results and potential power loss. Empirical p-values obtained via resampling methods are computationally costly for highly significant p-values and the results can be conservative due to the discrete nature of resampling tests. Based on the observation that only the individuals containing minor alleles contribute to the score statistics, we develop an efficient resampling method for single and multiple variant score-based tests that can adjust for covariates. Our method can improve computational efficiency >1000-fold over conventional resampling for low MAC variant sets. We ameliorate the conservativeness of results through the use of mid-p-values. Using the estimated minimum achievable p-value for each test, we calibrate QQ plots and provide an effective number of tests. In analysis of a case-control study with deep exome sequence, we demonstrate that our methods are both well calibrated and also reduce computation time significantly compared with resampling methods.

Integration of genome-wide association study (GWAS) signals with expression quantitative trait loci (eQTL) studies enables identification of candidate genes. However, evaluating whether nearby signals may share causal variants, termed colocalization, is affected by the presence of allelic heterogeneity, different variants at the same locus impacting the same phenotype. We previously identified eQTL in subcutaneous adipose tissue from 770 participants in the Metabolic Syndrome in Men (METSIM) study and detected 15 eQTL signals that colocalized with GWAS signals for waist-hip ratio adjusted for body mass index (WHRadjBMI) from the Genetic Investigation of Anthropometric Traits consortium. Here, we reevaluated evidence of colocalization using two approaches, conditional analysis and the Bayesian test COLOC, and show that providing COLOC with approximate conditional summary statistics at multi-signal GWAS loci can reconcile disagreements in colocalization classification between the two tests. Next, we performed conditional analysis on the METSIM subcutaneous adipose tissue data to identify conditionally distinct or secondary eQTL signals. We used the two approaches to test for colocalization with WHRadjBMI GWAS signals and evaluated the differences in colocalization classification between the two tests. Through these analyses, we identified four GWAS signals colocalized with secondary eQTL signals for FAM13A, SSR3, GRB14 and FMO1. Thus, at loci with multiple eQTL and/or GWAS signals, analyzing each signal independently enabled additional candidate genes to be identified.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

This study examined perceived substance use, opioid knowledge, and barriers to Black people accessing treatment for substance and opioid use disorder (SUD/OUD).Thirty-nine participants completed the community survey and The Brief Opioid Overdose Knowledge questionnaire. Qualitative interviews were conducted with 18 stakeholders and 9 people with SUD/OUD.Out of 39 participants, <50 % knew where to refer someone for treatment and fewer knew where to access naloxone. Majority of the stakeholders and people with SUD/OUD reported stigma as a treatment barrier.Studies related to provider anti-stigma trainings and psychoeducation for Black people living in the rural South are warranted.

The aim of the Finland-United States Investigation of NIDDM Genetics (FUSION) study is to identify genes that predispose to type 2 diabetes or are responsible for variability in diabetes-related traits via a positional cloning and positional candidate gene approach. In a previously published genome-wide scan of 478 Finnish affected sibling pair (ASP) families (FUSION 1), the strongest linkage results were on chromosomes 20 and 11. We now report a second genome-wide scan using an independent set of 242 Finnish ASP families (FUSION 2), a detailed analysis of the combined set of 737 FUSION 1 + 2 families (495 updated FUSION 1 families), and fine mapping of the regions of chromosomes 11 and 20. The strongest FUSION 2 linkage results were on chromosomes 6 (maximum logarithm of odds score [MLS] = 2.30 at 95 cM) and 14 (MLS = 1.80 at 57 cM). For the combined FUSION 1 + 2 families, three results were particularly notable: chromosome 11 (MLS = 2.98 at 82 cM), chromosome 14 (MLS = 2.74 at 58 cM), and chromosome 6 (MLS = 2.66 at 96 cM). We obtained smaller FUSION 1 + 2 MLSs on chromosomes X (MLS = 1.27 at 152 cM) and 20p (MLS = 1.21 at 20 cM). Among the 10 regions that showed nominally significant evidence for linkage in FUSION 1, four (on chromosomes 6, 11, 14, and X) also showed evidence for linkage in FUSION 2 and stronger evidence for linkage in the combined FUSION 1 + 2 sample.

Mitochondria play an integral role in ATP production in cells and are involved in glucose metabolism and insulin secretion, suggesting that variants in the mitochondrial genome may contribute to diabetes susceptibility. In a study of Finnish families ascertained for type 2 diabetes mellitus (T2DM), we genotyped single nucleotide polymorphisms (SNPs) based on phylogenetic networks. These SNPs defined eight major haplogroups and subdivided groups H and U, which are common in Finns. We evaluated association with both diabetes disease status and up to 14 diabetes-related traits for 762 cases, 402 non-diabetic controls, and 465 offspring of genotyped females. Haplogroup J showed a trend toward association with T2DM affected status (OR 1.69, P=0.056) that became slightly more significant after excluding cases with affected fathers (OR 1.77, P=0.045). We also genotyped non-haplogroup-tagging SNPs previously reported to show evidence for association with diabetes or related traits. Our data support previous evidence for association of T16189C with reduced ponderal index at birth and also show evidence for association with reduced birthweight but not with diabetes status. Given the multiple tests performed and the significance levels obtained, this study suggests that mitochondrial genome variants may play at most a modest role in glucose metabolism in the Finnish population. Furthermore, our data do not support a reported maternal inheritance pattern of T2DM but instead show a strong effect of recall bias.

Abstract Motivation: Finding the potential functional significance of SNPs is a major bottleneck in understanding genome-wide SNP scanning results, as the related functional data are distributed across many different databases. The SNP Function Portal is designed to be a clearing house for all public domain SNP functional annotation data, as well as in-house functional annotations derived from different data sources. It currently contains SNP functional annotations in six major categories including genomic elements, transcription regulation, protein function, pathway, disease and population genetics. Besides extensive SNP functional annotations, the SNP Function Portal includes a powerful search engine that accepts different types of genetic markers as input and identifies all genetically related SNPs based on the HapMap Phase II data as well as the relationship of different markers to known genes. As a result, our system allows users to identify the potential biological impact of genetic markers and complex relationships among genetic markers and genes, and it greatly facilitates knowledge discovery in genome-wide SNP scanning experiments. Availability: Contact: mengf@umich.edu

Abstract The pattern and nature of linkage disequilibrium in the human genome is being studied and catalogued as part of the International HapMap Project [:2003 Nature 426:789–796]. A key goal of the HapMap Project is to enable identification of tag single nucleotide polymorphisms (SNPs) that capture a substantial portion of common human genetic variability while requiring only a small fraction of SNPs to be genotyped [International HapMap Consortium, 2005: Nature 437:1299–1320]. In the current study, we examined the effectiveness of using the CEU HapMap database to select tag SNPs for a Finnish sample. We selected SNPs in a 17.9‐Mb region of chromosome 14 based on pairwise linkage disequilibrium (r 2 ) estimates from the HapMap CEU sample, and genotyped 956 of these SNPs in 1,425 Finnish individuals. An excess of SNPs showed significantly different allele frequencies between the HapMap CEU and the Finnish samples, consistent with population‐specific differences. However, we observed strong correlations between the two samples for estimates of allele frequencies, r 2 values, and haplotype frequencies. Our results demonstrate that the HapMap CEU samples provide an adequate basis for tag SNP selection in Finnish individuals, without the need to create a map specifically for the Finnish population, and suggest that the four‐population HapMap data will provide useful information for tag SNP selection beyond the specific populations from which they were sampled. Genet. Epidemiol . 2006. © 2005 Wiley‐Liss, Inc.

Rare fully penetrant mutations in AKT2 are an established cause of monogenic disorders of glucose metabolism. Recently, a novel partial loss-of-function AKT2 coding variant (p.Pro50Thr) was identified that is nearly specific to Finns (frequency 1.1%), with the low-frequency allele associated with an increase in fasting plasma insulin level and risk of type 2 diabetes. The effects of the p.Pro50Thr AKT2 variant (p.P50T/AKT2) on insulin-stimulated glucose uptake (GU) in the whole body and in different tissues have not previously been investigated. We identified carriers (N = 20) and matched noncarriers (N = 25) for this allele in the population-based Metabolic Syndrome in Men (METSIM)study and invited these individuals back for positron emission tomography study with [18F]-fluorodeoxyglucose during euglycemic hyperinsulinemia. When we compared p.P50T/AKT2 carriers to noncarriers, we found a 39.4% reduction in whole-body GU (P = 0.006) and a 55.6% increase in the rate of endogenous glucose production (P = 0.038). We found significant reductions in GU in multiple tissues-skeletal muscle (36.4%), liver (16.1%), brown adipose (29.7%), and bone marrow (32.9%)-and increases of 16.8-19.1% in seven tested brain regions. These data demonstrate that the p.P50T substitution of AKT2 influences insulin-mediated GU in multiple insulin-sensitive tissues and may explain, at least in part, the increased risk of type 2 diabetes in p.P50T/AKT2 carriers.

Advancing ability to predict who is likely to develop depression holds great potential in reducing the disease burden. Here, we use the predictable and large increase in depression with physician training stress to identify predictors of depression. Applying the major depressive disorder polygenic risk score (MDD-PRS) derived from the most recent Psychiatric Genomics Consortium–UK Biobank–23andMe genome-wide association study to 5,227 training physicians, we found that MDD-PRS predicted depression under training stress (β = 0.095, P = 4.7 × 10−16) and that MDD-PRS was more strongly associated with depression under stress than at baseline (MDD-PRS × stress interaction β = 0.036, P = 0.005). Further, known risk factors accounted for substantially less of the association between MDD-PRS and depression when under stress than at baseline, suggesting that MDD-PRS adds unique predictive power in depression prediction. Finally, we found that low MDD-PRS may have particular use in identifying individuals with high resilience. Together, these findings suggest that MDD-PRS holds promise in furthering our ability to predict vulnerability and resilience under stress. Using physician stress as a model stressor, Fang et al. demonstrate that the polygenic risk score for major depressive disorder is a stronger predictor of depression under stress than under baseline conditions and may be particularly useful for identifying resilience.

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Summary Aims To investigate association and linkage between DNA sequence variants in the aldose reductase (AR) gene on chromosome 7q35 and diabetic nephropathy (DN) in Type 1 diabetes mellitus. Methods By sequencing the promoter region and 10 exons in eight DN cases and eight controls, a frequent bi‐allelic polymorphism (C‐106T) was discovered. This polymorphism and the known 5′ALR2 dinucleotide repeat polymorphism were genotyped in unrelated cases with advanced nephropathy ( n = 221) and unrelated controls with normoalbuminuria ( n = 193). For a family based study, 166 case‐trios (case and both parents) and 83 control‐trios (control and both parents) were also genotyped. Results In the case–control study, carriers of the Z‐2 allele of the 5′ALR2 polymorphism had a significantly higher risk of DN than non‐carriers (odds ratios: 1.6 for heterozygotes and 2.1 for homozygotes, P &lt; 0.05 for each). The same was true for carriers of the T allele of the C‐106T polymorphism (odds ratios: 1.6 for heterozygotes and 1.9 for homozygotes, P &lt; 0.05 for each). Moreover, the haplotype carrying both risk alleles was in excess in DN cases. In the family study, transmission of risk alleles from heterozygous parents was consistent with the case–control study, excess transmission in case‐trios and deficient in control‐trios. Conclusions Association between DN and two DNA sequence variants in the promoter region of the AR gene implicates the polyol pathway in the development of kidney complications in Type 1 diabetes mellitus. Further examination of the molecular mechanisms underlying these findings may provide insight into the pathogenesis of DN.

We have identified B10, a plasma membrane protein previously defined by a monoclonal antibody, as an alkaline phosphodiesterase I (APDE) expressed in the plasma membrane of rat hepatocytes and enterocytes, with a restricted apical distribution. B10 complementary DNA (cDNA) was cloned from a rat intestinal library screened with a polyclonal antibody directed to the hepatic protein. Two distinct B10 clones with an open reading frame of 2,625 bp were obtained that differed only by 12 bases in the coding region. One B10 clone had a single base difference with gp130RB13-6 cDNA, which was recently cloned in rat fetal brain. B10/gp130RB13-6 had 50% identity at the amino acid level with the plasma cell antigen PC-1, an APDE cloned in the mouse and in human. Anti-B10 antibodies immunoprecipitated 34% of the APDE activity in liver plasma membranes and over 95% of the APDE activity in intestinal cells. Most of the remaining activity in hepatocytes (44%) could be immunoprecipitated by antibodies directed to PC-1. APDE activity immunoprecipitated with anti-B10 antibodies was found in the apical rat liver plasma membrane fractions on a sucrose gradient whereas most of the remaining APDE activity was associated with the basolateral fractions, which contained PC-1. By immunofluorescence, B10 was localized to the apical surfaces of hepatocytes and enterocytes whereas PC-1 was present on the basolateral surfaces of hepatocytes. B10/gp130RB13-6 and rat PC-1 are a unique example of distinct molecules having similar enzymatic activity but different apical/basolateral location, and possibly different functions.

To fine map genes, investigators often test for disease-marker association in chromosomal regions with evidence for linkage. Given a marker allele tentatively associated with disease, one would ask if this allele, or one in linkage disequilibrium (LD) with it, could account in part for the observed linkage signal. This question can be addressed by determining if families selected on the basis of the presence of the tentatively associated allele show stronger evidence of linkage as measured by increased allele sharing identical by descent (IBD) by affected family members. However, common selection strategies can be biased for or against linkage in the marker region, even given no disease-marker association. We define unbiased selection schemes and extend the definition to allow weighted selection on the basis of all genotyped family members. For affected-sibship data, we describe three genotype-based weight variables, corresponding to dominant, recessive, and additive models. We then introduce a test for association of a family weight variable with excess IBD sharing. This test allows us to determine if the linkage signal in a region can be attributed in part to the presence of a marker allele, either because of direct involvement in disease etiology or because of LD with a predisposing genetic variant. For samples of 500 affected sib pairs, the tests are powerful in detection of genotype-IBD sharing association, even for disease models with sib relative risk as low as lambda S=1.1, or when evidence for linkage is absent because of sampling variation. This makes our method a new tool for detecting linkage as well as association, especially in regions harboring a candidate gene. We have implemented these methods in the software package GIST (Genotype-IBD Sharing Test).

Recent studies from this laboratory have demonstrated that a form of hexokinase characteristic of rapidly growing, highly glycolytic tumor cells is bound to an outer mitochondrial membrane receptor complex containing a Mr 35,000 pore protein (D. M. Parry and P. L. Pedersen, J. Biol. Chem., 258: 10904-10912, 1983; R. A. Nakashima, et al., Biochemistry, 25: 1015-1021, 1986). In new studies reported here the specificity of this receptor complex for binding hexokinase is defined, and a purification scheme is described which leads to a homogeneous and bindable form of the tumor hexokinase. In the AS-30D hepatoma, hexokinase activity is elevated more than 100-fold relative to liver tissue. The relative increase in hexokinase activity is 8 times greater than that of any other glycolytic enzyme. Hexokinase is the only glycolytic enzyme of AS-30D cells to exhibit a mitochondrial/cytoplasmic specific activity ratio greater than 1, showing a 3.5-fold elevation in the mitochondrial fraction. Purification of hexokinase is accomplished by preferential solubilization of the mitochondrial bound enzyme with glucose-6-phosphate, followed by high-performance liquid chromatography on gel permeation and anion exchange columns. The final fraction has a specific activity of 144 units per mg of protein, with a Km for glucose of 0.13 mM and for ATP of 1.4 mM. The purified tumor enzyme migrates as a single species upon sodium dodecyl sulfate: polyacrylamide gel electrophoresis with an apparent molecular weight of 98,000. Significantly, the purified tumor enzyme retains its activity for mitochondrial binding. Additional results derived from chromatographic, polyclonal antibody, and amino acid analysis studies indicate that the predominant rat hepatoma hexokinase species is related most closely to isozymic form(s) of the enzyme commonly referred to as type II, and least related to the liver type IV isozyme (glucokinase).

Bipolar-related subphenotypes that cluster within families may help identify subsets of patients that are more genetically homogeneous. Environmental or assessment factors that segregate by family may influence estimates of familiality. We aimed to determine familiality of subphenotypes of bipolar disorder (BP), accounting for effects of age, sex, diagnosis, and site/wave of ascertainment. We studied 589 sibships with 1416 siblings affected with bipolar I (BPI), schizoaffective disorder, bipolar type (SAB), bipolar II (BPII), or recurrent unipolar depression (RUDD). Sibships were from families with > or =2 BPI cases collected by the NIMH Bipolar Genetics Initiative (NIMHBGI). Rapid cycling showed the strongest evidence for familiality [odds ratio (OR) (95%CI) = 2.02 (1.43, 2.85), P = 6.0 x 10(-5)] in a model including age, sex, diagnosis, and site/wave of ascertainment. Additional significantly familial traits were comorbid alcohol abuse/dependence (P = 2 x 10(-4)) and comorbid panic disorder (P = 8 x 10(-3)), as well as psychosis, suicidal thoughts, and rapid mood switching (P = 6 x 10(-3) - 0.03). Omission of the effect of site/wave of ascertainment from the model inflated the significance level of the apparent familial association of almost all subphenotypes from one to four orders of magnitude. We have found evidence of familiality for subphenotypes of BP. In multicenter samples, familiality may be overestimated if variability in diagnosis of subphenotypes between site/wave of ascertainment is not considered.

There are strong associations between measures of inflammation and type 2 diabetes, but the causal directions of these associations are not known. We tested the hypothesis that common gene variants known to alter circulating levels of inflammatory proteins, or known to alter autoimmune-related disease risk, influence type 2 diabetes risk. We selected 46 variants: (1) eight variants known to alter circulating levels of inflammatory proteins, including those in the IL18, IL1RN, IL6R, MIF, PAI1 (also known as SERPINE1) and CRP genes; and (2) 38 variants known to predispose to autoimmune diseases, including type 1 diabetes. We tested the associations of these variants with type 2 diabetes using a meta-analysis of 4,107 cases and 5,187 controls from the Wellcome Trust Case Control Consortium, the Diabetes Genetics Initiative, and the Finland–United States Investigation of NIDDM studies. We followed up associated variants (p < 0.01) in a further set of 3,125 cases and 3,596 controls from the UK. We found no evidence that inflammatory or autoimmune disease variants are associated with type 2 diabetes (at p ≤ 0.01). The OR observed between the variant altering IL-18 levels, rs2250417, and type 2 diabetes (OR 1.00 [95% CI 0.99–1.03]), is much lower than that expected given (1) the effect of the variant on IL-18 levels (0.28 SDs per allele); and (2) estimates, based on other studies, of the correlation between IL-18 levels and type 2 diabetes risk (approximate OR 1.15 [95% CI 1.09–1.21] per 0.28 SD increase in IL-18 levels). Our study provided no evidence that variants known to alter measures of inflammation, autoimmune or inflammatory disease risk, including type 1 diabetes, alter type 2 diabetes risk.

Genome-wide association studies (GWAS) have identified numerous loci associated with various complex traits for which the underlying susceptibility gene(s) remain unknown. In a GWAS for high-density lipoprotein-cholesterol (HDL-C) level, one strongly associated locus contains at least two biologically compelling candidates, methylmalonic aciduria cblB type (MMAB) and mevalonate kinase (MVK). To detect evidence of cis-acting regulation at this locus, we measured relative allelic expression of transcribed SNPs in five genes using human hepatocyte samples heterozygous for the transcribed SNP. If an HDL-C-associated SNP allele differentially regulates mRNA level in cis, samples heterozygous both for a transcribed SNP and an HDL-C-associated SNP should display allelic expression imbalance (AEI) of the transcribed SNP. We designed statistical tests to detect AEI in a comprehensive set of linkage disequilibrium (LD) scenarios between the transcribed SNP and an HDL-C-associated SNP (rs7298565) in phase unknown samples. We observed significant AEI of 22% in MMAB (P = 1.4 × 10−13, transcribed SNP rs11067231), and the allele associated with lower HDL-C level was associated with greater MMAB transcript level. The same rs7298565 allele was also associated with higher MMAB mRNA level (P = 0.0081) and higher MMAB protein level (P = 0.0020). In contrast, MVK, UBE3B, KCTD10 and ACACB did not show significant AEI (P ≥ 0.05). These data suggest MMAB is the most likely gene influencing HDL-C levels at this locus and demonstrate that measuring AEI at loci containing more than one candidate gene can prioritize genes for functional studies.

We describe a reference panel of 64,976 human haplotypes at 39,235,157 SNPs constructed using whole genome sequence data from 20 studies of predominantly European ancestry. Using this resource leads to accurate genotype imputation at minor allele frequencies as low as 0.1%, a large increase in the number of SNPs tested in association studies and can help to discover and refine causal loci. We describe remote server resources that allow researchers to carry out imputation and phasing consistently and efficiently.

High cancerous inhibitor of PP2A (CIP2A) protein levels at diagnosis of chronic myeloid leukaemia (CML) are predictive of disease progression in imatinib-treated patients. It is not known whether this is true in patients treated with second generation tyrosine kinase inhibitors (2G TKI) from diagnosis, and whether 2G TKIs modulate the CIP2A pathway. Here, we show that patients with high diagnostic CIP2A levels who receive a 2G TKI do not progress, unlike those treated with imatinib (P=<0.0001). 2G TKIs induce more potent suppression of CIP2A and c-Myc than imatinib. The transcription factor E2F1 is elevated in high CIP2A patients and following 1 month of in vivo treatment 2G TKIs suppress E2F1 and reduce CIP2A; these effects are not seen with imatinib. Silencing of CIP2A, c-Myc or E2F1 in K562 cells or CML CD34+ cells reactivates PP2A leading to BCR-ABL suppression. CIP2A increases proliferation and this is only reduced by 2G TKIs. Patients with high CIP2A levels should be offered 2G TKI treatment in preference to imatinib. 2G TKIs disrupt the CIP2A/c-Myc/E2F1 positive feedback loop, leading to lower disease progression risk. The data supports the view that CIP2A inhibits PP2Ac, stabilising E2F1, creating a CIP2A/c-Myc/E2F1 positive feedback loop, which imatinib cannot overcome.

Significance Contributions of rare variants to common and complex traits such as type 2 diabetes (T2D) are difficult to measure. This paper describes our results from deep whole-genome analysis of large Mexican-American pedigrees to understand the role of rare-sequence variations in T2D and related traits through enriched allele counts in pedigrees. Our study design was well-powered to detect association of rare variants if rare variants with large effects collectively accounted for large portions of risk variability, but our results did not identify such variants in this sample. We further quantified the contributions of common and rare variants in gene expression profiles and concluded that rare expression quantitative trait loci explain a substantive, but minor, portion of expression heritability.

Abstract Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3–2.8, one-sided p = 6.0 × 10 −4 ), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Finnish samples have been extensively utilized in studying single-gene disorders, where the founder effect has clearly aided in discovery, and more recently in genome-wide association studies of complex traits, where the founder effect has had less obvious impacts. As the field starts to explore rare variants' contribution to polygenic traits, it is of great importance to characterize and confirm the Finnish founder effect in sequencing data and to assess its implications for rare-variant association studies. Here, we employ forward simulation, guided by empirical deep resequencing data, to model the genetic architecture of quantitative polygenic traits in both the general European and the Finnish populations simultaneously. We demonstrate that power of rare-variant association tests is higher in the Finnish population, especially when variants' phenotypic effects are tightly coupled with fitness effects and therefore reflect a greater contribution of rarer variants. SKAT-O, variable-threshold tests, and single-variant tests are more powerful than other rare-variant methods in the Finnish population across a range of genetic models. We also compare the relative power and efficiency of exome array genotyping to those of high-coverage exome sequencing. At a fixed cost, less expensive genotyping strategies have far greater power than sequencing; in a fixed number of samples, however, genotyping arrays miss a substantial portion of genetic signals detected in sequencing, even in the Finnish founder population. As genetic studies probe sequence variation at greater depth in more diverse populations, our simulation approach provides a framework for evaluating various study designs for gene discovery.

Cancerous inhibitor of protein phosphatase 2A (CIP2A) is a predictive biomarker of disease progression in many malignancies, including imatinib-treated chronic myeloid leukemia (CML). Although high CIP2A levels correlate with disease progression in CML, the underlying molecular mechanisms remain elusive. In a screen of diagnostic chronic phase samples from patients with high and low CIP2A protein levels, high CIP2A levels correlate with an antiapoptotic phenotype, characterized by downregulation of proapoptotic BCL-2 family members, including BIM, PUMA and HRK, and upregulation of the antiapoptotic protein BCL-XL. These results suggest that the poor prognosis of patients with high CIP2A levels is due to an antiapoptotic phenotype. Disrupting this antiapoptotic phenotype by inhibition of BCL-XL via RNA interference or A-1331852, a novel, potent and BCL-XL-selective inhibitor, resulted in extensive apoptosis either alone or in combination with imatinib, dasatinib or nilotinib, both in cell lines and in primary CD34+ cells from patients with high levels of CIP2A. These results demonstrate that BCL-XL is the major antiapoptotic survival protein and may be a novel therapeutic target in CML.

We have determined the half-lives and amounts per hepatocyte of the polymeric IgA receptor (pIgA-R) and four rat hepatocyte plasma membrane proteins and subsequently have predicted their rates of synthesis and possible routes of degradation. Using in vivo pulse-chase metabolic labeling with L-[35S]cysteine, we found that the pIgA-R had an apparent half-life of 1.1 h. Additional metabolic labeling experiments showed that CE9, HA4, and HA321 had apparent half-lives of 4-5 days, and dipeptidyl peptidase IV had an apparent half-life of 9 days. To quantify the amount of each protein per hepatocyte, homogenates and a standard curve of purified protein were compared by immunoblotting. We found that these proteins were present at 1-8 x 10(6) molecules/hepatocyte. The calculated rate of synthesis for pIgA-R was 1.6 x 10(6) molecules/hepatocyte/h, whereas the others were synthesized at much lower rates (0.9-5 x 10(4) molecules/hepatocyte/h). Using immunoblot analysis, we found that pIgA-R was released into bile at a rate of 30%/h (700%/day), whereas dipeptidyl peptidase IV and HA4 were released at a rate of 2-3%/day. While the majority of the loss of pIgA-R from hepatocytes occurred by release into the bile, less than 30% of the degradation of dipeptidyl peptidase IV and HA4 could be accounted for by this pathway, suggesting that the remaining molecules must be retrieved from the apical surface before degradation.

We have evaluated the relationship between systolic blood pressure (SBP) and age, body mass index (BMI), waist circumference, sodium to potassium ratio (Na/K), and tobacco use in an urban African population. We conducted a random, population-based, cross-sectional survey of people 25 years and older in Marondera, Zimbabwe, with over-sampling in older age groups (n = 775), using a method comparable to that used in International Collaborative Study on Hypertension in Blacks (ICSHIB). The age-adjusted prevalences of hypertension in Marondera (SBP >/=140/DBP >/=90/antihypertensive medication) were 30% for women and 21% for men. The average BMI was 26.3 kg/m2 for women and 21.4 kg/m2 for men. The prevalence of hypertension had a steep association with age and in women ranged from 15% (25-34 years) to 63% (55 years and over) and in men from 9% to 47%. No tobacco use in women and greater Na/K ratio in spot urines in men were significantly associated with an increased SBP. In both men and women the levels of hypertension and SBP were strongly positively associated with BMI, although the relationship appeared to plateau in women with a BMI greater than >/=25 kg/m2. At a given BMI, men and women had similar SBPs and prevalences of hypertension. There is a very high prevalence of hypertension among urban Zimbabweans, particularly among women. Under the assumption the studies are comparable, the prevalence of hypertension in Zimbabwean women (41%) and men (26%) after age adjustment to the ICSHIB populations, appeared higher than almost all of the ICSHIB populations, including those with higher average body mass indexes. Journal of Human Hypertension (2000) 14, 65-73.

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N = 38,580) informative for adult waist circumference (WC) and waist-hip ratio (WHR).We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height.Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P = 1.9610 211 ) and MSRA (WC, P = 8.9610 29 ).A third locus, near LYPLAL1, was associated with WHR in women only (P = 2.6610 28 ).The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution.By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.

OBJECTIVE Linkage of the chromosome 1q21–25 region to type 2 diabetes has been demonstrated in multiple ethnic groups. We performed common variant fine-mapping across a 23-Mb interval in a multiethnic sample to search for variants responsible for this linkage signal. RESEARCH DESIGN AND METHODS In all, 5,290 single nucleotide polymorphisms (SNPs) were successfully genotyped in 3,179 type 2 diabetes case and control subjects from eight populations with evidence of 1q linkage. Samples were ascertained using strategies designed to enhance power to detect variants causal for 1q linkage. After imputation, we estimate ∼80% coverage of common variation across the region (r 2 &amp;gt; 0.8, Europeans). Association signals of interest were evaluated through in silico replication and de novo genotyping in ∼8,500 case subjects and 12,400 control subjects. RESULTS Association mapping of the 23-Mb region identified two strong signals, both of which were restricted to the subset of European-descent samples. The first mapped to the NOS1AP (CAPON) gene region (lead SNP: rs7538490, odds ratio 1.38 [95% CI 1.21–1.57], P = 1.4 × 10−6, in 999 case subjects and 1,190 control subjects); the second mapped within an extensive region of linkage disequilibrium that includes the ASH1L and PKLR genes (lead SNP: rs11264371, odds ratio 1.48 [1.18–1.76], P = 1.0 × 10−5, under a dominant model). However, there was no evidence for association at either signal on replication, and, across all data (&amp;gt;24,000 subjects), there was no indication that these variants were causally related to type 2 diabetes status. CONCLUSIONS Detailed fine-mapping of the 23-Mb region of replicated linkage has failed to identify common variant signals contributing to the observed signal. Future studies should focus on identification of causal alleles of lower frequency and higher penetrance.