A genome‐wide association study for reading and language abilities in two population cohorts

Candidate genes have been identified for both reading and language, but most of the heritable variance in these traits remains unexplained. Here, we report a genome‐wide association meta‐analysis of two large cohorts: population samples of Australian twins and siblings aged 12–25 years (n = 1177 from 538 families), and a younger cohort of children of the UK Avon Longitudinal Study of Parents and their Children (aged 8 and 9 years; maximum n = 5472). Suggestive association was indicated for reading measures and non‐word repetition (NWR), with the greatest support found for single nucleotide polymorphisms (SNPs) in the pseudogene, ABCC13 (P = 7.34 × 10^?8), and the gene, DAZAP1 (P = 1.32 × 10^?6). Gene‐based analyses showed significant association (P < 2.8 × 10^?6) for reading and spelling with genes CD2L1, CDC2L2 and RCAN3 in two loci on chromosome 1. Some support was found for the same SNPs having effects on both reading skill and NWR, which is compatible with behavior genetic evidence for influences of reading acquisition on phonological‐task performance. The results implicate novel candidates for study in additional cohorts for reading and language abilities.     

Genome‐wide association analyses of child genotype effects and parent‐of‐origin effects in specific language impairment

Specific language impairment (SLI) is a neurodevelopmental disorder that affects linguistic abilities when development is otherwise normal. We report the results of a genome‐wide association study of SLI which included parent‐of‐origin effects and child genotype effects and used 278 families of language‐impaired children. The child genotype effects analysis did not identify significant associations. We found genome‐wide significant paternal parent‐of‐origin effects on chromosome 14q12 (P = 3.74 × 10^?8) and suggestive maternal parent‐of‐origin effects on chromosome 5p13 (P = 1.16 × 10^?7). A subsequent targeted association of six single‐nucleotide‐polymorphisms (SNPs) on chromosome 5 in 313 language‐impaired individuals and their mothers from the ALSPAC cohort replicated the maternal effects, albeit in the opposite direction (P = 0.001); as fathers' genotypes were not available in the ALSPAC study, the replication analysis did not include paternal parent‐of‐origin effects. The paternally‐associated SNP on chromosome 14 yields a non‐synonymous coding change within the NOP9 gene. This gene encodes an RNA‐binding protein that has been reported to be significantly dysregulated in individuals with schizophrenia. The region of maternal association on chromosome 5 falls between the PTGER4 and DAB2 genes, in a region previously implicated in autism and ADHD. The top SNP in this association locus is a potential expression QTL of ARHGEF19 (also called WGEF) on chromosome 1. Members of this protein family have been implicated in intellectual disability. In summary, this study implicates parent‐of‐origin effects in language impairment, and adds an interesting new dimension to the emerging picture of shared genetic etiology across various neurodevelopmental disorders .     

Obesity‐susceptibility loci and the tails of the pediatric BMI distribution

Objective: To determine whether previously identified adult obesity susceptibility loci were associated uniformly with childhood BMI across the BMI distribution. Design and Methods: Children were recruited through the Children's Hospital of Philadelphia (n = 7,225). Associations between the following loci and BMI were assessed using quantile regression: FTO (rs3751812), MC4R (rs12970134), TMEM18 (rs2867125), BDNF (rs6265), TNNI3K (rs1514175), NRXN3 (rs10146997), SEC16B (rs10913469), and GNPDA2 (rs13130484). BMI z‐score (age and gender adjusted) was modeled as the dependent variable, and genotype risk score (sum of risk alleles carried at the 8 loci) was modeled as the independent variable. Results: Each additional increase in genotype risk score was associated with an increase in BMI z‐score at the 5th, 15th, 25th, 50th, 75th, 85th, and 95th BMI z‐score percentiles by 0.04 (±0.02, P = 0.08), 0.07 (±0.01, P = 9.58 × 10^?7), 0.07 (±0.01, P = 1.10 × 10^?8), 0.09 (±0.01, P = 3.13 × 10^?22), 0.11 (±0.01, P = 1.35 × 10^?25), 0.11 (±0.01, P = 1.98 × 10^?20), and 0.06 (±0.01, P = 2.44 × 10^?6), respectively. Each additional increase in genotype risk score was associated with an increase in mean BMI z‐score by 0.08 (±0.01, P = 4.27 × 10^?20). Conclusion: Obesity risk alleles were more strongly associated with increases in BMI z‐score at the upper tail compared to the lower tail of the distribution.     

