Variant-to-gene mapping followed by cross-species genetic screening identifies GPI-anchor biosynthesis as a regulator of sleep

Circadian clock gene polymorphisms implicated in human pathologies

Circadian rhythms, ~24 h cycles of physiological and behavioral processes, can be synchronized by external signals (e.g., light) and persist even in their absence. Consequently, dysregulation of circadian rhythms adversely affects the well-being of the organism. This timekeeping system is generated and sustained by a genetically encoded endogenous mechanism composed of interlocking transcriptional/translational feedback loops that generate rhythmic expression of core clock genes. Genome-wide association studies (GWAS) and forward genetic studies show that SNPs in clock genes influence gene regulation and correlate with the risk of developing various conditions. We discuss genetic variations in core clock genes that are associated with various phenotypes, their implications for human health, and stress the need for thorough studies in this domain of circadian regulation.

Cross-species evidence for a developmental origin of adult hypersomnia with loss of synaptic adhesion molecules beat-Ia/CADM2

Idiopathic hypersomnia (IH) is a poorly-understood sleep disorder characterized by excessive daytime sleepiness despite normal nighttime sleep. Combining human genomics with behavioral and mechanistic studies in fish and flies, we uncover a role for beat-Ia/CADM2 , synaptic adhesion molecules of the immunoglobulin superfamily, in excessive sleepiness. Neuronal knockdown of Drosophila beat-Ia results in sleepy flies and loss of the vertebrate ortholog of beat-Ia , CADM2 , results in sleepy fish. We delineate a developmental function for beat-Ia in synaptic elaboration of neuropeptide F (NPF) neurites projecting to the suboesophageal zone (SEZ) of the fly brain. Brain connectome and experimental evidence demonstrate these NPF outputs synapse onto a subpopulation of SEZ GABAergic neurons to stabilize arousal. NPF is the Drosophila homolog of vertebrate neuropeptide Y (NPY), and an NPY receptor agonist restores sleep to normal levels in zebrafish lacking CADM2 . These findings point towards NPY modulation as a treatment target for human hypersomnia.

Characterization of non-coding variants associated with transcription factor binding through ATAC-seq-defined footprint QTLs in liver

Non-coding variants discovered by genome-wide association studies (GWAS) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint". Here, we sought to identify variants associated with TF-binding, or "footprint quantitative trait loci" (fpQTLs) in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 693 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 80% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding site disruption in disease and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

Implicating type 2 diabetes effector genes in relevant metabolic cellular models using promoter-focused Capture-C

AIMS/HYPOTHESIS

Genome-wide association studies (GWAS) have identified hundreds of type 2 diabetes loci, with the vast majority of signals located in non-coding regions; as a consequence, it remains largely unclear which 'effector' genes these variants influence. Determining these effector genes has been hampered by the relatively challenging cellular settings in which they are hypothesised to confer their effects.

METHODS

To implicate such effector genes, we elected to generate and integrate high-resolution promoter-focused Capture-C, assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq datasets to characterise chromatin and expression profiles in multiple cell lines relevant to type 2 diabetes for subsequent functional follow-up analyses: EndoC-BH1 (pancreatic beta cell), HepG2 (hepatocyte) and Simpson-Golabi-Behmel syndrome (SGBS; adipocyte).

RESULTS

The subsequent variant-to-gene analysis implicated 810 candidate effector genes at 370 type 2 diabetes risk loci. Using partitioned linkage disequilibrium score regression, we observed enrichment for type 2 diabetes and fasting glucose GWAS loci in promoter-connected putative cis-regulatory elements in EndoC-BH1 cells as well as fasting insulin GWAS loci in SGBS cells. Moreover, as a proof of principle, when we knocked down expression of the SMCO4 gene in EndoC-BH1 cells, we observed a statistically significant increase in insulin secretion.

CONCLUSIONS/INTERPRETATION

These results provide a resource for comparing tissue-specific data in tractable cellular models as opposed to relatively challenging primary cell settings.

DATA AVAILABILITY

Raw and processed next-generation sequencing data for EndoC-BH1, HepG2, SGBS_undiff and SGBS_diff cells are deposited in GEO under the Superseries accession GSE262484. Promoter-focused Capture-C data are deposited under accession GSE262496. Hi-C data are deposited under accession GSE262481. Bulk ATAC-seq data are deposited under accession GSE262479. Bulk RNA-seq data are deposited under accession GSE262480.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

3D chromatin-based variant-to-gene maps across 57 human cell types reveal the cellular and genetic architecture of autoimmune disease susceptibility

