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.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic human embryonic stem cell-derived ß-like cells.

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on ß cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on ß cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with ß cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.

Modeling islet enhancers using deep learning identifies candidate causal variants at loci associated with T2D and glycemic traits.

Genetic association studies have identified hundreds of independent signals associated with type 2 diabetes (T2D) and related traits. Despite these successes, the identification of specific causal variants underlying a genetic association signal remains challenging. In this study, we describe a deep learning (DL) method to analyze the impact of sequence variants on enhancers. Focusing on pancreatic islets, a T2D relevant tissue, we show that our model learns islet-specific transcription factor (TF) regulatory patterns and can be used to prioritize candidate causal variants. At 101 genetic signals associated with T2D and related glycemic traits where multiple variants occur in linkage disequilibrium, our method nominates a single causal variant for each association signal, including three variants previously shown to alter reporter activity in islet-relevant cell types. For another signal associated with blood glucose levels, we biochemically test all candidate causal variants from statistical fine-mapping using a pancreatic islet beta cell line and show biochemical evidence of allelic effects on TF binding for the model-prioritized variant. To aid in future research, we publicly distribute our model and islet enhancer perturbation scores across ~67 million genetic variants. We anticipate that DL methods like the one presented in this study will enhance the prioritization of candidate causal variants for functional studies.

Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 hours.

Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycemia, beta cell glucotoxicity, and ultimately type 2 diabetes (T2D). In this study, we sought to explore the effects of hyperglycemia on human pancreatic islet (HPI) gene expression by exposing HPIs from two donors to low (2.8mM) and high (15.0mM) glucose concentrations over 24 hours, assaying the transcriptome at seven time points using single-cell RNA sequencing (scRNA-seq). We modeled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Across all cell types, we identified 1,528 genes associated with time, 1,185 genes associated with glucose exposure, and 845 genes associated with interaction effects between time and glucose. We clustered differentially expressed genes across cell types and found 347 modules of genes with similar expression patterns across time and glucose conditions, including two beta cell modules enriched in genes associated with T2D. Finally, by integrating genomic features from this study and genetic summary statistics for T2D and related traits, we nominate 363 candidate effector genes that may underlie genetic associations for T2D and related traits.

Bone dysplasia in Hutchinson-Gilford progeria syndrome is associated with dysregulated differentiation and function of bone cell populations.

Hutchinson-Gilford progeria syndrome (HGPS) is a premature aging disorder affecting tissues of mesenchymal origin. Most individuals with HGPS harbor a de novo c.1824C > T (p.G608G) mutation in the gene encoding lamin A (LMNA), which activates a cryptic splice donor site resulting in production of the toxic "progerin" protein. Clinical manifestations include growth deficiency, lipodystrophy, sclerotic dermis, cardiovascular defects, and bone dysplasia. Here we utilized the LmnaG609G knock-in (KI) mouse model of HGPS to further define mechanisms of bone loss associated with normal and premature aging disorders. Newborn skeletal staining of KI mice revealed altered rib cage shape and spinal curvature, and delayed calvarial mineralization with increased craniofacial and mandibular cartilage content. MicroCT analysis and mechanical testing of adult femurs indicated increased fragility associated with reduced bone mass, recapitulating the progressive bone deterioration that occurs in HGPS patients. We investigated mechanisms of bone loss in KI mice at the cellular level in bone cell populations. Formation of wild-type and KI osteoclasts from marrow-derived precursors was inhibited by KI osteoblast-conditioned media in vitro, suggesting a secreted factor(s) responsible for decreased osteoclasts on KI trabecular surfaces in vivo. Cultured KI osteoblasts exhibited abnormal differentiation characterized by reduced deposition and mineralization of extracellular matrix with increased lipid accumulation compared to wild-type, providing a mechanism for altered bone formation. Furthermore, quantitative analyses of KI transcripts confirmed upregulation of adipogenic genes both in vitro and in vivo. Thus, osteoblast phenotypic plasticity, inflammation and altered cellular cross-talk contribute to abnormal bone formation in HGPS mice.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived ß-like cells.

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined ß-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-ß cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

Human pancreatic islet microRNAs implicated in diabetes and related traits by large-scale genetic analysis.

