Differential DNA methylation of steatosis and non-alcoholic fatty liver disease in adolescence.

BACKGROUND AND AIMS: Epigenetic modifications are associated with hepatic fat accumulation and non-alcoholic fatty liver disease (NAFLD). However, few epigenetic modifications directly implicated in such processes have been identified during adolescence, a critical developmental window where physiological changes could influence future disease trajectory. To investigate the association between DNA methylation and NAFLD in adolescence, we undertook discovery and validation of novel methylation marks, alongside replication of previously reported marks. APPROACH AND RESULTS: We performed a DNA methylation epigenome-wide association study (EWAS) on DNA from whole blood from 707 Raine Study adolescents phenotyped for steatosis score and NAFLD by ultrasound at age 17. Next, we performed pyrosequencing validation of loci within the most 100 strongly associated differentially methylated CpG sites (dmCpGs) for which ≥ 2 probes per gene remained significant across four statistical models with a nominal p value < 0.007. EWAS identified dmCpGs related to three genes (ANK1, MIR10a, PTPRN2) that met our criteria for pyrosequencing. Of the dmCpGs and surrounding loci that were pyrosequenced (ANK1 n = 6, MIR10a n = 7, PTPRN2 n = 3), three dmCpGs in ANK1 and two in MIR10a were significantly associated with NAFLD in adolescence. After adjustment for waist circumference only dmCpGs in ANK1 remained significant. These ANK1 CpGs were also associated with γ-glutamyl transferase and alanine aminotransferase concentrations. Three of twenty-two differentially methylated dmCpGs previously associated with adult NAFLD were associated with NAFLD in adolescence (all adjusted p < 2.3 × 10-3). CONCLUSIONS: We identified novel DNA methylation loci associated with NAFLD and serum liver biochemistry markers during adolescence, implicating putative dmCpG/gene regulatory pathways and providing insights for future mechanistic studies.

Adiposity associated DNA methylation signatures in adolescents are related to leptin and perinatal factors.

Epigenetics links perinatal influences with later obesity. We identifed differentially methylated CpG (dmCpG) loci measured at 17 years associated with concurrent adiposity measures and examined whether these were associated with hsCRP, adipokines, and early life environmental factors. Genome-wide DNA methylation from 1192 Raine Study participants at 17 years, identified 29 dmCpGs (Bonferroni corrected p < 1.06E-07) associated with body mass index (BMI), 10 with waist circumference (WC) and 9 with subcutaneous fat thickness. DmCpGs within Ras Association (RalGDS/AF-6), Pleckstrin Homology Domains 1 (RAPH1), Musashi RNA-Binding Protein 2 (MSI2), and solute carrier family 25 member 10 (SLC25A10) are associated with both BMI and WC. Validation by pyrosequencing confirmed these associations and showed that MSI2 , SLC25A10 , and RAPH1 methylation was positively associated with serum leptin. These were  also associated with the early environment; MSI2 methylation (ß = 0.81, p = 0.0004) was associated with pregnancy maternal smoking, SLC25A10 (CpG2 ß = 0.12, p = 0.002) with pre- and early pregnancy BMI, and RAPH1 (ß = -1.49, p = 0.036) with gestational weight gain. Adjusting for perinatal factors, methylation of the dmCpGs within MSI2, RAPH1, and SLC25A10 independently predicted BMI, accounting for 24% of variance. MSI2 methylation was additionally associated with BMI over time (17 years old ß = 0.026, p = 0.0025; 20 years old ß = 0.027, p = 0.0029) and between generations (mother ß = 0.044, p = 7.5e-04). Overall findings suggest that DNA methylation in MSI2, RAPH1, and SLC25A10 in blood may be robust markers, mediating through early life factors.

Ratio of Klebsiella/Bifidobacterium in early life correlates with later development of paediatric allergy.

