ABSTRACT
Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments (P<0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders.
ABSTRACT
BACKGROUND: Autoimmune cytopenias (AICs) regularly occur in profoundly IgG-deficient patients with common variable immunodeficiency (CVID). The isotypes, antigenic targets, and origin(s) of their disease-causing autoantibodies are unclear. OBJECTIVE: We sought to determine reactivity, clonality, and provenance of AIC-associated IgM autoantibodies in patients with CVID. METHODS: We used glycan arrays, patient erythrocytes, and platelets to determine targets of CVID IgM autoantibodies. Glycan-binding profiles were used to identify autoreactive clones across B-cell subsets, specifically circulating marginal zone (MZ) B cells, for sorting and IGH sequencing. The locations, transcriptomes, and responses of tonsillar MZ B cells to different TH- cell subsets were determined by confocal microscopy, RNA-sequencing, and cocultures, respectively. RESULTS: Autoreactive IgM coated erythrocytes and platelets from many CVID patients with AICs (CVID+AIC). On glycan arrays, CVID+AIC plasma IgM narrowly recognized erythrocytic i antigens and platelet i-related antigens and failed to bind hundreds of pathogen- and tumor-associated carbohydrates. Polyclonal i antigen-recognizing B-cell receptors were highly enriched among CVID+AIC circulating MZ B cells. Within tonsillar tissues, MZ B cells secreted copious IgM when activated by the combination of IL-10 and IL-21 or when cultured with IL-10/IL-21-secreting FOXP3-CD25hi T follicular helper (Tfh) cells. In lymph nodes from immunocompetent controls, MZ B cells, plentiful FOXP3+ regulatory T cells, and rare FOXP3-CD25+ cells that represented likely CD25hi Tfh cells all localized outside of germinal centers. In CVID+AIC lymph nodes, cellular positions were similar but CD25hi Tfh cells greatly outnumbered regulatory cells. CONCLUSIONS: Our findings indicate that glycan-reactive IgM autoantibodies produced outside of germinal centers may contribute to the autoimmune pathogenesis of CVID.
Subject(s)
Autoantibodies , Blood Platelets , Common Variable Immunodeficiency , Erythrocytes , Immunoglobulin M , Polysaccharides , Humans , Immunoglobulin M/immunology , Immunoglobulin M/blood , Erythrocytes/immunology , Common Variable Immunodeficiency/immunology , Polysaccharides/immunology , Blood Platelets/immunology , Autoantibodies/immunology , Autoantibodies/blood , Male , Female , B-Lymphocyte Subsets/immunology , AdultABSTRACT
The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.
Subject(s)
3' Untranslated Regions , Apoptosis Regulatory Proteins , Membrane Proteins , Pediatric Obesity , Child , Humans , 3' Untranslated Regions/genetics , Alleles , Cell Differentiation/genetics , Chromosomes, Human, Pair 12/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study , Human Embryonic Stem Cells/metabolism , Neurons/metabolism , Pediatric Obesity/genetics , Polymorphism, Single Nucleotide/genetics , Membrane Proteins/genetics , Apoptosis Regulatory Proteins/geneticsABSTRACT
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.
Subject(s)
Genome-Wide Association Study , Polymorphism, Single Nucleotide , Sleep Initiation and Maintenance Disorders , Polymorphism, Single Nucleotide/genetics , Humans , Sleep Initiation and Maintenance Disorders/genetics , Induced Pluripotent Stem Cells , Animals , Neural Stem Cells , Linkage Disequilibrium/geneticsABSTRACT
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.
ABSTRACT
Genome wide association study (GWAS)-implicated bone mineral density (BMD) signals have been shown to localize in cis-regulatory regions of distal effector genes using 3D genomic methods. Detailed characterization of such genes can reveal novel causal genes for BMD determination. Here, we elected to characterize the "DNM3" locus on chr1q24, where the long non-coding RNA DNM3OS and the embedded microRNA MIR199A2 (miR-199a-5p) are implicated as effector genes contacted by the region harboring variation in linkage disequilibrium with BMD-associated sentinel single nucleotide polymorphism, rs12041600. During osteoblast differentiation of human mesenchymal stem/progenitor cells (hMSC), miR-199a-5p expression was temporally decreased and correlated with the induction of osteoblastic transcription factors RUNX2 and Osterix. Functional relevance of miR-199a-5p downregulation in osteoblastogenesis was investigated by introducing miR-199a-5p mimic into hMSC. Cells overexpressing miR-199a-5p depicted a cobblestone-like morphological change and failed to produce BMP2-dependent extracellular matrix mineralization. Mechanistically, a miR-199a-5p mimic modified hMSC propagated normal SMAD1/5/9 signaling and expressed osteoblastic transcription factors RUNX2 and Osterix but depicted pronounced upregulation of SOX9 and enhanced expression of essential chondrogenic genes ACAN, COMP, and COL10A1. Mineralization defects, morphological changes, and enhanced chondrogenic gene expression associated with miR-199a-5p mimic over-expression were restored with miR-199a-5p inhibitor suggesting specificity of miR-199a-5p in chondrogenic fate specification. The expression of both the DNM3OS and miR-199a-5p temporally increased and correlated with hMSC chondrogenic differentiation. Although miR-199a-5p overexpression failed to further enhance chondrogenesis, blocking miR-199a-5p activity significantly reduced chondrogenic pellet size, extracellular matrix deposition, and chondrogenic gene expression. Taken together, our results indicate that oscillating miR-199a-5p levels dictate hMSC osteoblast or chondrocyte terminal fate. Our study highlights a functional role of miR-199a-5p as a BMD effector gene at the DNM3 BMD GWAS locus, where patients with cis-regulatory genetic variation which increases miR-199a-5p expression could lead to reduced osteoblast activity.
ABSTRACT
Ikaros is a transcriptional factor required for conventional T cell development, differentiation, and anergy. While the related factors Helios and Eos have defined roles in regulatory T cells (Treg), a role for Ikaros has not been established. To determine the function of Ikaros in the Treg lineage, we generated mice with Treg-specific deletion of the Ikaros gene (Ikzf1). We find that Ikaros cooperates with Foxp3 to establish a major portion of the Treg epigenome and transcriptome. Ikaros-deficient Treg exhibit Th1-like gene expression with abnormal production of IL-2, IFNg, TNFa, and factors involved in Wnt and Notch signaling. While Ikzf1-Treg-cko mice do not develop spontaneous autoimmunity, Ikaros-deficient Treg are unable to control conventional T cell-mediated immune pathology in response to TCR and inflammatory stimuli in models of IBD and organ transplantation. These studies establish Ikaros as a core factor required in Treg for tolerance and the control of inflammatory immune responses.
Subject(s)
Forkhead Transcription Factors , Gene Expression Regulation , Ikaros Transcription Factor , T-Lymphocytes, Regulatory , Animals , Ikaros Transcription Factor/metabolism , Ikaros Transcription Factor/genetics , T-Lymphocytes, Regulatory/immunology , T-Lymphocytes, Regulatory/metabolism , Forkhead Transcription Factors/metabolism , Forkhead Transcription Factors/genetics , Mice , Mice, KnockoutABSTRACT
BACKGROUND: Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms. METHODS: We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts. FINDINGS: We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro. INTERPRETATION: Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting. FUNDING: This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.
Subject(s)
Carpal Tunnel Syndrome , Humans , Carpal Tunnel Syndrome/genetics , Genome-Wide Association Study , Muscle, Skeletal , Chromosome Mapping , Chromatin/genetics , DNA-Binding Proteins/geneticsABSTRACT
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.
ABSTRACT
The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.
ABSTRACT
BACKGROUND AND AIMS: A key histopathological feature of inflammatory bowel disease is damage to the mucosa, including breakdown of the epithelial barrier. Human enteroids and colonoids are a critical bench-to-bedside tool for studying the epithelium in inflammatory bowel disease. The goal of the current study was to define transcriptional differences in healthy versus diseased subjects that are sustained in enteroids and colonoids, including from disease-spared tissue. METHODS: Biopsies and matching enteroid or colonoid cultures from pediatric patients with ileal Crohn disease (N = 6) and control subjects (N = 17) were subjected to RNA sequencing followed by bioinformatic and machine learning analyses. Late passage enteroids were exposed to cytokines to assess durable transcriptional differences. RESULTS: We observed substantial overlap of pathways upregulated in Crohn disease in enteroids and ileal biopsies, as well as colonoids and rectal biopsies. KEGG pathways for cytokine-cytokine receptor interaction, chemokine signaling, protein export, and Toll-like receptor signaling were upregulated in both ileal and rectal biopsies, as well as enteroids and colonoids. In vitro cytokine exposure reactivated genes previously increased in biopsies. Machine learning predicted biopsy location (100% accuracy) and donor disease status (83% accuracy). A random forest classifier generated using ileal enteroids identified rectal colonoids from ileal Crohn disease subjects with 80% accuracy. CONCLUSION: We confirmed transcriptional profiles of Crohn disease biopsies are expressed in enteroids and colonoids. Furthermore, transcriptomic data from disease-spared rectal tissue can identify patients with ileal Crohn disease. Our data support the use of patient enteroids and colonoids as critical translational tools for the study of inflammatory bowel disease.
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
Leukemia , Neuroblastoma , RNA, Long Noncoding , Adult , Humans , Child , RNA, Long Noncoding/genetics , RNA, Long Noncoding/metabolism , Gene Expression Profiling , Neuroblastoma/genetics , Leukemia/genetics , Genomics , Gene Regulatory Networks , Gene Expression Regulation, NeoplasticABSTRACT
Hypoxia induced intestinal barrier injury, microbial translocation, and local/systemic inflammation may contribute to high-altitude associated gastrointestinal complications or symptoms of acute mountain sickness (AMS). Therefore, we tested the hypothesis that six-hours of hypobaric hypoxia increases circulating markers of intestinal barrier injury and inflammation. A secondary aim was to determine if the changes in these markers were different between those with and without AMS. Thirteen participants were exposed to six hours of hypobaric hypoxia, simulating an altitude of 4572â m. Participants completed two 30-minute bouts of exercise during the early hours of hypoxic exposure to mimic typical activity required by those at high altitude. Pre- and post-exposure blood samples were assessed for circulating markers of intestinal barrier injury and inflammation. Data below are presented as mean ± standard deviation or median [interquartile range]. Intestinal fatty acid binding protein (Δ251 [103-410] pgâ¢mL-1; p = 0.002, d = 0.32), lipopolysaccharide binding protein (Δ2 ± 2.4 µgâ¢mL-1; p = 0.011; d = 0.48), tumor necrosis factor-α (Δ10.2 [3-42.2] pgâ¢mL-1; p = 0.005; d = 0.25), interleukin-1ß (Δ1.5 [0-6.7] pgâ¢mL-1 p = 0.042; d = 0.18), and interleukin-1 receptor agonist (Δ3.4 [0.4-5.2] pgâ¢mL-1p = 0.002; d = 0.23) increased from pre- to post-hypoxia. Six of the 13 participants developed AMS; however, the pre- to post-hypoxia changes for each marker were not different between those with and without AMS (p > 0.05 for all indices). These data provide evidence that high altitude exposures can lead to intestinal barrier injury, which may be an important consideration for mountaineers, military personnel, wildland firefighters, and athletes who travel to high altitudes to perform physical work or exercise.
Subject(s)
Altitude Sickness , Physical Exertion , Humans , Hypoxia , Altitude Sickness/complications , Altitude Sickness/diagnosis , Altitude Sickness/metabolism , Altitude , InflammationABSTRACT
Genome-wide association studies (GWAS) in humans have identified loci robustly associated with several heritable diseases or traits, yet little is known about the functional roles of the underlying causal variants in regulating sleep duration or quality. We applied an ATAC-seq/promoter focused Capture C strategy in human iPSC-derived neural progenitors to carry out a "variant-to-gene" mapping campaign that identified 88 candidate sleep effector genes connected to relevant GWAS signals. To functionally validate the role of the implicated effector genes in sleep regulation, we performed a neuron-specific RNA interference screen in the fruit fly, Drosophila melanogaster, followed by validation in zebrafish. This approach identified a number of genes that regulate sleep including a critical role for glycosylphosphatidylinositol (GPI)-anchor biosynthesis. These results provide the first physical variant-to-gene mapping of human sleep genes followed by a model organism-based prioritization, revealing a conserved role for GPI-anchor biosynthesis in sleep regulation.
Subject(s)
Drosophila melanogaster , Glycosylphosphatidylinositols , Animals , Humans , Glycosylphosphatidylinositols/genetics , Drosophila melanogaster/genetics , Genome-Wide Association Study/methods , Zebrafish/genetics , Chromosome Mapping , Genetic Testing , Sleep/geneticsABSTRACT
The use of masks in public settings and when around people has been recommended to limit the spread of Coronavirus disease 2019 (COVID-19) by major public health agencies. Several different types of masks classified as either medical- or non-medical grade are commonly used among the public. However, concerns with difficulty breathing, re-breathing exhaled carbon dioxide, a decrease in arterial oxygen saturation, and a decrease in exercise performance have been raised regarding the use of mask during exercise. We review the current knowledge related to the effect of different masks during exercise on cardiorespiratory, metabolic, thermoregulatory, and perceptual responses. As such, the current literature seems to suggest that there are minimal changes to cardiovascular, metabolic, and no changes to thermoregulatory parameters with facemask use. However, differences in ventilatory parameters have been reported with submaximal and maximal intensity exercise to volitional fatigue. Literature on perceptual responses to exercise indicate an impact on ratings of perceived exertion, dyspnea, and overall discomfort dependent on mask use as well as exercise intensity. In conclusion, data from the current literature suggests a minimal impact on physiological, perceptual, and thermoregulatory responses dependent on the type of mask used during exercise.
Subject(s)
COVID-19 , Masks , Humans , Exercise/physiology , Respiration , Oximetry , DyspneaABSTRACT
Bellovary, Bryanne N., Andrew D. Wells, Zachary J. Fennel, Jeremy B. Ducharme, Jonathan M. Houck, Trevor J. Mayschak, Ann L. Gibson, Scott N. Drum, and Christine M. Mermier. Could orthostatic stress responses predict acute mountain sickness susceptibility before high altitude travel? A pilot study. High Alt Med Biol. 24:19-26, 2023. Purpose: This study assessed head-up tilt (HUT) responses in relation to acute mountain sickness (AMS)-susceptibility during hypoxic exposure. Materials and Methods: Fifteen participants completed three lab visits: (1) protocol familiarization and cycle maximal oxygen consumption (VO2max) test; (2) HUT test consisting of supine rest for 20 minutes followed by 70° tilting for ≤40 minutes; and (3) 6 hours of hypobaric hypoxic exposure (4,572 m) where participants performed two 30-minute cycling bouts separated by 1 hour at a 50% VO2max workload within the first 3 hours and rested when not exercising. During HUT, systolic blood pressure (SBP), diastolic blood pressure, heart rate (HR), and variability (blood pressure variability [BPV] and HR variability [HRV]) were measured continuously. The AMS scores were determined after 6 hours of exposure. Correlations determined relationships between HUT cardiovascular responses and AMS scores. Repeated-measures analysis of variance (ANOVA) assessed differences between those with and without AMS symptoms during HUT. Results: Higher AMS scores correlated with greater change in SBP variability (r = 0.52, p = 0.048) and blunted changes in HRV (root mean square of successive differences between normal heartbeats r = 0.81, p = 0.001, percentage of adjacent normal sinus intervals that differ by more than 50 milliseconds [pNN50] r = 0.87, p < 0.001) during HUT. A pNN50 interaction (p = 0.02) suggested elevated cardiac sympathetic activity at baseline and a blunted increase in cardiac sympathetic influence throughout HUT in those with AMS (pNN50 baseline: AMS = 26.2% ± 15.3%, no AMS = 51.0% ± 13.5%; first 3 minutes into HUT: AMS = 17.2% ± 19.1%, no AMS = 17.1% ± 10.9%; end of HUT: AMS = 6.2% ± 9.1%, no AMS 11.0% ± 10.0%). Conclusions: The results suggest autonomic responses via HUT differ in AMS-susceptible individuals. Changes in HRV and BPV during HUT may be a promising predictive measurement for AMS-susceptibility, but further research is needed for confirmation.