The methodological and ethical concerns of genetic studies of same-sex sexual behavior.

Same-sex sexual behavior has long interested genetics researchers in part because, while there is evidence of heritability, the trait as typically defined is associated with fewer offspring. Investigations of this phenomenon began in the 1990s with linkage studies and continue today with the advent of genome-wide association studies. As this body of research grows, so does critical scientific and ethical review of it. Here, we provide a targeted overview of existing genetics studies on same-sex sexual behavior, highlight the ethical and scientific considerations of this nascent field, and provide recommendations developed by the authors to enhance social and ethical responsibility.

Human vascularized macrophage-islet organoids to model immune-mediated pancreatic ß cell pyroptosis upon viral infection.

There is a paucity of human models to study immune-mediated host damage. Here, we utilized the GeoMx spatial multi-omics platform to analyze immune cell changes in COVID-19 pancreatic autopsy samples, revealing an accumulation of proinflammatory macrophages. Single-cell RNA sequencing (scRNA-seq) analysis of human islets exposed to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or coxsackievirus B4 (CVB4) viruses identified activation of proinflammatory macrophages and ß cell pyroptosis. To distinguish viral versus proinflammatory-macrophage-mediated ß cell pyroptosis, we developed human pluripotent stem cell (hPSC)-derived vascularized macrophage-islet (VMI) organoids. VMI organoids exhibited enhanced marker expression and function in both ß cells and endothelial cells compared with separately cultured cells. Notably, proinflammatory macrophages within VMI organoids induced ß cell pyroptosis. Mechanistic investigations highlighted TNFSF12-TNFRSF12A involvement in proinflammatory-macrophage-mediated ß cell pyroptosis. This study established hPSC-derived VMI organoids as a valuable tool for studying immune-cell-mediated host damage and uncovered the mechanism of ß cell damage during viral exposure.

Human Vascularized Macrophage-Islet Organoids to Model Immune-Mediated Pancreatic ß cell Pyroptosis upon Viral Infection.

There is a paucity of human models to study immune-mediated host damage. Here, we utilized the GeoMx spatial multi-omics platform to analyze immune cell changes in COVID-19 pancreatic autopsy samples, revealing an accumulation of proinflammatory macrophages. Single cell RNA-seq analysis of human islets exposed to SARS-CoV-2 or Coxsackievirus B4 (CVB4) viruses identified activation of proinflammatory macrophages and ß cell pyroptosis. To distinguish viral versus proinflammatory macrophage-mediated ß cell pyroptosis, we developed human pluripotent stem cell (hPSC)-derived vascularized macrophage-islet (VMI) organoids. VMI organoids exhibited enhanced marker expression and function in both ß cells and endothelial cells compared to separately cultured cells. Notably, proinflammatory macrophages within VMI organoids induced ß cell pyroptosis. Mechanistic investigations highlighted TNFSF12-TNFRSF12A involvement in proinflammatory macrophage-mediated ß cell pyroptosis. This study established hPSC-derived VMI organoids as a valuable tool for studying immune cell-mediated host damage and uncovered mechanism of ß cell damage during viral exposure.

CNTools: A computational toolbox for cellular neighborhood analysis from multiplexed images.

Recent studies show that cellular neighborhoods play an important role in evolving biological events such as cancer and diabetes. Therefore, it is critical to accurately and efficiently identify cellular neighborhoods from spatially-resolved single-cell transcriptomic data or single-cell resolution tissue imaging data. In this work, we develop CNTools, a computational toolbox for end-to-end cellular neighborhood analysis on annotated cell images, comprising both the identification and analysis steps. It includes state-of-the-art cellular neighborhood identification methods and post-identification smoothing techniques, with our newly proposed Cellular Neighbor Embedding (CNE) method and Naive Smoothing technique, as well as several established downstream analysis approaches. We applied CNTools on three real-world CODEX datasets and evaluated identification methods with smoothing techniques quantitatively and qualitatively. It shows that CNE with Naive Smoothing overall outperformed other methods and revealed more convincing biological insights. We also provided suggestions on how to choose proper identification methods and smoothing techniques according to input data.

Untangling the genetics of beta cell dysfunction and death in type 1 diabetes.

BACKGROUND: Type 1 diabetes (T1D) is a complex multi-system disease which arises from both environmental and genetic factors, resulting in the destruction of insulin-producing pancreatic beta cells. Over the past two decades, human genetic studies have provided new insight into the etiology of T1D, including an appreciation for the role of beta cells in their own demise. SCOPE OF REVIEW: Here, we outline models supported by human genetic data for the role of beta cell dysfunction and death in T1D. We highlight the importance of strong evidence linking T1D genetic associations to bona fide candidate genes for mechanistic and therapeutic consideration. To guide rigorous interpretation of genetic associations, we describe molecular profiling approaches, genomic resources, and disease models that may be used to construct variant-to-gene links and to investigate candidate genes and their role in T1D. MAJOR CONCLUSIONS: We profile advances in understanding the genetic causes of beta cell dysfunction and death at individual T1D risk loci. We discuss how genetic risk prediction models can be used to address disease heterogeneity. Further, we present areas where investment will be critical for the future use of genetics to address open questions in the development of new treatment and prevention strategies for T1D.

