Çiğdem Köroğlu

Objective: Rare cases of monogenic obesity, which may respond to specific therapeutics, can remain undetected in populations in which polygenic obesity is prevalent. This study examined rare DNA variation in established monogenic obesity genes within a community using whole‐exome sequence data from 6803 longitudinally studied individuals.Methods: Exome data across 15 monogenic obesity genes were analyzed for nonsynonymous variants observed in any child with a maximum BMI z score > 2 (N = 279) but not observed in a child with a maximum BMI z score ≤ 0 (n = 1542) or that occurred in adults in the top 5th percentile of BMI (n = 263) but not in adults below the median BMI (n = 2629). Variants were then functionally analyzed using luciferase assays.Results: The comparisons between cases of obesity and controls identified eight missense variants in six genes: DYRK1B, KSR2, MC4R, NTRK2, PCSK1, and SIM1. Among these, MC4R p.A303P and p.R165G were previously shown to impair MC4R function. Functional analyses of the remaining six variants suggest that KSR2 p.I402F and p.T193I and NTRK2 p.S249Y alter protein function.Conclusions: In addition to MC4R, rare missense variants in KSR2 and NTRK2 may potentially explain the severe obesity observed for the carriers.

There is a collective push to diversify human genetic studies by including underrepresented populations. However, analyzing DNA sequence reads involves the initial step of aligning the reads to the GRCh38/hg38 reference genome which is inadequate for non-European ancestries. In this study, using long-read sequencing technology, we constructed de novo genome assemblies from two indigenous Americans from Arizona (IAZ). Each assembly included ∼17 Mb of DNA sequence not present [nonreference sequence (NRS)] in hg38, which consists mostly of repeat elements. Forty NRSs totaling 240 kb were uniquely anchored to the hg38 primary assembly generating a modified hg38-NRS reference genome. DNA sequence alignment and variant calling were then conducted with whole-genome sequencing (WGS) sequencing data from 387 IAZ using both the hg38 and modified hg38-NRS reference maps. Variant calling with the hg38-NRS map identified ∼50,000 single-nucleotide variants present in at least 5% of the WGS samples which were not detected with the hg38 reference map. We also directly assessed the NRSs positioned within genes. Seventeen NRSs anchored to regions including an identical 187 bp NRS found in both de novo assemblies. The NRS is located in HCN2 79 bp downstream of Exon 3 and contains several putative transcriptional regulatory elements. Genotyping of the HCN2-NRS revealed that the insertion is enriched in IAZ (minor allele frequency = 0.45) compared to other reference populations tested. This study shows that inclusion of population-specific NRSs can dramatically change the variant profile in an underrepresented ethnic groups and thereby lead to the discovery of previously missed common variations.

There is a collective push to diversify human genetic studies by including underrepresented populations. However, analyzing DNA sequence reads involves the initial step of aligning the reads to the GRCh38/hg38 reference genome which is inadequate for non-European ancestries. To help address this issue, we created a modified hg38 reference map using de novo sequence assemblies from Indigenous Americans living in Arizona (IAZ). Using HiFi SMRT long-read sequencing technology, we generated de novo genome assemblies for one female and one male IAZ individual. Each assembly included ∼17 Mb of DNA sequence not present (non-reference sequence; NRS) in hg38, which consists mostly of repeat elements. Forty NRSs totaling 240 kb were uniquely anchored to the hg38 primary assembly generating a modified hg38-NRS reference genome. DNA sequence alignment and variant calling were then conducted with WGS sequencing data from 387 IAZ cohorts using both the hg38 and modified hg38-NRS reference maps. Variant calling with the hg38-NRS map identified ∼50,000 single nucleotide variants present in at least 5% of the WGS samples which were not detected with the hg38 reference map. We also directly assessed the NRSs positioned within genes. Seventeen NRSs anchored to regions including an identical 187 bp NRS found in both de novo assemblies. The NRS is located in HCN2 79 bp downstream of exon 3 and contains several putative transcriptional regulatory elements. Genotyping of the HCN2 -NRS revealed that the insertion is enriched in IAZ (MAF = 0.45) compared to Caucasians (MAF = 0.15) and African Americans (MAF = 0.03). This study shows that inclusion of population-specific NRSs can dramatically change the variant profile in an under-represented ethnic groups and thereby lead to the discovery of previously missed common variations. AUTHOR SUMMARY: GRCh38/hg38 reference genome has been the standard reference for large-scale human genetics studies. However, it does not adequately represent sequences of non-European ancestry. In this study, using long-read sequencing technology, we constructed de novo sequence assemblies from two Indigenous Americans from Arizona. We then compared the de novo assemblies to the hg38 reference genome to identify non-reference sequences (NRSs). We integrated these NRSs into our whole-genome sequencing (WGS) variant calling pipeline to improve read alignment and variant detection. We also directly assessed the NRSs positioned within genes. Inclusion of population-specific NRSs dramatically changed the variant profile of our study group with under-represented ethnicity, revealing common variation not detected by our previous population-level WGS and genotyping studies.

