Gene-based copy number variation study reveals a microdeletion at 12q24 that influences height in the Korean population

Association of coffee consumption with type 2 diabetes and glycemic traits: a Mendelian randomization study

BACKGROUND/OBJECTIVES

Habitual coffee consumption was inversely associated with type 2 diabetes (T2D) and hyperglycemia in observational studies, but the causality of the association remains uncertain. This study tested a causal association of genetically predicted coffee consumption with T2D using the Mendelian randomization (MR) method.

SUBJECTS/METHODS

We used five single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs) associated with habitual coffee consumption in a previous genome-wide association study among Koreans. We analyzed the associations between IVs and T2D, fasting blood glucose (FBG), 2h-postprandial glucose (2h-PG), and glycated haemoglobin (HbA1C) levels. The MR results were further evaluated by standard sensitivity tests for possible pleiotropism.

RESULTS

MR analysis revealed that increased genetically predicted coffee consumption was associated with a reduced prevalence of T2D; ORs per one-unit increment of log-transformed cup per day of coffee consumption ranged from 0.75 (0.62-0.90) for the weighted mode-based method to 0.79 (0.62-0.99) for Wald ratio estimator. We also used the inverse-variance-weighted method, weighted median-based method, MR-Egger method, and MR-PRESSO method. Similarly, genetically predicted coffee consumption was inversely associated with FBG and 2h-PG levels but not with HbA1c. Sensitivity measures gave similar results without evidence of pleiotropy.

CONCLUSIONS

A genetic predisposition to habitual coffee consumption was inversely associated with T2D prevalence and lower levels of FBG and 2h-PG profiles. Our study warrants further exploration.

Bitter Taste Receptor TAS2R38 Genetic Variation (rs10246939), Dietary Nutrient Intake, and Bio-Clinical Parameters in Koreans

Differential bitterness perception associated with genetic polymorphism in the bitter taste receptor gene taste 2 receptor member 38 (TAS2R38) may influence an individual's food preferences, nutrition consumption, and eventually chronic nutrition-related disorders including cardiovascular disease. Therefore, the effect of genetic variations on nutritional intake and clinical markers needs to be elaborated for health and disease prevention. In this study, we conducted sex-stratified analysis to examine the association between genetic variant TAS2R38 rs10246939 A > G with daily nutritional intake, blood pressure, and lipid parameters in Korean adults (males = 1,311 and females = 2,191). We used the data from the Multi Rural Communities Cohort, Korean Genome and Epidemiology Study. Findings suggested that the genetic variant TAS2R38 rs10246939 was associated with dietary intake of micronutrients including calcium (adjusted p = 0.007), phosphorous (adjusted p = 0.016), potassium (adjusted p = 0.022), vitamin C (adjusted p = 0.009), and vitamin E (adjusted p = 0.005) in females. However, this genetic variant did not influence blood glucose, lipid profile parameters, and other blood pressure markers. These may suggest that this genetic variation is associated with nutritional intake, but its clinical effect was not found. More studies are needed to explore whether TAS2R38 genotype may be a potential predictive marker for the risk of metabolic diseases via modulation of dietary intake.

TAS1R2 sweet taste receptor genetic variation and dietary intake in Korean females

Taste receptor type 1, member 2 (TAS1R2) controls the oral sensing of sweetness. Genetic variations in TAS1R2 have been shown to be associated with differential sweetness intensity and varying carbohydrate intake levels among individuals. This study examined whether rs7534618 A > C in TAS1R2 is associated with dietary behavior and energy nutrient intake in Korean females. A cross-sectional design utilizing data from the Multi-Rural Communities Cohort Study, which was a nationwide epidemiological research project in Korea, was applied in this study. In total, 2198 females were analyzed to evaluate the differences in macronutrient intake levels and intake of carbohydrate-rich and sweet-tasting foods between the rs7534618 genotypes. The findings suggest that individuals with the CC minor genotype tended to have lower carbohydrate but higher fat intake than subjects with the A* genotype (p = 0.035 and p = 0.042, respectively). Subjects with the CC genotype also exhibited less intake of total grains but greater intake of bread than those with the A* genotype (p = 0.017 and p = 0.006, respectively). However, these observed associations were statistically modest (false discovery rate adjusted p > 0.05). In conclusion, TAS1R2 rs7534618 is not a decisive genetic modifier of nutrition and dietary intake in Korean females. However, given the paucity of studies, these putative associations between the TAS1R variation and dietary intake may be referred for further sensory genetic studies in Koreans.

