TFTenricher: a python toolbox for annotation enrichment analysis of transcription factor target genes

Background Transcription factors (TFs) are the upstream regulators that orchestrate gene expression, and therefore a centrepiece in bioinformatics studies. While a core strategy to understand the biological context of genes and proteins includes annotation enrichment analysis, such as Gene Ontology term enrichment, these methods are not well suited for analysing groups of TFs. This is particularly true since such methods do not aim to include downstream processes, and given a set of TFs, the expected top ontologies would revolve around transcription processes. Results We present the TFTenricher, a Python toolbox that focuses specifically at identifying gene ontology terms, cellular pathways, and diseases that are over-represented among genes downstream of user-defined sets of human TFs. We evaluated the inference of downstream gene targets with respect to false positive annotations, and found an inference based on co-expression to best predict downstream processes. Based on these downstream genes, the TFTenricher uses some of the most common databases for gene functionalities, including GO, KEGG and Reactome, to calculate functional enrichments. By applying the TFTenricher to differential expression of TFs in 21 diseases, we found significant terms associated with disease mechanism, while the gene set enrichment analysis on the same dataset predominantly identified processes related to transcription. Conclusions and availability The TFTenricher package enables users to search for biological context in any set of TFs and their downstream genes. The TFTenricher is available as a Python 3 toolbox at https://github.com/rasma774/Tftenricher, under a GNU GPL license and with minimal dependencies.

Early white matter pathology in the fornix of the limbic system in Huntington disease

Huntington disease (HD) is a fatal neurodegenerative disorder caused by an expanded CAG repeat in the huntingtin (HTT) gene. The typical motor symptoms have been associated with basal ganglia pathology. However, psychiatric and cognitive symptoms often precede the motor component and may be due to changes in the limbic system. Recent work has indicated pathology in the hypothalamus in HD but other parts of the limbic system have not been extensively studied. Emerging evidence suggests that changes in HD also include white matter pathology. Here we investigated if the main white matter tract of the limbic system, the fornix, is affected in HD. We demonstrate that the fornix is 34% smaller already in prodromal HD and 41% smaller in manifest HD compared to controls using volumetric analyses of MRI of the IMAGE-HD study. In post-mortem fornix tissue from HD cases, we confirm the smaller fornix volume in HD which is accompanied by signs of myelin breakdown and reduced levels of the transcription factor myelin regulating factor but detect no loss of oligodendrocytes. Further analyses using RNA-sequencing demonstrate downregulation of oligodendrocyte identity markers in the fornix of HD cases. Analysis of differentially expressed genes based on transcription-factor/target-gene interactions also revealed enrichment for binding sites of SUZ12 and EZH2, components of the Polycomb Repressive Complex 2, as well as RE1 Regulation Transcription Factor. Taken together, our data show that there is early white matter pathology of the fornix in the limbic system in HD likely due to a combination of reduction in oligodendrocyte genes and myelin break down.

ACE2 protein expression within isogenic cell lines is heterogeneous and associated with distinct transcriptomes

The membrane protein angiotensin-converting enzyme 2 (ACE2) is a physiologic regulator of the renin-angiotensin system and the cellular receptor for the SARS-CoV-2 virus. Prior studies of ACE2 expression have primarily focused on mRNA abundance, with investigation at the protein level limited by uncertain specificity of commercial ACE2 antibodies. Here, we report our development of a sensitive and specific flow cytometry-based assay for cellular ACE2 protein abundance. Application of this approach to multiple cell lines revealed an unexpected degree of cellular heterogeneity, with detectable ACE2 protein in only a subset of cells in each isogenic population. This heterogeneity was mediated at the mRNA level by transcripts predominantly initiated from the ACE2 proximal promoter. ACE2 expression was heritable but not fixed over multiple generations of daughter cells, with gradual drift toward the original heterogeneous background. RNA-seq profiling identified distinct transcriptomes of ACE2-expressing relative cells to non-expressing cells, with enrichment in functionally related genes and transcription factor target sets. Our findings provide a validated approach for the specific detection of ACE2 protein at the surface of single cells, support an epigenetic mechanism of ACE2 gene regulation, and identify specific pathways associated with ACE2 expression in HuH7 cells.

