Genome-wide scan for serum ghrelin detects linkage on chromosome 1p36 in Hispanic children: results from the Viva La Familia study

Genetic variants in ALDH1L1 and GLDC influence the serine-to-glycine ratio in Hispanic children

BACKGROUND

Glycine is a proteogenic amino acid that is required for numerous metabolic pathways, including purine, creatine, heme, and glutathione biosynthesis. Glycine formation from serine, catalyzed by serine hydroxy methyltransferase, is the major source of this amino acid in humans. Our previous studies in a mouse model have shown a crucial role for the 10-formyltetrahydrofolate dehydrogenase enzyme in serine-to-glycine conversion.

OBJECTIVES

We sought to determine the genomic influence on the serine-glycine ratio in 803 Hispanic children from 319 families of the Viva La Familia cohort.

METHODS

We performed a genome-wide association analysis for plasma serine, glycine, and the serine-glycine ratio in Sequential Oligogenic Linkage Analysis Routines while accounting for relationships among family members.

RESULTS

All 3 parameters were significantly heritable (h2 = 0.22-0.78; P < 0.004). The strongest associations for the serine-glycine ratio were with single nucleotide polymorphisms (SNPs) in aldehyde dehydrogenase 1 family member L1 (ALDH1L1) and glycine decarboxylase (GLDC) and for glycine with GLDC (P < 3.5 × 10-8; effect sizes, 0.03-0.07). No significant associations were found for serine. We also conducted a targeted genetic analysis with ALDH1L1 exonic SNPs and found significant associations between the serine-glycine ratio and rs2886059 (β = 0.68; SE, 0.25; P = 0.006) and rs3796191 (β = 0.25; SE, 0.08; P = 0.003) and between glycine and rs3796191 (β = -0.08; SE, 0.02; P = 0.0004). These exonic SNPs were further associated with metabolic disease risk factors, mainly adiposity measures (P < 0.006). Significant genetic and phenotypic correlations were found for glycine and the serine-glycine ratio with metabolic disease risk factors, including adiposity, insulin sensitivity, and inflammation-related phenotypes [estimate of genetic correlation = -0.37 to 0.35 (P < 0.03); estimate of phenotypic correlation = -0.19 to 0.13 (P < 0.006)]. The significant genetic correlations indicate shared genetic effects among glycine, the serine-glycine ratio, and adiposity and insulin sensitivity phenotypes.

CONCLUSIONS

Our study suggests that ALDH1L1 and GLDC SNPs influence the serine-to-glycine ratio and metabolic disease risk.

Genome-wide analyses of multiple obesity-related cytokines and hormones informs biology of cardiometabolic traits

BACKGROUND

A complex set of perturbations occur in cytokines and hormones in the etiopathogenesis of obesity and related cardiometabolic conditions such as type 2 diabetes (T2D). Evidence for the genetic regulation of these cytokines and hormones is limited, particularly in African-ancestry populations. In order to improve our understanding of the biology of cardiometabolic traits, we investigated the genetic architecture of a large panel of obesity- related cytokines and hormones among Africans with replication analyses in African Americans.

METHODS

We performed genome-wide association studies (GWAS) in 4432 continental Africans, enrolled from Ghana, Kenya, and Nigeria as part of the Africa America Diabetes Mellitus (AADM) study, for 13 obesity-related cytokines and hormones, including adipsin, glucose-dependent insulinotropic peptide (GIP), glucagon-like peptide-1 (GLP-1), interleukin-1 receptor antagonist (IL1-RA), interleukin-6 (IL-6), interleukin-10 (IL-10), leptin, plasminogen activator inhibitor-1 (PAI-1), resistin, visfatin, insulin, glucagon, and ghrelin. Exact and local replication analyses were conducted in African Americans (n = 7990). The effects of sex, body mass index (BMI), and T2D on results were investigated through stratified analyses.

RESULTS

GWAS identified 39 significant (P value < 5 × 10-8) loci across all 13 traits. Notably, 14 loci were African-ancestry specific. In this first GWAS for adipsin and ghrelin, we detected 13 and 4 genome-wide significant loci respectively. Stratified analyses by sex, BMI, and T2D showed a strong effect of these variables on detected loci. Eight novel loci were successfully replicated: adipsin (3), GIP (1), GLP-1 (1), and insulin (3). Annotation of these loci revealed promising links between these adipocytokines and cardiometabolic outcomes as illustrated by rs201751833 for adipsin and blood pressure and locus rs759790 for insulin level and T2D in lean individuals.

CONCLUSIONS