Genetic risk converges on regulatory networks mediating early type 2 diabetes.

Type 2 diabetes mellitus (T2D), a major cause of worldwide morbidity and mortality, is characterized by dysfunction of insulin-producing pancreatic islet ß cells1,2. T2D genome-wide association studies (GWAS) have identified hundreds of signals in non-coding and ß cell regulatory genomic regions, but deciphering their biological mechanisms remains challenging3-5. Here, to identify early disease-driving events, we performed traditional and multiplexed pancreatic tissue imaging, sorted-islet cell transcriptomics and islet functional analysis of early-stage T2D and control donors. By integrating diverse modalities, we show that early-stage T2D is characterized by ß cell-intrinsic defects that can be proportioned into gene regulatory modules with enrichment in signals of genetic risk. After identifying the ß cell hub gene and transcription factor RFX6 within one such module, we demonstrated multiple layers of genetic risk that converge on an RFX6-mediated network to reduce insulin secretion by ß cells. RFX6 perturbation in primary human islet cells alters ß cell chromatin architecture at regions enriched for T2D GWAS signals, and population-scale genetic analyses causally link genetically predicted reduced RFX6 expression with increased T2D risk. Understanding the molecular mechanisms of complex, systemic diseases necessitates integration of signals from multiple molecules, cells, organs and individuals, and thus we anticipate that this approach will be a useful template to identify and validate key regulatory networks and master hub genes for other diseases or traits using GWAS data.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits.

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

Using a modular massively parallel reporter assay to discover context-specific regulatory grammars in type 2 diabetes.

Recent genome-wide association studies have established that most complex disease-associated loci are found in noncoding regions where defining their function is nontrivial. In this study, we leverage a modular massively parallel reporter assay (MPRA) to uncover sequence features linked to context-specific regulatory activity. We screened enhancer activity across a panel of 198-bp fragments spanning over 10k type 2 diabetes- and metabolic trait-associated variants in the 832/13 rat insulinoma cell line, a relevant model of pancreatic beta cells. We explored these fragments' context sensitivity by comparing their activities when placed up-or downstream of a reporter gene, and in combination with either a synthetic housekeeping promoter (SCP1) or a more biologically relevant promoter corresponding to the human insulin gene ( INS ). We identified clear effects of MPRA construct design on measured fragment enhancer activity. Specifically, a subset of fragments (n = 702/11,656) displayed positional bias, evenly distributed across up- and downstream preference. A separate set of fragments exhibited promoter bias (n = 698/11,656), mostly towards the cell-specific INS promoter (73.4%). To identify sequence features associated with promoter preference, we used Lasso regression with 562 genomic annotations and discovered that fragments with INS promoter-biased activity are enriched for HNF1 motifs. HNF1 family transcription factors are key regulators of glucose metabolism disrupted in maturity onset diabetes of the young (MODY), suggesting genetic convergence between rare coding variants that cause MODY and common T2D-associated regulatory variants. We designed a follow-up MPRA containing HNF1 motif-enriched fragments and observed several instances where deletion or mutation of HNF1 motifs disrupted the INS promoter-biased enhancer activity, specifically in the beta cell model but not in a skeletal muscle cell line, another diabetes-relevant cell type. Together, our study suggests that cell-specific regulatory activity is partially influenced by enhancer-promoter compatibility and indicates that careful attention should be paid when designing MPRA libraries to capture context-specific regulatory processes at disease-associated genetic signals.

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis.

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

Three-dimensional chromatin re-organization during muscle stem cell aging.

Age-related skeletal muscle atrophy or sarcopenia is a significant societal problem that is becoming amplified as the world's population continues to increase. The regeneration of damaged skeletal muscle is mediated by muscle stem cells, but in old age muscle stem cells become functionally attenuated. The molecular mechanisms that govern muscle stem cell aging encompass changes across multiple regulatory layers and are integrated by the three-dimensional organization of the genome. To quantitatively understand how hierarchical chromatin architecture changes during muscle stem cell aging, we generated 3D chromatin conformation maps (Hi-C) and integrated these datasets with multi-omic (chromatin accessibility and transcriptome) profiles from bulk populations and single cells. We observed that muscle stem cells display static behavior at global scales of chromatin organization during aging and extensive rewiring of local contacts at finer scales that were associated with variations in transcription factor binding and aberrant gene expression. These data provide insights into genome topology as a regulator of molecular function in stem cell aging.

FixItFelix: improving genomic analysis by fixing reference errors.

