Common Variants of FTO Are Associated with Childhood Obesity in a Cross-Sectional Study of 3,126 Urban Indian Children

Background

FTO variants are robustly associated with obesity and related traits in many population and shown to have variable impact during life course. Although studies have shown association of FTO variants with adiposity in adult Indian, its association in Indian children is yet to be confirmed.

Methods

Here we examined association of FTO variants (rs9939609 and rs8050136) with obesity and related anthropometric and biochemical traits in 3,126 Indian children (aged 11–17 years) including 2,230 normal-weight and 896 over-weight/obese children. We also compared effects observed in the present study with that observed in previous studies on South Asian adults and children of other ethnic groups.

Results

The variant rs9939609 showed significant association with risk of obesity [OR = 1.21, P = 2.5×10⁻³] and its measures BMI, weight, waist circumference and hip circumference [β range = 0.11 to 0.14 Z-score units; P range = 1.3×10⁻⁴ to 1.6×10⁻⁷] in children. The observed effect sizes in Indian children were similar to those reported for European children. Variant rs9939609 explained 0.88% of BMI variance in Indian children. The effect sizes of rs9939609 on BMI and WC were ∼2 fold higher in children than adults. Interestingly rs9939609 was also associated with serum levels of thyroid stimulating hormone (TSH) [β = 0.10 Z-score, P = 5.8×10⁻³]. The other variant rs8050136 was in strong linkage disequilibrium with rs9939609 (r² = 0.97) and provided similar association results.

Conclusion

The study provides first report of association of FTO variants with obesity and related anthropometric traits in Indian children with higher impact in children compared to adults. We also demonstrated association of FTO variant with serum levels of TSH, indicating putative influence of FTO in hypothalamic-pituitary-thyroid axis.     

Targeting 160 Candidate Genes for Blood Pressure Regulation with a Genome-Wide Genotyping Array

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region±10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10⁻³) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15×10⁻⁵), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5×10⁻⁵) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20×10⁻⁵). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10⁻³) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.     

Contribution of Common PCSK1 Genetic Variants to Obesity in 8,359 Subjects from Multi-Ethnic American Population

Common PCSK1 variants (notably rs6232 and rs6235) have been shown to be associated with obesity in European, Asian and Mexican populations. To determine whether common PCSK1 variants contribute to obesity in American population, we conducted association analyses in 8,359 subjects using two multi-ethnic American studies: the Coronary Artery Risk Development in Young Adults (CARDIA) study and the Multi-Ethnic Study of Atherosclerosis (MESA). By evaluating the contribution of rs6232 and rs6235 in each ethnic group, we found that in European-American subjects from CARDIA, only rs6232 was associated with BMI (P = 0.006) and obesity (P = 0.018) but also increased the obesity incidence during the 20 years of follow-up (HR = 1.53 [1.07–2.19], P = 0.019). Alternatively, in African-American subjects from CARDIA, rs6235 was associated with BMI (P = 0.028) and obesity (P = 0.018). Further, by combining the two case-control ethnic groups from the CARDIA study in a meta-analysis, association between rs6235 and obesity risk remained significant (OR = 1.23 [1.05–1.45], P = 9.5×10⁻³). However, neither rs6232 nor rs6235 was associated with BMI or obesity in the MESA study. Interestingly, rs6232 was associated with BMI (P = 4.2×10⁻³) and obesity (P = 3.4×10⁻³) in the younger European-American group combining samples from the both studies [less than median age (53 years)], but not among the older age group (P = 0.756 and P = 0.935 for BMI and obesity, respectively). By combining all the case-control ethnic groups from CARDIA and MESA in a meta-analysis, we found no significant association for the both variants and obesity risk. Finally, by exploring the full PCSK1 locus, we observed that no variant remained significant after correction for multiple testing. These results indicate that common PCSK1 variants (notably rs6232 and rs6235) contribute modestly to obesity in multi-ethnic American population. Further, these results suggest that the association of rs6232 with obesity may be age-dependent in European-Americans. However, multiple replication studies in multi-ethnic American population are needed to confirm our findings.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

Common variants in CRP and LEPR influence high sensitivity C-reactive protein levels in North Indians

Background

High sensitivity C-reactive protein (hsCRP) levels are shown to be influenced by genetic variants in Europeans; however, little is explored in Indian population.

Methods