Our study identified genetic variants underlying variation in multiple adipocytokines, including the first loci for adipsin and ghrelin. We identified population differences in variants associated with adipocytokines and highlight the importance of stratification for discovery of loci. The high number of African-specific loci detected emphasizes the need for GWAS in African-ancestry populations, as these loci could not have been detected in other populations. Overall, our work contributes to the understanding of the biology linking adipocytokines to cardiometabolic traits.

Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children

CONTEXT

Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied.

OBJECTIVE

We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children.

METHODS

We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry.

RESULTS

Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7).

CONCLUSIONS

We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

Mendelian Inconsistent Signatures from 1314 Ancestrally Diverse Family Trios Distinguish Biological Variation from Sequencing Error

Next-generation sequencing enables advances in the clinical application of genomics by providing high-throughput detection of genomic variation. However, next-generation sequencing technologies, especially whole-genome sequencing (WGS), are often associated with a high false-positive rate. Trio-based WGS can contribute significantly towards improved quality control methods. Mendelian-inconsistent calls (MIC) in parent–child trios are commonly attributed to erroneous sequencing calls, as the true de novo mutation rate is extremely low compared with MIC incidence. Here, we analyzed WGS data from 1314 mother, father, and child trios across ethnically diverse populations with the goal of characterizing MIC. Genotype calls in a trio can be used to assign different signatures to MIC. MIC occur more frequently within repeats but show varying distribution and error mechanisms across repeat types. MIC are enriched within poly-A/T runs in short interspersed nuclear elements. Alignability scores, allele balance, and relative parental read depth vary among MIC signatures and these differences should be considered when designing filters for MIC reduction. MIC cluster in germline deletions and these MIC also segregate with population. Our results provide a basis for making decisions on how each MIC type should be evaluated before discarding them as errors or including them in alternative applications. With the reduction of sequencing cost, family trio whole genome and exome analysis are being performed more routinely in clinical practice. We provide a reference that can be used for annotating MIC with their frequencies in a larger population to aid in the filtering of candidate de novo mutations.

The Relevance of Genomic Signatures at Adhesion GPCR Loci in Humans

Adhesion G protein-coupled receptors (aGPCRs) have a long evolutionary history dating back to very basal unicellular eukaryotes. Almost every vertebrate is equipped with a set of different aGPCRs. Genomic sequence data of several hundred extinct and extant species allows for reconstruction of aGPCR phylogeny in vertebrates and non-vertebrates in general but also provides a detailed view into the recent evolutionary history of human aGPCRs. Mining these sequence sources with bioinformatic tools can unveil many facets of formerly unappreciated aGPCR functions. In this review, we extracted such information from the literature and open public sources and provide insights into the history of aGPCR in humans. This includes comprehensive analyses of signatures of selection, variability of human aGPCR genes, and quantitative traits at human aGPCR loci. As indicated by a large number of genome-wide genotype-phenotype association studies, variations in aGPCR contribute to specific human phenotypes. Our survey demonstrates that aGPCRs are significantly involved in adaptation processes, phenotype variations, and diseases in humans.

Appetite-related peptides in childhood and adolescence: role of ghrelin, PYY, and GLP-1

During childhood and adolescence, a number of factors, including age, puberty, sex, race, and body composition, may contribute to differences in satiety, food intake, and appetite-related peptides. These peptides include the orexigenic peptide ghrelin and anorexigenic gut peptides peptide YY (PYY) and glucagon-like peptide-1 (GLP-1). For example, lower fasting ghrelin levels, lower postprandial ghrelin suppression, and blunted PYY and GLP-1 responses to food intake could contribute to a dysregulation of appetite in already obese children and adolescents. Whereas, changes in these peptides observed during puberty could facilitate growth. A greater understanding of the major moderating factors of appetite-related peptides in the pediatric population is essential to improve interpretation of study findings and for effective tailoring of strategies targeting appetite control to individuals. While more studies are needed, there is some evidence to suggest that exercise-based lifestyle interventions could be a potential therapeutic strategy to improve appetite-peptide profiles in overweight and obese children and adolescents. The aim of this review is (i) to discuss the potential moderating factors of ghrelin, PYY, and GLP-1, including age and puberty, sex, race and body composition; and (ii) to examine the effects of exercise interventions on these appetite-related gut peptides in children and adolescents.

Significant genotype by diet (G × D) interaction effects on cardiometabolic responses to a pedigree-wide, dietary challenge in vervet monkeys (Chlorocebus aethiops sabaeus)

