Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites.

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Directional dominance on stature and cognition in diverse human populations.

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

GWAS identifies population-specific new regulatory variants in FUT6 associated with plasma B12 concentrations in Indians.

Vitamin B12 is an important cofactor in one-carbon metabolism whose dysregulation is associated with various clinical conditions. Indians have a high prevalence of B12 deficiency but little is known about the genetic determinants of circulating B12 concentrations in Indians. We performed a genome-wide association study in 1001 healthy participants in the Pune Maternal Nutrition Study (PMNS), replication studies in 3418 individuals from other Indian cohorts and by meta-analysis identified new variants, rs3760775 (P = 1.2 × 10-23) and rs78060698 (P = 8.3 × 10-17) in FUT6 to be associated with circulating B12 concentrations. Although in-silico analysis replicated both variants in Europeans, differences in the effect allele frequency, effect size and the linkage disequilibrium structure of credible set variants with the reported variants suggest population-specific characteristics in this region. We replicated previously reported variants rs602662, rs601338 in FUT2, rs3760776, rs708686 in FUT6, rs34324219 in TCN1 (all P < 5 × 10-8), rs1131603 in TCN2 (P = 3.4 × 10-5), rs12780845 in CUBN (P = 3.0 × 10-3) and rs2270655 in MMAA (P = 2.0 × 10-3). Circulating B12 concentrations in the PMNS and Parthenon study showed a significant decline with increasing age (P < 0.001), however, the genetic contribution to B12 concentrations remained constant. Luciferase reporter and electrophoretic-mobility shift assay for the FUT6 variant rs78060698 using HepG2 cell line demonstrated strong allele-specific promoter and enhancer activity and differential binding of HNF4α, a key regulator of expression of various fucosyltransferases. Hence, the rs78060698 variant, through regulation of fucosylation may control intestinal host-microbial interaction which could influence B12 concentrations. Our results suggest that in addition to established genetic variants, population-specific variants are important in determining plasma B12 concentrations.

DNA methylation at the suppressor of cytokine signaling 3 (SOCS3) gene influences height in childhood.

Human height is strongly influenced by genetics but the contribution of modifiable epigenetic factors is under-explored, particularly in low and middle-income countries (LMIC). We investigate links between blood DNA methylation and child height in four LMIC cohorts (n = 1927) and identify a robust association at three CpGs in the suppressor of cytokine signaling 3 (SOCS3) gene which replicates in a high-income country cohort (n = 879). SOCS3 methylation (SOCS3m)-height associations are independent of genetic effects. Mendelian randomization analysis confirms a causal effect of SOCS3m on height. In longitudinal analysis, SOCS3m explains a maximum 9.5% of height variance in mid-childhood while the variance explained by height polygenic risk score increases from birth to 21 years. Children's SOCS3m is associated with prenatal maternal folate and socio-economic status. In-vitro characterization confirms a regulatory effect of SOCS3m on gene expression. Our findings suggest epigenetic modifications may play an important role in driving child height in LMIC.

Babies of South Asian and European Ancestry Show Similar Associations With Genetic Risk Score for Birth Weight Despite the Smaller Size of South Asian Newborns.

Size at birth is known to be influenced by various fetal and maternal factors, including genetic effects. South Asians have a high burden of low birth weight and cardiometabolic diseases, yet studies of common genetic variations underpinning these phenotypes are lacking. We generated independent, weighted fetal genetic scores (fGSs) and maternal genetic scores (mGSs) from 196 birth weight-associated variants identified in Europeans and conducted an association analysis with various fetal birth parameters and anthropometric and cardiometabolic traits measured at different follow-up stages (5-6-year intervals) from seven Indian and Bangladeshi cohorts of South Asian ancestry. The results from these cohorts were compared with South Asians in UK Biobank and the Exeter Family Study of Childhood Health, a European ancestry cohort. Birth weight increased by 50.7 g and 33.6 g per SD of fGS (P = 9.1 × 10-11) and mGS (P = 0.003), respectively, in South Asians. A relatively weaker mGS effect compared with Europeans indicates possible different intrauterine exposures between Europeans and South Asians. Birth weight was strongly associated with body size in both childhood and adolescence (P = 3 × 10-5 to 1.9 × 10-51); however, fGS was associated with body size in childhood only (P < 0.01) and with head circumference, fasting glucose, and triglycerides in adults (P < 0.01). The substantially smaller newborn size in South Asians with comparable fetal genetic effect to Europeans on birth weight suggests a significant role of factors related to fetal growth that were not captured by the present genetic scores. These factors may include different environmental exposures, maternal body size, health and nutritional status, etc. Persistent influence of genetic loci on size at birth and adult metabolic syndrome in our study supports a common genetic mechanism that partly explains associations between early development and later cardiometabolic health in various populations, despite marked differences in phenotypic and environmental factors in South Asians.

Intrauterine Programming of Diabetes and Adiposity.

The prevalence of diabetes and adiposity has increased at an alarming rate and together they contribute to the rise in morbidity and mortality worldwide. Genetic studies till date have succeeded in explaining only a proportion of heritability, while a major component remains unexplained. Early life determinants of future risk of these diseases are likely contributors to the missing heritability and thus have a significant potential in disease prevention. Epidemiological and animal studies show the importance of intrauterine and early postnatal environment in programming of the fetus to adverse metabolic outcomes and support the notion of Developmental Origins of Health and Disease (DOHaD). Emerging evidence highlights the role of epigenetic mechanisms in mediating effects of environmental exposures, which in certain instances may exhibit intergenerational transmission even in the absence of exposure. In this article, we will discuss the complexity of diabetes and increased adiposity and mechanisms of programming of these adverse metabolic conditions.

Association of common genetic variants with lipid traits in the Indian population.

Genome-wide association studies (GWAS) have been instrumental in identifying novel genetic variants associated with altered plasma lipid levels. However, these quantitative trait loci have not been tested in the Indian population, where there is a poorly understood and growing burden of cardiometabolic disorders. We present the association of six single nucleotide polymorphisms in 1671 sib pairs (3342 subjects) with four lipid traits: total cholesterol, triglycerides, high density lipoprotein cholesterol (HDL-C) and low density lipoprotein cholesterol (LDL-C). We also investigated the interaction effects of gender, location, fat intake and physical activity. Each copy of the risk allele of rs964184 at APOA1 was associated with 1.06 mmol/l increase in triglycerides (SE = 0.049; p = 0.006), rs3764261 at CETP with 1.02 mmol/l increase in both total cholesterol (SE = 0.042; p = 0.017) and HDL-C (SE = 0.041; p = 0.008), rs646776 at CELSR2-PSRC1-SORT1 with 0.96 mmol/l decrease in cholesterol (SE = 0.043; p = 0.0003) and 0.15 mmol/l decrease in LDL-C levels (SE = 0.043; p = 0.0003) and rs2954029 at TRIB1 with 1.02 mmol/l increase in HDL-C (SE = 0.039; p = 0.047). A combined risk score of APOA1 and CETP loci predicted an increase of 1.25 mmol/l in HDL-C level (SE = 0.312; p = 0.0007). Urban location and sex had strong interaction effects on the genetic association of most of the studied loci with lipid traits. To conclude, we validated four genetic variants (identified by GWAS in western populations) associated with lipid traits in the Indian population. The interaction effects found here may explain the sex-specific differences in lipid levels and their heritability. Urbanization appears to influence the nature of the association with GWAS lipid loci in this population. However, these findings will require replication in other Indian populations.