Search

BackgroundNon-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver.MethodsTo identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis.FindingsWe discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels.InterpretationOur results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease.FundingThe work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

The treatment of brain abscess induced by carbapenem-resistant Pseudomonas aeruginosa (CRPA) is a clinical challenge around the world. Apart from novel β-lactam/β-lactamase inhibitors and polymyxins, there are few sufficiently powerful antibiotics that are effective against CRPA-induced infections. Considering the blood-brain barrier factor, there are even fewer drugs that can be used to treat intracranial CRPA-induced infections. In this article, we reported a case of CRPA-induced brain abscess that was successfully treated with intravenous ceftazidime/avibactam and intrathecal colistimethate sodium in a child after intracranial tumor resection.

The prevalence of non-alcoholic fatty liver disease (NAFLD), now also known as metabolic dysfunction-associated fatty liver disease (MAFLD), is rapidly increasing worldwide due to the ongoing obesity epidemic. However, currently the NALFD diagnosis requires non-readily available imaging technologies or liver biopsy, which has drastically limited the sample sizes of NAFLD studies and hampered the discovery of its genetic component. Here we utilized the large UK Biobank (UKB) to accurately estimate the NAFLD status in UKB based on common serum traits and anthropometric measures. Scoring all individuals in UKB for NAFLD risk resulted in 28,396 NAFLD cases and 108,652 healthy individuals at a >90% confidence level. Using this imputed NAFLD status to perform the largest NAFLD genome-wide association study (GWAS) to date, we identified 94 independent (R2 < 0.2) NAFLD GWAS loci, of which 90 have not been identified before; built a polygenic risk score (PRS) model to predict the genetic risk of NAFLD; and used the GWAS variants of imputed NAFLD for a tissue-aware Mendelian randomization analysis that discovered a significant causal effect of NAFLD on coronary artery disease (CAD). In summary, we accurately estimated the NAFLD status in UKB using common serum traits and anthropometric measures, which empowered us to identify 90 GWAS NAFLD loci, build NAFLD PRS, and discover a significant causal effect of NAFLD on CAD.

BackgroundObesity predisposes individuals to multiple cardiometabolic disorders, including type 2 diabetes (T2D). As body mass index (BMI) cannot reliably differentiate fat from lean mass, the metabolically detrimental abdominal obesity has been estimated using waist-hip ratio (WHR). Waist-hip ratio adjusted for body mass index (WHRadjBMI) in turn is a well-established sex-specific marker for abdominal fat and adiposity, and a predictor of adverse metabolic outcomes, such as T2D. However, the underlying genes and regulatory mechanisms orchestrating the sex differences in obesity and body fat distribution in humans are not well understood.MethodsWe searched for genetic master regulators of WHRadjBMI by employing integrative genomics approaches on human subcutaneous adipose RNA sequencing (RNA-seq) data (n ~ 1400) and WHRadjBMI GWAS data (n ~ 700,000) from the WHRadjBMI GWAS cohorts and the UK Biobank (UKB), using co-expression network, transcriptome-wide association study (TWAS), and polygenic risk score (PRS) approaches. Finally, we functionally verified our genomic results using gene knockdown experiments in a human primary cell type that is critical for adipose tissue function.ResultsHere, we identified an adipose gene co-expression network that contains 35 obesity GWAS genes and explains a significant amount of polygenic risk for abdominal obesity and T2D in the UKB (n = 392,551) in a sex-dependent way. We showed that this network is preserved in the adipose tissue data from the Finnish Kuopio Obesity Study and Mexican Obesity Study. The network is controlled by a novel adipose master transcription factor (TF), TBX15, a WHRadjBMI GWAS gene that regulates the network in trans. Knockdown of TBX15 in human primary preadipocytes resulted in changes in expression of 130 network genes, including the key adipose TFs, PPARG and KLF15, which were significantly impacted (FDR < 0.05), thus functionally verifying the trans regulatory effect of TBX15 on the WHRadjBMI co-expression network.ConclusionsOur study discovers a novel key function for the TBX15 TF in trans regulating an adipose co-expression network of 347 adipose, mitochondrial, and metabolically important genes, including PPARG, KLF15, PPARA, ADIPOQ, and 35 obesity GWAS genes. Thus, based on our converging genomic, transcriptional, and functional evidence, we interpret the role of TBX15 to be a main transcriptional regulator in the adipose tissue and discover its importance in human abdominal obesity.

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

