Amanda J. Bennett

A genome-wide association study identifies protein quantitative trait loci (pQTLs)

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts – cis effects, and elsewhere in the genome – trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8×10−57), CCL4L1 (p = 3.9×10−21), IL18 (p = 6.8×10−13), LPA (p = 4.4×10−10), GGT1 (p = 1.5×10−7), SHBG (p = 3.1×10−7), CRP (p = 6.4×10−6) and IL1RN (p = 7.3×10−6) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8×10−40), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways.

Meeting the demand for crop production : the challenge of yield decline in crops grown in short rotations

There is a trend world‐wide to grow crops in short rotation or in monoculture, particularly in conventional agriculture. This practice is becoming more prevalent due to a range of factors including economic market trends, technological advances, government incentives, and retailer and consumer demands. Land‐use intensity will have to increase further in future in order to meet the demands of growing crops for both bioenergy and food production, and long rotations may not be considered viable or practical. However, evidence indicates that crops grown in short rotations or monoculture often suffer from yield decline compared to those grown in longer rotations or for the first time. Numerous factors have been hypothesised as contributing to yield decline, including biotic factors such as plant pathogens, deleterious rhizosphere microorganisms, mycorrhizas acting as pathogens, and allelopathy or autotoxicity of the crop, as well as abiotic factors such as land management practices and nutrient availability. In many cases, soil microorganisms have been implicated either directly or indirectly in yield decline. Although individual factors may be responsible for yield decline in some cases, it is more likely that combinations of factors interact to cause the problem. However, evidence confirming the precise role of these various factors is often lacking in field studies due to the complex nature of cropping systems and the numerous interactions that take place within them. Despite long‐term knowledge of the yield‐decline phenomenon, there are few tools to counteract it apart from reverting to longer crop rotations or break crops. Alternative cropping and management practices such as double‐cropping or inter‐cropping, tillage and organic amendments may prove valuable for combating some of the negative effects seen when crops are grown in short rotation. Plant breeding continues to be important, although this does require a specific breeding target to be identified. This review identifies gaps in our understanding of yield decline, particularly with respect to the complex interactions occurring between the different components of agro‐ecosystems, which may well influence food security in the 21st Century.

Association between common variation at the FTO locus and changes in body mass index from infancy to late childhood: the complex nature of genetic association through growth and development.

An age-dependent association between variation at the FTO locus and BMI in children has been suggested. We meta-analyzed associations between the FTO locus (rs9939609) and BMI in samples, aged from early infancy to 13 years, from 8 cohorts of European ancestry. We found a positive association between additional minor (A) alleles and BMI from 5.5 years onwards, but an inverse association below age 2.5 years. Modelling median BMI curves for each genotype using the LMS method, we found that carriers of minor alleles showed lower BMI in infancy, earlier adiposity rebound (AR), and higher BMI later in childhood. Differences by allele were consistent with two independent processes: earlier AR equivalent to accelerating developmental age by 2.37% (95% CI 1.87, 2.87, p = 10−20) per A allele and a positive age by genotype interaction such that BMI increased faster with age (p = 10−23). We also fitted a linear mixed effects model to relate genotype to the BMI curve inflection points adiposity peak (AP) in infancy and AR. Carriage of two minor alleles at rs9939609 was associated with lower BMI at AP (−0.40% (95% CI: −0.74, −0.06), p = 0.02), higher BMI at AR (0.93% (95% CI: 0.22, 1.64), p = 0.01), and earlier AR (−4.72% (−5.81, −3.63), p = 10−17), supporting cross-sectional results. Overall, we confirm the expected association between variation at rs9939609 and BMI in childhood, but only after an inverse association between the same variant and BMI in infancy. Patterns are consistent with a shift on the developmental scale, which is reflected in association with the timing of AR rather than just a global increase in BMI. Results provide important information about longitudinal gene effects and about the role of FTO in adiposity. The associated shifts in developmental timing have clinical importance with respect to known relationships between AR and both later-life BMI and metabolic disease risk.

Type 2 diabetes risk alleles are associated with reduced size at birth.

OBJECTIVE Low birth weight is associated with an increased risk of type 2 diabetes. The mechanisms underlying this association are unknown and may represent intrauterine programming or two phenotypes of one genotype. The fetal insulin hypothesis proposes that common genetic variants that reduce insulin secretion or action may predispose to type 2 diabetes and also reduce birth weight, since insulin is a key fetal growth factor. We tested whether common genetic variants that predispose to type 2 diabetes also reduce birth weight. RESEARCH DESIGN AND METHODS We genotyped single-nucleotide polymorphisms (SNPs) at five recently identified type 2 diabetes loci (CDKAL1, CDKN2A/B, HHEX-IDE, IGF2BP2, and SLC30A8) in 7,986 mothers and 19,200 offspring from four studies of white Europeans. We tested the association between maternal or fetal genotype at each locus and birth weight of the offspring. RESULTS We found that type 2 diabetes risk alleles at the CDKAL1 and HHEX-IDE loci were associated with reduced birth weight when inherited by the fetus (21 g [95% CI 11–31], P = 2 × 10−5, and 14 g [4–23], P = 0.004, lower birth weight per risk allele, respectively). The 4% of offspring carrying four risk alleles at these two loci were 80 g (95% CI 39–120) lighter at birth than the 8% carrying none (Ptrend = 5 × 10−7). There were no associations between birth weight and fetal genotypes at the three other loci or maternal genotypes at any locus. CONCLUSIONS Our results are in keeping with the fetal insulin hypothesis and provide robust evidence that common disease-associated variants can alter size at birth directly through the fetal genotype.