Several studies have reported that intestinal microbial colonisation patterns differ between non-allergic and allergic infants. However, the microbial signature underlying the pathogenesis of allergies remains unclear. We aim to gain insight into the development of the intestinal microbiota of healthy infants and infants who develop allergy in early life, and identify potential microbiota biomarkers of later allergic disease. Using a case-control design in a Chinese sub-cohort of a Singaporean birth cohort (GUSTO), we utilised 16S rRNA gene sequencing to assess intestinal microbial composition and diversity of 21 allergic and 18 healthy infants at 3 weeks, 3 months and 6 months of age, and correlated the microbiota with allergy at ages 18 and 36 months. Pronounced differences in intestinal microbiota composition between allergic and healthy infants were observed at 3 months of age. The intestine of healthy infants was colonised with higher abundance of commensal Bifidobacterium. Conversely, Klebsiella, an opportunistic pathogen, was significantly enriched in the allergic infants. Interestingly, infants with a high Klebsiella/Bifidobacterium (K/B) ratio (above the population median K/B ratio) at age 3 months had an odds ratio of developing allergy by 3 years of age of 9.00 (95% confidence interval 1.46-55.50) compared to those with low K/B ratio. This study demonstrated a relationship between the ratio of genera Klebsiella and Bifidobacterium during early infancy and development of paediatric allergy in childhood. Our study postulates that an elevated K/B ratio in early infancy could be a potential indicator of an increased risk of allergy development. This line of research might enable future intervention strategies in early life to prevent or treat allergy. Our study provides new insights into microbial signatures associated with childhood allergy, in particular, suggests that an elevated K/B ratio could be a potential early-life microbiota biomarker of allergic disease.

Does the maternal vaginal microbiota play a role in seeding the microbiota of neonatal gut and nose?

The acquisition and early maturation of infant microbiota is not well understood despite its likely influence on later health. We investigated the contribution of the maternal microbiota to the microbiota of infant gut and nose in the context of mode of delivery and feeding. Using 16S rRNA sequencing and specific qPCR, we profiled microbiota of 42 mother-infant pairs from the GUSTO birth cohort, at body sites including maternal vagina, rectum and skin; and infant stool and nose. In our study, overlap between maternal vaginal microbiota and infant faecal microbiota was minimal, while the similarity between maternal rectal microbiota and infant microbiota was more pronounced. However, an infant's nasal and gut microbiota were no more similar to that of its own mother, than to that of unrelated mothers. These findings were independent of delivery mode. We conclude that the transfer of maternal vaginal microbes play a minor role in seeding infant stool microbiota. Transfer of maternal rectal microbiota could play a larger role in seeding infant stool microbiota, but approaches other than the generally used analyses of community similarity measures are likely to be needed to quantify bacterial transmission. We confirmed the clear difference between microbiota of infants born by Caesarean section compared to vaginally delivered infants and the impact of feeding mode on infant gut microbiota. Only vaginally delivered, fully breastfed infants had gut microbiota dominated by Bifidobacteria. Our data suggest that reduced transfer of maternal vaginal microbial is not the main mechanism underlying the differential infant microbiota composition associated with Caesarean delivery. The sources of a large proportion of infant microbiota could not be identified in maternal microbiota, and the sources of seeding of infant gut and nasal microbiota remain to be elucidated.

Low birth weight associates with hippocampal gene expression.

Birth weight predicts the lifetime risk for psychopathology suggesting that the quality of fetal development influences the predisposition for mental disorders. The connectivity and synaptic network of the hippocampus are implicated in depression, schizophrenia and anxiety. We thus examined the underlying molecular adaptations in the hippocampus as a function of the fetal conditions associated with low birth weight. We used tissues from the non-human primate, Macaca fascicularis, to identify changes in hippocampal gene expression early in postnatal development associated with naturally occurring low compared with normal birth weight. Microarrays were used to analyze gene expression and DNA methylation in the hippocampus of five low- and five normal-birth weight neonates. Real-time PCR was employed to validate differentially expressed genes. Birth weight associated with altered global transcription in the hippocampus. Hierarchical clustering of gene expression profiles from 24,154 probe sets grouped all samples except one by their birth weight status. Differentially expressed genes were enriched in biological processes associated with neuronal projection, positive regulation of transcription and apoptosis. About 4% of the genes with differential expression co-varied with DNA methylation levels. The data suggest that low birth weight is closely associated with hippocampal gene expression with a small epigenetic underpinning by DNA methylation in neonates. The data also provide a potential molecular basis for the developmental origin of an enhanced risk for mental disorders.