Objective: This study aimed to identify genetic variants enriched in Southwest American Indian (SWAI) individuals that associate with BMI.Methods: Whole genome sequencing data (n = 296) were used to identify potentially functional variants that are common in SWAI individuals (minor allele frequency ≥10%) but rare in other ethnic groups (minor allele frequency < 0.1%). Enriched variants were tested for association with BMI in 5,870 SWAI individuals. One variant was studied using a luciferase reporter, and haplotypes that included this variant were analyzed for association with various measures of obesity (n = 917‐5,870), 24‐hour energy expenditure (24‐h EE; n = 419), and skeletal muscle biopsy expression data (n = 207).Results: A 5′ untranslated region variant in cytochrome b5 type A (CYB5A), rs548402150, met the enrichment criteria and associated with increased BMI (β = 2%, p = 0.004). Functionally, rs548402150 decreased luciferase expression by 30% (p = 0.003) and correlated with decreased skeletal muscle CYB5A expression (β = −0.5 SD, p = 0.0008). Combining rs548402150 with two splicing quantitative trait loci in CYB5A identified a haplotype carried almost exclusively in SWAI individuals that associated with increased BMI (β = 3%, p = 0.0003) and decreased CYB5A expression, whereas the most common haplotype in all ethnic groups associated with lower BMI and percentage of body fatness, increased 24‐h EE, and increased CYB5A expression.Conclusions: Further studies on the effects of CYB5A on 24‐h EE and BMI may provide insights into obesity‐related physiology.

Objective: In an ongoing effort to identify the genetic variation that contributes to obesity in American Indians, known Bardet–Biedl syndrome (BBS) genes were analyzed for an effect on BMI and leptin signaling.Methods: Potentially deleterious variants (Combined Annotation Dependent Depletion score > 20) in BBS genes were identified in whole‐exome sequence data from 6,851 American Indians informative for BMI. Common variants (detected in ≥ 10 individuals) were analyzed for association with BMI; rare variants (detected in < 10 individuals) were analyzed for mean BMI of carriers. Functional assessment of variants’ effect on signal transducer and activator of transcription 3 (STAT3) activity was performed in vitro.Results: One common variant, rs59252892 (Thr549Ile) in BBS9, was associated with BMI (P = 0.0008, β = 25% increase per risk allele). Among rare variants for which carriers had severe obesity (mean BMI > 40 kg/m²), four were in BBS9. In vitro analysis of BBS9 found the Ile allele at Thr549Ile had a 20% increase in STAT3 activity compared with the Thr allele (P = 0.01). Western blot analysis showed the Ile allele had a 15% increase in STAT3 phosphorylation (P = 0.006). Comparable functional results were observed with Ser545Gly and Val209Leu but not Leu665Phe and Lys810Glu.Conclusions: Potentially functional variants in BBS genes in American Indians are reported. However, functional evidence supporting a causal role for BBS9 in obesity is inconclusive.

Background/objectives Nighttime eating (NE) behavior has a genetic component and predicts weight gain. We hypothesized that some genetic variants, which affect NE would also show an effect on body mass index (BMI). We aimed to determine which known BMI variants associate with NE in Southwestern American Indians (SWAIs), who are at elevated risk for obesity. Methods Known BMI variants from the GIANT-UK Biobank meta-analysis (N = 700,000) were analysed in SWAIs characterized for NE during an inpatient 3-day protocol. Variants were analysed for association with NE using whole-genome sequence data from 50 SWAIs (23 cases and 27 controls) and selected variants were genotyped in an additional 32 SWAIs (13 NE cases and 19 controls). Variants associated with NE in a meta-analysis of the two SWAI samples were further analysed for association with nightly caloric intake and functionality in hypothalamus, pituitary, and adrenal tissues. Results Variants were identified where the allele that associated with increased BMI in the GIANT-UK Biobank meta-analysis (P ≤ 1 × 10⁻⁸) also had a P < 0.01 for increased NE in the SWAI meta-analysis. These variants were captured by six tagSNPs. Comparison of the nightly calorie intake by genotype and eQTL data from relevant tissues highlighted rs3753612 upstream of HCRTR1. Conclusions Our strategy led to the HCRTR1 locus, which has previously been linked to sleep regulation and feeding. Although this is an intriguing candidate gene for NE, further studies in larger samples and different populations are required to validate the role of HCRTR1 in NE.

Amyotrophic lateral sclerosis (ALS) is a motor neuron disease eventually leading to death from respiratory failure. Recessive inheritance is very rare. Here, we describe the clinical findings in a consanguineous family with five men afflicted with recessive ALS and the identification of the homozygous mutation responsible for the disorder. The onset of the disease ranged from 12 to 35 years of age, with variable disease progressions. We performed clinical investigations including metabolic and paraneoplastic screening, cranial and cervical imaging, and electrophysiology. We mapped the disease gene to 9p21.1-p12 with a LOD score of 5.2 via linkage mapping using genotype data for single-nucleotide polymorphism markers and performed exome sequence analysis to identify the disease-causing gene variant. We also Sanger sequenced all coding sequences of SIGMAR1, a gene reported as responsible for juvenile ALS in a family. We did not find any mutation in SIGMAR1. Instead, we identified a novel homozygous missense mutation p.(His705Arg) in GNE which was predicted as damaging by online tools. GNE has been associated with inclusion body myopathy and is expressed in many tissues. We propose that the GNE mutation underlies the pathology in the family.