Skip to main content

BRIEF RESEARCH REPORT article

Front. Genet., 07 March 2022
Sec. Genetics of Common and Rare Diseases

Targeted Exome Sequencing of Genes Involved in Rare CNVs in Early-Onset Severe Obesity

Petra Loid,,
Petra Loid1,2,3*Minna Pekkinen,,Minna Pekkinen1,2,3Taina MustilaTaina Mustila4Pivi TossavainenPäivi Tossavainen5Heli Viljakainen,Heli Viljakainen2,6Anna Lindstrand,Anna Lindstrand7,8Outi Mkitie,,,,Outi Mäkitie1,2,3,7,8
  • 1Children's Hospital, University of Helsinki and Helsinki University Hospital, Helsinki, Finland
  • 2Folkhälsan Research Center, Genetics Research Program, Helsinki, Finland
  • 3Research Program for Clinical and Molecular Metabolism, Faculty of Medicine, University of Helsinki, Helsinki, Finland
  • 4City of Turku Wellfare Services, Diabetes Care, Turku, Finland
  • 5Department of Pediatrics, PEDEGO Research Unit, Medical Research Center, Oulu University Hospital, University of Oulu, Oulu, Finland
  • 6Faculty of Medicine, University of Helsinki, Helsinki, Finland
  • 7Department of Molecular Medicine and Surgery, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
  • 8Department of Clinical Genetics, Karolinska University Hospital, Stockholm, Sweden

Context: Rare copy number variants (CNVs) have been associated with the development of severe obesity. However, the potential disease-causing contribution of individual genes within the region of CNVs is often not known.

Objective: Screening of rare variants in genes involved in CNVs in Finnish patients with severe early-onset obesity to find candidate genes linked to severe obesity.

Methods: Custom-made targeted exome sequencing panel to search for rare (minor allele frequency <0.1%) variants in genes affected by previously identified CNVs in 92 subjects (median age 14 years) with early-onset severe obesity (median body mass index (BMI) Z-score + 4.0).

Results: We identified thirteen rare heterozygous variants of unknown significance in eleven subjects in twelve of the CNV genes. Two rare missense variants (p.Pro405Arg and p.Tyr232Cys) were found in SORCS1, a gene highly expressed in the brain and previously linked to diabetes risk. Four rare variants were in genes in the proximal 16p11.2 region (a frameshift variant in TAOK2 and missense variants in SEZ6L2, ALDOA and KIF22) and three rare missense variants were in genes in the 22q11.21 region (AIFM3, ARVCF and KLHL22).

Conclusion: We report several rare variants in CNV genes in subjects with childhood obesity. However, the role of the individual genes in the previously identified rare CNVs to development of obesity remains uncertain. More studies are needed to understand the potential role of the specific genes within obesity associated CNVs.

Introduction