The current version of the human reference genome, GRCh38, contains a number of errors including 1.2 Mbp of falsely duplicated and 8.04 Mbp of collapsed regions. These errors impact the variant calling of 33 protein-coding genes, including 12 with medical relevance. Here, we present FixItFelix, an efficient remapping approach, together with a modified version of the GRCh38 reference genome that improves the subsequent analysis across these genes within minutes for an existing alignment file while maintaining the same coordinates. We showcase these improvements over multi-ethnic control samples, demonstrating improvements for population variant calling as well as eQTL studies.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation.

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits.

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

Genetic effects on liver chromatin accessibility identify disease regulatory variants.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures.

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Multi-omic profiling of pituitary thyrotropic cells and progenitors.

BACKGROUND: The pituitary gland is a neuroendocrine organ containing diverse cell types specialized in secreting hormones that regulate physiology. Pituitary thyrotropes produce thyroid-stimulating hormone (TSH), a critical factor for growth and maintenance of metabolism. The transcription factors POU1F1 and GATA2 have been implicated in thyrotrope fate, but the transcriptomic and epigenomic landscapes of these neuroendocrine cells have not been characterized. The goal of this work was to discover transcriptional regulatory elements that drive thyrotrope fate. RESULTS: We identified the transcription factors and epigenomic changes in chromatin that are associated with differentiation of POU1F1-expressing progenitors into thyrotropes using cell lines that represent an undifferentiated Pou1f1 lineage progenitor (GHF-T1) and a committed thyrotrope line that produces TSH (TαT1). We compared RNA-seq, ATAC-seq, histone modification (H3K27Ac, H3K4Me1, and H3K27Me3), and POU1F1 binding in these cell lines. POU1F1 binding sites are commonly associated with bZIP transcription factor consensus binding sites in GHF-T1 cells and Helix-Turn-Helix (HTH) or basic Helix-Loop-Helix (bHLH) factors in TαT1 cells, suggesting that these classes of transcription factors may recruit or cooperate with POU1F1 binding at unique sites. We validated enhancer function of novel elements we mapped near Cga, Pitx1, Gata2, and Tshb by transfection in TαT1 cells. Finally, we confirmed that an enhancer element near Tshb can drive expression in thyrotropes of transgenic mice, and we demonstrate that GATA2 enhances Tshb expression through this element. CONCLUSION: These results extend the ENCODE multi-omic profiling approach to the pituitary gland, which should be valuable for understanding pituitary development and disease pathogenesis.

Chromatin information content landscapes inform transcription factor and DNA interactions.

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

CHD7 promotes neural progenitor differentiation in embryonic stem cells via altered chromatin accessibility and nascent gene expression.

CHARGE syndrome, a rare multiple congenital anomaly condition, is caused by haploinsufficiency of the chromatin remodeling protein gene CHD7 (Chromodomain helicase DNA binding protein 7). Brain abnormalities and intellectual disability are commonly observed in individuals with CHARGE, and neuronal differentiation is reduced in CHARGE patient-derived iPSCs and conditional knockout mouse brains. However, the mechanisms of CHD7 function in nervous system development are not well understood. In this study, we asked whether CHD7 promotes gene transcription in neural progenitor cells via changes in chromatin accessibility. We used Chd7 null embryonic stem cells (ESCs) derived from Chd7 mutant mouse blastocysts as a tool to investigate roles of CHD7 in neuronal and glial differentiation. Loss of Chd7 significantly reduced neuronal and glial differentiation. Sholl analysis showed that loss of Chd7 impaired neuronal complexity and neurite length in differentiated neurons. Genome-wide studies demonstrated that loss of Chd7 leads to modified chromatin accessibility (ATAC-seq) and differential nascent expression (Bru-Seq) of neural-specific genes. These results suggest that CHD7 acts preferentially to alter chromatin accessibility of key genes during the transition of NPCs to neurons to promote differentiation. Our results form a basis for understanding the cell stage-specific roles for CHD7-mediated chromatin remodeling during cell lineage acquisition.

BAF60a deficiency uncouples chromatin accessibility and cold sensitivity from white fat browning.

Brown and beige fat share a remarkably similar transcriptional program that supports fuel oxidation and thermogenesis. The chromatin-remodeling machinery that governs genome accessibility and renders adipocytes poised for thermogenic activation remains elusive. Here we show that BAF60a, a subunit of the SWI/SNF chromatin-remodeling complexes, serves an indispensable role in cold-induced thermogenesis in brown fat. BAF60a maintains chromatin accessibility at PPARγ and EBF2 binding sites for key thermogenic genes. Surprisingly, fat-specific BAF60a inactivation triggers more pronounced cold-induced browning of inguinal white adipose tissue that is linked to induction of MC2R, a receptor for the pituitary hormone ACTH. Elevated MC2R expression sensitizes adipocytes and BAF60a-deficient adipose tissue to thermogenic activation in response to ACTH stimulation. These observations reveal an unexpected dichotomous role of BAF60a-mediated chromatin remodeling in transcriptional control of brown and beige gene programs and illustrate a pituitary-adipose signaling axis in the control of thermogenesis.