Person:

Liang, Liming

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.661028), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6610213; SOX6, p = 6.4610210) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.

Dyslipidemia has been associated with type 2 diabetes, but it remains unclear whether dyslipidemia plays a causal role in type 2 diabetes. We aimed to examine the association between the genetic predisposition to dyslipdemia and type 2 diabetes risk. The current study included 2,447 patients with type 2 diabetes and 3,052 control participants of European ancestry from the Nurses’ Health Study and the Health Professionals Follow-up Study. Genetic predisposition to dyslipidemia was estimated by three genotype scores of lipids (LDL cholesterol, HDL cholesterol, and triglycerides) on the basis of the established loci for blood lipids. Linear relation analysis indicated that the HDL cholesterol and triglyceride genotype scores, but not the LDL cholesterol genotype score, were linearly related to elevated type 2 diabetes risk. Each point of the HDL cholesterol and triglyceride genotype scores was associated with a 3% (odds ratio [OR] 1.03 [95% CI 1.01–1.04]) and a 2% (1.02 [1.00–1.04]) increased risk of developing type 2 diabetes, respectively. The ORs were 1.39 (1.17–1.65) and 1.19 (1.01–1.41) for type 2 diabetes by comparing extreme quartiles of the HDL cholesterol genotype score and triglyceride genotype score, respectively. In conclusion, genetic predisposition to low HDL cholesterol or high triglycerides is related to elevated type 2 diabetes risk.

To pinpoint the exact location of the etiological variant/s present at 1q21.1 harboring FCRL1-5 and CD5L genes, we carried out a refined association study in the entire FCRL region in 1,536 patients with Graves’ disease (GD) and 1,516 sex-matched controls by imputation analysis, logistic regression, and cis-eQTL analysis. Among 516 SNPs with P<0.05 in the initial GWAS scan, the strongest signals associated with GD and correlated to FCRL3 expression were located at a cluster of SNPs including rs7528684 and rs3761959. And the allele-specific effects for rs3761959 and rs7528684 on FCRL3 expression level revealed that the risk alleles A of rs3761959 and C of rs7528684 were correlated with the elevated expression level of FCRL3 whether in PBMCs or its subsets, especially in (CD19^+) B cells and (CD8^+) T subsets. Next, the combined analysis with 5,300 GD cases and 4,916 control individuals confirmed FCRL3 was a susceptibility gene of GD in Chinese Han populations, and rs3761959 and rs7528684 met the genome-wide association significance level ((P_{combined}) = 2.27×(10^{−12}) and 7.11×(10^{−13}), respectively). Moreover, the haplotypes with the risk allele A of rs3761959 and risk allele C of rs7528684 were associated with GD risk. Finally, our epigenetic analysis suggested the disease-associated C allele of rs7528684 increased affinity for NF-KB transcription factor. Above data indicated that FCRL3 gene and its proxy SNP rs7528684 may be involved in the pathogenesis of GD by excessive inhibiting B cell receptor signaling and the impairment of suppressing function of Tregs.

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new genome-wide association studies (GWAS) and one prior scan, totaling 3,857 cases and 7,666 controls of European ancestry, with additional genotyping of nine promising SNPs in 1,359 cases and 4,557 controls. In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; (P=2.33x10^{-21})), rs2523607 at 6p21.33 (HLA-B; (2.40x10^{-10})), rs79480871 at 2p23.3 (NCOA1; (P=4.23x10^{-8})), and two independent SNPs, rs13255292 and rs4733601, at 8q24.21 (PVT1; (P=9.98x10^{-13}) and (P=3.63x10^{-11}), respectively). These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL.

Context Anxiety disorders are common, with a lifetime prevalence of 20% in the U.S., and are responsible for substantial burdens of disability, missed work days and health care utilization. To date, no causal genetic variants have been identified for anxiety, anxiety disorders, or related traits. Objective: To investigate whether a phobic anxiety symptom score was associated with 3 alternative polygenic risk scores, derived from external genome-wide association studies of anxiety, an internally estimated agnostic polygenic score, or previously identified candidate genes. Design: Longitudinal follow-up study. Using linear and logistic regression we investigated whether phobic anxiety was associated with polygenic risk scores derived from internal, leave-one out genome-wide association studies, from 31 candidate genes, and from out-of-sample genome-wide association weights previously shown to predict depression and anxiety in another cohort. Setting and Participants: Study participants (n = 11,127) were individuals from the Nurses' Health Study and Health Professionals Follow-up Study. Main Outcome Measure: Anxiety symptoms were assessed via the 8-item phobic anxiety scale of the Crown Crisp Index at two time points, from which a continuous phenotype score was derived. Results: We found no genome-wide significant associations with phobic anxiety. Phobic anxiety was also not associated with a polygenic risk score derived from the genome-wide association study beta weights using liberal p-value thresholds; with a previously published genome-wide polygenic score; or with a candidate gene risk score based on 31 genes previously hypothesized to predict anxiety. Conclusion: There is a substantial gap between twin-study heritability estimates of anxiety disorders ranging between 20–40% and heritability explained by genome-wide association results. New approaches such as improved genome imputations, application of gene expression and biological pathways information, and incorporating social or environmental modifiers of genetic risks may be necessary to identify significant genetic predictors of anxiety.