Fetal lung and placental methylation is associated with in utero nicotine exposure

In utero smoke exposure has been shown to have detrimental effects on lung function and to be associated with persistent wheezing and asthma in children. One potential mechanism of IUS effects could be alterations in DNA methylation, which may have life-long implications. The goal of this study was to examine the association between DNA methylation and nicotine exposure in fetal lung and placental tissue in early development; nicotine exposure in this analysis represents a likely surrogate for in-utero smoke. We performed an epigenome-wide analysis of DNA methylation in fetal lung tissue (n = 85, 41 smoke exposed (48%), 44 controls) and the corresponding placental tissue samples (n = 80, 39 smoke exposed (49%), 41 controls) using the Illumina HumanMethylation450 BeadChip array. Differential methylation analyses were conducted to evaluate the variation associated with nicotine exposure. The most significant CpG sites in the fetal lung analysis mapped to the PKP3 (P = 2.94 × 10^?03), ANKRD33B (P = 3.12 × 10^?03), CNTD2 (P = 4.9 × 10^?03) and DPP10 (P = 5.43 × 10^?03) genes. In the placental methylome, the most significant CpG sites mapped to the GTF2H2C and GTF2H2D genes (P = 2.87 × 10^?06 ? 3.48 × 10^?05). One hundred and one unique CpG sites with P-values < 0.05 were concordant between lung and placental tissue analyses. Gene Set Enrichment Analysis demonstrated enrichment of specific disorders, such as asthma and immune disorders. Our findings demonstrate an association between in utero nicotine exposure and variable DNA methylation in fetal lung and placental tissues, suggesting a role for DNA methylation variation in the fetal origins of chronic diseases.     

The DYX2 locus and neurochemical signaling genes contribute to speech sound disorder and related neurocognitive domains

A major milestone of child development is the acquisition and use of speech and language. Communication disorders, including speech sound disorder (SSD), can impair a child's academic, social and behavioral development. Speech sound disorder is a complex, polygenic trait with a substantial genetic component. However, specific genes that contribute to SSD remain largely unknown. To identify associated genes, we assessed the association of the DYX2 dyslexia risk locus and markers in neurochemical signaling genes (e.g., nicotinic and dopaminergic) with SSD and related endophenotypes. We first performed separate primary associations in two independent samples – Cleveland SSD (210 affected and 257 unaffected individuals in 127 families) and Denver SSD (113 affected individuals and 106 unaffected individuals in 85 families) – and then combined results by meta‐analysis. DYX2 markers, specifically those in the 3′ untranslated region of DCDC2 (P = 1.43 × 10^?4), showed the strongest associations with phonological awareness. We also observed suggestive associations of dopaminergic‐related genes ANKK1 (P = 1.02 × 10^?2) and DRD2 (P = 9.22 × 10^?3) and nicotinic‐related genes CHRNA3 (P = 2.51 × 10^?3) and BDNF (P = 8.14 × 10^?3) with case–control status and articulation. Our results further implicate variation in putative regulatory regions in the DYX2 locus, particularly in DCDC2, influencing language and cognitive traits. The results also support previous studies implicating variation in dopaminergic and nicotinic neural signaling influencing human communication and cognitive development. Our findings expand the literature showing genetic factors (e.g., DYX2) contributing to multiple related, yet distinct neurocognitive domains (e.g., dyslexia, language impairment, and SSD). How these factors interactively yield different neurocognitive and language‐related outcomes remains to be elucidated.     

Effects and mechanisms of ghrelin on cardiac microvascular endothelial cells in rats