Genetic Determinants of Height Growth Assessed Longitudinally from Infancy to Adulthood in the Northern Finland Birth Cohort 1966

Recent genome-wide association (GWA) studies have identified dozens of common variants associated with adult height. However, it is unknown how these variants influence height growth during childhood. We derived peak height velocity in infancy (PHV1) and puberty (PHV2) and timing of pubertal height growth spurt from parametric growth curves fitted to longitudinal height growth data to test their association with known height variants. The study consisted of N = 3,538 singletons from the prospective Northern Finland Birth Cohort 1966 with genotype data and frequent height measurements (on average 20 measurements per person) from 0–20 years. Twenty-six of the 48 variants tested associated with adult height (p<0.05, adjusted for sex and principal components) in this sample, all in the same direction as in previous GWA scans. Seven SNPs in or near the genes HHIP, DLEU7, UQCC, SF3B4/SV2A, LCORL, and HIST1H1D associated with PHV1 and five SNPs in or near SOCS2, SF3B4/SV2A, C17orf67, CABLES1, and DOT1L with PHV2 (p<0.05). We formally tested variants for interaction with age (infancy versus puberty) and found biologically meaningful evidence for an age-dependent effect for the SNP in SOCS2 (p = 0.0030) and for the SNP in HHIP (p = 0.045). We did not have similar prior evidence for the association between height variants and timing of pubertal height growth spurt as we had for PHVs, and none of the associations were statistically significant after correction for multiple testing. The fact that in this sample, less than half of the variants associated with adult height had a measurable effect on PHV1 or PHV2 is likely to reflect limited power to detect these associations in this dataset. Our study is the first genetic association analysis on longitudinal height growth in a prospective cohort from birth to adulthood and gives grounding for future research on the genetic regulation of human height during different periods of growth.

Type 2 Diabetes TCF7L2 Risk Genotypes Alter Birth Weight: A Study of 24,053 Individuals

The role of genes in normal birth-weight variation is poorly understood, and it has been suggested that the genetic component of fetal growth is small. Type 2 diabetes genes may influence birth weight through maternal genotype, by increasing maternal glycemia in pregnancy, or through fetal genotype, by altering fetal insulin secretion. We aimed to assess the role of the recently described type 2 diabetes gene TCF7L2 in birth weight. We genotyped the polymorphism rs7903146 in 15,709 individuals whose birth weight was available from six studies and in 8,344 mothers from three studies. Each fetal copy of the predisposing allele was associated with an 18-g (95% confidence interval [CI] 7-29 g) increase in birth weight (P=.001) and each maternal copy with a 30-g (95% CI 15-45 g) increase in offspring birth weight (P=2.8x10-5). Stratification by fetal genotype suggested that the association was driven by maternal genotype (31-g [95% CI 9-48 g] increase per allele; corrected P=.003). Analysis of diabetes-related traits in 10,314 nondiabetic individuals suggested the most likely mechanism is that the risk allele reduces maternal insulin secretion (disposition index reduced by ~0.15 standard deviation; P=1x10-4), which results in increased maternal glycemia in pregnancy and hence increased offspring birth weight. We combined information with the other common variant known to alter fetal growth, the -30G-->A polymorphism of glucokinase (rs1799884). The 4% of offspring born to mothers carrying three or four risk alleles were 119 g (95% CI 62-172 g) heavier than were the 32% born to mothers with none (for overall trend, P=2x10-7), comparable to the impact of maternal smoking during pregnancy. In conclusion, we have identified the first type 2 diabetes-susceptibility allele to be reproducibly associated with birth weight. Common gene variants can substantially influence normal birth-weight variation.

Total Zinc Intake May Modify the Glucose-Raising Effect of a Zinc Transporter (SLC30A8) Variant: A 14-Cohort Meta-analysis