Copy number variants (CNVs) are structural variants that include duplications or deletions of genomic regions ranging from 50 bp to several megabases (Zarrei et al., 2015). These structural variants can arise by different mutational mechanisms such as replication, repair of double-strand breaks and recombination (Gu et al., 2008; Carvalho and Lupski, 2016). CNVs contribute to human diversity but also to several diseases, including neurodevelopmental disorders and severe obesity (Wang et al., 2010; Wheeler et al., 2013; D'Angelo et al., 2018; Saeed et al., 2020).

Deletion in 16p11.2 region is the most established CNV linked to early-onset severe obesity (Bochukova et al., 2010; Walters et al., 2010). CNVs also contribute to syndromic obesity (Vuillaume et al., 2014; Pettersson et al., 2017; D’Angelo et al., 2018; Selvanayagam et al., 2018; Saeed et al., 2020; Madeo et al., 2021). The identification of rare CNVs in obesity has led to the discovery of several novel candidate genes. However, in many cases the disease-causing gene within the region of CNVs is unclear. Indeed, unraveling the pathogenicity of CNVs remains a major challenge (Nowakowska, 2017). Pathogenic CNVs are usually large and contain genes with high evolutionary copy number conservation. It has been proposed that dosage sensitivity of individual genes within the region of CNVs is a frequent cause of CNV pathogenicity. Haploinsufficiency is commonly used as a model of gene dosage sensitivity and predictions of haploinsufficiency have been used to predict the impact of CNVs (Huang et al., 2010; Lek et al., 2016; Rice and McLysaght, 2017).

We previously identified rare CNVs enriched in subjects with severe childhood obesity (Pettersson et al., 2017). In a case-control study of 90 obese subjects and 67 normal-weight controls, using array comparative genomic hybridization, we identified rare CNVs in 17 obese subjects. However, the pathogenicity of several of the identified CNVs remained unclear and the potential phenotypic contribution of the genes within the CNVs is not known. We hypothesized that rare variants in the genes within the CNVs might contribute to the development of obesity. We designed a custom-made targeted exome sequencing panel, including all genes located in the regions of the rare CNVs identified in our cohort of patients with childhood obesity, in search for rare variants and novel candidate genes.

Methods

Study Subjects

This study included 92 children and adolescents with severe early-onset obesity, defined as height-adjusted weight >60% before 10 years of age according to Finnish growth standards (Saari et al., 2011). The study subjects were recruited through pediatric endocrine clinics at Helsinki University Hospital, Seinäjoki Central Hospital and Oulu University Hospital during years 2011–2017. We have previously investigated the study subjects by targeted exome sequencing for rare variants in 24 genes related to the leptin-melanocortin pathway or hypothalamic development (Loid et al., 2020). Altogether 80% of the patients in this study have also been screened for rare CNVs by array comparative genomic hybridization (a-CGH). In our previous study we identified rare CNVs in 17 obese subjects (Pettersson et al., 2017) and 14 of these 17 subjects with rare CNVs participated in the present study. The study was approved by Research Ethics Committees of Hospital District of Helsinki and Uusimaa, the Pirkanmaa Hospital District and the Northern Ostrobothnia Hospital District. Informed written consents were obtained from all the participants or their parents/guardians. Clinical and growth data were collected from hospital records and anthropometric measurements were obtained during a study visit. Age-and gender specific BMI Z-scores were derived based on the World Health Organization reference values (www.who.int/childgrowth/standards).

Targeted Exome Sequencing

We designed the probes for targeted exome sequencing using SeqCap EZ Choice Library and Nimble Design (Roche NimbleGen, United States). DNA capture and sequencing were performed at Oxford Genomics Centre. Sequence reads were aligned to reference genome GRCh37. In this study, we included 117 genes involved in previously identified rare CNVs (Pettersson et al., 2017). The genes are listed in Supplementary Table S1. We filtered variants using Varaft (version 2.17) (Desvignes et al., 2018). We selected coding and splice site variants with an allele frequency of <0.001 in the Genome Aggregation Database (gnomAD) database (http://gnomad.broadinstitute.org), the 1000 Genomes Project (http://www.internationalgenome.org), and the Sequencing Initiative Suomi database (SISu) (http://www.sisuproject.fi). In a second analysis, we also evaluated coding and splice site variants with an allele frequency of <0.005 to explore whether change of allele frequency threshold would add any potentially interesting discoveries, considering that severe obesity is not an extremely rare phenotype. In silico predictions were performed using Sorting Intolerant From Tolerant (SIFT) (Kumar P. et al., 2009), Polymorphism Phenotyping version 2 (PolyPhen2) (http://genetics.bwh.harvard.edu/pph2/), MutationTaster2 (Schwarz et al., 2014), Combined Annotation Dependent Depletion (CADD) scores (https://cadd.gs.washington.edu/), Mendelian Clinically Applicable Pathogenicity (M-CAP) (http://bejerano.stanford.edu/mcap/), Rare Exome Variant Ensemble Learner (REVEL) score (Ioannidis et al., 2016) and splice AI (https://spliceailookup.broadinstitute.org). We selected rare variants predicted to be pathogenic by several of the abovementioned prediction tools including a CADD score >20. Integrative Genomics Viewer (IGV) was used for visualization of the variants. The identified variants were classified according to the American Collage of Medical Genetics (ACMG) guidelines (Richards et al., 2015) using VarSome tool (Kopanos et al., 2019). Parental and siblings’ samples was obtained from two available families with the identified rare variants and Sanger sequencing was performed to determine inheritance pattern.

Results

Characteristics of Study Subjects

The cohort in this study comprised 92 subjects (51% males) with median age of 13.7 years (interquartile range (IQR) 10.6–16.8 years) and median BMI Z-score of +4.0 (IQR + 3.4 to + 4.9). The study subjects fulfilled our inclusion criteria (height-adjusted weight >60%) at median age of 6.0 years (IQR 4.5–7.0 years). Of the 68 school-aged study subjects for whom information on learning difficulties was available, ten (15%) had learning difficulties. Two study subjects had developmental delay based on previous clinical investigations. In these two patients Prader-Willi syndrome had been ruled out. Ten participants had insulin resistance based on oral glucose tolerance test.

Targeted Exome Sequencing

On average 99.3% of the bases in the coding exons were covered by at least 100 reads and the mean read depth of coverage was 750X. We identified 13 rare heterozygous variants in 11 subjects with early-onset obesity. Four of them had learning difficulties and one had developmental delay with autistic features. The variants with allele frequency of <0.001 were in 12 of the genes in the panel, namely AIFM3, ALDOA, ARVCF, BRAF, CHD1L, KIF22, KLHL22, MCTP2, SEZ6L2, SORCS1, TAOK2 and XRN1. Table 1 presents the rare variants with allele frequencies in gnomAD and in silico prediction values. For the rare splicing variant in CHDL1 we also assessed the SpliceAI score, which was 0.990 (donor loss). In one of the genes (SORCS1) we identified two rare variants; p.Pro405Arg and p.Tyr232Cys. Figure 1 presents the segregation of the SORCS1 variants in the two families. In both families the variant was inherited from a parent with overweight. However, also other family members without the variant were overweight or obese.

TABLE 1
www.frontiersin.org

TABLE 1. The identified rare variants with minor allele frequencies (MAF) in gnomAD and in silico prediction values.

FIGURE 1
www.frontiersin.org

FIGURE 1. Pedigrees of the families with missense variants in SORCS1. Square = male, circle = female, +/− heterozygote, −/− wildtype. Black square/circle = obesity, gray square/circle = overweight.

Table 2 presents the cytogenic location of the genes with the rare variants. Four of the genes (ALDOA, KIF22, SEZ6L2 and TAOK2) were in the 16p11.2 region and three of the genes (AIFM3, ARVCF, KLHL22) were in the 22q11.21 region. We further evaluated the gnomAD constrained metrics for the genes that were involved in the rare deletions that we previously identified. Table 3 presents the predicted constraint metrics. Three of the genes (BRAF, TAOK2 and XRN1) had a pLI = 1, indicating dosage sensitivity to loss-of-function variants.

TABLE 2
www.frontiersin.org

TABLE 2. The cytogenic location for the 12 genes with identified rare variants.

TABLE 3
www.frontiersin.org

TABLE 3. Predicted constrained metrics for the genes involved in deletions.

We also explored coding and splice variants with an allele frequency of <0.005. These variants were in nine genes; ARVCF, C16orf54, C22orf39, ESF1, GNB1L, KLHL22, MCTP2, MVP and TRMT2A (Supplementary Table S2). We did not find any potential pathogenic/likely pathogenic variants or candidate genes in this second analysis.

Furthermore, we evaluated the subjects in our study with the rare CNVs to look for additional rare variants in the specific genes corresponding to the regions of the deletion since a heterozygous deletion can also affect the function of a gene by unmasking a rare pathogenic variant on the other allele following recessive inheritance. No rare variants were found in these genes, thus excluding biallelic presence of damaging variants.

Discussion

CNVs are causing several diseases. However, in many cases it is challenging to determine the pathogenicity and clinical impact of the CNVs and to pinpoint the potential contribution of specific genes to the phenotype. In addition, the clinical relevance of rare CNVs can be difficult to determine because of incomplete penetrance and variable phenotypic expression. We performed targeted exome sequencing to identify rare variants in 117 genes located in previously identified CNV regions in subjects with severe early-onset obesity (Pettersson et al., 2017) in search for potential candidate genes. We found rare heterozygous variants with allele frequencies <0.1% in twelve of the CNV genes.

The 16p11.2 deletion is the most common CNV linked to severe obesity. However, the mechanism by which the proximal 16p11.2 deletion causes the phenotype is poorly understood and none of the genes in this region have so far been identified as a solely causative gene for the complex phenotype associated with proximal 16p11.2 CNVs. The patients with 16p11.2 deletions often present with various neurodevelopmental and neurobehavioral conditions, including autism spectrum disorders (ASD) and speech and language problems (Zufferey et al., 2012; Chung et al., 2021). The link between obesity and neurodevelopmental disorders may be due to common genes involved in the neurodevelopment that affect pathways regulating energy balance. Recently, new approaches for fine-mapping the genes in 16p11.2 region have been developed using GWAS and biobank data. Vysotskiy et al. performed an association analysis between imputed gene expression levels and BMI and they found four genes in the 16p11.2 region (TMEM219, SPN, TAOK2 and INO80E) associated with BMI (Vysotskiy et al., 2021). We identified rare heterozygous variants in four genes in the proximal 16p11.2 region; a frameshift variant in TAOK2 and heterozygous missense variants in SEZ6L2, ALDOA and KIF22. These four genes are potential candidate genes for the phenotype in patients with 16p11.2 CNVs. TAOK2 has a pLI of 1, suggesting rare loss-of-function variants may be deleterious (Lek et al., 2016). TAOK2 has been implicated in neurodevelopmental disorders and TAOK2 mutations found in patients with ASD showed impaired dendrite and synaptic development; the BMIs in the affected patients were not reported (Richter et al., 2019).TAOK2 heterozygous and knockout mice show gene brain morphological and behavioral abnormalities (Richter et al., 2019). SEZ6L2 has also been proposed as a candidate gene for ASD because of high expression in the brain (Kumar R. A. et al., 2009; Konyukh et al., 2011). Subcutaneous adipose tissue expression of ALDOA has been correlated with weight regain after dietary intervention (Roumans et al., 2017). ALDOA encodes an enzyme involved in glycolysis and the upregulation of glucose intake may be related to the risk of weight regain (Roumans et al., 2017). Furthermore, analysis of human fetal gene expression data found that KIF22 and ALDOA are significantly enriched in progenitors compared to post-mitotic cells making them candidates for having specific roles in neurogenesis in the developing human fetal cerebral cortex. It has been suggested that neurogenesis is disrupted in 16p11.2 CNVs (Morson et al., 2021). The participants in our study with rare variants in ALDOA, KIF22, TAOK2 and SEZ6L2 did not present with ASD or neurodevelopmental disorders (Table 1). However, the subjects with the missense variants in KIF22 and SEZ6L2 had learning difficulties. To our knowledge, no variants in these four genes have been reported in monogenic obesity. We consider the contribution of the identified rare variants to development of obesity unclear.

Two of the CNVs (16p11.2 deletion and 22q11.21 duplication) that we previously identified in our obese cohort include over 20 genes. In large CNVs, several genes within the locus may contribute to the phenotype. Disruption of multiple genes may be responsible for the complex phenotype seen in patients with 16p11.2 deletions, explaining the difficulty of identifying a single causative gene in the 16p11.2 region. Jensen et al. proposed that haploinsufficiency of certain genes may be altered by haploinsufficiency of other genes in the same CNV region and that these genes interact with each other via shared biological pathways (Jensen and Girirajan, 2019).

Interestingly, we identified heterozygous missense variants in SORCS1 in two subjects. Another study investigating genetic variants in 72 individuals with obesity also found two rare missense SORCS1 variants of unknown significance (Ginete et al., 2021). SORCS1 is highly expressed in the central nervous system and encodes a vacuolar protein sorting 10 domain-containing receptor that binds neuropeptides. SORSC1 has been associated with insulin secretion and diabetes risk (Goodarzi et al., 2007) and loss of both Sorcs1 and Sorcs3 in mutant mice caused increased adiposity and enhanced food intake, but not increased body weight (Subkhangulova et al., 2018). The participants in our study with the rare variants in SORSC1 did not present with insulin resistance or diabetes. In each family the variants were inherited from a parent with overweight, but overweight/obesity was present also in family members without the variant and therefore the contribution of these variants to the phenotype remains uncertain. More studies are needed to investigate the molecular mechanisms of SORCS1 and the possible contribution of this gene to early-onset obesity.

Another gene of interest is XRN1 located in 3q23 region. XRN1 is highly constrained for both loss-of-function and missense variants. This gene encodes 5′-3′ Exoribonuclease 1, which is involved in replication-dependent histone mRNA degradation. No XRN1 variants in humans have been associated with obesity, but a recent study investigating forebrain-specific Xrn1 knockout mice found that lack of Xrn1 in neurons led to obesity, hyperphagia, leptin resistance and hyperglycemia. The researchers proposed that Xrn1 in the hypothalamus is important for the mRNA degradation in regulating gene expression in the leptin signaling pathway (Takaoka et al., 2021). In our study, we identified a heterozygous missense variant in XRN1 of uncertain significance in a participant with isolated obesity. Further studies are required to explore the clinical relevance of this variant and biological function of XRN1.

Our study was limited by the lack of a control group of normal-weight subjects. However, we compared the allele frequencies of the variants we found in this study with the Finnish population in gnomAD and SISu project, with genotype data of more than 10,000 Finns. Because of the methodological limitation of targeted exome sequencing we may have missed variants in non-coding parts of the selected genes and variants in novel genes that were not included in the panel but may contribute to obesity. We may also have missed more common variants potentially contributing to the phenotype. Other limitations of this study were unavailability of parental DNA samples in part of the cohort and the lack of functional evaluation of the genetic variants. Further studies are needed to explore the possible role of the identified variants in the development of obesity.

In conclusion, the genes and mechanisms mediating the obesity phenotype in the previously identified CNVs remains uncertain and more studies are needed to provide insight to the potential contribution of the rare CNVs and the specific genes to severe obesity.

Data Availability Statement

The datasets presented in this article are not readily available because Data cannot be shared publicly because the data consists sensitive patient data. More specifically the data consists of individual clinical data and individual genotypes for young children. Data are available from the Helsinki University Hospital’s Institutional Data Access/Ethics Committee for researchers who meet the criteria for access to confidential data. Requests to access the datasets should be directed to Outi Mäkitie MD, Ph.D., outi.makitie@helsinki.fi.

Ethics Statement

The studies involving human participants were reviewed and approved by the Research Ethics Committees of Hospital District of Helsinki and Uusimaa, the Pirkanmaa Hospital District and the Northern Ostrobothnia Hospital District. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

PL, TM, HV, MP, and OM: study design. PL, TM, HV, and PT: data collection. PL, MP, and OM: data analysis and data interpretation. PL, MP, and OM: drafting of manuscript. All authors contributed to manuscript revision and approved the final version.

Funding

Academy of Finland; Sigrid Jusélius Foundation; Foundation for Pediatric Research; Folkhälsan Research Foundation; Päivikki ja Sakari Sohlberg Foundation; Stiftelsen Dorothea Olivia, Karl Walter och Jarl Walter Perkléns minne; Finska Läkaresällskapet; Swedish Research Council; Novo Nordisk Foundation; University of Helsinki through the Doctoral Program in Clinical Research; Helsinki University Hospital research funds.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank all the study subjects and their family members for participating in the study. We thank RN Päivi Turunen for help with data collection.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.839349/full#supplementary-material

References

Bochukova, E. G., Huang, N., Keogh, J., Henning, E., Purmann, C., Blaszczyk, K., et al. (2010). Large, Rare Chromosomal Deletions Associated with Severe Early-Onset Obesity. Nature 463 (7281), 666–670. doi:10.1038/nature08689

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, C. M. B., and Lupski, J. R. (2016). Mechanisms Underlying Structural Variant Formation in Genomic Disorders. Nat. Rev. Genet. 17 (4), 224–238. doi:10.1038/nrg.2015.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, W. K., Roberts, T. P., Sherr, E. H., Snyder, L. G., and Spiro, J. E. (2021). 16p11.2 Deletion Syndrome. Curr. Opin. Genet. Develop. 68, 49–56. doi:10.1016/j.gde.2021.01.011

CrossRef Full Text | Google Scholar

D’Angelo, C. S., Varela, M. C., de Castro, C. I. E., Otto, P. A., Perez, A. B. A., Lourenço, C. M., et al. (2018). Chromosomal Microarray Analysis in the Genetic Evaluation of 279 Patients with Syndromic Obesity. Mol. Cytogenet. 11, 14. doi:10.1186/s13039-018-0363-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Desvignes, J.-P., Bartoli, M., Delague, V., Krahn, M., Miltgen, M., Béroud, C., et al. (2018). VarAFT: a Variant Annotation and Filtration System for Human Next Generation Sequencing Data. Nucleic Acids Res. 46 (W1), W545–w553. doi:10.1093/nar/gky471

PubMed Abstract | CrossRef Full Text | Google Scholar

Ginete, C., Serrasqueiro, B., Silva-Nunes, J., Veiga, L., and Brito, M. (2021). Identification of Genetic Variants in 65 Obesity Related Genes in a Cohort of Portuguese Obese Individuals. Genes 12 (4), 603. doi:10.3390/genes12040603

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodarzi, M. O., Lehman, D. M., Taylor, K. D., Guo, X., Cui, J., Quiñones, M. J., et al. (2007). SORCS1: a Novel Human Type 2 Diabetes Susceptibility Gene Suggested by the Mouse. Diabetes 56 (7), 1922–1929. doi:10.2337/db06-1677

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, W., Zhang, F., and Lupski, J. R. (2008). Mechanisms for Human Genomic Rearrangements. Pathogenetics 1 (1), 4. doi:10.1186/1755-8417-1-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, N., Lee, I., Marcotte, E. M., and Hurles, M. E. (2010). Characterising and Predicting Haploinsufficiency in the Human Genome. Plos Genet. 6 (10), e1001154. doi:10.1371/journal.pgen.1001154

PubMed Abstract | CrossRef Full Text | Google Scholar

Ioannidis, N. M., Rothstein, J. H., Pejaver, V., Middha, S., McDonnell, S. K., Baheti, S., et al. (2016). REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants. Am. J. Hum. Genet. 99 (4), 877–885. doi:10.1016/j.ajhg.2016.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, M., and Girirajan, S. (2019). An Interaction-Based Model for Neuropsychiatric Features of Copy-Number Variants. Plos Genet. 15 (1), e1007879. doi:10.1371/journal.pgen.1007879

PubMed Abstract | CrossRef Full Text | Google Scholar

Konyukh, M., Delorme, R., Chaste, P., Leblond, C., Lemière, N., Nygren, G., et al. (2011). Variations of the Candidate SEZ6L2 Gene on Chromosome 16p11.2 in Patients with Autism Spectrum Disorders and in Human Populations. PLoS One 6 (3), e17289. doi:10.1371/journal.pone.0017289

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopanos, C., Tsiolkas, V., Kouris, A., Chapple, C. E., Albarca Aguilera, M., Meyer, R., et al. (2019). VarSome: the Human Genomic Variant Search Engine. Bioinformatics 35 (11), 1978–1980. doi:10.1093/bioinformatics/bty897

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, P., Henikoff, S., and Ng, P. C. (2009a). Predicting the Effects of Coding Non-synonymous Variants on Protein Function Using the SIFT Algorithm. Nat. Protoc. 4 (7), 1073–1081. doi:10.1038/nprot.2009.86

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, R. A., Marshall, C. R., Badner, J. A., Babatz, T. D., Mukamel, Z., Aldinger, K. A., et al. (2009b). Association and Mutation Analyses of 16p11.2 Autism Candidate Genes. PLoS One 4 (2), e4582. doi:10.1371/journal.pone.0004582

PubMed Abstract | CrossRef Full Text | Google Scholar

Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., et al. (2016). Analysis of Protein-Coding Genetic Variation in 60,706 Humans. Nature 536 (7616), 285–291. doi:10.1038/nature19057

PubMed Abstract | CrossRef Full Text | Google Scholar

Loid, P., Mustila, T., Mäkitie, R. E., Viljakainen, H., Kämpe, A., Tossavainen, P., et al. (2020). Rare Variants in Genes Linked to Appetite Control and Hypothalamic Development in Early-Onset Severe Obesity. Front. Endocrinol. 11, 81. doi:10.3389/fendo.2020.00081

PubMed Abstract | CrossRef Full Text | Google Scholar

Madeo, S. F., Stanghellini, I., Predieri, B., Ciancia, S., Leo, F., Bruzzi, P., et al. (2021). Copy Number Variation Analysis Increases the Number of Candidate Loci Associated with Pediatric Obesity. Horm. Res. Paediatr. 94, 251–262. doi:10.1159/000519299

PubMed Abstract | CrossRef Full Text | Google Scholar

Morson, S., Yang, Y., Price, D. J., and Pratt, T. (2021). Expression of Genes in the 16p11.2 Locus during Development of the Human Fetal Cerebral Cortex. Cereb. Cortex 31 (9), 4038–4052. doi:10.1093/cercor/bhab067

PubMed Abstract | CrossRef Full Text | Google Scholar

Nowakowska, B. (2017). Clinical Interpretation of Copy Number Variants in the Human Genome. J. Appl. Genet. 58 (4), 449–457. doi:10.1007/s13353-017-0407-4

CrossRef Full Text | Google Scholar

Pettersson, M., Viljakainen, H., Loid, P., Mustila, T., Pekkinen, M., Armenio, M., et al. (2017). Copy Number Variants Are Enriched in Individuals with Early-Onset Obesity and Highlight Novel Pathogenic Pathways. J. Clin. Endocrinol. Metab. 102 (8), 3029–3039. doi:10.1210/jc.2017-00565

PubMed Abstract | CrossRef Full Text | Google Scholar

Rice, A. M., and McLysaght, A. (2017). Dosage Sensitivity Is a Major Determinant of Human Copy Number Variant Pathogenicity. Nat. Commun. 8, 14366. doi:10.1038/ncomms14366

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and Guidelines for the Interpretation of Sequence Variants: a Joint Consensus Recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 17 (5), 405–424. doi:10.1038/gim.2015.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, M., Murtaza, N., Scharrenberg, R., White, S. H., Johanns, O., Walker, S., et al. (2019). Altered TAOK2 Activity Causes Autism-Related Neurodevelopmental and Cognitive Abnormalities through RhoA Signaling. Mol. Psychiatry 24 (9), 1329–1350. doi:10.1038/s41380-018-0025-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Roumans, N. J. T., Vink, R. G., Bouwman, F. G., Fazelzadeh, P., van Baak, M. A., and Mariman, E. C. M. (2017). Weight Loss-Induced Cellular Stress in Subcutaneous Adipose Tissue and the Risk for Weight Regain in Overweight and Obese Adults. Int. J. Obes. 41 (6), 894–901. doi:10.1038/ijo.2016.221

CrossRef Full Text | Google Scholar

Saari, A., Sankilampi, U., Hannila, M.-L., Kiviniemi, V., Kesseli, K., and Dunkel, L. (2011). New Finnish Growth References for Children and Adolescents Aged 0 to 20 years: Length/height-For-Age, Weight-For-Length/height, and Body Mass index-for-age. Ann. Med. 43 (3), 235–248. doi:10.3109/07853890.2010.515603

PubMed Abstract | CrossRef Full Text | Google Scholar

Saeed, S., Arslan, M., Manzoor, J., Din, S. M., Janjua, Q. M., Ayesha, H., et al. (2020). Genetic Causes of Severe Childhood Obesity: A Remarkably High Prevalence in an Inbred Population of Pakistan. Diabetes 69 (7), 1424–1438. doi:10.2337/db19-1238

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwarz, J. M., Cooper, D. N., Schuelke, M., and Seelow, D. (2014). MutationTaster2: Mutation Prediction for the Deep-Sequencing Age. Nat. Methods 11 (4), 361–362. doi:10.1038/nmeth.2890

PubMed Abstract | CrossRef Full Text | Google Scholar

Selvanayagam, T., Walker, S., Gazzellone, M. J., Kellam, B., Cytrynbaum, C., Stavropoulos, D. J., et al. (2018). Genome-wide Copy Number Variation Analysis Identifies Novel Candidate Loci Associated with Pediatric Obesity. Eur. J. Hum. Genet. 26 (11), 1588–1596. doi:10.1038/s41431-018-0189-0

CrossRef Full Text | Google Scholar

Subkhangulova, A., Malik, A. R., Hermey, G., Popp, O., Dittmar, G., Rathjen, T., et al. (2018). SORCS 1 and SORCS 3 Control Energy Balance and Orexigenic Peptide Production. EMBO Rep. 19 (4), e44810. doi:10.15252/embr.201744810

PubMed Abstract | CrossRef Full Text | Google Scholar

Takaoka, S., Yanagiya, A., Mohamed, H. M. A., Higa, R., Abe, T., Inoue, K. I., et al. (2021). Neuronal XRN1 is Required for Maintenance of Whole-Body Metabolic Homeostasis. iScience 24 (10), 103151. doi:10.1016/j.isci.2021.103151

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuillaume, M.-L., Naudion, S., Banneau, G., Diene, G., Cartault, A., Cailley, D., et al. (2014). New Candidate Loci Identified by Array-CGH in a Cohort of 100 Children Presenting with Syndromic Obesity. Am. J. Med. Genet. 164 (8), 1965–1975. doi:10.1002/ajmg.a.36587

CrossRef Full Text | Google Scholar

Vysotskiy, M., Zhong, X., Miller-Fleming, T. W., Zhou, D., Cox, N. J., and Weiss, L. A. (2021). Integration of Genetic, Transcriptomic, and Clinical Data Provides Insight into 16p11.2 and 22q11.2 CNV Genes. Genome Med. 13 (1), 172. doi:10.1186/s13073-021-00972-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Walters, R. G., Jacquemont, S., Valsesia, A., de Smith, A. J., Martinet, D., Andersson, J., et al. (2010). A New Highly Penetrant Form of Obesity Due to Deletions on Chromosome 16p11.2. Nature 463 (7281), 671–675. doi:10.1038/nature08727

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Li, W.-D., Glessner, J. T., Grant, S. F. A., Hakonarson, H., and Price, R. A. (2010). Large Copy-Number Variations Are Enriched in Cases with Moderate to Extreme Obesity. Diabetes 59 (10), 2690–2694. doi:10.2337/db10-0192

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheeler, E., Huang, N., Bochukova, E. G., Keogh, J. M., Lindsay, S., Garg, S., et al. (2013). Genome-wide SNP and CNV Analysis Identifies Common and Low-Frequency Variants Associated with Severe Early-Onset Obesity. Nat. Genet. 45 (5), 513–517. doi:10.1038/ng.2607

PubMed Abstract | CrossRef Full Text | Google Scholar

Zarrei, M., MacDonald, J. R., Merico, D., and Scherer, S. W. (2015). A Copy Number Variation Map of the Human Genome. Nat. Rev. Genet. 16 (3), 172–183. doi:10.1038/nrg3871

PubMed Abstract | CrossRef Full Text | Google Scholar

Zufferey, F., Sherr, E. H., Beckmann, N. D., Hanson, E., Maillard, A. M., Hippolyte, L., et al. (2012). A 600 Kb Deletion Syndrome at 16p11.2 Leads to Energy Imbalance and Neuropsychiatric Disorders. J. Med. Genet. 49 (10), 660–668. doi:10.1136/jmedgenet-2012-101203

CrossRef Full Text | Google Scholar

Keywords: childhood obesity, gene panel, rare variants, candidate genes, copy number variants

Citation: Loid P, Pekkinen M, Mustila T, Tossavainen P, Viljakainen H, Lindstrand A and Mäkitie O (2022) Targeted Exome Sequencing of Genes Involved in Rare CNVs in Early-Onset Severe Obesity. Front. Genet. 13:839349. doi: 10.3389/fgene.2022.839349

Received: 19 December 2021; Accepted: 07 February 2022;
Published: 07 March 2022.

Edited by:

Tommaso Pippucci, Policlinico Sant’Orsola-Malpighi, Italy

Reviewed by:

Olivia J. Veatch, University of Kansas Medical Center, United States
Edoardo Errichiello, University of Pavia, Italy

Copyright © 2022 Loid, Pekkinen, Mustila, Tossavainen, Viljakainen, Lindstrand and Mäkitie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Petra Loid, petra.loid@helsinki.fi

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.