Genetic studies have identified ≥240 loci associated with the risk of type 2 diabetes (T2D), yet most of these loci lie in non-coding regions, masking the underlying molecular mechanisms. Recent studies investigating mRNA expression in human pancreatic islets have yielded important insights into the molecular drivers of normal islet function and T2D pathophysiology. However, similar studies investigating microRNA (miRNA) expression remain limited. Here, we present data from 63 individuals, the largest sequencing-based analysis of miRNA expression in human islets to date. We characterized the genetic regulation of miRNA expression by decomposing the expression of highly heritable miRNAs into cis- and trans-acting genetic components and mapping cis-acting loci associated with miRNA expression [miRNA-expression quantitative trait loci (eQTLs)]. We found i) 84 heritable miRNAs, primarily regulated by trans-acting genetic effects, and ii) 5 miRNA-eQTLs. We also used several different strategies to identify T2D-associated miRNAs. First, we colocalized miRNA-eQTLs with genetic loci associated with T2D and multiple glycemic traits, identifying one miRNA, miR-1908, that shares genetic signals for blood glucose and glycated hemoglobin (HbA1c). Next, we intersected miRNA seed regions and predicted target sites with credible set SNPs associated with T2D and glycemic traits and found 32 miRNAs that may have altered binding and function due to disrupted seed regions. Finally, we performed differential expression analysis and identified 14 miRNAs associated with T2D status-including miR-187-3p, miR-21-5p, miR-668, and miR-199b-5p-and 4 miRNAs associated with a polygenic score for HbA1c levels-miR-216a, miR-25, miR-30a-3p, and miR-30a-5p.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19.

BACKGROUND: COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain. SUBJECTS/METHODS: In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry. RESULTS: Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression. CONCLUSIONS: Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Protein expression of the gp78 E3 ligase predicts poor breast cancer outcome based on race.

Women of African ancestry suffer higher rates of breast cancer mortality compared with all other groups in the United States. Though the precise reasons for these disparities remain unclear, many recent studies have implicated a role for differences in tumor biology. Using an epitope-validated antibody against the endoplasmic reticulum-associated E3 ligase, gp78, we show that elevated levels of gp78 in patient breast cancer cells predict poor survival. Moreover, high levels of gp78 are associated with poor outcomes in both ER+ and ER- tumors, and breast cancers expressing elevated amounts of gp78 protein are enriched in gene expression pathways that influence cell cycle, metabolism, receptor-mediated signaling, and cell stress response pathways. In multivariate analysis adjusted for subtype and grade, gp78 protein is an independent predictor of poor outcomes in women of African ancestry. Furthermore, gene expression signatures, derived from patients stratified by gp78 protein expression, are strong predictors of recurrence and pathological complete response in retrospective clinical trial data and share many common features with gene sets previously identified to be overrepresented in breast cancers based on race. These findings implicate a prominent role for gp78 in tumor progression and offer insights into our understanding of racial differences in breast cancer outcomes.

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.

Kaiso (ZBTB33) subcellular partitioning functionally links LC3A/B, the tumor microenvironment, and breast cancer survival.

The use of digital pathology for the histomorphologic profiling of pathological specimens is expanding the precision and specificity of quantitative tissue analysis at an unprecedented scale; thus, enabling the discovery of new and functionally relevant histological features of both predictive and prognostic significance. In this study, we apply quantitative automated image processing and computational methods to profile the subcellular distribution of the multi-functional transcriptional regulator, Kaiso (ZBTB33), in the tumors of a large racially diverse breast cancer cohort from a designated health disparities region in the United States. Multiplex multivariate analysis of the association of Kaiso's subcellular distribution with other breast cancer biomarkers reveals novel functional and predictive linkages between Kaiso and the autophagy-related proteins, LC3A/B, that are associated with features of the tumor immune microenvironment, survival, and race. These findings identify effective modalities of Kaiso biomarker assessment and uncover unanticipated insights into Kaiso's role in breast cancer progression.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition.

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.

Single-cell transcriptomics from human pancreatic islets: sample preparation matters.

Single-cell RNA sequencing (scRNA-seq) of human primary tissues is a rapidly emerging tool for investigating human health and disease at the molecular level. However, optimal processing of solid tissues presents a number of technical and logistical challenges, especially for tissues that are only available at autopsy, which includes pancreatic islets, a tissue that is highly relevant to diabetes. To assess the possible effects of different sample preparation protocols on fresh islet samples, we performed a detailed comparison of scRNA-seq data generated with islets isolated from a human donor but processed according to four treatment strategies, including fixation and cryopreservation. We found significant and reproducible differences in the proportion of cell types identified, and more minor effects on cell-specific patterns of gene expression. Fresh islets from a second donor confirmed gene expression signatures of alpha and beta subclusters. These findings may well apply to other tissues, emphasizing the need for careful consideration when choosing processing methods, comparing results between different studies, and/or interpreting data in the context of multiple cell types from preserved tissue.

