Clarence K. Zhang

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohns disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohns disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10−16, odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

A common variant in the patatin-like phospholipase 3 gene (PNPLA3) is associated with fatty liver disease in obese children and adolescents†‡

The genetic factors associated with susceptibility to nonalcoholic fatty liver disease (NAFLD) in pediatric obesity remain largely unknown. Recently, a nonsynonymous single‐nucleotide polymorphism (rs738409), in the patatin‐like phospholipase 3 gene (PNPLA3) has been associated with hepatic steatosis in adults. In a multiethnic group of 85 obese youths, we genotyped the PNLPA3 single‐nucleotide polymorphism, measured hepatic fat content by magnetic resonance imaging and insulin sensitivity by the insulin clamp. Because PNPLA3 might affect adipogenesis/lipogenesis, we explored the putative association with the distribution of adipose cell size and the expression of some adipogenic/lipogenic genes in a subset of subjects who underwent a subcutaneous fat biopsy. Steatosis was present in 41% of Caucasians, 23% of African Americans, and 66% of Hispanics. The frequency of PNPLA3(rs738409) G allele was 0.324 in Caucasians, 0.183 in African Americans, and 0.483 in Hispanics. The prevalence of the G allele was higher in subjects showing hepatic steatosis. Surprisingly, subjects carrying the G allele showed comparable hepatic glucose production rates, peripheral glucose disposal rate, and glycerol turnover as the CC homozygotes. Carriers of the G allele showed smaller adipocytes than those with CC genotype (P = 0.005). Although the expression of PNPLA3, PNPLA2, PPARγ2(peroxisome proliferator‐activated receptor gamma 2), SREBP1c(sterol regulatory element binding protein 1c), and ACACA(acetyl coenzyme A carboxylase) was not different between genotypes, carriers of the G allele showed lower leptin (LEP)(P = 0.03) and sirtuin 1 (SIRT1) expression (P = 0.04). Conclusion: A common variant of the PNPLA3 gene confers susceptibility to hepatic steatosis in obese youths without increasing the level of hepatic and peripheral insulin resistance. The rs738409 PNPLA3 G allele is associated with morphological changes in adipocyte cell size. (HEPATOLOGY 2010.)

A Genome-Wide Association Study on Obesity and Obesity-Related Traits

Large-scale genome-wide association studies (GWAS) have identified many loci associated with body mass index (BMI), but few studies focused on obesity as a binary trait. Here we report the results of a GWAS and candidate SNP genotyping study of obesity, including extremely obese cases and never overweight controls as well as families segregating extreme obesity and thinness. We first performed a GWAS on 520 cases (BMI>35 kg/m2) and 540 control subjects (BMI<25 kg/m2), on measures of obesity and obesity-related traits. We subsequently followed up obesity-associated signals by genotyping the top ∼500 SNPs from GWAS in the combined sample of cases, controls and family members totaling 2,256 individuals. For the binary trait of obesity, we found 16 genome-wide significant signals within the FTO gene (strongest signal at rs17817449, P = 2.5×10−12). We next examined obesity-related quantitative traits (such as total body weight, waist circumference and waist to hip ratio), and detected genome-wide significant signals between waist to hip ratio and NRXN3 (rs11624704, P = 2.67×10−9), previously associated with body weight and fat distribution. Our study demonstrated how a relatively small sample ascertained through extreme phenotypes can detect genuine associations in a GWAS.

Variant in the glucokinase regulatory protein (GCKR) gene is associated with fatty liver in obese children and adolescents

