Cytosolic nucleic acid sensing and mitochondrial transcriptomic changes as early triggers of metabolic disease in db/db mice.

Animal models of diabetes, such as db/db mice, are a useful tool for deciphering the genetic background of molecular changes at the initial stages of disease development. Our goal was to find early transcriptomic changes in three tissues involved in metabolism regulation in db/db mice: adipose tissue, muscle tissue and liver tissue. Nine animals (three per time point) were studied. Tissues were collected at 8, 12 and 16 weeks of age. Transcriptome-wide analysis was performed using mRNA-seq. Libraries were sequenced on NextSeq (Illumina). Differential expression (DE) analysis was performed with edgeR. The analysis of the gene expression profile shared by all three tissues revealed eight upregulated genes (Irf7, Sp100, Neb, Stat2, Oas2, Rtp4, H2-T24 and Oasl2) as early as between 8 and 12 weeks of age. The most pronounced differences were found in liver tissue: nine DE genes between 8 and 12 weeks of age (Irf7, Ly6a, Ly6g6d, H2-Dma, Pld4, Ly86, Fcer1g, Ly6e and Idi1) and five between 12 and 16 weeks of age (Irf7, Plac8, Ifi44, Xaf1 and Ly6a) (adj. p-value < 0.05). The mitochondrial transcriptomic profile also changed with time: we found two downregulated genes in mice between 8 and 12 weeks old (Ckmt2 and Cox6a2) and five DE genes between 12 and 16 weeks of age (Mavs, Tomm40L, Mtfp1, Ckmt2 and Cox6a2). The KEGG pathway analysis showed significant enrichment in pathways related to the autoimmune response and cytosolic DNA sensing. Our results suggest an important involvement of the immunological response, mainly cytosolic nucleic acid sensing, and mitochondrial signalling in the early stages of diabetes and obesity.

MicroRNA composition of plasma extracellular vesicles: a harbinger of late cardiotoxicity of doxorubicin.

BACKGROUND: The use of doxorubicin is associated with an increased risk of acute and long-term cardiomyopathy. Despite the constantly growing number of cancer survivors, little is known about the transcriptional mechanisms which progress in the time leading to a severe cardiac outcome. It is also unclear whether long-term transcriptomic alterations related to doxorubicin use are similar to transcriptomic patterns present in patients suffering from other cardiomyopathies. METHODS: We have sequenced miRNA from total plasma and extracellular vesicles (EVs) from 66 acute lymphoblastic leukemia (ALL) survivors and 61 healthy controls (254 samples in total). We then analyzed processes regulated by differentially expressed circulating miRNAs and cross-validated results with the data of patients with clinically manifested cardiomyopathies. RESULTS: We found that especially miRNAs contained within EVs may be informative in terms of cardiomyopathy development and may regulate pathways related to neurotrophin signaling, transforming growth factor beta (TGFß) or epidermal growth factor receptors (ErbB). We identified vesicular miR-144-3p and miR-423-3p as the most variable between groups and significantly correlated with echocardiographic parameters and, respectively, for plasma: let-7g-5p and miR-16-2-3p. Moreover, vesicular miR-144-3p correlates with the highest number of echocardiographic parameters and is differentially expressed in the circulation of patients with dilated cardiomyopathy. We also found that distribution of particular miRNAs between of plasma and EVs (proportion between compartments) e.g., miR-184 in ALL, is altered, suggesting changes within secretory and miRNA sorting mechanisms. CONCLUSIONS: Our results show that transcriptomic changes resulting from doxorubicin induced myocardial injury are reflected in circulating miRNA levels and precede development of the late onset cardiomyopathy phenotype. Among miRNAs related to cardiac function, we found vesicular miR-144-3p and miR-423-3p, as well as let-7g-5p and miR-16-2-3p contained in the total plasma. Selection of source for such studies (plasma or EVs) is of critical importance, as distribution of some miRNA between plasma and EVs is altered in ALL survivors, in comparison to healthy people, which suggests that doxorubicin-induced changes include miRNA sorting and export to extracellular space.

Variants influencing age at diagnosis of HNF1A-MODY.