Nutrient composition of a diet (D) has been shown to interact with genetic predispositions (G) to affect various lipid phenotypes. Our aim in this study was to confirm G × D interaction and determine whether the interaction extends to other cardiometabolic risk factors such as glycemic measures and body weight. Subjects were vervet monkeys (Chlorocebus aethiops sabaeus; n = 309) from a multigenerational pedigreed colony initially fed with a plant-based diet, standard primate diet (18% calories from protein, 13% from fat, and 69% from carbohydrates), and subsequently challenged for 8 weeks with a diet modeled on the typical American diet (18% calories from protein, 35% from fat, and 47% from carbohydrates). Our results showed that although exposure to the challenge diet did not result in significant changes in weight, most lipid and glycemic biomarkers moved in an adverse direction (P < 0.01). Quantitative genetic analyses showed that cardiometabolic phenotypes were significantly heritable under both dietary conditions (P < 0.05), and there was significant evidence of G × D interaction for these phenotypes. We observed significant differences in the additive genetic variances for most lipid phenotypes (P < 10(-4) ), indicating that the magnitude of genetic effects varies by diet. Furthermore, genetic correlations between diets differed significantly from 1 with respect to insulin, body weight, and some lipid phenotypes (P < 0.01). This implied that distinct genetic effects are involved in the regulation of these phenotypes under the two dietary conditions. These G × D effects confirm and extend previous observations in baboons (Papio sp.) and suggest that mimicking the typical human nutritional environment can reveal genetic influences that might not be observed in animals consuming standard, plant-based diets.

Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population

Genetic variants responsible for susceptibility to obesity and its comorbidities among Hispanic children have not been identified. The VIVA LA FAMILIA Study was designed to genetically map childhood obesity and associated biological processes in the Hispanic population. A genome-wide association study (GWAS) entailed genotyping 1.1 million single nucleotide polymorphisms (SNPs) using the Illumina Infinium technology in 815 children. Measured genotype analysis was performed between genetic markers and obesity-related traits i.e., anthropometry, body composition, growth, metabolites, hormones, inflammation, diet, energy expenditure, substrate utilization and physical activity. Identified genome-wide significant loci: 1) corroborated genes implicated in other studies (MTNR1B, ZNF259/APOA5, XPA/FOXE1 (TTF-2), DARC, CCR3, ABO); 2) localized novel genes in plausible biological pathways (PCSK2, ARHGAP11A, CHRNA3); and 3) revealed novel genes with unknown function in obesity pathogenesis (MATK, COL4A1). Salient findings include a nonsynonymous SNP (rs1056513) in INADL (p = 1.2E-07) for weight; an intronic variant in MTNR1B associated with fasting glucose (p = 3.7E-08); variants in the APOA5-ZNF259 region associated with triglycerides (p = 2.5-4.8E-08); an intronic variant in PCSK2 associated with total antioxidants (p = 7.6E-08); a block of 23 SNPs in XPA/FOXE1 (TTF-2) associated with serum TSH (p = 5.5E-08 to 1.0E-09); a nonsynonymous SNP (p = 1.3E-21), an intronic SNP (p = 3.6E-13) in DARC identified for MCP-1; an intronic variant in ARHGAP11A associated with sleep duration (p = 5.0E-08); and, after adjusting for body weight, variants in MATK for total energy expenditure (p = 2.7E-08) and in CHRNA3 for sleeping energy expenditure (p = 6.0E-08). Unprecedented phenotyping and high-density SNP genotyping enabled localization of novel genetic loci associated with the pathophysiology of childhood obesity.

A QTL for genotype by sex interaction for anthropometric measurements in Alaskan Eskimos (GOCADAN Study) on chromosome 19q12-13

Variation in anthropometric measurements due to sexual dimorphism can be the result of genotype by sex interactions (G×S). The purpose of this study was to examine the sex-specific genetic architecture in anthropometric measurements in Alaskan Eskimos from the Genetics of Coronary Artery Disease in Alaska Natives (GOCADAN) study. Maximum likelihood-based variance components decomposition methods, implemented in SOLAR, were used for G×S analyses. Anthropometric measurements included BMI, waist circumference (WC), waist/height ratio, percent body fat (%BF), and subscapular and triceps skinfolds. Except for WC, mean values of all phenotypes were significantly different in men and women (P < 0.05). All anthropometric measures were significantly heritable (P < 0.001). In a preliminary analysis not allowing for G×S interaction, evidence of linkage was detected between markers D19S414 and D19S220 on chromosome 19 for WC (logarithm of odds (lod) = 3.5), %BF (lod = 1.7), BMI (lod = 2.4), waist/height ratio (lod = 2.5), subscapular (lod = 2.1), and triceps skinfolds (lod = 1.9). In subsequent analyses which allowed for G×S interaction, linkage was again found between these traits and the same two markers on chromosome 19 with significantly improved lod scores for: WC (lod = 4.5), %BF (lod = 3.8), BMI (lod = 3.5), waist/height ratio (lod = 3.2), subscapular (lod = 3.0), and triceps skinfolds (lod = 2.9). These results support the evidence of a G×S interaction in the expression of genetic effects resulting in sexual dimorphism in anthropometric phenotypes and identify the chromosome 19q12-13 region as important for adiposity-related traits in Alaskan Eskimos.

