Epigenetic alterations in patients with anorexia nervosa-a systematic review.

Anorexia nervosa (AN) is a complex metabolic and psychological disorder that is influenced by both heritable genetic components and environmental factors. Exposure to various environmental influences can lead to epigenetically induced changes in gene expression. Epigenetic research in AN is still in its infancy, and studies to date are limited in determining clear, valid links to disease onset and progression are limited. Therefore, the aim of this systematic review was to compile and critically evaluate the available results of epigenetic studies specifically in AN and to provide recommendations for future studies. In accordance with the PRISMA guidelines, a systematic literature search was performed in three different databases (PubMed, Embase, and Web of Science) through May 2023. Twenty-three original papers or conference abstracts on epigenetic studies in AN were collected. Epigenome-wide association studies (EWASs), which analyze DNA methylation across the genome in patients with AN and identify potential disease-relevant changes in promoter/regulatory regions of genes, are the most promising for future research. To date, five EWASs on AN have been published, suggesting a potential reversibility of malnutrition-induced epigenetic changes once patients recover. Hence, determining differential DNA methylation levels could serve as a biomarker for disease status or early diagnosis and might be involved in disease progression or chronification. For future research, EWASs with a larger sample size, longitudinal study design and uniform methods should be performed to contribute to the understanding of the pathophysiology of AN, the development of individual interventions and a better prognosis for affected patients.

Unexpected identification of obesity-associated mutations in LEP and MC4R genes in patients with anorexia nervosa.

Mutations leading to a reduced or loss of function in genes of the leptin-melanocortin system confer a risk for monogenic forms of obesity. Yet, gain of function variants in the melanocortin-4-receptor (MC4R) gene predispose to a lower BMI. In individuals with reduced body weight, we thus expected mutations leading to an enhanced function in the respective genes, like leptin (LEP) and MC4R. Therefore, we have Sanger sequenced the coding regions of LEP and MC4R in 462 female patients with anorexia nervosa (AN), and 445 healthy-lean controls. In total, we have observed four and eight variants in LEP and MC4R, respectively. Previous studies showed different functional in vitro effects for the detected frameshift and non-synonymous variants: (1) LEP: reduced/loss of function (p.Val94Met), (2) MC4R: gain of function (p.Val103Ile, p.Ile251Leu), reduced or loss of function (p.Thr112Met, p.Ser127Leu, p.Leu211fsX) and without functional in vitro data (p.Val50Leut). In LEP, the variant p.Val94Met was detected in one patient with AN. For MC4R variants, one patient with AN carried the frameshift variant p.Leu211fsX. One patient with AN was heterozygous for two variants at the MC4R (p.Val103Ile and p.Ser127Leu). All other functionally relevant variants were detected in similar frequencies in patients with AN and lean individuals.

Evaluation of the MC3R gene pertaining to body weight and height regulation and puberty development.

Recent studies reported an impact of the melanocortin 3 receptor (MC3R) on the regulation of body weight, linear growth and puberty timing. Previously, allele p.44Ile of a frequent non-synonymous variant (NSV) p.Val44Ile was reported to be associated with decreased lean body mass (LBM) and later puberty in both sexes. We Sanger sequenced the coding region of MC3R in 185 children or adolescents with short normal stature (SNS) or 258 individuals with severe obesity, and 192 healthy-lean individuals. Eleven variants (six NSVs) were identified. In-silico analyses ensued. Three rare loss-of-function (LoF) variants (p.Phe45Ser, p.Arg220Ser and p.Ile298Ser) were only found in severely obese individuals. One novel highly conserved NSV (p.Ala214Val), predicted to increase protein stability, was detected in a single lean female. In the individuals with SNS, we observed deviation from Hardy-Weinberg Equilibrium (HWE) (p = 0.012) for p.Val44Ile (MAF = 11.62%). Homozygous p.44Ile carriers with SNS had an increased BMI, but this effect did not remain significant after Bonferroni correction. In line with previous findings, the detected LoF NSVs may suggest that dysfunction in MC3R is associated with decreased body height, obesity and delayed puberty.

Lipocalin 2 - mutation screen and serum levels in patients with anorexia nervosa or obesity and in lean individuals.