OBJECTIVE Many genetic variants have been associated with glucose homeostasis and type 2 diabetes in genome-wide association studies. Zinc is an essential micronutrient that is important for β-cell function and glucose homeostasis. We tested the hypothesis that zinc intake could influence the glucose-raising effect of specific variants. RESEARCH DESIGN AND METHODS We conducted a 14-cohort meta-analysis to assess the interaction of 20 genetic variants known to be related to glycemic traits and zinc metabolism with dietary zinc intake (food sources) and a 5-cohort meta-analysis to assess the interaction with total zinc intake (food sources and supplements) on fasting glucose levels among individuals of European ancestry without diabetes. RESULTS We observed a significant association of total zinc intake with lower fasting glucose levels (β-coefficient ± SE per 1 mg/day of zinc intake: −0.0012 ± 0.0003 mmol/L, summary P value = 0.0003), while the association of dietary zinc intake was not significant. We identified a nominally significant interaction between total zinc intake and the SLC30A8 rs11558471 variant on fasting glucose levels (β-coefficient ± SE per A allele for 1 mg/day of greater total zinc intake: −0.0017 ± 0.0006 mmol/L, summary interaction P value = 0.005); this result suggests a stronger inverse association between total zinc intake and fasting glucose in individuals carrying the glucose-raising A allele compared with individuals who do not carry it. None of the other interaction tests were statistically significant. CONCLUSIONS Our results suggest that higher total zinc intake may attenuate the glucose-raising effect of the rs11558471 SLC30A8 (zinc transporter) variant. Our findings also support evidence for the association of higher total zinc intake with lower fasting glucose levels.

Transcript Expression Data from Human Islets Links Regulatory Signals from Genome-Wide Association Studies for Type 2 Diabetes and Glycemic Traits to Their Downstream Effectors.

The intersection of genome-wide association analyses with physiological and functional data indicates that variants regulating islet gene transcription influence type 2 diabetes (T2D) predisposition and glucose homeostasis. However, the specific genes through which these regulatory variants act remain poorly characterized. We generated expression quantitative trait locus (eQTL) data in 118 human islet samples using RNA-sequencing and high-density genotyping. We identified fourteen loci at which cis-exon-eQTL signals overlapped active islet chromatin signatures and were coincident with established T2D and/or glycemic trait associations. ‎At some, these data provide an experimental link between GWAS signals and biological candidates, such as DGKB and ADCY5. At others, the cis-signals implicate genes with no prior connection to islet biology, including WARS and ZMIZ1. At the ZMIZ1 locus, we show that perturbation of ZMIZ1 expression in human islets and beta-cells influences exocytosis and insulin secretion, highlighting a novel role for ZMIZ1 in the maintenance of glucose homeostasis. Together, these findings provide a significant advance in the mechanistic insights of T2D and glycemic trait association loci.

Mutations in HNF1A Result in Marked Alterations of Plasma Glycan Profile

A recent genome-wide association study identified hepatocyte nuclear factor 1-α (HNF1A) as a key regulator of fucosylation. We hypothesized that loss-of-function HNF1A mutations causal for maturity-onset diabetes of the young (MODY) would display altered fucosylation of N-linked glycans on plasma proteins and that glycan biomarkers could improve the efficiency of a diagnosis of HNF1A-MODY. In a pilot comparison of 33 subjects with HNF1A-MODY and 41 subjects with type 2 diabetes, 15 of 29 glycan measurements differed between the two groups. The DG9-glycan index, which is the ratio of fucosylated to nonfucosylated triantennary glycans, provided optimum discrimination in the pilot study and was examined further among additional subjects with HNF1A-MODY (n = 188), glucokinase (GCK)-MODY (n = 118), hepatocyte nuclear factor 4-α (HNF4A)-MODY (n = 40), type 1 diabetes (n = 98), type 2 diabetes (n = 167), and nondiabetic controls (n = 98). The DG9-glycan index was markedly lower in HNF1A-MODY than in controls or other diabetes subtypes, offered good discrimination between HNF1A-MODY and both type 1 and type 2 diabetes (C statistic ≥0.90), and enabled us to detect three previously undetected HNF1A mutations in patients with diabetes. In conclusion, glycan profiles are altered substantially in HNF1A-MODY, and the DG9-glycan index has potential clinical value as a diagnostic biomarker of HNF1A dysfunction.

Variation in the glucose transporter gene SLC2A2 is associated with glycemic response to metformin

Metformin is the first-line antidiabetic drug with over 100 million users worldwide, yet its mechanism of action remains unclear. Here the Metformin Genetics (MetGen) Consortium reports a three-stage genome-wide association study (GWAS), consisting of 13,123 participants of different ancestries. The C allele of rs8192675 in the intron of SLC2A2, which encodes the facilitated glucose transporter GLUT2, was associated with a 0.17% (P = 6.6 × 10−14) greater metformin-induced reduction in hemoglobin A1c (HbA1c) in 10,577 participants of European ancestry. rs8192675 was the top cis expression quantitative trait locus (cis-eQTL) for SLC2A2 in 1,226 human liver samples, suggesting a key role for hepatic GLUT2 in regulation of metformin action. Among obese individuals, C-allele homozygotes at rs8192675 had a 0.33% (3.6 mmol/mol) greater absolute HbA1c reduction than T-allele homozygotes. This was about half the effect seen with the addition of a DPP-4 inhibitor, and equated to a dose difference of 550 mg of metformin, suggesting rs8192675 as a potential biomarker for stratified medicine.