Racial Differences in the Association Between Luminal Master Regulator Gene Expression Levels and Breast Cancer Survival.

PURPOSE: Compared with their European American (EA) counterparts, African American (AA) women are more likely to die from breast cancer in the United States. This disparity is greatest in hormone receptor-positive subtypes. Here we uncover biological factors underlying this disparity by comparing functional expression and prognostic significance of master transcriptional regulators of luminal differentiation. EXPERIMENTAL DESIGN: Data and biospecimens from 262 AA and 293 EA patients diagnosed with breast cancer from 2001 to 2010 at a major medical center were analyzed by IHC for functional biomarkers of luminal differentiation, including estrogen receptor (ESR1) and its pioneer factors, FOXA1 and GATA3. Integrated comparison of protein levels with network-level gene expression analysis uncovered predictive correlations with race and survival. RESULTS: Univariate or multivariate HRs for overall survival, estimated from digital IHC scoring of nuclear antigen, show distinct differences in the magnitude and significance of these biomarkers to predict survival based on race: ESR1 [EA HR = 0.47; 95% confidence interval (CI), 0.31-0.72 and AA HR = 0.77; 95% CI, 0.48-1.18]; FOXA1 (EA HR = 0.38; 95% CI, 0.23-0.63 and AA HR = 0.53; 95% CI, 0.31-0.88), and GATA3 (EA HR = 0.36; 95% CI, 0.23-0.56; AA HR = 0.57; CI, 0.56-1.4). In addition, we identify genes in the downstream regulons of these biomarkers highly correlated with race and survival. CONCLUSIONS: Even within clinically homogeneous tumor groups, regulatory networks that drive mammary luminal differentiation reveal race-specific differences in their association with clinical outcome. Understanding these biomarkers and their downstream regulons will elucidate the intrinsic mechanisms that drive racial disparities in breast cancer survival.

Epigenetic re-wiring of breast cancer by pharmacological targeting of C-terminal binding protein.

The C-terminal binding protein (CtBP) is an NADH-dependent dimeric family of nuclear proteins that scaffold interactions between transcriptional regulators and chromatin-modifying complexes. Its association with poor survival in several cancers implicates CtBP as a promising target for pharmacological intervention. We employed computer-assisted drug design to search for CtBP inhibitors, using quantitative structure-activity relationship (QSAR) modeling and docking. Functional screening of these drugs identified 4 compounds with low toxicity and high water solubility. Micro molar concentrations of these CtBP inhibitors produces significant de-repression of epigenetically silenced pro-epithelial genes, preferentially in the triple-negative breast cancer cell line MDA-MB-231. This epigenetic reprogramming occurs through eviction of CtBP from gene promoters; disrupted recruitment of chromatin-modifying protein complexes containing LSD1, and HDAC1; and re-wiring of activating histone marks at targeted genes. In functional assays, CtBP inhibition disrupts CtBP dimerization, decreases cell migration, abolishes cellular invasion, and improves DNA repair. Combinatorial use of CtBP inhibitors with the LSD1 inhibitor pargyline has synergistic influence. Finally, integrated correlation of gene expression in breast cancer patients with nuclear levels of CtBP1 and LSD1, reveals new potential therapeutic vulnerabilities. These findings implicate a broad role for this class of compounds in strategies for epigenetically targeted therapeutic intervention.

Leptospire genomic diversity revealed by microarray-based comparative genomic hybridization.

Comparative genomic hybridization was used to compare genetic diversity of five strains of Leptospira (Leptospira interrogans serovars Bratislava, Canicola, and Hebdomadis and Leptospira kirschneri serovars Cynopteri and Grippotyphosa). The array was designed based on two available sequenced Leptospira reference genomes, those of L. interrogans serovar Copenhageni and L. interrogans serovar Lai. A comparison of genetic contents showed that L. interrogans serovar Bratislava was closest to the reference genomes while L. kirschneri serovar Grippotyphosa had the least similarity to the reference genomes. Cluster analysis indicated that L. interrogans serovars Bratislava and Hebdomadis clustered together first, followed by L. interrogans serovar Canicola, before the two L. kirschneri strains. Confirmed/potential virulence factors identified in previous research were also detected in the tested strains.

Generation and validation of a Shewanella oneidensis MR-1 clone set for protein expression and phage display.