MC3R gene polymorphisms are associated with early childhood adiposity gain and infant appetite in an Asian population.

BACKGROUND: Polymorphic variants within human melanocortin-3 receptor gene (MC3R) gene have been associated with obesity. However, its influence on infancy and early childhood adiposity has not been reported before. OBJECTIVES: We assessed associations between genotype at polymorphic sites within MC3R with early childhood adiposity and interaction with early childhood appetitive traits. METHODS: We studied 1090 singletons in an Asian mother-offspring cohort genotyped for MC3R and in a subgroup (n = 422) who had completed Child Eating Behaviour Questionnaires (CEBQ) at 12 months. Children were followed from birth to 48 months, and up to 10 measurements of body mass index and five measures of triceps and subscapular skin-folds were obtained. RESULTS: Independent of potential confounders, each additional MC3R minor allele copy was associated with greater body mass index standard deviation score [B{95% confidence interval}: 0.004 units/month {0.001,0.007}; p = 0.007], triceps [0.009 mm/month {0.001,0.02}; p = 0.021] and subscapular skin-fold [0.008 mm/month {0.002,0.01}; p = 0.011] gain velocity in the first 48 months. Each additional MC3R minor allele copy was also associated with increased odds of overweight [odds ratio {95% confidence interval}: 1.48{1.17-1.88}] and obesity [1.58{1.10-2.28}] in the first 48 months. Every additional copy of MC3R minor allele was positively associated with 'slowness-in-eating' appetitive trait [0.24{0.06,0.39}, p = 0.006]; however, the relationship between 'slowness-in-eating' with adiposity gain was not statistically significant. CONCLUSIONS: Our findings support the role of MC3R genetic variants in adiposity gain during early childhood.

Genome-wide methylation analysis identifies differentially methylated CpG loci associated with severe obesity in childhood.

Childhood obesity is a major public health issue. Here we investigated whether differential DNA methylation was associated with childhood obesity. We studied DNA methylation profiles in whole blood from 78 obese children (mean BMI Z-score: 2.6) and 71 age- and sex-matched controls (mean BMI Z-score: 0.1). DNA samples from obese and control groups were pooled and analyzed using the Infinium HumanMethylation450 BeadChip array. Comparison of the methylation profiles between obese and control subjects revealed 129 differentially methylated CpG (DMCpG) loci associated with 80 unique genes that had a greater than 10% difference in methylation (P-value < 0.05). The top pathways enriched among the DMCpGs included developmental processes, immune system regulation, regulation of cell signaling, and small GTPase-mediated signal transduction. The associations between the methylation of selected DMCpGs with childhood obesity were validated using sodium bisulfite pyrosequencing across loci within the FYN, PIWIL4, and TAOK3 genes in individual subjects. Three CpG loci within FYN were hypermethylated in obese individuals (all P < 0.01), while obesity was associated with lower methylation of CpG loci within PIWIL4 (P = 0.003) and TAOK3 (P = 0.001). After building logistic regression models, we determined that a 1% increase in methylation in TAOK3, multiplicatively decreased the odds of being obese by 0.91 (95% CI: 0.86 - 0.97), and an increase of 1% methylation in FYN CpG3, multiplicatively increased the odds of being obese by 1.03 (95% CI: 0.99 - 1.07). In conclusion, these findings provide evidence that childhood obesity is associated with specific DNA methylation changes in whole blood, which may have utility as biomarkers of obesity risk.

Measuring epigenetics as the mediator of gene/environment interactions in DOHaD.

Analysis of DNA methylation data in epigenome-wide association studies provides many bioinformatics and statistical challenges. Not least of these, are the non-independence of individual DNA methylation marks from each other, from genotype and from technical sources of variation. In this review we discuss DNA methylation data from the Infinium450K array and processing methodologies to reduce technical variation. We describe recent approaches to harness the concordance of neighbouring DNA methylation values to improve power in association studies. We also describe how the non-independence of genotype and DNA methylation has been used to infer causality (in the case of Mendelian randomization approaches); suggest the mediating effect of DNA methylation in linking intergenic single nucleotide polymorphisms, identified in genome-wide association studies, to phenotype; and to uncover the widespread influence of gene and environment interactions on methylation levels.