Context: The bone-derived adipokine lipocalin-2 is relevant for body weight regulation by stimulating the leptin-melanocortin pathway. Objective: We aimed to (i) detect variants in the lipocalin-2 gene (LCN2) which are relevant for body weight regulation and/or anorexia nervosa (AN); (ii) describe and characterize the impact of LCN2 and MC4R variants on circulating lipocalin-2 level. Methods: Sanger sequencing of the coding region of LCN2 in 284 children and adolescents with severe obesity or 287 patients with anorexia nervosa. In-silico analyses to evaluate functional implications of detected LCN2 variants. TaqMan assays for rare non-synonymous variants (NSVs) in additional independent study groups. Serum levels of lipocalin-2 were measured by ELISA in 35 females with NSVs in either LCN2 or MC4R, and 33 matched controls without NSVs in the two genes. Results: Fourteen LCN2-variants (five NSVs) were detected. LCN2-p.Leu6Pro and p.Gly9Val located in the highly conserved signal peptide region may induce functional consequences. The secondary structure change of lipocalin-2 due to LCN2-p.Val89Ile may decrease solubility and results in a low lipocalin-2 level in a heterozygotes carrier (female recovered from AN). Lean individuals had lower lipocalin-2 levels compared to patients with obesity (p = 0.033). Conclusion: Lipocalin-2 levels are positively associated with body mass index (BMI). Single LCN2-variants might have a profound effect on lipocalin-2 levels.

Prevalence estimates of putatively pathogenic leptin variants in the gnomAD database.

Homozygosity for pathogenic variants in the leptin gene leads to congenital leptin deficiency causing severe early-onset obesity. This monogenic form of obesity has mainly been detected in patients from consanguineous families. Prevalence estimates for the general population using the Exome Aggregation Consortium (ExAC) database reported a low frequency of leptin mutations. One in approximately 15 million individuals will be homozygous for a deleterious leptin variant. With the present study, we aimed to extend these findings utilizing the augmented Genome Aggregation Database (gnomAD) v2.1.1 including more than 140,000 samples. In total, 68 non-synonymous and 7 loss-of-function leptin variants were deposited in gnomAD. By predicting functional implications with the help of in silico tools, like SIFT, PolyPhen2 and MutationTaster2021, the prevalence of hetero- and homozygosity for putatively pathogenic variants (n = 32; pathogenic prediction by at least two tools) in the leptin gene were calculated. Across all populations, the estimated prevalence for heterozygosity for functionally relevant variants was approximately 1:2,100 and 1:17,830,000 for homozygosity. This prevalence deviated between the individual populations. Accordingly, people from East Asia and individuals of mixed ethnicities ('Others') were at greater risk to carry a possibly damaging leptin variant. Generally, this study emphasises the scarcity of pathogenic leptin variants in the general population with varying prevalence for distinct study groups.

Evidence for correlations between BMI-associated SNPs and circRNAs.

Circular RNAs (circRNAs) are regulators of processes like adipogenesis. Their expression can be modulated by SNPs. We analysed links between BMI-associated SNPs and circRNAs. First, we detected an enrichment of BMI-associated SNPs on circRNA genomic loci in comparison to non-significant variants. Analysis of sex-stratified GWAS data revealed that circRNA genomic loci encompassed more genome-wide significant BMI-SNPs in females than in males. To explore whether the enrichment is restricted to BMI, we investigated nine additional GWAS studies. We showed an enrichment of trait-associated SNPs in circRNAs for four analysed phenotypes (body height, chronic kidney disease, anorexia nervosa and autism spectrum disorder). To analyse the influence of BMI-affecting SNPs on circRNA levels in vitro, we examined rs4752856 located on hsa_circ_0022025. The analysis of heterozygous individuals revealed an increased level of circRNA derived from the BMI-increasing SNP allele. We conclude that genetic variation may affect the BMI partly through circRNAs.

PTBP2 - a gene with relevance for both Anorexia nervosa and body weight regulation.

Genetic factors are relevant for both eating disorders and body weight regulation. A recent genome-wide association study (GWAS) for anorexia nervosa (AN) detected eight genome-wide significant chromosomal loci. One of these loci, rs10747478, was also genome-wide and significantly associated with body mass index (BMI). The nearest coding gene is the Polypyrimidine Tract Binding Protein 2 gene (PTBP2). To detect mutations in PTBP2, Sanger sequencing of the coding region was performed in 192 female patients with AN (acute or recovered) and 191 children or adolescents with (extreme) obesity. Twenty-five variants were identified. Twenty-three of these were predicted to be pathogenic or functionally relevant in at least one in silico tool. Two novel synonymous variants (p.Ala77Ala and p.Asp195Asp), one intronic SNP (rs188987764), and the intronic deletion (rs561340981) located in the highly conserved region of PTBP2 may have functional consequences. Ten of 20 genes interacting with PTBP2 were studied for their impact on body weight regulation based on either previous functional studies or GWAS hits for body weight or BMI. In a GWAS for BMI (Pulit et al. 2018), the number of genome-wide significant associations at the PTBP2 locus was different between males (60 variants) and females (two variants, one of these also significant in males). More than 65% of these 61 variants showed differences in the effect size pertaining to BMI between sexes (absolute value of Z-score >2, two-sided p < 0.05). One LD block overlapping 5'UTR and all coding regions of PTBP2 comprises 56 significant variants in males. The analysis based on sex-stratified BMI GWAS summary statistics implies that PTBP2 may have a more pronounced effect on body weight regulation in males than in females.