Recently, the single nucleotide polymorphism (SNP) identified as rs1260326, in the glucokinase regulatory protein (GCKR), was associated with hypertriglyceridemia in adults. Because accumulation of triglycerides in hepatocytes represents the hallmark of steatosis, we aimed to investigate whether this variant might be associated with fatty liver (hepatic fat content, HFF%). Moreover, because recently rs738409 in the PNPLA3 and rs2854116 in the APOC3 were associated with fatty liver, we explored how the GCKR SNP and these two variants jointly influence hepatosteatosis. We studied 455 obese children and adolescents (181 Caucasians, 139 African Americans, and 135 Hispanics). All underwent an oral glucose tolerance test and fasting lipoprotein subclasses measurement by proton nuclear magnetic resonance. A subset of 142 children underwent a fast gradient magnetic resonance imaging to measure the HFF%. The rs1260326 was associated with elevated triglycerides (Caucasians P = 0.00014; African Americans P = 0.00417), large very low‐density lipoprotein (VLDL) (Caucasians P = 0.001; African Americans, P = 0.03), and with fatty liver (Caucasians P = 0.034; African Americans P = 0.00002; and Hispanics P = 0.016). The PNPLA3, but not the APOC3 rs2854116 SNP, was associated with fatty liver but not with triglyceride levels. There was a joint effect between the PNPLA3 and GCKR SNPs, explaining 32% of HFF% variance in Caucasians (P = 0.00161), 39.0% in African Americans (P = 0.00000496), and 15% in Hispanics (P = 0.00342). Conclusion: The rs1260326 in GCKR is associated with hepatic fat accumulation along with large VLDL and triglyceride levels. GCKR and PNPLA3 act together to convey susceptibility to fatty liver in obese youths. (Hepatology 2012)

Genome-Wide Association Study of Alcohol Dependence Implicates KIAA0040 on Chromosome 1q

Previous studies using SAGE (the Study of Addiction: Genetics and Environment) and COGA (the Collaborative Study on the Genetics of Alcoholism) genome-wide association study (GWAS) data sets reported several risk loci for alcohol dependence (AD), which have not yet been well replicated independently or confirmed by functional studies. We combined these two data sets, now publicly available, to increase the study power, in order to identify replicable, functional, and significant risk regions for AD. A total of 4116 subjects (1409 European-American (EA) cases with AD, 1518 EA controls, 681 African-American (AA) cases, and 508 AA controls) underwent association analysis. An additional 443 subjects underwent expression quantitative trait locus (eQTL) analysis. Genome-wide association analysis was performed in EAs to identify significant risk genes. All available markers in the genome-wide significant risk genes were tested in AAs for associations with AD, and in six HapMap populations and two European samples for associations with gene expression levels. We identified a unique genome-wide significant gene—KIAA0040—that was enriched with many replicable risk SNPs for AD, all of which had significant cis-acting regulatory effects. The distributions of −log(p) values for SNP-disease and SNP-expression associations for all markers in the TNN–KIAA0040 region were consistent across EAs, AAs, and five HapMap populations (0.369⩽r⩽0.824; 2.8 × 10−9⩽p⩽0.032). The most significant SNPs in these populations were in high LD, concentrating in KIAA0040. Finally, expression of KIAA0040 was significantly (1.2 × 10−11⩽p⩽1.5 × 10−6) associated with the expression of numerous genes in the neurotransmitter systems or metabolic pathways previously associated with AD. We concluded that KIAA0040 might harbor a causal variant for AD and thus might directly contribute to risk for this disorder. KIAA0040 might also contribute to the risk of AD via neurotransmitter systems or metabolic pathways that have previously been implicated in the pathophysiology of AD. Alternatively, KIAA0040 might regulate the risk via some interactions with flanking genes TNN and TNR. TNN is involved in neurite outgrowth and cell migration in hippocampal explants, and TNR is an extracellular matrix protein expressed primarily in the central nervous system.

A Comparison of Methods for Clustering 16S rRNA Sequences into OTUs