Atopic dermatitis (AD) is the most common dermatological disease of childhood. Many children with AD have asthma and AD shares regions of genetic linkage with psoriasis, another chronic inflammatory skin disease. We present here a genome-wide association study (GWAS) of childhood-onset AD in 1563 European cases with known asthma status and 4054 European controls. Using Illumina genotyping followed by imputation, we generated 268 034 consensus genotypes and in excess of 2 million single nucleotide polymorphisms (SNPs) for analysis. Association signals were assessed for replication in a second panel of 2286 European cases and 3160 European controls. Four loci achieved genome-wide significance for AD and replicated consistently across all cohorts. These included the epidermal differentiation complex (EDC) on chromosome 1, the genomic region proximal to LRRC32 on chromosome 11, the RAD50/IL13 locus on chromosome 5 and the major histocompatibility complex (MHC) on chromosome 6; reflecting action of classical HLA alleles. We observed variation in the contribution towards co-morbid asthma for these regions of association. We further explored the genetic relationship between AD, asthma and psoriasis by examining previously identified susceptibility SNPs for these diseases. We found considerable overlap between AD and psoriasis together with variable coincidence between allergic rhinitis (AR) and asthma. Our results indicate that the pathogenesis of AD incorporates immune and epidermal barrier defects with combinations of specific and overlapping effects at individual loci.

While research has suggested that being married may confer a health advantage, few studies to date have investigated the role of marital status in the development of type 2 diabetes. We examined whether men who are not married have increased risk of incident type 2 diabetes in the Health Professionals Follow-up Study. Men (n = 41,378) who were free of T2D in 1986, were followed for ≤22 years with biennial reports of T2D, marital status and covariates. Cox proportional hazard models were used to compare risk of incident T2D by marital status (married vs unmarried and married vs never married, divorced/separated, or widowed). There were 2,952 cases of incident T2D. Compared to married men, unmarried men had a 16% higher risk of developing T2D (95%CI:1.04,1.30), adjusting for age, family history of diabetes, ethnicity, lifestyle and body mass index (BMI). Relative risks (RR) for developing T2D differed for divorced/separated (1.09 [95%CI: 0.94,1.27]), widowed (1.29 [95%CI:1.06,1.57]), and never married (1.17 [95%CI:0.91,1.52]) after adjusting for age, family history of diabetes and ethnicity. Adjusting for lifestyle and BMI, the RR for T2D associated with widowhood was no longer significant (RR:1.16 [95%CI:0.95,1.41]). When allowing for a 2-year lag period between marital status and disease, RRs of T2D for widowers were augmented and borderline significant (RR:1.24 [95%CI:1.00,1.54]) after full adjustment. In conclusion, not being married, and more specifically, widowhood was more consistently associated with an increased risk of type 2 diabetes in men and this may be mediated, in part, through unfavorable changes in lifestyle, diet and adiposity.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10−8). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

OBJECTIVE Genetic variants near IRS1 are associated with features of the metabolic syndrome (MetS). We examined whether genetic variants near IRS1 might modulate the effects of diets varying in fat content on the MetS status in a 2-year weight-loss trial. RESEARCH DESIGN AND METHODS Two variants near IRS1, rs1522813 and rs2943641, were genotyped in 738 overweight/obese adults (age 60 ± 9 years; BMI 32.7 ± 3.9 kg/m2) randomly assigned to one of four weight-loss diets (a deficit of 750 kcal/day of caloric intake from baseline) varying in macronutrient contents for 2 years. We compared MetS status of high-fat (40% of caloric intake; n = 370) and low-fat (20% caloric intake; n = 368) diet groups differentiated by genotypes (rs1522813 A-allele carriers and noncarriers and rs2943641T-allele carriers and noncarriers). RESULTS Among rs1522813 A-allele carriers, the reversion rates of the MetS were higher in the high-fat diet group than those in the low-fat diet group over the 2-year intervention (P = 0.002), while no significant difference between diet groups was observed among noncarriers (P = 0.27). The genetic modulation on dietary effect was independent of weight changes. The odds ratio (OR) for the 2-year reversion of the MetS was 2.88 (95% CI 1.25–6.67) comparing the high-fat and low-fat diets among rs1522813 A-allele carriers, while the corresponding OR was 0.83 (0.36–1.92) in noncarriers. The variant rs2943641 was not observed to modulate dietary effects on the MetS status. CONCLUSIONS Our data suggest that high-fat weight-loss diets might be more effective in the management of the MetS compared with low-fat diets among individuals with the A-allele of the rs1522813 variant near IRS1.

A previous genome-wide association study in asthma revealed putative associations that merit further investigation. In this study, the genome-wide significant associations of SNPs at the 5% false discovery rate were examined in independent groups of severe asthmatics. The panel consisted of 397 severe asthmatic adults, 116 severe asthmatic children, and a collection of 207 family-trios with an asthmatic proband. Three SNPs in the PDCD4 gene (rs6585018:G>A, rs1322997:C>A, and rs34104444:G>A) were significantly associated with severe childhood asthma (P values: 0.003, 0.002, 0.004) and total immunoglobulin E (IgE) levels (P values: 0.034, 0.041, 0.052). In an independent group of 234 asthmatic children and 652 controls, PDCD4 SNPs rs1407696:T>G and rs11195360:T>C were associated with total IgE levels (P values: 0.006, 0.014). In silico analysis of PDCD4 locus showed that rs6585018:G>A had the potential to affect MYB transcription factor binding, shown to act as a PDCD4-transcription inducer. Electromobility shift assays and reporter assays revealed that rs6585018:G>A alters MYB binding thereby influencing the expression of PDCD4. SNPs within MYB itself confer susceptibility to eosinophilia and asthma. Our association between a variant MYB binding site in PDCD4 and the severest form of childhood asthma therefore suggests that PDCD4 is a novel molecule of importance to asthmatic inflammatory responses.