Ghrelin is thought to directly exert a protective effect on the cardiovascular system, specifically by promoting vascular endothelial cell function. Our study demonstrates the ability of ghrelin to promote rat CMEC (cardiac microvascular endothelial cell) proliferation, migration and NO (nitric oxide) secretion. CMECs were isolated from left ventricle of adult male Sprague—Dawley rat by enzyme digestion and maintained in endothelial cell medium. Dil‐ac‐LDL (1,1′‐dioctadecyl‐3,3,3′,3′‐ tetramethylindocarbocyanine‐labelled acetylated low‐density lipoprotein) intake assays were used to identify CMECs. Cells were split into five groups and treated with varying concentrations of ghrelin as follows: one control non‐treated group; three ghrelin dosage groups (1×10^?9, 1×10^?8, 1×10^?7 mol/l) and one ghrelin+PI3K inhibitor group (1×10^?7 mol/l ghrelin+20 μmol/l LY294002). After 24 h treatment, cell proliferation capability was measured by MTT [3‐(4,5‐dimethylthiazol‐2‐yl)‐2,5‐diphenyl‐2H‐tetrazolium bromide] assay and Western blot for PCNA (proliferating cell nuclear antigen) protein expression. Migration of CMECs was detected by transwell assays, and NO secretion of CMECs was measured via nitrate reduction. Protein expression of AKT and phosphorylated AKT in CMECs was measured by Western blot after exposure to various concentrations of ghrelin and the PI3K inhibitor LY294002. Our results indicate that ghrelin significantly enhanced cell growth at concentrations of 10^?8 mol/l (0.271±0.041 compared with 0.199±0.021, P=0.03) and 10^?7 mol/l (0.296±0.039 compared with 0.199±0.021, P<0.01). However, addition of the PI3K/AKT inhibitor LY294002 inhibited the ghrelin‐mediated enhancement in cell proliferation (0.227±0.042 compared with 0.199±0.021, P=0.15). At a concentration between 10^?8 and 10^?7 mol/l, ghrelin caused a significant increase in the number of migrated cells compared with the control group (126±9 compared with 98±7, P=0.02; 142±6 compared with 98±7, P<0.01), whereas no such change could be observed in the presence of 20 μmol/l of the PI3K/Akt inhibitor LY294002 (103±7 compared with 98±7, P=0.32). Ghrelin treatment significantly enhanced NO production in a dose‐dependent fashion compared with the untreated control group [(39.93±2.12) μmol/l compared with (30.27±2.71) μmol/l, P=0.02; (56.80±1.98) μmol/l compared with (30.27±2.71) μmol/l, P<0.01]. However, pretreatment with 20 μmol/l LY294002 inhibited the ghrelin‐stimulated increase in NO secretion [(28.97±1.64) μmol/l compared with (30.27±2.71) μmol/l, P=0.37]. In summary, we have found that ghrelin treatment promotes the proliferation, migration and NO secretion of CMECs through activation of PI3K/AKT signalling pathway.     

ALDH2 gene G487A polymorphism and coronary artery disease: a meta‐analysis including 5644 participants

Several studies indicate the mitochondrial Aldehyde Dehydrogenase‐2 (ALDH2) gene G487A polymorphism may be correlated with coronary artery disease (CAD) susceptibility, but a clear consensus has yet to be reached. To elucidate the relationship between the ALDH2 gene G487A polymorphism and CAD within the Chinese population, a meta‐analysis of 5644 subjects from nine individual studies was performed. Pooled odds ratios (ORs) and their corresponding 95% confidence intervals were assessed using random or fixed‐effect models depending the heterogeneity existence or not. Our meta‐analysis found a significant association between ALDH2 gene G487A polymorphism and CAD in the Chinese population under allele (OR: 1.830, 95% CI: 1.560–2.140, P = 1.36 × 10^?13), recessive (OR: 1.920, 95% CI: 1.530–2.390, P = 1.20 × 10^?8), dominant (OR: 1.593, 95% CI: 1.336–1.900, P = 2.22 × 10^?7), homozygous (OR: 2.280, 95% CI: 1.810–2.870, P = 3.17 × 10^?12) and heterozygous genetic models (OR: 3.330, 95% CI: 2.070–5.370, P = 7.81 × 10^?7). The positive correlation between the ALDH2 gene G487A polymorphism and CAD makes the mutation a strong candidate as a genetic risk marker for CAD. Through further analysis, we also found that A allele carriers of ALDH2 gene G487A polymorphism may be particularly susceptible to CAD.     