Coffee Consumption, Genetic Polymorphisms, and the Risk of Type 2 Diabetes Mellitus: A Pooled Analysis of Four Prospective Cohort Studies

The association between coffee consumption and the risk of type 2 diabetes may vary by genetic variants. Our study addresses the question of whether the incidence of type 2 diabetes is related to the consumption of coffee and whether this relationship is modified by polymorphisms related to type 2 diabetes. We performed a pooled analysis of four Korean prospective studies that included 71,527 participants; median follow-up periods ranged between 2 and 13 years. All participants had completed a validated food-frequency questionnaire (FFQ) at baseline. The odds ratios (ORs) and 95% confidence intervals (CIs) for type 2 diabetes were calculated using logistic regression models. The ORs were combined using a fixed or random effects model depending on the heterogeneity across the studies. Compared with 0 to <0.5 cups/day of coffee consumption, the OR for type 2 diabetes was 0.89 (95% CI: 0.80-0.98, p for trend = 0.01) for ≥3 cups/day of coffee consumption. We did not observe significant interactions by five single nucleotide polymorphisms (SNPs) related to type 2 diabetes (CDKAL1 rs7756992, CDKN2A/B rs10811661, KCNJ11 rs5215, KCNQ1 rs163184, and PEPD rs3786897) in the association between coffee and the risk of type 2 diabetes. We found that coffee consumption was inversely associated with the risk of type 2 diabetes.

Associations of ORMDL1 gene copy number variations with growth traits in four Chinese sheep breeds

Copy number variations (CNVs) are gains and losses of genomic sequence of more than 50 bp between two individuals of a species. Also, CNV is considered to be one of the main elements affecting the phenotypic diversity and evolutionary adaptation of animals. ORMDL sphingolipid biosynthesis regulator 1 (ORMDL1) is a protein-coding gene associated with diseases and development. In our study, the polymorphism of ORMDL1 gene copy numbers in four Chinese sheep breeds (abbreviated CK, HU, STH, and LTH) was detected. In addition, we analyzed the transcriptional expression level of ORMDL1 gene in different tissues of sheep and examined the association of ORMDL1 CNV with growth traits. The statistical analysis revealed that ORMDL1 CNV was remarkably correlated with body height, heart girth, and circumference of cannon bone in HU sheep (P<0.05), and there are significant effects on body weight, body height, body length, chest depth, and height of hip cross in STH sheep (P<0.05). In conclusion, our results provide a basis for the relationship between CNV of ORMDL1 gene and sheep growth traits, suggesting that ORMDL1 CNV may be considered a promising marker for the molecular breeding of Chinese sheep.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

Variation in the TAS2R38 Bitterness Receptor Gene Was Associated with Food Consumption and Obesity Risk in Koreans

Bitterness-sensing protein taste receptor type-2 member 38 (TAS2R38, T2R38) mediates taste perception and various physiological responses, including energy- and adiposity-related mechanisms. This study examined whether the genetic variant rs10246939 C > T in TAS2R38 was associated with food intake and body size as well as obesity risk. Data from the Korean Multi-Rural Communities Cohort study (1338 males and 2229 females) were analyzed to obtain the intake of six food groups, alcohol consumption, smoking status and anthropometric measurements, including height, weight, waist and hip circumference, and body mass index (BMI), according to the rs10246939 genotype. Findings suggested that females with the TT genotype consumed more fruit (adjusted p = 0.025) and had significantly higher body weights (adjusted p = 0.046) and BMIs (adjusted p = 0.003) than individuals with other genotypes. Having the TT genotype also increased the risk of obesity by 1.75-fold (95% confidence interval: 1.31–2.36) in females. The genetic variation had a minimal influence on the males’ dietary intake, but tended to increase the adiposity risk. In conclusion, TAS2R38 rs10246939 variation was associated with Koreans’ dietary intake and increased their risk of obesity. Although more detailed statistical analyses in the larger cohort are required, current study suggested that, as a genetic predictive marker, TAS2R38 bitterness receptor variations may have a large implication in obesity prevention and treatment.

The role of GPCRs in bone diseases and dysfunctions