Autophagy-Related Genes in Atherosclerosis

Background. Atherosclerosis (AS) is a common chronic vascular inflammatory disease and one of the main causes of cardiovascular/cerebrovascular diseases (CVDs). Autophagy-related genes (ARGs) play a crucial part in pathophysiological processes of AS. However, the expression profile of ARGs has rarely been adopted to explore the relationship between autophagy and AS. Therefore, using the expression profile of ARGs to explore the relationship between autophagy and AS may provide new insights for the treatment of CVDs. Methods. The differentially expressed ARGs of the GSE57691 dataset were obtained from the Human Autophagy Database (HADb) and the Gene Expression Omnibus (GEO) database, and the GSE57691 dataset contains 9 aortic atheroma tissues and 10 normal aortic tissues. The differentially expressed ARGs of the GSE57691 dataset were analyzed by protein-protein interaction (PPI), gene ontology analysis (GO), and Kyoto Encyclopedia of Genes and Genomes analysis (KEGG) and were chosen to explore related miRNAs/transcriptional factors. Results. The GSE57691 dataset had a total of 41 differentially expressed ARGs. The GO analysis results revealed that ARGs were mainly enriched in autophagy, autophagosome, and protein serine/threonine kinase activity. KEGG analysis results showed that ARGs were mainly enriched in autophagy-animal and longevity regulating signaling pathways. Expressions of ATG5, MAP1LC3B, MAPK3, MAPK8, and RB1CC1 were regarded as focus in the PPI regulatory networks. Furthermore, 11 related miRNAs and 6 related transcription factors were obtained by miRNAs/transcription factor target network analysis. Conclusions. Autophagy and ARGs may play a vital role in regulating the pathophysiology of AS.

Phenotypic Robustness of Epidermal Stem Cell Number in C. elegans Is Modulated by the Activity of the Conserved N-acetyltransferase nath-10/NAT10

Individual cells and organisms experience perturbations from internal and external sources, yet manage to buffer these to produce consistent phenotypes, a property known as robustness. While phenotypic robustness has often been examined in unicellular organisms, it has not been sufficiently studied in multicellular animals. Here, we investigate phenotypic robustness in Caenorhabditis elegans seam cells. Seam cells are stem cell-like epithelial cells along the lateral edges of the animal, which go through asymmetric and symmetric divisions contributing cells to the hypodermis and neurons, while replenishing the stem cell reservoir. The terminal number of seam cells is almost invariant in the wild-type population, allowing the investigation of how developmental precision is achieved. We report here that a loss-of-function mutation in the highly conserved N-acetyltransferase nath-10/NAT10 increases seam cell number variance in the isogenic population. RNA-seq analysis revealed increased levels of mRNA transcript variability in nath-10 mutant populations, which may have an impact on the phenotypic variability observed. Furthermore, we found disruption of Wnt signaling upon perturbing nath-10 function, as evidenced by changes in POP-1/TCF nuclear distribution and ectopic activation of its GATA transcription factor target egl-18. These results highlight that NATH-10/NAT-10 can influence phenotypic variability partly through modulation of the Wnt signaling pathway.

ACE2 protein expression within isogenic cell lines is heterogeneous and associated with distinct transcriptomes