A comprehensive gene collection for S. oneidensis was constructed using the lambda recombinase (Gateway) cloning system. A total of 3584 individual ORFs (85%) have been successfully cloned into the entry plasmids. To validate the use of the clone set, three sets of ORFs were examined within three different destination vectors constructed in this study. Success rates for heterologous protein expression of S. oneidensis His- or His/GST-tagged proteins in E. coli were approximately 70%. The ArcA and NarP transcription factor proteins were tested in an in vitro binding assay to demonstrate that functional proteins can be successfully produced using the clone set. Further functional validation of the clone set was obtained from phage display experiments in which a phage encoding thioredoxin was successfully isolated from a pool of 80 different clones after three rounds of biopanning using immobilized anti-thioredoxin antibody as a target. This clone set complements existing genomic (e.g., whole-genome microarray) and other proteomic tools (e.g., mass spectrometry-based proteomic analysis), and facilitates a wide variety of integrated studies, including protein expression, purification, and functional analyses of proteins both in vivo and in vitro.

Influence of bicarbonate, sulfate, and electron donors on biological reduction of uranium and microbial community composition.

A microcosm study was performed to investigate the effect of ethanol and acetate on uranium(VI) biological reduction and microbial community changes under various geochemical conditions. Each microcosm contained an uranium-contaminated sediment (up to 2.8 g U/kg) suspended in buffer with bicarbonate at concentrations of either 1 or 40 mM and sulfate at either 1.1 or 3.2 mM. Ethanol or acetate was used as an electron donor. Results indicate that ethanol yielded in significantly higher U(VI) reduction rates than acetate. A low bicarbonate concentration (1 mM) was favored for U(VI) bioreduction to occur in sediments, but high concentrations of bicarbonate (40 mM) and sulfate (3.2 mM) decreased the reduction rates of U(VI). Microbial communities were dominated by species from the Geothrix genus and Proteobacteria phylum in all microcosms. However, species in the Geobacteraceae family capable of reducing U(VI) were significantly enriched by ethanol and acetate in low-bicarbonate buffer. Ethanol increased the population of unclassified Desulfuromonales, while acetate increased the population of Desulfovibrio. Additionally, species in the Geobacteraceae family were not enriched in high-bicarbonate buffer, but the Geothrix and the unclassified Betaproteobacteria species were enriched. This study concludes that ethanol could be a better electron donor than acetate for reducing U(VI) under given experimental conditions, and electron donor and groundwater geochemistry alter microbial communities responsible for U(VI) reduction.

Diversity of functional genes for methanotrophs in sediments associated with gas hydrates and hydrocarbon seeps in the Gulf of Mexico.

Methanotrophs are ubiquitous in soil, fresh water and the open ocean, but have not been well characterized in deep-sea hydrocarbon seeps and gas hydrates, where methane is unusually abundant. Here we report the presence of new functional genes for the aerobic oxidation of methane by methanotrophs in marine sediments associated with gas hydrates and hydrocarbon seeps in the Gulf of Mexico. Samples were collected from two hydrate locations (GC185 and GC234): one hydrocarbon-seep location at a brine pool (GC233) and one background-marine location about 1.2 miles north of the brine pool (NBP). Community DNA was extracted from each location to establish clone libraries for the pmoA functional gene using a PCR-based cloning approach. Three hundred and ninety clones were screened by sequencing and 46 operational taxonomic units were obtained. Eight operational taxonomic units were present in every sample; one of them was predominant and accounted for 22.8-25.3% of each clone library. Principal-component analysis indicated that samples GC185 and GC234 were closely related and, along with GC233, were significantly different from NBP. These results indicate that methanotrophic communities may be similarly impacted by hydrocarbons at the gas-hydrate and seep sites, and can be distinguished from methanotrophic communities in the normal marine sediment. Furthermore, cluster analysis showed that 84.8% of operational taxonomic units from all samples formed distinct clusters, which could not be grouped with any published pmoA sequences, indicating that a considerable number of novel methanotrophic species may exist in the Gulf of Mexico.

Transcriptome of a Nitrosomonas europaea mutant with a disrupted nitrite reductase gene (nirK).

Global gene expression was compared between the Nitrosomonas europaea wild type and a nitrite reductase-deficient mutant using a genomic microarray. Forty-one genes were differentially regulated between the wild type and the nirK mutant, including the nirK operon, genes for cytochrome c oxidase, and seven iron uptake genes. Relationships of differentially regulated genes to the nirK mutant phenotype are discussed.