Profile of common prostate cancer risk variants in an unscreened Romanian population

To find sequence variants affecting prostate cancer (PCA) susceptibility in an unscreened Romanian population we use a genome‐wide association study (GWAS). The study population included 990 unrelated pathologically confirmed PCA cases and 1034 male controls. DNA was genotyped using Illumina SNP arrays, and 24.295.558 variants were imputed using the 1000 Genomes data set. An association test was performed between the imputed markers and PCA. A systematic literature review for variants associated with PCA risk identified 115 unique variants that were tested in the Romanian sample set. Thirty of the previously reported SNPs replicated (P‐value < 0.05), with the strongest associations observed at: 8q24.21, 11q13.3, 6q25.3, 5p15.33, 22q13.2, 17q12 and 3q13.2. The replicated variants showing the most significant association in Romania are rs1016343 at 8q24.21 (P = 2.2 × 10^?4), rs7929962 at 11q13.3 (P = 2.7 × 10^?4) and rs9364554 at 6q25.2 (P = 4.7 × 10^?4). None of the variants tested in the Romanian GWAS reached genome‐wide significance (P‐value <5 × 10^?8) but 807 markers had P‐values <1 × 10^?4. Here, we report the results of the first GWAS of PCA performed in a Romanian population. Our study provides evidence that a substantial fraction of previously validated PCA variants associate with risk in this unscreened Romanian population.     

Changes in satellite‐derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006

Monitoring changes in vegetation growth has been the subject of considerable research during the past several decades, because of the important role of vegetation in regulating the terrestrial carbon cycle and the climate system. In this study, we combined datasets of satellite‐derived Normalized Difference Vegetation Index (NDVI) and climatic factors to analyze spatio‐temporal patterns of changes in vegetation growth and their linkage with changes in temperature and precipitation in temperate and boreal regions of Eurasia (> 23.5°N) from 1982 to 2006. At the continental scale, although a statistically significant positive trend of average growing season NDVI is observed (0.5 × 10^?3 year^?1, P = 0.03) during the entire study period, there are two distinct periods with opposite trends in growing season NDVI. Growing season NDVI has first significantly increased from 1982 to 1997 (1.8 × 10^?3 year^?1, P < 0.001), and then decreased from 1997 to 2006 (?1.3 × 10^?3 year^?1, P = 0.055). This reversal in the growing season NDVI trends over Eurasia are largely contributed by spring and summer NDVI changes. Both spring and summer NDVI significantly increased from 1982 to 1997 (2.1 × 10^?3 year^?1, P = 0.01; 1.6 × 10^?3 year^?1P < 0.001, respectively), but then decreased from 1997 to 2006, particularly summer NDVI which may be related to the remarkable decrease in summer precipitation (?2.7 mm yr^?1, P = 0.009). Further spatial analyses supports the idea that the vegetation greening trend in spring and summer that occurred during the earlier study period 1982–1997 was either stalled or reversed during the following study period 1997–2006. But the turning point of vegetation NDVI is found to vary across different regions.     

Variant in the 3′ Region of the IκBα Gene Associated With Insulin Resistance in Hispanic Americans: The IRAS Family Study

The IKKβ/NF‐κB pathway is known to play an important role in inflammatory response and has also recently been implicated in the process of insulin resistance. We hypothesized that one or more variants in the IκBα gene (NFKBIA) or surrounding untranslated regions would be associated with insulin sensitivity (S_I) in Hispanic‐American families. We tested for association between 25 single‐nucleotide polymorphisms (SNPs) in and near NFKBIA and S_I in 981 individuals in 90 Hispanic‐American families from the Insulin Resistance Atherosclerosis (IRAS) Family Study. SNP rs1951276 in the 3′ flanking region of NFKBIA was associated with S_I in the San Antonio (SA) sample after adjusting for age, gender, and admixture (uncorrected P = 1.69 × 10^?5; conservative Bonferroni correction P = 3.38 × 10^?4). Subjects with at least one A allele for NFKBIA rs1951276 had ～29% lower S_I compared to individuals homozygous for the G allele in the SA sample. Although not statistically significant, the effect was in the same direction in the San Luis Valley (SLV) sample alone (P = 0.348) and was significant in the combined SA and SLV samples (P = 5.37 × 10^?4; presence of A allele associated with ～20% lower S_I). In SA, when adjusted for subcutaneous adipose tissue area (SAT, cm²), the association was modestly attenuated (P = 1.25 × 10^?3), but the association remained highly significant after adjustment for visceral adipose tissue area (VAT, cm²; P = 4.41 × 10^?6). These results provide corroborating evidence that the NF‐κB/IKKβ pathway may mediate obesity‐induced insulin resistance in humans.     