A novel locus for arterial hypertension on chromosome 1p36 maps to a metabolic syndrome trait cluster in the Sorbs, a Slavic population isolate in Germany

OBJECTIVE

Genome-wide linkage studies and genome-wide association studies have not as yet identified major genes contributing to primary hypertension in the general population. This state-of-affairs suggests considerable heterogeneity with small contributing effects for primary hypertension, or other complex genetic traits, in outbred populations. Isolated populations, as recent data from Iceland and French Canada suggest, could offer a solution to this problem.

METHODS

We studied a Slavic isolate in Germany, the Sorbs, and genotyped 1040 polymorphic microsatellite markers in 87 multigeneration families.

RESULTS

Our genome-wide linkage scan revealed a locus on chromosome 1p36.13 at D1S3669-D1S2826 (40.95 cM Marshfield coordinates; logarithm of the odds = 3.45, nominal P = 0.00003) that reached genome-wide significance (P = 0.004), indicating the increased power in isolated populations. The chromosome 1 locus maps to a region in which traits such as diabetes, hyperlipidemia, obesity and BMI cluster.

CONCLUSION

Our results suggest that this locus contributes to the metabolic syndrome, and that further attention in this and other populations is warranted.

Agreement among type 2 diabetes linkage studies but a poor correlation with results from genome-wide association studies

AIMS/HYPOTHESIS

Little of the genetic basis for type 2 diabetes has been explained, despite numerous genetic linkage studies and the discovery of multiple genes in genome-wide association (GWA) studies. To begin to resolve the genetic component of this disease, we searched for sites at which genetic results had been corroborated in different studies, in the expectation that replication among studies should direct us to the genomic locations of causative genes with more confidence than the results of individual studies.

METHODS

We have mapped the physical location of results from 83 linkage reports (for type 2 diabetes and diabetes precursor quantitative traits [QTs, e.g. plasma insulin levels]) and recent large GWA reports (for type 2 diabetes) onto the same human genome sequence to identify replicated results in diabetes genetic 'hot spots'.

RESULTS

Genetic linkage has been found at least ten times at 18 different locations, and at least five times in 56 locations. All replication clusters contained study populations from more than one ethnic background and most contained results for both diabetes and QTs. There is no close relationship between the GWA results and linkage clusters, and the nine best replication clusters have no nearby GWA result.

CONCLUSIONS/INTERPRETATION

Many of the genes for type 2 diabetes remain unidentified. This analysis identifies the broad location of yet to be identified genes on 6q, 1q, 18p, 2q, 20q, 17pq, 8p, 19q and 9q. The discrepancy between the linkage and GWA studies may be explained by the presence of multiple, uncommon, mildly deleterious polymorphisms scattered throughout the regulatory and coding regions of genes for type 2 diabetes.

A quantitative trait locus on chromosome 13q affects fasting glucose levels in Hispanic children

OBJECTIVE

The prevalence of childhood obesity has increased dramatically in the United States. Early presentation of type 2 diabetes has been observed in children and adolescents, especially in the Hispanic population. The genetic contribution of glucose homeostasis related to childhood obesity is poorly understood. The objective of this study was to localize quantitative trait loci influencing fasting serum glucose levels in Hispanic children participating in the Viva La Familia Study.

DESIGN

Subjects were 1030 children ascertained through an overweight child from 319 Hispanic families. Fasting serum glucose levels were measured enzymatically, and genetic linkage analyses were conducted using SOLAR software.

RESULTS

Fasting glucose was heritable, with a heritability of 0.62 +/- 0.08 (P < 0.01). Genome-wide scan mapped fasting serum glucose to markers D13S158-D13S173 on chromosome 13q (LOD score of 4.6). A strong positional candidate gene is insulin receptor substrate 2, regulator of glucose homeostasis and a candidate gene for obesity. This region was reported previously to be linked to obesity- and diabetes-related phenotypes.

CONCLUSIONS

A quantitative trait locus on chromosome 13q contributes to the variation in fasting serum glucose levels in Hispanic children at high risk for obesity.