Knowledge of epigenetically regulated biomarkers linked to obesity development is still scarce. Improving molecular understanding of the involved factors and pathways would improve obesity phenotype characterization and reveal potentially relevant targets for obesity intervention. The Illumina Infinium HumanMethylation450 BeadChip was used in a leucocyte epigenome-wide association study (EWAS) to quantify differential DNA methylation in 60 lean compared with 60 obese young women. Replication was done in monozygotic twins discordant for obesity. At adolescence and adulthood, the two weight groups differed significantly in obesity-related traits and metabolic risk factors. Differential hypomethylation was overrepresented in obese compared to lean women. In the adjusted model, the EWAS revealed 10 differentially methylated CpG sites linked to 8 gene loci COX6A1P2/FGD2, SBNO2, TEX41, RPS6KA2, IGHE/IGHG1/IGHD, DMAP1, SOCS3, and SETBP1- and an enhancer locus at chromosome 2 (2p25.1). The sites linked to TEX41, IGHE/IGHG1/IGHD, DMAP1, and SETBP1 were novel findings, while COX6A1P/FGD2, SBNO2, RPS6KA2, and SOCS3 had been identified previously with concordant direction of effects. RPS6KA2, DMAP1, and SETBP1 were replicated in the BMI-discordant monozygotic twin cohort using the FDR of 5%.

Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene (IGH) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2, and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 (p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations (p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.