A portion of the genetic basis for many common autoimmune disorders has been uncovered by genome-wide association studies (GWAS), but GWAS do not reveal causal variants, effector genes, or the cell types impacted by disease-associated variation. We have generated 3D genomic datasets consisting of promoter-focused Capture-C, Hi-C, ATAC-seq, and RNA-seq and integrated these data with GWAS of 16 autoimmune traits to physically map disease-associated variants to the effector genes they likely regulate in 57 human cell types. These 3D maps of gene cis-regulatory architecture are highly powered to identify the cell types most likely impacted by disease-associated genetic variation compared to 1D genomic features, and tend to implicate different effector genes than eQTL approaches in the same cell types. Most of the variants implicated by these cis-regulatory architectures are highly trait-specific, but nearly half of the target genes connected to these variants are shared across multiple autoimmune disorders in multiple cell types, suggesting a high level of genetic diversity and complexity among autoimmune diseases that nonetheless converge at the level of target gene and cell type. Substantial effector gene sharing led to the common enrichment of similar biological networks across disease and cell types. However, trait-specific pathways representing potential areas for disease-specific intervention were identified. To test this, we pharmacologically validated squalene synthase, a cholesterol biosynthetic enzyme encoded by the FDFT1 gene implicated by our approach in MS and SLE, as a novel immunomodulatory drug target controlling inflammatory cytokine production by human T cells. These data represent a comprehensive resource for basic discovery of gene cis-regulatory mechanisms, and the analyses reported reveal mechanisms by which autoimmune-associated variants act to regulate gene expression, function, and pathology across multiple, distinct tissues and cell types.

The Genetics of Human Sleep and Sleep Disorders

Healthy sleep is vital for humans to achieve optimal health and longevity. Poor sleep and sleep disorders are strongly associated with increased morbidity and mortality. However, the importance of good sleep continues to be underrecognized. Mechanisms regulating sleep and its functions in humans remain mostly unclear even after decades of dedicated research. Advancements in gene sequencing techniques and computational methodologies have paved the way for various genetic analysis approaches, which have provided some insights into human sleep genetics. This review summarizes our current knowledge of the genetic basis underlying human sleep traits and sleep disorders. We also highlight the use of animal models to validate genetic findings from human sleep studies and discuss potential molecular mechanisms and signaling pathways involved in the regulation of human sleep.

Perturbation of the insomnia WDR90 genome-wide association studies locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome-wide association studies (GWAS) have identified loci for sleep-related traits, they do not directly uncover the underlying causal variants and corresponding effector genes. The majority of such variants reside in non-coding regions and are therefore presumed to impact cis-regulatory elements. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated phosphatidyl inositol glycan (PIG)-Q as a functionally relevant gene at the insomnia "WDR90" GWAS locus. However, importantly that effort did not characterize the corresponding underlying causal variant. Specifically, our previous 3D genomic datasets nominated a shortlist of three neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium within an intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. We sought to investigate the influence of these SNPs collectively and then individually on PIG-Q modulation to pinpoint the causal "regulatory" variant. Starting with gross level perturbation, deletion of the entire region in NPCs via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from individual luciferase reporter assays for each SNP in iPSCs revealed that the region with the rs3752495 risk allele (RA) induced a ~2.5-fold increase in luciferase expression. Importantly, rs3752495 also exhibited an allele-specific effect, with the RA increasing the luciferase expression by ~2-fold versus the non-RA. In conclusion, our variant-to-function approach and in vitro validation implicate rs3752495 as a causal insomnia variant embedded within WDR90 while modulating the expression of the distally located PIG-Q.

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

A conserved role for frizzled in sleep architecture

Previous studies of natural variants in Drosophila melanogaster implicated the Wnt signaling receptor frizzled in sleep. Given that the Wnt signaling pathway is highly conserved across species, we hypothesized that frizzled class receptor 1 (Fzd1), the murine homolog of frizzled, would also have a role in sleep. Using a CRISPR transgenic approach, we removed most of the Fzd1 coding region from C57BL/6N mice. We used a video assay to measure sleep characteristics in Fzd1-deficient mice. As Wnt signaling is known to affect visuospatial memory, we also examined the impact of the deletion on learning and memory using the novel object recognition (NOR) paradigm. Fzd1-deficient mice had altered sleep compared to littermate controls. The mice did not respond differently to the NOR paradigm compared to controls but did display anxiety-like behavior. Our strategy demonstrates that the study of natural variation in Drosophila sleep translates into candidate genes for sleep in vertebrate species such as the mouse.

Integrative Genomic Analyses Identify LncRNA Regulatory Networks across Pediatric Leukemias and Solid Tumors

Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies.

SIGNIFICANCE

Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks.

Perturbation of the insomnia WDR90 GWAS locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome wide association studies (GWAS) have been crucial for the identification of loci associated with sleep traits and disorders, the method itself does not directly uncover the underlying causal variants and corresponding effector genes. The overwhelming majority of such variants reside in non-coding regions and are therefore presumed to impact the activity of cis-regulatory elements, such as enhancers. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated PIG-Q as a functionally relevant gene at the insomnia 'WDR90' locus. However, importantly that effort did not characterize the corresponding underlying causal variant at this GWAS signal. Specifically, our genome-wide ATAC-seq and high-resolution promoter-focused Capture C datasets generated in this cell setting brought our attention to a shortlist of three tightly neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium in a candidate intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. The objective of this study was to investigate the influence of the proxy SNPs collectively and then individually on PIG-Q modulation and to pinpoint the causal "regulatory" variant among the three SNPs. Starting at a gross level perturbation, deletion of the entire region harboring all three SNPs in human iPSC-derived neural progenitor cells via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from more refined individual luciferase reporter assays for each of the three SNPs in iPSCs revealed that the intronic region with the rs3752495 risk allele induced a ~2.5-fold increase in luciferase expression (n=10). Importantly, rs3752495 also exhibited an allele specific effect, with the risk allele increasing the luciferase expression by ~2-fold compared to the non-risk allele. In conclusion, our variant-to-function approach and subsequent in vitro validation implicates rs3752495 as a causal insomnia risk variant embedded at the WDR90-PIG-Q locus.