Recent studies of 16S rRNA sequences through next-generation sequencing have revolutionized our understanding of the microbial community composition and structure. One common approach in using these data to explore the genetic diversity in a microbial community is to cluster the 16S rRNA sequences into Operational Taxonomic Units (OTUs) based on sequence similarities. The inferred OTUs can then be used to estimate species, diversity, composition, and richness. Although a number of methods have been developed and commonly used to cluster the sequences into OTUs, relatively little guidance is available on their relative performance and the choice of key parameters for each method. In this study, we conducted a comprehensive evaluation of ten existing OTU inference methods. We found that the appropriate dissimilarity value for defining distinct OTUs is not only related with a specific method but also related with the sample complexity. For data sets with low complexity, all the algorithms need a higher dissimilarity threshold to define OTUs. Some methods, such as, CROP and SLP, are more robust to the specific choice of the threshold than other methods, especially for shorter reads. For high-complexity data sets, hierarchical cluster methods need a more strict dissimilarity threshold to define OTUs because the commonly used dissimilarity threshold of 3% often leads to an under-estimation of the number of OTUs. In general, hierarchical clustering methods perform better at lower dissimilarity thresholds. Our results show that sequence abundance plays an important role in OTU inference. We conclude that care is needed to choose both a threshold for dissimilarity and abundance for OTU inference.

A Novel, Functional and Replicable Risk Gene Region for Alcohol Dependence Identified by Genome-Wide Association Study

Several genome-wide association studies (GWASs) reported tens of risk genes for alcohol dependence, but most of them have not been replicated or confirmed by functional studies. The present study used a GWAS to search for novel, functional and replicable risk gene regions for alcohol dependence. Associations of all top-ranked SNPs identified in a discovery sample of 681 African-American (AA) cases with alcohol dependence and 508 AA controls were retested in a primary replication sample of 1,409 European-American (EA) cases and 1,518 EA controls. The replicable associations were then subjected to secondary replication in a sample of 6,438 Australian family subjects. A functional expression quantitative trait locus (eQTL) analysis of these replicable risk SNPs was followed-up in order to explore their cis-acting regulatory effects on gene expression. We found that within a 90 Mb region around PHF3-PTP4A1 locus in AAs, a linkage disequilibrium (LD) block in PHF3-PTP4A1 formed the only peak associated with alcohol dependence at p<10−4. Within this block, 30 SNPs associated with alcohol dependence in AAs (1.6×10−5≤p≤0.050) were replicated in EAs (1.3×10−3≤p≤0.038), and 18 of them were also replicated in Australians (1.8×10−3≤p≤0.048). Most of these risk SNPs had strong cis-acting regulatory effects on PHF3-PTP4A1 mRNA expression across three HapMap samples. The distributions of −log(p) values for association and functional signals throughout this LD block were highly consistent across AAs, EAs, Australians and three HapMap samples. We conclude that the PHF3-PTP4A1 region appears to harbor a causal locus for alcohol dependence, and proteins encoded by PHF3 and/or PTP4A1 might play a functional role in the disorder.

The influence of depression on quality of life in patients with inflammatory bowel disease.

Background: Inflammatory bowel disease is a chronic inflammatory disorder of the gastrointestinal tract that significantly impacts the health-related quality of life (HR-QOL). A decreased HR-QOL has been demonstrated in patients with active disease compared with patients in remission. In this cross-sectional study, we examined the role of depression and disease activity as independent factors in predicting patient’s HR-QOL. Methods: Hundred and five patients with either Crohn’s disease (CD) or ulcerative colitis (UC) were enrolled. Disease activity was evaluated using Crohn’s Disease Activity Index or Seo’s Activity Index. Depressive symptoms were evaluated using Beck’s Depression Inventory-II and Beck’s Depression Inventory for Primary Care (BDI-PC). HR-QOL was evaluated using the Short Inflammatory Bowel Disease Questionnaire. Simple and multiple regressions were performed on quality of life score with demographic and clinical variables as predictors. Results: The prevalence of depression in our study population is 25%. In patients with both CD and UC, depression is the most significant predictor to a poor HR-QOL (in CD, P = 8.22 × 10−6; in UC, P = 2.02 × 10−6). HR-QOL is weakly affected by disease activity (in CD, P = 0.110; in UC, P = 0.00492). In CD, biological use displays positive effect on HR-QOL (P = 0.00780). In total, the proportion of variance explained by all predictors is 61% for CD and 53% for UC, whereas the depression alone explains 44% and 36%. Conclusions: Our study demonstrates the importance of depression toward the quality of life in patients with inflammatory bowel disease. The diagnosis of depression should be actively sought out and treated in outpatient inflammatory bowel disease practices.