β‐Cell Function and Insulin Sensitivity in Adolescents From an OGTT

Given the increase in the incidence of insulin resistance, obesity, and type 2 diabetes in children and adolescents, it would be of paramount importance to assess quantitative indices of insulin secretion and action during a physiological perturbation, such as a meal or an oral glucose‐tolerance test (OGTT). A minimal model method is proposed to measure quantitative indices of insulin secretion and action in adolescents from an oral test. A 7 h, 21‐sample OGTT was performed in 11 adolescents. The C‐peptide minimal model was identified on C‐peptide and glucose data to quantify indices of β‐cell function: static φ_s and dynamic φ_d responsivity to glucose from which total responsivity φ was also measured. The glucose minimal model was identified on glucose and insulin data to estimate insulin sensitivity, S_I, which was compared to a reference measure, S_I^ref, provided by a tracer method. Disposition indices, which adjust insulin secretion for insulin action, were then calculated. Indices of β‐cell function were φ_s = 51.35 ± 8.89 × 10^?9min^?1, φ_d = 1,392 ± 258 × 10^?9, and φ = 82.09 ± 17.70 × 10^?9min^?1. Insulin sensitivity was S_I = 14.19 ± 2.73 × 10^?4, not significantly different from S_I^ref = 14.96 ± 3.04 × 10^?4 dl/kg·min per µU/ml, and well correlated: r = 0.98, P < 0.0001, thus indicating that S_I can be accurately measured from an oral test. Disposition indices were DI_s = 1,040 ± 201 × 10^?14 dl/kg/min² per pmol/l, DI_d = 33,178 ± 10,720 × 10^?14 dl/kg/min per pmol/l, DI = 1,844 ± 522 × 10^?14 dl/kg/min² per pmol/l. Virtually the same minimal model assessment was obtained with a reduced 3 h, 9‐sample protocol. OGTT interpreted with C‐peptide and glucose minimal model has the potential to provide novel insight regarding the regulation of glucose metabolism in adolescents, and to evaluate the effect of obesity and interventions such as diet and exercise.     

Flow injection chemiluminescence determination of 6‐mercaptopurine based on a new system of potassium permanganate–thioacetamide–sodium hexametaphosphate

A novel chemiluminescence method for the determination of 6‐mercaptopurine was established based on 6‐mercaptopurine inhibition of the chemiluminescence emission of potassium permanganate–thioacetamide–sodium hexametaphosphate system. The peak height was proportional to log 6‐mercaptopurine concentration in the range 7.0 × 10^?10 to 1.0 × 10^?7 g/mL and the detection limit was 1.9 × 10^?11 g/mL (S/N = 3). The relative standard deviation was 1.5% for the determination of 8.0 × 10^?8 g/mL 6‐mercaptopurine (n = 11). The proposed sensor was successfully applied to the analysis of 6‐mercaptopurine in human serum samples. Copyright © 2009 John Wiley & Sons, Ltd.     

A genome‐wide association study of bipolar disorder and comorbid migraine