The superfamily of G protein-coupled receptors (GPCRs) contains immense structural and functional diversity and mediates a myriad of biological processes upon activation by various extracellular signals. Critical roles of GPCRs have been established in bone development, remodeling, and disease. Multiple human GPCR mutations impair bone development or metabolism, resulting in osteopathologies. Here we summarize the disease phenotypes and dysfunctions caused by GPCR gene mutations in humans as well as by deletion in animals. To date, 92 receptors (5 glutamate family, 67 rhodopsin family, 5 adhesion, 4 frizzled/taste2 family, 5 secretin family, and 6 other 7TM receptors) have been associated with bone diseases and dysfunctions (36 in humans and 72 in animals). By analyzing data from these 92 GPCRs, we found that mutation or deletion of different individual GPCRs could induce similar bone diseases or dysfunctions, and the same individual GPCR mutation or deletion could induce different bone diseases or dysfunctions in different populations or animal models. Data from human diseases or dysfunctions identified 19 genes whose mutation was associated with human BMD: 9 genes each for human height and osteoporosis; 4 genes each for human osteoarthritis (OA) and fracture risk; and 2 genes each for adolescent idiopathic scoliosis (AIS), periodontitis, osteosarcoma growth, and tooth development. Reports from gene knockout animals found 40 GPCRs whose deficiency reduced bone mass, while deficiency of 22 GPCRs increased bone mass and BMD; deficiency of 8 GPCRs reduced body length, while 5 mice had reduced femur size upon GPCR deletion. Furthermore, deficiency in 6 GPCRs induced osteoporosis; 4 induced osteoarthritis; 3 delayed fracture healing; 3 reduced arthritis severity; and reduced bone strength, increased bone strength, and increased cortical thickness were each observed in 2 GPCR-deficiency models. The ever-expanding number of GPCR mutation-associated diseases warrants accelerated molecular analysis, population studies, and investigation of phenotype correlation with SNPs to elucidate GPCR function in human diseases.

Analysis of copy number variation in the NEDD4L gene potentially implicated in body height in the Japanese population

Recently it has been recognized that a considerable number of copy number variations (CNVs) are associated with diseases and other complex human traits. In our previous study, we developed a simple quantitative real-time PCR (Q-PCR) method for analysis of CNV copy number, which had the advantage of obviating the need for reference DNA with a known copy number. Using DNA samples obtained from 231 Japanese individuals, we applied this method for analyzing the copy number of a candidate CNV associated with body height, located in the neural precursor cell expressed, developmentally down-regulated 4-like, E3 ubiquitin protein ligase (NEDD4L) gene. In addition, the appropriateness of the results was evaluated and confirmed by quantification of amplicons with an Agilent 2100 Bioanalyzer. The NEDD4L gene encodes a member of the Nedd4 family of HECT domain E3 ubiquitin ligases. The target CNV located in the intron has been found to be significantly associated with height variation in Chinese. However, it remains unknown whether such an association exists in other populations, including Japanese. Analysis of the correlations between copy number and body height using ANOVA revealed no statistically significant correlations in Japanese.

The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits

We introduce the design and implementation of a new array, the Korea Biobank Array (referred to as KoreanChip), optimized for the Korean population and demonstrate findings from GWAS of blood biochemical traits. KoreanChip comprised >833,000 markers including >247,000 rare-frequency or functional variants estimated from >2,500 sequencing data in Koreans. Of the 833 K markers, 208 K functional markers were directly genotyped. Particularly, >89 K markers were presented in East Asians. KoreanChip achieved higher imputation performance owing to the excellent genomic coverage of 95.38% for common and 73.65% for low-frequency variants. From GWAS (Genome-wide association study) using 6,949 individuals, 28 associations were successfully recapitulated. Moreover, 9 missense variants were newly identified, of which we identified new associations between a common population-specific missense variant, rs671 (p.Glu457Lys) of ALDH2, and two traits including aspartate aminotransferase (P = 5.20 × 10⁻¹³) and alanine aminotransferase (P = 4.98 × 10⁻⁸). Furthermore, two novel missense variants of GPT with rare frequency in East Asians but extreme rarity in other populations were associated with alanine aminotransferase (rs200088103; p.Arg133Trp, P = 2.02 × 10⁻⁹ and rs748547625; p.Arg143Cys, P = 1.41 × 10⁻⁶). These variants were successfully replicated in 6,000 individuals (P = 5.30 × 10⁻⁸ and P = 1.24 × 10⁻⁶). GWAS results suggest the promising utility of KoreanChip with a substantial number of damaging variants to identify new population-specific disease-associated rare/functional variants.

Simple screening method for copy number variations associated with physical features