Biological rhythms and task allocation in ant colonies

Task allocation in ant colonies, mediated by social interactions, regulates which individuals perform each task, and when they are active, in response to the current situation. Many tasks are performed in a daily temporal pattern. An ant's biological clock depends on patterns of gene expression that are regulated using a negative feedback loop which is synchronized to earth's rotation by external cues. An individual's biological clock can shift in response to social cues, and this plasticity contributes to task switching. Daily rhythms in individual ant behavior combine, via interactions within and across task groups, to adjust the collective behavior of colonies. Further work is needed to elucidate how the social cues that lead to task switching influence the molecular mechanisms that generate clock outputs associated with each task, and to investigate the evolution of temporal patterns in task allocation in relation to ecological factors.

Pleiotropic effects of a high confidence Autism Spectrum Disorder gene, arid1b, on zebrafish sleep

Sleep fulfills critical functions in neurodevelopment, such as promoting synaptic plasticity, neuronal wiring, and brain connectivity which are critical phenomena in Autism Spectrum Disorder (ASD) pathophysiology. Sleep disturbance, specifically insomnia, accompanies ASD and is associated with more severe core symptoms (e.g., social impairment). It is possible that focusing on identifying effective ways to treat sleep problems can help alleviate other ASD-related symptoms. A body of evidence indicates shared mechanisms and neurobiological substrates between sleep and ASD and investigation of these may inform therapeutic effects of improving sleep at both behavioral and molecular levels. In this study, we tested if sleep and social behavior were different in a zebrafish model with the arid1b gene mutated compared to controls. This gene was selected for study as expert curations conducted for the Simons Foundation for Autism Research Institute (SFARI) Gene database define it is as a 'high confidence' ASD gene (i.e., clearly implicated) encoding a chromatin remodeling protein. Homozygous arid1b mutants displayed increased arousability and light sleep compared to their heterozygous and wild type counterparts, based on testing a mechano-acoustic stimulus presenting different vibration frequencies of increasing intensity to detect sleep depth. In addition, decreased social preference was observed in arid1b heterozygous and homozygous mutant zebrafish. The behavioral phenotypes reported in our study are in line with findings from mouse models and human studies and demonstrate the utility of zebrafish as a vertebrate model system with high throughput phenotyping in the investigation of changes in sleep in models relevant to ASD. Furthermore, we demonstrate the importance of including assessments of arousal threshold when studying sleep using in vivo models.

Thalamic neurons drive distinct forms of motor asymmetry that are conserved in teleost and dependent on visual evolution

Brain laterality is a prominent feature in Bilateria, where neural functions are favored in a single brain hemisphere. These hemispheric specializations are thought to improve behavioral performance and are commonly observed as sensory or motor asymmetries, such as handedness in humans. Despite its prevalence, our understanding of the neural and molecular substrates instructing functional lateralization is limited. Moreover, how functional lateralization is selected for or modulated throughout evolution is poorly understood. While comparative approaches offer a powerful tool for addressing this question, a major obstacle has been the lack of a conserved asymmetric behavior in genetically tractable organisms. Previously, we described a robust motor asymmetry in larval zebrafish. Following the loss of illumination, individuals show a persistent turning bias that is associated with search pattern behavior with underlying functional lateralization in the thalamus. This behavior permits a simple yet robust assay that can be used to address fundamental principles underlying lateralization in the brain across taxa. Here, we take a comparative approach and show that motor asymmetry is conserved across diverse larval teleost species, which have diverged over the past 200 million years. Using a combination of transgenic tools, ablation, and enucleation, we show that teleosts exhibit two distinct forms of motor asymmetry, vision-dependent and - independent. These asymmetries are directionally uncorrelated, yet dependent on the same subset of thalamic neurons. Lastly, we leverage Astyanax sighted and blind morphs, which show that fish with evolutionarily derived blindness lack both retinal-dependent and -independent motor asymmetries, while their sighted surface conspecifics retained both forms. Our data implicate that overlapping sensory systems and neuronal substrates drive functional lateralization in a vertebrate brain that are likely targets for selective modulation during evolution.

Variant-to-gene, gene-to-trait: finding new sleep genes through GWAS

Distilling insomnia genome-wide association study (GWAS) variants, Palermo and colleagues identified several genes that participate in sleep regulation in two different model organisms. This workflow sets off an innovative strategy to extract biological relevance from large human genomic databases.