BACKGROUND: HNF1A-MODY is a monogenic form of diabetes caused by variants in the HNF1A gene. Different HNF1A variants are associated with differences in age of disease onset, but other factors are postulated to influence this trait. Here, we searched for genetic variants influencing age of HNF1A-MODY onset. METHODS: Blood samples from 843 HNF1A-MODY patients from Czech Republic, France, Poland, Slovakia, the UK and the US were collected. A validation set consisted of 121 patients from the US. We conducted a genome-wide association study in 843 HNF1A-MODY patients. Samples were genotyped using Illumina Human Core arrays. The core analysis was performed using the GENESIS package in R statistical software. Kinship coefficients were estimated with the KING and PC-Relate algorithms. In the linear mixed model, we accounted for year of birth, sex, and location of the HNF1A causative variant. RESULTS: A suggestive association with age of disease onset was observed for rs2305198 (p = 2.09E-07) and rs7079157 (p = 3.96E-06) in the HK1 gene, rs2637248 in the LRMDA gene (p = 2.44E-05), and intergenic variant rs2825115 (p = 2.04E-05). Variant rs2637248 reached nominal significance (p = 0.019), while rs7079157 (p = 0.058) and rs2825115 (p = 0.068) showed suggestive association with age at diabetes onset in the validation set. CONCLUSIONS: rs2637248 in the LRMDA gene is associated with age at diabetes onset in HNF1A-MODY patients.

The transcriptome-wide association search for genes and genetic variants which associate with BMI and gestational weight gain in women with type 1 diabetes.

BACKGROUND: Clinical data suggest that BMI and gestational weight gain (GWG) are strongly interconnected phenotypes; however, the genetic basis of the latter is rather unclear. Here we aim to find genes and genetic variants which influence BMI and/or GWG. METHODS: We have genotyped 316 type 1 diabetics using Illumina Infinium Omni Express Exome-8 v1.4 arrays. The GIANT, ARIC and T2D-GENES summary statistics were used for TWAS (performed with PrediXcan) in adipose tissue. Next, the analysis of association of imputed expression with BMI in the general and diabetic cohorts (Analysis 1 and 2) or GWG (Analysis 3 and 4) was performed, followed by variant association analysis (1 Mb around identified loci) with the mentioned phenotypes. RESULTS: In Analysis 1 we have found 175 BMI associated genes and 19 variants (p < 10-4) which influenced GWG, with the strongest association for rs11465293 in CCL24 (p = 3.18E-06). Analysis 2, with diabetes included in the model, led to discovery of 1812 BMI associated loci and 207 variants (p < 10-4) influencing GWG, with the strongest association for rs9690213 in PODXL (p = 9.86E-07). In Analysis 3, among 648 GWG associated loci, 2091 variants were associated with BMI (FDR < 0.05). In Analysis 4, 7 variants in GWG associated loci influenced BMI in the ARIC cohort. CONCLUSIONS: Here, we have shown that loci influencing BMI might have an impact on GWG and GWG associated loci might influence BMI, both in the general and T1DM cohorts. The results suggest that both phenotypes are related to insulin signaling, glucose homeostasis, mitochondrial metabolism, ubiquitinoylation and inflammatory responses.

Dysregulation of Transcription Factor Activity During Formation of Cancer-Associated Fibroblasts.

The reciprocal interactions between cancer cells and the quiescent fibroblasts leading to the activation of cancer-associated fibroblasts (CAFs) serve an important role in cancer progression. Here, we investigated the activation of transcription factors (TFs) in prostate fibroblasts (WPMY cell line) co-cultured with normal prostate or tumorous cells (RWPE1 and RWPE2 cell lines, respectively). After indirect co-cultures, we performed mRNA-seq and predicted TF activity using mRNA expression profiles with the Systems EPigenomics Inference of Regulatory Activity (SEPIRA) package and the GTEx and mRNA-seq data of 483 cultured fibroblasts. The initial differential expression analysis between time points and experimental conditions showed that co-culture with normal epithelial cells mainly promotes an inflammatory response in fibroblasts, whereas with the cancerous epithelial, it stimulates transformation by changing the expression of the genes associated with microfilaments. TF activity analysis revealed only one positively regulated TF in the RWPE1 co-culture alone, while we observed dysregulation of 45 TFs (7 decreased activity and 38 increased activity) uniquely in co-culture with RWPE2. Pathway analysis showed that these 45 dysregulated TFs in fibroblasts co-cultured with RWPE2 cells may be associated with the RUNX1 and PTEN pathways. Moreover, we showed that observed dysregulation could be associated with FER1L4 expression. We conclude that phenotypic changes in fibroblast responses to co-culturing with cancer epithelium result from orchestrated dysregulation of signaling pathways that favor their transformation and motility rather than proinflammatory status. This dysregulation can be observed both at the TF and transcriptome levels.