Single-nucleus RNA sequencing (snRNA-seq) measures gene expression in individual nuclei instead of cells, allowing for unbiased cell type characterization in solid tissues. We observe that snRNA-seq is commonly subject to contamination by high amounts of ambient RNA, which can lead to biased downstream analyses, such as identification of spurious cell types if overlooked. We present a novel approach to quantify contamination and filter droplets in snRNA-seq experiments, called Debris Identification using Expectation Maximization (DIEM). Our likelihood-based approach models the gene expression distribution of debris and cell types, which are estimated using EM. We evaluated DIEM using three snRNA-seq data sets: (1) human differentiating preadipocytes in vitro, (2) fresh mouse brain tissue, and (3) human frozen adipose tissue (AT) from six individuals. All three data sets showed evidence of extranuclear RNA contamination, and we observed that existing methods fail to account for contaminated droplets and led to spurious cell types. When compared to filtering using these state of the art methods, DIEM better removed droplets containing high levels of extranuclear RNA and led to higher quality clusters. Although DIEM was designed for snRNA-seq, our clustering strategy also successfully filtered single-cell RNA-seq data. To conclude, our novel method DIEM removes debris-contaminated droplets from single-cell-based data fast and effectively, leading to cleaner downstream analysis. Our code is freely available for use at https://github.com/marcalva/diem.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

We present Bisque, a tool for estimating cell type proportions in bulk expression. Bisque implements a regression-based approach that utilizes single-cell RNA-seq (scRNA-seq) or single-nucleus RNA-seq (snRNA-seq) data to generate a reference expression profile and learn gene-specific bulk expression transformations to robustly decompose RNA-seq data. These transformations significantly improve decomposition performance compared to existing methods when there is significant technical variation in the generation of the reference profile and observed bulk expression. Importantly, compared to existing methods, our approach is extremely efficient, making it suitable for the analysis of large genomic datasets that are becoming ubiquitous. When applied to subcutaneous adipose and dorsolateral prefrontal cortex expression datasets with both bulk RNA-seq and snRNA-seq data, Bisque replicates previously reported associations between cell type proportions and measured phenotypes across abundant and rare cell types. We further propose an additional mode of operation that merely requires a set of known marker genes.

Identifying gene-environment interactions (GxEs) contributing to human cardiometabolic disorders is challenging. Here we apply a reverse GxE candidate search by deriving candidate variants from promoter-enhancer interactions that respond to dietary fatty acid challenge through altered chromatin accessibility in human primary adipocytes. We then test all variants residing in the lipid-responsive open chromatin sites within adipocyte promoter-enhancer contacts for interaction effects between the genotype and dietary saturated fat intake on body mass index (BMI) in the UK Biobank. We discover 14 novel GxE variants in 12 lipid-responsive promoters, including well-known lipid genes (LIPE, CARM1, and PLIN2) and novel genes, such as LDB3, for which we provide further functional and integrative genomics evidence. We further identify 24 GxE variants in enhancers, totaling 38 new GxE variants for BMI in the UK Biobank, demonstrating that molecular genomics data produced in physiologically relevant contexts can discover new functional GxE mechanisms in humans.

Chromosomal abnormality is one of the important causes of dysplasia in children. However, due to regional and ethnic differences, the reported rates of chromosomal abnormalities in patients with dysplasia vary greatly. Moreover, the clinical manifestations in children with rare chromosomal diseases were heterogeneous. So, we retrospectively analyzed the karyotype results of 436 children with dysplasia and conducted a detailed analysis of rare chromosomal diseases. The results showed that chromosomal abnormalities were present in 181 of 436 cases. Intellectual disability, dysmorphology, congenital malformations, the disorder of sexual development, and short stature were the main five clinical symptoms in children with chromosomal abnormalities. Moreover, 136 cases of Trisomy 21 (Tri21) were detected, of which 130 were standard Tri21, 5 were robertsonian Tri21, and 1 was chimera type. In addition, 16 cases of rare abnormal karyotype, including complex Tri21, complex Turner syndrome, 4p-syndrome, 18q-syndrome, and 5p-syndrome, were also detected. In summary, chromosome abnormality is one of the important causes of dysplasia in children. Furthermore, prenatal screening and diagnosis could play a great significance in preventing dysplasia in children. In addition, the retrospective analysis of rare cases is valuable for clinical diagnosis and risk assessment of recurrence.