The membrane protein angiotensin-converting enzyme 2 (ACE2) is a physiologic regulator of the renin-angiotensin system and the cellular receptor for the SARS-CoV-2 virus. Prior studies of ACE2 expression have primarily focused on mRNA abundance, with investigation at the protein level limited by uncertain specificity of commercial ACE2 antibodies. Here, we report our development of a sensitive and specific flow cytometry-based assay for cellular ACE2 protein abundance. Application of this approach to multiple cell lines revealed an unexpected degree of cellular heterogeneity, with detectable ACE2 protein in only a subset of cells in each isogenic population. This heterogeneity was mediated at the mRNA level by transcripts predominantly initiated from the ACE2 proximal promoter. ACE2 expression was heritable but not fixed over multiple generations of daughter cells, with gradual drift toward the original heterogeneous background. RNA-seq profiling identified distinct transcriptomes of ACE2-expressing relative cells to non-expressing cells, with enrichment in functionally related genes and transcription factor target sets. Our findings provide a validated approach for the specific detection of ACE2 protein at the surface of single cells, support an epigenetic mechanism ACE2 gene regulation, and identify specific pathways associated with ACE2 expression in HuH7 cells.

Identification of Pivotal MicroRNAs and Target Genes Associated with Persistent Atrial Fibrillation Based on Bioinformatics Analysis

Atrial fibrillation (AF) is one of the most common supraventricular arrhythmias worldwide. However, the specific molecular mechanism underlying AF remains unclear. Our study is aimed at identifying pivotal microRNAs (miRNAs) and targeting genes associated with persistent AF (pAF) using bioinformatics analysis. Three gene expression array datasets (GSE31821, GSE41177, and GSE79768) and an miRNA expression array dataset (GSE68475) associated with pAF were downloaded. Differentially expressed genes (DEGs) were identified using the LIMMA package, and differentially expressed miRNAs (DEMs) were screened from GSE68475. Target genes for DEMs were predicted using the miRTarBase database, and intersections between these target genes and DEGs were selected for further analysis, including the generation of protein–protein interaction (PPI) network, miRNA–transcription factor–target regulatory network, and drug–gene network. A total of 264 DEGs and 40 DEMs were identified between the pAF and control groups. Functional and pathway enrichment analyses of up- and downregulated DEGs were performed. The common genes (CGs) were primarily enriched in the phosphoinositide 3-kinase- (PI3K-) protein kinase B (Akt) signaling pathway, negative regulation of cell division, and response to hypoxia. The PPI network, miRNA–transcription factor–target regulatory network, and drug–gene network were constructed using Cytoscape. The present study revealed several novel miRNAs and genes involved in pAF. We speculated that miR-4298, miR-3125, miR-4306, and miR-671-5p could represent significant miRNAs that act on the target gene superoxide dismutase 2 (SOD2) during the development of pAF and may serve as essential biomarkers for pAF diagnosis and treatment. Moreover, MYC might function in pAF pathogenesis through the PI3K–Akt signaling pathway.

Identification of bacterial virulence factors based on an integration of experimental and in silico transcription factor target discovery

In silico transcription factors binding sites analysis in Dickeya dadantii genome matches published experimental results and predicts new regulators for genes inducing soft-rot disease in plants.

Tissue-specific transcription factor target identification in the Caenorhabditis elegans epidermis using targeted DamID

SummaryTranscription factors are key orchestrators of development in multicellular animals and display complex patterns of expression, as well as tissue-specific binding to targets. However, our ability to map transcription factor-target interactions in specific tissues of intact animals remains limited. We introduce here targeted DamID (TaDa) as a method to identify transcription factor targets with tissue-specific resolution in C. elegans. We focus on the epidermis as a paradigm and demonstrate that TaDa circumvents problems with Dam-associated toxicity and allows reproducible identification of putative targets. Using a combination of TaDa and single-molecule FISH (smFISH), we refine the positions of LIN-22 and NHR-25 within the epidermal gene network. We reveal direct links between these two factors and the cell differentiation programme, as well as the Wnt signalling pathway. Our results illustrate how TaDa and smFISH can be used to dissect the architecture of tissue-specific gene regulatory networks.HighlightsTaDa circumvents Dam-associated toxicity by keeping levels of Dam expression low.TaDa allows the recovery of tissue-specific methylation profiles representing TF binding.Methylation signal is enriched in regulatory regions of the genome.LIN-22 and NHR-25 targets reveal a link to cell differentiation and Wnt signalling.