Herein, we comprehensively evaluated association of all previously reported genetic determinants of hsCRP levels, including 18 cis (proximal to CRP gene) and 73 trans-acting (distal to CRP gene) variants in 4,200 North Indians of Indo-European ethnicity. First, we evaluated association of 91 variants from 12 candidate loci with hsCRP levels in 2,115 North Indians (1,042 non-diabetic subjects and 1,073 patients with type 2 diabetes). Then, cis and trans-acting variants contributing maximally to hsCRP level variation were further replicated in an independent 2,085 North Indians (1,047 patients with type 2 diabetes and 1,038 non-diabetic subjects).

Results

We found association of 12 variants from CRP, LEPR, IL1A, IL6, and IL6R with hsCRP levels in non-diabetic subjects. However, only rs3093059-CRP [β = 0.33, P = 9.6×10⁻⁵] and the haplotype harboring rs3093059 risk allele [β = 0.32 µg/mL, P = 1.4×10⁻⁴/P_perm = 9.0×10⁻⁴] retained significance after correcting for multiple testing. The cis-acting variant rs3093059-CRP had maximum contribution to the variance in hsCRP levels (1.14%). Among, trans-acting variants, rs1892534-LEPR was observed to contribute maximally to hsCRP level variance (0.59%). Associations of rs3093059-CRP and rs1892534-LEPR were confirmed by replication and attained higher significance after meta-analysis [β_meta = 0.26/0.22; P_meta = 4.3×10⁻⁷/7.4×10⁻³ and β_meta = −0.15/−0.12; P_meta = 2.0×10⁻⁶/1.6×10⁻⁶ for rs3093059 and rs1892534, respectively in non-diabetic subjects and all subjects taken together].

Conclusion

In conclusion, we identified rs3093059 in CRP and rs1892534 in LEPR as major cis and trans-acting contributor respectively, to the variance in hsCRP levels in North Indian population.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Ancestry-Shift Refinement Mapping of the C6orf97-ESR1 Breast Cancer Susceptibility Locus

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor α (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case∶control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2×10⁻³), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9×10⁻⁴) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9×10⁻⁷), was without significant heterogeneity between ancestries (P_het = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.     

Rare Genomic Structural Variants in Complex Disease: Lessons from the Replication of Associations with Obesity

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the ‘missing heritability’. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10⁻⁴ (95% confidence interval [9.6×10⁻⁵–3.1×10⁻⁴]); accounts overall for 0.5% [0.19%–0.82%] of severe childhood obesity cases (P = 3.8×10⁻¹⁰; odds ratio = 25.0 [9.9–60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m⁻² [1.8–10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.     

Identification of Novel Genetic Loci Associated with Thyroid Peroxidase Antibodies and Clinical Thyroid Disease

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto''s thyroiditis), as well as autoimmune hyperthyroidism (Graves'' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10⁻⁸) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68–2.81, P = 8.1×10⁻⁸), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26–1.82, P = 2.9×10⁻⁶), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66–0.89, P = 6.5×10⁻⁴). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves'' disease (OR: 1.37, 95% CI 1.22–1.54, P = 1.2×10⁻⁷ and OR: 1.25, 95% CI 1.12–1.39, P = 6.2×10⁻⁵). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18–2.10, P = 1.9×10⁻³). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.     

The Effects of Genetic Variation in FTO rs9939609 on Obesity and Dietary Preferences in Chinese Han Children and Adolescents

The association of the rs9939609 single nucleotide polymorphism in FTO gene with obesity has been extensively investigated in studies of populations of European, African, and Asian ancestry. However, inconsistent results have been reported in Asian populations, and the relationship of FTO variation and dietary behaviors has only rarely been examined in Chinese children and adolescents. The aim of this study was to assess the association of rs9939609 with obesity and dietary preferences in childhood in a Chinese population. Epidemiological data including dietary preferences were collected in interviews using survey questionnaires, and rs9939609 genotype was determined by real-time PCR. The associations of rs9939609 genotypes with obesity and dietary preferences were analyzed by multivariate logistic regression using both additive and dominant models. The results showed that subjects with a TA or AA genotype had an increased risk of obesity compared with the TT participants; the odds ratios (ORs) were 1.47 (95% CI: 1.25–1.71, P = 1.73×10⁻⁶), and 3.32 (95% CI: 2.01–5.47, P = 2.68×10⁻⁶), respectively. After adjusting for age and gender, body mass index, waist circumference, hip circumference, systolic blood pressure, diastolic blood pressure, fasting blood glucose, triglycerides, and low-density lipoprotein cholesterol were higher, and high-density lipoprotein cholesterol was lower in TA and AA participants than in those with the TT genotype. After additionally controlling for body mass index, the association remained significant only for systolic blood pressure (P = 0.005). Compared with TT participants, those with the AA genotype were more likely to prefer a meat-based diet (OR = 2.81, 95% CI: 1.52–5.21). The combined OR for obesity in participants with TA/AA genotypes and preference for a meat-based diet was 4.04 (95% CI: 2.8–5.81) compared with the TT participants who preferred a plant-based diet. These findings indicate the genetic variation of rs9939609 is associated with obesity and dietary preferences in Chinese children and adolescents.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

A Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

The role of inflammatory pathway genetic variation on maternal metabolic phenotypes during pregnancy

Background

Since mediators of inflammation are associated with insulin resistance, and the risk of developing diabetes mellitus and gestational diabetes, we hypothesized that genetic variation in members of the inflammatory gene pathway impact glucose levels and related phenotypes in pregnancy. We evaluated this hypothesis by testing for association between genetic variants in 31 inflammatory pathway genes in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) cohort, a large multiethnic multicenter study designed to address the impact of glycemia less than overt diabetes on pregnancy outcome.

Results

Fasting, 1-hour, and 2-hour glucose, fasting and 1-hour C-peptide, and HbA1c levels were measured in blood samples obtained from HAPO participants during an oral glucose tolerance test at 24-32 weeks gestation. We tested for association between 458 SNPs mapping to 31 genes in the inflammatory pathway and metabolic phenotypes in 3836 European ancestry and 1713 Thai pregnant women. The strongest evidence for association was observed with TNF alpha and HbA1c (rs1052248; 0.04% increase per allele C; p-value = 4.4×10⁻⁵), RETN and fasting plasma glucose (rs1423096; 0.7 mg/dl decrease per allele A; p-value = 1.1×10⁻⁴), IL8 and 1 hr plasma glucose (rs2886920; 2.6 mg/dl decrease per allele T; p-value = 1.3×10⁻⁴), ADIPOR2 and fasting C-peptide (rs2041139; 0.55 ug/L decrease per allele A; p-value = 1.4×10⁻⁴), LEPR and 1-hour C-peptide (rs1171278; 0.62 ug/L decrease per allele T; p-value = 2.4×10⁻⁴), and IL6 and 1-hour plasma glucose (rs6954897; −2.29 mg/dl decrease per allele G, p-value = 4.3×10⁻⁴).

Conclusions

Based on the genes surveyed in this study the inflammatory pathway is unlikely to have a strong impact on maternal metabolic phenotypes in pregnancy although variation in individual members of the pathway (e.g. RETN, IL8, ADIPOR2, LEPR, IL6, and TNF alpha,) may contribute to metabolic phenotypes in pregnant women.     

Gene × Physical Activity Interactions in Obesity: Combined Analysis of 111,421 Individuals of European Ancestry

Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication studies are lacking. Therefore, we tested whether the aggregate effect of these loci is diminished in adults of European ancestry reporting high levels of physical activity. Twelve obesity-susceptibility loci were genotyped or imputed in 111,421 participants. A genetic risk score (GRS) was calculated by summing the BMI-associated alleles of each genetic variant. Physical activity was assessed using self-administered questionnaires. Multiplicative interactions between the GRS and physical activity on BMI were tested in linear and logistic regression models in each cohort, with adjustment for age, age², sex, study center (for multicenter studies), and the marginal terms for physical activity and the GRS. These results were combined using meta-analysis weighted by cohort sample size. The meta-analysis yielded a statistically significant GRS × physical activity interaction effect estimate (P_interaction = 0.015). However, a statistically significant interaction effect was only apparent in North American cohorts (n = 39,810, P_interaction = 0.014 vs. n = 71,611, P_interaction = 0.275 for Europeans). In secondary analyses, both the FTO rs1121980 (P_interaction = 0.003) and the SEC16B rs10913469 (P_interaction = 0.025) variants showed evidence of SNP × physical activity interactions. This meta-analysis of 111,421 individuals provides further support for an interaction between physical activity and a GRS in obesity disposition, although these findings hinge on the inclusion of cohorts from North America, indicating that these results are either population-specific or non-causal.     

Distribution and Medical Impact of Loss-of-Function Variants in the Finnish Founder Population

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5–5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻⁸) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹⁷). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻⁴), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.