Antibodies Enhance the Suppressive Activity of Extracellular Vesicles in Mouse Delayed-Type Hypersensitivity.

Previously, we showed that mouse delayed-type hypersensitivity (DTH) can be antigen-specifically downregulated by suppressor T cell-derived miRNA-150 carried by extracellular vesicles (EVs) that target antigen-presenting macrophages. However, the exact mechanism of the suppressive action of miRNA-150-targeted macrophages on effector T cells remained unclear, and our current studies aimed to investigate it. By employing the DTH mouse model, we showed that effector T cells were inhibited by macrophage-released EVs in a miRNA-150-dependent manner. This effect was enhanced by the pre-incubation of EVs with antigen-specific antibodies. Their specific binding to MHC class II-expressing EVs was proved in flow cytometry and ELISA-based experiments. Furthermore, by the use of nanoparticle tracking analysis and transmission electron microscopy, we found that the incubation of macrophage-released EVs with antigen-specific antibodies resulted in EVs' aggregation, which significantly enhanced their suppressive activity in vivo. Nowadays, it is increasingly evident that EVs play an exceptional role in intercellular communication and selective cargo transfer, and thus are considered promising candidates for therapeutic usage. However, EVs appear to be less effective than their parental cells. In this context, our current studies provide evidence that antigen-specific antibodies can be easily used for increasing EVs' biological activity, which has great therapeutic potential.

Application of information theoretical approaches to assess diversity and similarity in single-cell transcriptomics.

Single-cell transcriptomics offers a powerful way to reveal the heterogeneity of individual cells. To date, many information theoretical approaches have been proposed to assess diversity and similarity, and characterize the latent heterogeneity in transcriptome data. Diversity implies gene expression variations and can facilitate the identification of signature genes; while, similarity unravels co-expression patterns for cell type clustering. In this review, we summarized 16 measures of information theory used for evaluating diversity and similarity in single-cell transcriptomic data, provide references and shed light on selected theoretical properties when there is a need to select proper measurements in general cases. We further provide an R package assembling discussed approaches to improve the researchers own single-cell transcriptome study. At last, we prospected further applications of diversity and similarity measures in support of depicting heterogeneity in single-cell multi-omics data.

Mitochondrial GWAS and association of nuclear - mitochondrial epistasis with BMI in T1DM patients.

BACKGROUND: BMI is a strong indicator of complications from type I diabetes, especially under intensive treatment. METHODS: We have genotyped 435 type 1 diabetics using Illumina Infinium Omni Express Exome-8 v1.4 arrays and performed mitoGWAS on BMI. We identified additive interactions between mitochondrial and nuclear variants in genes associated with mitochondrial functioning MitoCarta2.0 and confirmed and refined the results on external cohorts: the Framingham Heart Study (FHS) and GTEx data. Linear mixed model analysis was performed using the GENESIS package in R/Bioconductor. RESULTS: We find a borderline significant association between the mitochondrial variant rs28357980, localized to MT-ND2, and BMI (ß = - 0.69, p = 0.056). This BMI association was confirmed on 1889 patients from FHS cohort (ß = - 0.312, p = 0.047). Next, we searched for additive interactions between mitochondrial and nuclear variants. MT-ND2 variants interacted with variants in the genes SIRT3, ATP5B, CYCS, TFB2M and POLRMT. TFB2M is a mitochondrial transcription factor and together with TFAM creates a transcription promoter complex for the mitochondrial polymerase POLRMT. We have found an interaction between rs3021088 in MT-ND2 and rs6701836 in TFB2M leading to BMI decrease (inter_pval = 0.0241), while interaction of rs3021088 in MT-ND2 and rs41542013 in POLRMT led to BMI increase (inter_pval = 0.0004). The influence of these interactions on BMI was confirmed in external cohorts. CONCLUSIONS: Here, we have shown that variants in the mitochondrial genome as well as additive interactions between mitochondrial and nuclear SNPs influence BMI in T1DM and general cohorts.