A new, improved and generalizable approach for the analysis of biological data generated by -omic platforms.

The principles embodied by the Developmental Origins of Health and Disease (DOHaD) view of 'life history' trajectory are increasingly underpinned by biological data arising from molecular-based epigenomic and transcriptomic studies. Although a number of 'omic' platforms are now routinely and widely used in biology and medicine, data generation is frequently confounded by a frequency distribution in the measurement error (an inherent feature of the chemistry and physics of the measurement process), which adversely affect the accuracy of estimation and thus, the inference of relationships to other biological measures such as phenotype. Based on empirical derived data, we have previously derived a probability density function to capture such errors and thus improve the confidence of estimation and inference based on such data. Here we use published open source data sets to calculate parameter values relevant to the most widely used epigenomic and transcriptomic technologies Then by using our own data sets, we illustrate the benefits of this approach by specific application, to measurement of DNA methylation in this instance, in cases where levels of methylation at specific genomic sites represents either (1) a response variable or (2) an independent variable. Further, we extend this formulation to consideration of the 'bivariate' case, in which the co-dependency of methylation levels at two distinct genomic sites is tested for biological significance. These tools not only allow greater accuracy of measurement and improved confidence of functional inference, but in the case of epigenomic data at least, also reveal otherwise cryptic information.

Molecular cloning and characterisation of a novel GABAB-related G-protein coupled receptor.

Using a homology-based bioinformatics approach we have analysed human genomic sequence and identified the human and rodent orthologues of a novel putative seven transmembrane G protein coupled receptor, termed GABA(BL). The amino acid sequence homology of these cDNAs compared to GABA(B1) and GABA(B2) led us to postulate that GABA(BL) was a putative novel GABA(B) receptor subunit. The C-terminal sequence of GABA(BL) contained a putative coiled-coil domain, di-leucine and several RXR(R) ER retention motifs, all of which have been shown to be critical in GABA(B) receptor subunit function. In addition, the distribution of GABA(BL) in the central nervous system was reminiscent of that of the other known GABA(B) subunits. However, we were unable to detect receptor function in response to any GABA(B) ligands when GABA(BL) was expressed in isolation or in the presence of either GABA(B1) or GABA(B2). Therefore, if GABA(BL) is indeed a GABA(B) receptor subunit, its partner is a potentially novel receptor subunit or chaperone protein which has yet to be identified.

Molecular adaptation of alanine:glyoxylate aminotransferase targeting in primates.

The intermediary metabolic enzyme alanine:glyoxylate aminotransferase (AGT) is targeted to different organelles (mitochondria and/or peroxisomes) in different species. Possibly under the influence of dietary selection pressure, the subcellular distribution of AGT has changed on at least eight occasions during the evolution of mammals. AGT targeting is dependent on the variable use of two alternative transcription and translation initiation sites which determine whether or not the region encoding the N-terminal mitochondrial targeting sequence is contained within the open reading frame. In the present study, we sequenced the 5' region of the AGT gene, including both ancestral translation start sites, for 11 anthropoid primates and compared the results with data already available for two others. We show that while the more 3' of the two translation start sites is maintained in all species, the more 5' site has been lost in six species (five of seven catarrhines and one of six platyrrhines). In addition, the remaining two catarrhines, which have maintained the 5' translation start site, are predicted to have lost mitochondrial targeting by a different mechanism, possibly loss of the more 5' transcription start site. Analysis of the relative frequencies of nonsynonymous and synonymous mutations in the region encoding the extant or ancestral mitochondrial targeting sequences led us to suggest that there has been recent strong positive selection pressure to lose, or decrease the efficiency of, mitochondrial AGT targeting in several anthropoid lineages, and that the loss of mitochondrial targeting in this group of mammals is likely to have occurred on at least four, and possibly five, separate occasions.