Although recent studies provide evidence for a common genetic basis between complex traits and Mendelian disorders, a thorough quantification of their overlap in a phenotype-specific manner remains elusive. Here, we have quantified the overlap of genes identified through large-scale genome-wide association studies (GWASs) for 62 complex traits and diseases with genes containing mutations known to cause 20 broad categories of Mendelian disorders. We identified a significant enrichment of genes linked to phenotypically matched Mendelian disorders in GWAS gene sets; of the total 1,240 comparisons, a higher proportion of phenotypically matched or related pairs (n = 50 of 92 [54%]) than phenotypically unmatched pairs (n = 27 of 1,148 [2%]) demonstrated significant overlap, confirming a phenotype-specific enrichment pattern. Further, we observed elevated GWAS effect sizes near genes linked to phenotypically matched Mendelian disorders. Finally, we report examples of GWAS variants localized at the transcription start site or physically interacting with the promoters of genes linked to phenotypically matched Mendelian disorders. Our results are consistent with the hypothesis that genes that are disrupted in Mendelian disorders are dysregulated by non-coding variants in complex traits and demonstrate how leveraging findings from related Mendelian disorders and functional genomic datasets can prioritize genes that are putatively dysregulated by local and distal non-coding GWAS variants.

In the original version of this Article, Supplementary Table 10 contained incorrect primer sequences for the mobility shift assay for SNP rs4776984. These errors have now been fixed and the corrected version of the Supplementary Information PDF is available to download from the HTML version of the Article.

Increased adiposity is a hallmark of obesity and overweight, which affect 2.2 billion people world-wide. Understanding the genetic and molecular mechanisms that underlie obesity-related phenotypes can help to improve treatment options and drug development. Here we perform promoter Capture Hi-C in human adipocytes to investigate interactions between gene promoters and distal elements as a transcription-regulating mechanism contributing to these phenotypes. We find that promoter-interacting elements in human adipocytes are enriched for adipose-related transcription factor motifs, such as PPARG and CEBPB, and contribute to heritability of cis-regulated gene expression. We further intersect these data with published genome-wide association studies for BMI and BMI-related metabolic traits to identify the genes that are under genetic cis regulation in human adipocytes via chromosomal interactions. This integrative genomics approach identifies four cis-eQTL-eGene relationships associated with BMI or obesity-related traits, including rs4776984 and MAP2K5, which we further confirm by EMSA, and highlights 38 additional candidate genes.

C21H25EuN2O10, monoclinic, P21/n (no. 14), a = 14.1635(8) Å, b = 12.1381(6) Å, c = 14.4528(7) Å, β = 104.361(2)°, V = 2407.1(2) Å3, Z = 4, Rgt(F) = 0.0304, wRref(F2) = 0.0893, T = 296 K.

C16H11BF2N2, monoclinic, P21/c (no. 14), a = 9.0978(14) Å, b = 18.580(3) Å, c = 15.476(2) Å, β = 99.516(5)°, V = 2580.0(7) Å3, Z = 8, Rgt(F) = 0.0406, wRref(F2) = 0.1307, T = 193 K.