Recent studies of copy number variations (CNVs) associated with physical features, such as body mass index, body height or bone length, have suggested that such CNVs could serve as markers in forensic cases involving unidentified individuals. However, the process of cataloging CNVs has been slow because of the cumbersome nature and low reliability of the procedures involved. Here we describe a simple quantitative real-time PCR (Q-PCR) method for screening of medicolegally useful CNVs, which does not require reference DNA with known copy number. The first step is to prepare a chimeric plasmid vector including one copy each of the single-copy gene-specific sequence as the internal standard, and the target CNV-specific sequence. To assess the validity of this new method, we analyzed CNVs in the LTBP1 and ETV6 gene regions, both of which are candidate CNVs associated with body height. The PCR efficiencies for the single-copy (reference) gene and the target CNV were similar, indicating that quantitation was reliable. Furthermore, simulated analysis of the LTBP1 CNV using mock samples prepared by mixing vectors in varying proportions showed that this analytical method allowed correct determination of the LTBP1 copy number. These results demonstrated that our simple method has considerable potential for screening of trait-related CNVs that would be useful for forensic casework.

Identification of a Copy Number Variation on Chromosome 20q13.12 Associated with Osteoporotic Fractures in the Korean Population

Osteoporotic fractures (OFs) are critical hard outcomes of osteoporosis and are characterized by decreased bone strength induced by low bone density and microarchitectural deterioration in bone tissue. Most OFs cause acute pain, hospitalization, immobilization, and slow recovery in patients and are associated with increased mortality. A variety of genetic studies have suggested associations of genetic variants with the risk of OF. Genome-wide association studies have reported various single-nucleotide polymorphisms and copy number variations (CNVs) in European and Asian populations. To identify CNV regions associated with OF risk, we conducted a genome-wide CNV study in a Korean population. We performed logistic regression analyses in 1,537 Korean subjects (299 OF cases and 1,238 healthy controls) and identified a total of 8 CNV regions significantly associated with OF (p < 0.05). Then, one CNV region located on chromosome 20q13.12 was selected for experimental validation. The selected CNV region was experimentally validated by quantitative polymerase chain reaction. The CNV region of chromosome 20q13.12 is positioned upstream of a family of long non-coding RNAs, LINC01260. Our findings could provide new information on the genetic factors associated with the risk of OF.

GPR133 (ADGRD1), an adhesion G-protein-coupled receptor, is necessary for glioblastoma growth

Glioblastoma (GBM) is a deadly primary brain malignancy with extensive intratumoral hypoxia. Hypoxic regions of GBM contain stem-like cells and are associated with tumor growth and angiogenesis. The molecular mechanisms that regulate tumor growth in hypoxic conditions are incompletely understood. Here, we use primary human tumor biospecimens and cultures to identify GPR133 (ADGRD1), an orphan member of the adhesion family of G-protein-coupled receptors, as a critical regulator of the response to hypoxia and tumor growth in GBM. GPR133 is selectively expressed in CD133+ GBM stem cells (GSCs) and within the hypoxic areas of PPN in human biospecimens. GPR133 mRNA is transcriptionally upregulated by hypoxia in hypoxia-inducible factor 1α (Hif1α)-dependent manner. Genetic inhibition of GPR133 with short hairpin RNA reduces the prevalence of CD133+ GSCs, tumor cell proliferation and tumorsphere formation in vitro. Forskolin rescues the GPR133 knockdown phenotype, suggesting that GPR133 signaling is mediated by cAMP. Implantation of GBM cells with short hairpin RNA-mediated knockdown of GPR133 in the mouse brain markedly reduces tumor xenograft formation and increases host survival. Analysis of the TCGA data shows that GPR133 expression levels are inversely correlated with patient survival. These findings indicate that GPR133 is an important mediator of the hypoxic response in GBM and has significant protumorigenic functions. We propose that GPR133 represents a novel molecular target in GBM and possibly other malignancies where hypoxia is fundamental to pathogenesis.

International Union of Basic and Clinical Pharmacology. XCIV. Adhesion G protein-coupled receptors