Genome-wide association study of N370S homozygous Gaucher disease reveals the candidacy of CLN8 gene as a genetic modifier contributing to extreme phenotypic variation†

Mutations in GBA1 gene result in defective acid β‐glucosidase and the complex phenotype of Gaucher disease (GD) related to the accumulation of glucosylceramide‐laden macrophages. The phenotype is highly variable even among patients harboring identical GBA1 mutations. We hypothesize that modifier gene(s) underlie phenotypic diversity in GD and performed a GWAS study in Ashkenazi Jewish patients with type 1 GD (GD1), homozygous for N370S mutation. Patients were assigned to mild, moderate, or severe disease categories using composite disease severity scoring systems. Whole‐genome genotyping for >500,000 SNPs was performed to search for association signals using OQLS algorithm in 139 eligible patients. Several SNPs in linkage disequilibrium within the CLN8 gene locus were associated with the GD1 severity: SNP rs11986414 was associated with GD1 severity at P value 1.26 × 10−6. Compared to mild disease, risk allele A at rs11986414 conferred an odds ratio of 3.72 for moderate/severe disease. Loss of function mutations in CLN8 causes neuronal ceroid‐lipofuscinosis, but our results indicate that its increased expression may protect against severe GD1. In cultured skin fibroblasts, the relative expression of CLN8 was higher in mild GD compared to severely affected patients, in whom CLN8 risk alleles were overrepresented. In an in vitro cell model of GD, CLN8 expression was increased, which was further enhanced in the presence of bioactive substrate, glucosylsphingosine. Taken together, CLN8 is a candidate modifier gene for GD1 that may function as a protective sphingolipid sensor and/or in glycosphingolipid trafficking. Future studies should explore the role of CLN8 in pathophysiology of GD. Am. J. Hematol., 2012.

The role of macrophage migration inhibitory factor in autoimmune liver disease

The role of the cytokine, macrophage migration inhibitory factor (MIF), and its receptor, CD74, was assessed in autoimmune hepatitis (AIH) and primary biliary cirrhosis (PBC). Two MIF promoter polymorphisms, a functional −794 CATT5‐8 microsatellite repeat (rs5844572) and a −173 G/C single‐nucleotide polymorphism (rs755622), were analyzed in DNA samples from over 500 patients with AIH, PBC, and controls. We found a higher frequency of the proinflammatory and high‐expression −794 CATT7 allele in AIH, compared to PBC, whereas lower frequency was found in PBC, compared to both AIH and healthy controls. MIF and soluble MIF receptor (CD74) were measured by enzyme‐linked immunosorbent assay in 165 serum samples of AIH, PBC, and controls. Circulating serum and hepatic MIF expression was elevated in patients with AIH and PBC versus healthy controls. We also identified a truncated circulating form of the MIF receptor, CD74, that is released from hepatic stellate cells and that binds MIF, neutralizing its signal transduction activity. Significantly higher levels of CD74 were found in patients with PBC versus AIH and controls. Conclusions: These data suggest a distinct genetic and immunopathogenic basis for AIH and PBC at the MIF locus. Circulating MIF and MIF receptor profiles distinguish PBC from the more inflammatory phenotype of AIH and may play a role in pathogenesis and as biomarkers of these diseases. (Hepatology 2014;59:580–591)