ObjectiveWe recently identified a locus on chromosome 18q11.2 for high serum triglycerides in Mexicans. We hypothesize that the lead genome-wide association study single-nucleotide polymorphism rs9949617, or its linkage disequilibrium proxies, regulates 1 of the 5 genes in the triglyceride-associated region.Approach and resultsWe performed a linkage disequilibrium analysis and found 9 additional variants in linkage disequilibrium (r(2)>0.7) with the lead single-nucleotide polymorphism. To select the variants for functional analyses, we annotated the 10 variants using DNase I hypersensitive sites, transcription factor and chromatin states and identified rs17259126 as the lead candidate variant for functional in vitro validation. Using luciferase transcriptional reporter assay in liver HepG2 cells, we found that the G allele exhibits a significantly lower effect on transcription (P

C32H50Eu2N16O24, monoclinic, P21/c (no. 14), a = 7.685(5) Å, b = 13.482(5) Å, c = 23.240(5) Å, β = 93.578(5)°, V = 2403.2(19) Å3, Z = 2, Rgt(F) = 0.0235, wRref(F2) = 0.0577, T = 293 K.

C18H15BN2O2, monolinic, P21/n (no. 14), a = 13.61(2) Å, b = 7.389(11) Å, c = 16.34(2) Å, β = 113.91(2)°, V = 1502(4) Å3, Z = 4, R gt(F) = 0.0544, wR ref(F 2) = 0.1651, T = 273 K.

C17H13BN2O2, orthorhombic, Pbca (no. 61), a = 10.0820(14) Å, b = 13.1135(19) Å, c = 20.644(3) Å, V = 2729.4 Å3, Z = 8, Rgt(F) = 0.0373, wRref(F2) = 0.1005, T = 193 K.

C18H15BN2O2, monoclinic, P21/c (no. 14), a = 6.8522(9) Å, b = 8.4333(11) Å, c = 13.1908(18) Å, β = 95.880(5)°, V = 758.24(17) Å3, Z = 2, Rgt(F) = 0.0381, wRref(F2) = 0.1068, T = 193 K.

C18H17BN2O2, monoclinic, P21/c (no. 14), a = 12.5782(14) Å, b = 13.6724(14) Å, c = 8.7113(9) Å, β = 97.819(6)°, V = 1484.2(3) Å3, Z = 4, Rgt(F) = 0.0482, wRref(F2) = 0.1357, T = 193 K.

C16H12BBrN2, monoclinic, P21 (no. 4), a = 4.7061(3) Å, b = 10.4639(6) Å, c = 13.5675(8) Å, β = 97.512(2)°, V = 662.39(7) Å3, Z = 2, Rgt(F) = 0.0384, wRref(F2) = 0.0995, T = 189.98 K.

Context: Mitochondria are essential for cellular energy homeostasis, yet their role in subcutaneous adipose tissue (SAT) during different types of weight-loss interventions remains unknown.Objective: To investigate how SAT mitochondria change following diet-induced and bariatric surgery-induced weight-loss interventions in 4 independent weight-loss studies.Methods: The DiOGenes study is a European multicenter dietary intervention with an 8-week low caloric diet (LCD; 800 kcal/d; n = 261) and 6-month weight-maintenance (n = 121) period. The Kuopio Obesity Surgery study (KOBS) is a Roux-en-Y gastric bypass (RYGB) surgery study (n = 172) with a 1-year follow-up. We associated weight-loss percentage with global and 2210 mitochondria-related RNA transcripts in linear regression analysis adjusted for age and sex. We repeated these analyses in 2 studies. The Finnish CRYO study has a 6-week LCD (800-1000 kcal/d; n = 19) and a 10.5-month follow-up. The Swedish DEOSH study is a RYGB surgery study with a 2-year (n = 49) and 5-year (n = 37) follow-up.Results: Diet-induced weight loss led to a significant transcriptional downregulation of oxidative phosphorylation (DiOGenes; ingenuity pathway analysis [IPA] z-scores: -8.7 following LCD, -4.4 following weight maintenance; CRYO: IPA z-score: -5.6, all P < 0.001), while upregulation followed surgery-induced weight loss (KOBS: IPA z-score: 1.8, P < 0.001; in DEOSH: IPA z-scores: 4.0 following 2 years, 0.0 following 5 years). We confirmed an upregulated oxidative phosphorylation at the proteomics level following surgery (IPA z-score: 3.2, P < 0.001).Conclusions: Differentially regulated SAT mitochondria-related gene expressions suggest qualitative alterations between weight-loss interventions, providing insights into the potential molecular mechanistic targets for weight-loss success. [ABSTRACT FROM AUTHOR]