The Adhesion family forms a large branch of the pharmacologically important superfamily of G protein-coupled receptors (GPCRs). As Adhesion GPCRs increasingly receive attention from a wide spectrum of biomedical fields, the Adhesion GPCR Consortium, together with the International Union of Basic and Clinical Pharmacology Committee on Receptor Nomenclature and Drug Classification, proposes a unified nomenclature for Adhesion GPCRs. The new names have ADGR as common dominator followed by a letter and a number to denote each subfamily and subtype, respectively. The new names, with old and alternative names within parentheses, are: ADGRA1 (GPR123), ADGRA2 (GPR124), ADGRA3 (GPR125), ADGRB1 (BAI1), ADGRB2 (BAI2), ADGRB3 (BAI3), ADGRC1 (CELSR1), ADGRC2 (CELSR2), ADGRC3 (CELSR3), ADGRD1 (GPR133), ADGRD2 (GPR144), ADGRE1 (EMR1, F4/80), ADGRE2 (EMR2), ADGRE3 (EMR3), ADGRE4 (EMR4), ADGRE5 (CD97), ADGRF1 (GPR110), ADGRF2 (GPR111), ADGRF3 (GPR113), ADGRF4 (GPR115), ADGRF5 (GPR116, Ig-Hepta), ADGRG1 (GPR56), ADGRG2 (GPR64, HE6), ADGRG3 (GPR97), ADGRG4 (GPR112), ADGRG5 (GPR114), ADGRG6 (GPR126), ADGRG7 (GPR128), ADGRL1 (latrophilin-1, CIRL-1, CL1), ADGRL2 (latrophilin-2, CIRL-2, CL2), ADGRL3 (latrophilin-3, CIRL-3, CL3), ADGRL4 (ELTD1, ETL), and ADGRV1 (VLGR1, GPR98). This review covers all major biologic aspects of Adhesion GPCRs, including evolutionary origins, interaction partners, signaling, expression, physiologic functions, and therapeutic potential.

The genome-wide landscape of copy number variations in the MUSGEN study provides evidence for a founder effect in the isolated Finnish population

Here we characterized the genome-wide architecture of copy number variations (CNVs) in 286 healthy, unrelated Finnish individuals belonging to the MUSGEN study, where molecular background underlying musical aptitude and related traits are studied. By using Illumina HumanOmniExpress-12v.1.0 beadchip, we identified 5493 CNVs that were spread across 467 different cytogenetic regions, spanning a total size of 287.83 Mb (∼9.6% of the human genome). Merging the overlapping CNVs across samples resulted in 999 discrete copy number variable regions (CNVRs), of which ∼6.9% were putatively novel. The average number of CNVs per person was 20, whereas the average size of CNV per locus was 52.39 kb. Large CNVs (>1 Mb) were present in 4% of the samples. The proportion of homozygous deletions in this data set (∼12.4%) seemed to be higher when compared with three other populations. Interestingly, several CNVRs were significantly enriched in this sample set, whereas several others were totally depleted. For example, a CNVR at chr2p22.1 intersecting GALM was more common in this population (P=3.3706 × 10(-44)) than in African and other European populations. The enriched CNVRs, however, showed no significant association with music-related phenotypes. Moreover, the most common CNV locations in world's normal population cohorts (6q14.1, 11q11) were overrepresented in this population. Thus, the genome-wide CNV investigation in this Finnish sample set demonstrated features that are characteristic to isolated populations. Novel CNVRs and the functional implications of CNVs revealed in this study elucidate structural variation present in this population isolate, and may also serve as candidate gene loci for music-related traits.

Overweight condition and waist circumference and a candidate gene within the 12q24 locus

Aims

Obesity and obesity-associated phenotypes are linked to the chromosome12q24 locus, the non-insulin-dependent-diabetes 2 (NIDDM2) locus. The gene of proteasome modulator 9 (PSMD9) lies in the NIDDM2 region and is linked to type 2 diabetes (T2D), microvascular and macrovascular complications of T2D.

We aimed at studying whether the PSMD9 T2D risk single nucleotide polymorphisms (SNPs) IVS3+nt460, IVS3+nt437, and 197G are linked to obesity, overweight status and waist circumference in Italian T2D families.

Methods and results

We screened 200 Italians T2D siblings/families for PSMD9 variants. Using Merlin software, we performed non-parametric linkage analysis to test for linkage with obesity and overweight condition and variance component analysis to test for linkage with waist circumference in our T2D siblings/families dataset.

Our study shows that the PSMD9 SNPs IVS3+nt460, IVS3+nt437, and 197G are in linkage with overweight condition and waist circumference in Italians. The statistical power tests performed via simulations on real data confirm that the results are not due to random chance.

Conclusions

In summary, the linkage strategy using a homogeneous family/subject dataset can identify a gene contributing to a complex trait. PMSD9 may be at least one of the genes responsible for the linkage to obesity and obesity-associated phenotypes at the locus 12q24 in other populations.