Both migraine and bipolar affective disorder (BPAD) are complex phenotypes with significant genetic and nongenetic components. Epidemiological and clinical studies have showed a high degree of comorbidity between migraine and BPAD, and overlapping regions of linkage have been shown in numerous genome‐wide linkage studies. To identify susceptibility factors for the BPAD/migraine phenotype, we conducted a genome‐wide association study (GWAS) in 1001 cases with bipolar disorder collected through the NIMH Genetics Initiative for Bipolar Disorder and genotyped at 1 m single‐nucleotide polymorphisms (SNPs) as part of the Genetic Association Information Network (GAIN). We compared BPAD patients without any headache (n = 699) with BPAD patients with doctor diagnosed migraine (n = 56). The strongest evidence for association was found for several SNPs in a 317‐kb region encompassing the uncharacterized geneKIAA0564 {e.g. rs9566845 [OR = 4.98 (95% CI: 2.6–9.48), P = 7.7 × 10^?8] and rs9566867 (P = 8.2 × 10^?8)}. Although the level of signficance was significantly reduced when using the Fisher's exact test (as a result of the low count of cases with migraine), rs9566845 P = 1.4 × 10^?5 and rs9566867 P = 1.5 × 10^?5, this region remained the most prominent finding. Furthermore, marker rs9566845 was genotyped and found associated with migraine in an independent Norwegian sample of adult attention deficit hyperactivity disorder (ADHD) patients with and without comorbid migraine (n = 131 and n = 324, respectively), OR = 2.42 (1.18–4.97), P = 0.013. This is the first GWAS examining patients with bipolar disorder and comorbid migraine. These data suggest that genetic variants in the KIAA0564 gene region may predispose to migraine headaches in subgroups of patients with both BPAD and ADHD.     

A genome-wide association study identifies a new variant associated with word reading fluency in Chinese children

Reading disability exhibited defects in different cognitive domains, including word reading fluency, word reading accuracy, phonological awareness, rapid automatized naming and morphological awareness. To identify the genetic basis of Chinese reading disability, we conducted a genome-wide association study (GWAS) of the cognitive traits related to Chinese reading disability in 2284 unrelated Chinese children. Among the traits analyzed in the present GWAS, we detected one genome-wide significant association (p < 5 × 10⁻⁸) on word reading fluency for one SNP on 4p16.2, within EVC genes (rs6446395, p = 7.33 × 10⁻¹⁰). Rs6446395 also showed significant association with Chinese character reading accuracy (p = 2.95 × 10⁻⁴), phonological awareness (p = 7.11 × 10⁻³) and rapid automatized naming (p = 4.71 × 10⁻³), implying multiple effects of this variant. The eQTL data showed that rs6446395 affected EVC expression in the cerebellum. Gene-based analyses identified a gene (PRDM10) to be associated with word reading fluency at the genome-wide level. Our study discovered a new candidate susceptibility variant for reading ability and provided new insights into the genetics of developmental dyslexia in Chinese children.     

A genome‐wide search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Heschl's gyrus (HG) is a core region of the auditory cortex whose morphology is highly variable across individuals. This variability has been linked to sound perception ability in both speech and music domains. Previous studies show that variations in morphological features of HG, such as cortical surface area and thickness, are heritable. To identify genetic variants that affect HG morphology, we conducted a genome‐wide association scan (GWAS) meta‐analysis in 3054 healthy individuals using HG surface area and thickness as quantitative traits. None of the single nucleotide polymorphisms (SNPs) showed association P values that would survive correction for multiple testing over the genome. The most significant association was found between right HG area and SNP rs72932726 close to gene DCBLD2 (3q12.1; P = 2.77 × 10^?7). This SNP was also associated with other regions involved in speech processing. The SNP rs333332 within gene KALRN (3q21.2; P = 2.27 × 10^?6) and rs143000161 near gene COBLL1 (2q24.3; P = 2.40 × 10^?6) were associated with the area and thickness of left HG, respectively. Both genes are involved in the development of the nervous system. The SNP rs7062395 close to the X‐linked deafness gene POU3F4 was associated with right HG thickness (Xq21.1; P = 2.38 × 10^?6). This is the first molecular genetic analysis of variability in HG morphology.     

Meta‐analysis Added Power to Identify Variants in FTO Associated With Type 2 Diabetes and Obesity in the Asian Population

Several common variants in the intron 1 of FTO (fat mass and associated obesity) gene have been reliably associated with BMI and obesity in European populations. We analyzed two variants (rs9939609 and rs8050136) in 4,189 Chinese Han individuals and conducted a meta‐analysis of published studies in Asian population to investigate whether these variants are associated with type 2 diabetes (T2D) and obesity in Asian population. In this study, both the minor allele A of rs9939609 and the minor allele A of rs805136 were associated with increased risk of T2D, independent of measures of BMI; the odds ratios (ORs) per copy of the risk allele were 1.19 for rs9939609 (95% confidence interval (CI), 1.04–1.37; P = 0.01) and 1.22 for rs8050136 (95% CI, 1.07–1.40; P = 0.004) after adjusting for age, sex, and BMI. Our results also showed association with risk of obesity (rs9939609: OR = 1.39 (95% CI 1.04–1.85), P = 0.02; rs8050136: OR = 1.45 (95% CI 1.09–1.93), P = 0.01) but no association with overweight. These results were consistent with the pooled results from our meta‐analysis study (for diabetes, rs8050136, P = 1.3 × 10^?3; rs9939609, P = 9.8 × 10^?4; for obesity, rs8050136, P = 2.2 × 10^?7; rs9939609, P = 9.0 × 10^?9). Our findings indicate that the two variants (rs9939609 and rs8050136) in the FTO gene contribute to obesity and T2D in the Asian populations.     

Methylation of SOCS3 is inversely associated with metabolic syndrome in an epigenome-wide association study of obesity

Epigenetic mechanisms, including DNA methylation, mediate the interaction between gene and environment and may play an important role in the obesity epidemic. We assessed the relationship between DNA methylation and obesity in peripheral blood mononuclear cells (PBMCs) at 485,000 CpG sites across the genome in family members (8-90 y of age) using a discovery cohort (192 individuals) and a validation cohort (1,052 individuals) of Northern European ancestry. After Bonferroni-correction (P_α=0.05 = 1.31 × 10^?7) for genome-wide significance, we identified 3 loci, cg18181703 (SOCS3), cg04502490 (ZNF771), and cg02988947 (LIMD2), where methylation status was associated with body mass index percentile (BMI%), a clinical index for obesity in children, adolescents, and adults. These sites were also associated with multiple metabolic syndrome (MetS) traits, including central obesity, fat depots, insulin responsiveness, and plasma lipids. The SOCS3 methylation locus was also associated with the clinical definition of MetS. In the validation cohort, SOCS3 methylation status was found to be inversely associated with BMI% (P = 1.75 × 10^?6), waist to height ratio (P = 4.18 × 10^?7), triglycerides (P = 4.01 × 10^?4), and MetS (P = 4.01 × 10^?7), and positively correlated with HDL-c (P = 4.57 × 10^?8). Functional analysis in a sub cohort (333 individuals) demonstrated SOCS3 methylation and gene expression in PBMCs were inversely correlated (P = 2.93 × 10^?4) and expression of SOCS3 was positively correlated with status of MetS (P = 0.012). We conclude that epigenetic modulation of SOCS3, a gene involved in leptin and insulin signaling, may play an important role in obesity and MetS.     

Proteomics of the milk fat globule membrane from Camelus dromedarius

Camel milk has been widely characterized with regards to casein and whey proteins. However, in camelids, almost nothing is known about the milk fat globule membrane (MFGM), the membrane surrounding fat globules in milk. The purpose of this study was thus to identify MFGM proteins from Camelus dromedarius milk. Major MFGM proteins (namely, fatty acid synthase, xanthine oxidase, butyrophilin, lactadherin, and adipophilin) already evidenced in cow milk were identified in camel milk using MS. In addition, a 1D‐LC‐MS/MS approach led us to identify 322 functional groups of proteins associated with the camel MFGM. Dromedary MFGM proteins were then classified into functional categories using DAVID (the Database for Annotation, Visualization, and Integrated Discovery) bioinformatics resources. More than 50% of MFGM proteins from camel milk were found to be integral membrane proteins (mostly belonging to the plasma membrane), or proteins associated to the membrane. Enriched GO terms associated with MFGM proteins from camel milk were protein transport (p‐value = 1.73 × 10^?14), translation (p‐value = 1.08 × 10^?11), lipid biosynthetic process (p‐value = 6.72 × 10^?10), hexose metabolic process (p‐value = 1.89 × 10^?04), and actin cytoskeleton organization (p‐value = 2.72 × 10^?04). These findings will help to contribute to a better characterization of camel milk. Identified MFGM proteins from camel milk may also provide new insight into lipid droplet formation in the mammary epithelial cell.