Knockdown of Hyaluronan synthase 2 suppresses liver fibrosis in mice via induction of transcriptomic changes similar to 4MU treatment.

Hepatic fibrosis remains a significant clinical challenge due to ineffective treatments. 4-methylumbelliferone (4MU), a hyaluronic acid (HA) synthesis inhibitor, has proven safe in phase one clinical trials. In this study, we aimed to ameliorate liver fibrosis by inhibiting HA synthesis. We compared two groups of mice with CCl4-induced fibrosis, treated with 4-methylumbelliferone (4MU) and hyaluronan synthase 2 (HAS2) targeting siRNA (siHAS2). The administration of 4MU and siHAS2 significantly reduced collagen and HA deposition, as well as biochemical markers of hepatic damage induced by repeated CCl4 injections. The transcriptomic analysis revealed converging pathways associated with downstream HA signalling. 4MU- and siHAS2-treated fibrotic livers shared 405 upregulated and 628 downregulated genes. These genes were associated with xenobiotic and cholesterol metabolism, mitosis, endoplasmic reticulum stress, RNA processing, and myeloid cell migration. The functional annotation of differentially expressed genes (DEGs) in siHAS2-treated mice revealed attenuation of extracellular matrix-associated pathways. In comparison, in the 4MU-treated group, DEGs were related to lipid and bile metabolism pathways and cell cycle. These findings confirm that HAS2 is an important pharmacological target for suppressing hepatic fibrosis using siRNA.

A comparative analysis of genes differentially expressed between rete testis cells and Sertoli cells of the mouse testis.

The rete testis (RT) is a region of the mammalian testis that plays an important role in testicular physiology. The RT epithelium consists of cells sharing some well-known gene markers with supporting Sertoli cells (SCs). However, little is known about the differences in gene expression between these two cell populations. Here, we used fluorescence-activated cell sorting (FACS) to obtain pure cultures of neonatal RT cells and SCs and identified differentially expressed genes (DEGs) between these cell types using RNA sequencing (RNA-seq). We then compared our data with the RNA-seq data of other studies that examined RT cells and SCs of mice of different ages and generated a list of DEGs permanently upregulated in RT cells throughout testis development and in culture, which included 86 genes, and a list of 79 DEGs permanently upregulated in SCs. The analysis of studies on DMRT1 function revealed that nearly half of the permanent DEGs could be regulated by this SC upregulated transcription factor. We suggest that useful cell lineage markers and candidate genes for the specification of both RT cells and SCs may be present among these permanent DEGs.

Complex Analysis of Single-Cell RNA Sequencing Data.

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

4-Methylumbelliferone Targets Revealed by Public Data Analysis and Liver Transcriptome Sequencing.

4-methylumbelliferone (4MU) is a well-known hyaluronic acid synthesis inhibitor and an approved drug for the treatment of cholestasis. In animal models, 4MU decreases inflammation, reduces fibrosis, and lowers body weight, serum cholesterol, and insulin resistance. It also inhibits tumor progression and metastasis. The broad spectrum of effects suggests multiple and yet unknown targets of 4MU. Aiming at 4MU target deconvolution, we have analyzed publicly available data bases, including: 1. Small molecule library Bio Assay screening (PubChemBioAssay); 2. GO pathway databases screening; 3. Protein Atlas Database. We also performed comparative liver transcriptome analysis of mice on normal diet and mice fed with 4MU for two weeks. Potential targets of 4MU public data base analysis fall into two big groups, enzymes and transcription factors (TFs), including 13 members of the nuclear receptor superfamily regulating lipid and carbohydrate metabolism. Transcriptome analysis revealed changes in the expression of genes involved in bile acid metabolism, gluconeogenesis, and immune response. It was found that 4MU feeding decreased the accumulation of the glycogen granules in the liver. Thus, 4MU has multiple targets and can regulate cell metabolism by modulating signaling via nuclear receptors.

Mechanisms of CP190 Interaction with Architectural Proteins in Drosophila Melanogaster.

Most of the known Drosophila architectural proteins interact with an important cofactor, CP190, that contains three domains (BTB, M, and D) that are involved in protein-protein interactions. The highly conserved N-terminal CP190 BTB domain forms a stable homodimer that interacts with unstructured regions in the three best-characterized architectural proteins: dCTCF, Su(Hw), and Pita. Here, we identified two new CP190 partners, CG4730 and CG31365, that interact with the BTB domain. The CP190 BTB resembles the previously characterized human BCL6 BTB domain, which uses its hydrophobic groove to specifically associate with unstructured regions of several transcriptional repressors. Using GST pull-down and yeast two-hybrid assays, we demonstrated that mutations in the hydrophobic groove strongly affect the affinity of CP190 BTB for the architectural proteins. In the yeast two-hybrid assay, we found that architectural proteins use various mechanisms to improve the efficiency of interaction with CP190. Pita and Su(Hw) have two unstructured regions that appear to simultaneously interact with hydrophobic grooves in the BTB dimer. In dCTCF and CG31365, two adjacent regions interact simultaneously with the hydrophobic groove of the BTB and the M domain of CP190. Finally, CG4730 interacts with the BTB, M, and D domains of CP190 simultaneously. These results suggest that architectural proteins use different mechanisms to increase the efficiency of interaction with CP190.

Quantitative Analysis of the Dynamics of Maternal Gradients in the Early Drosophila Embryo.

Predetermination, formation, and maintenance of the primary morphogenetic gradient (bicoid, bcd) of the early Drosophila embryo involves many interrelated processes. Here we focus on the biological systems analysis of the bcd mRNA redistribution in an early embryo. The results of the quantitative analysis of experimental data, together with the results of their dynamic modeling, substantiate the role of active transport in the redistribution of the bcd mRNA. The role of the nonlinearity of degradation mechanisms in the mRNA pattern robustness is discussed.

Mechanism and functional role of the interaction between CP190 and the architectural protein Pita in Drosophila melanogaster.

BACKGROUND: Pita is required for Drosophila development and binds specifically to a long motif in active promoters and insulators. Pita belongs to the Drosophila family of zinc-finger architectural proteins, which also includes Su(Hw) and the conserved among higher eukaryotes CTCF. The architectural proteins maintain the active state of regulatory elements and the long-distance interactions between them. In particular, Pita is involved in the formation of several boundaries between regulatory domains that controlled the expression of three hox genes in the Bithorax complex (BX-C). The CP190 protein is recruited to chromatin through interaction with the architectural proteins. RESULTS: Using in vitro pull-down analysis, we precisely mapped two unstructured regions of Pita that interact with the BTB domain of CP190. Then we constructed transgenic lines expressing the Pita protein of the wild-type and mutant variants lacking CP190-interacting regions. We have demonstrated that CP190-interacting region of the Pita can maintain nucleosome-free open chromatin and is critical for Pita-mediated enhancer blocking activity in BX-C. At the same time, interaction with CP190 is not required for the in vivo function of the mutant Pita protein, which binds to the same regions of the genome as the wild-type protein. Unexpectedly, we found that CP190 was still associated with the most of genome regions bound by the mutant Pita protein, which suggested that other architectural proteins were continuing to recruit CP190 to these regions. CONCLUSIONS: The results directly demonstrate role of CP190 in insulation and support a model in which the regulatory elements are composed of combinations of binding sites that interact with several architectural proteins with similar functions.

Single-cell-resolution transcriptome map of human, chimpanzee, bonobo, and macaque brains.

Identification of gene expression traits unique to the human brain sheds light on the molecular mechanisms underlying human evolution. Here, we searched for uniquely human gene expression traits by analyzing 422 brain samples from humans, chimpanzees, bonobos, and macaques representing 33 anatomical regions, as well as 88,047 cell nuclei composing three of these regions. Among 33 regions, cerebral cortex areas, hypothalamus, and cerebellar gray and white matter evolved rapidly in humans. At the cellular level, astrocytes and oligodendrocyte progenitors displayed more differences in the human evolutionary lineage than the neurons. Comparison of the bulk tissue and single-nuclei sequencing revealed that conventional RNA sequencing did not detect up to two-thirds of cell-type-specific evolutionary differences.

Distinct Elements Confer the Blocking and Bypass Functions of the Bithorax Fab-8 Boundary.

Boundaries in the Drosophila bithorax complex (BX-C) enable the regulatory domains that drive parasegment-specific expression of the three Hox genes to function autonomously. The four regulatory domains (iab-5, iab-6, iab-7, and iab-8) that control the expression of the Abdominal-B (Abd-B) gene are located downstream of the transcription unit, and are delimited by the Mcp, Fab-6, Fab-7, and Fab-8 boundaries. These boundaries function to block cross talk between neighboring regulatory domains. In addition, three of the boundaries (Fab-6, Fab-7, and Fab-8) must also have bypass activity so that regulatory domains distal to the boundaries can contact the Abd-B promoter. In the studies reported here, we have undertaken a functional dissection of the Fab-8 boundary using a boundary-replacement strategy. Our studies indicate that the Fab-8 boundary has two separable subelements. The distal subelement blocks cross talk, but cannot support bypass. The proximal subelement has only minimal blocking activity but is able to mediate bypass. A large multiprotein complex, the LBC (large boundary complex), binds to sequences in the proximal subelement and contributes to its bypass activity. The same LBC complex has been implicated in the bypass activity of the Fab-7 boundary.

Complete reconstitution of bypass and blocking functions in a minimal artificial Fab-7 insulator from Drosophila bithorax complex.

Boundaries in the bithorax complex (BX-C) delimit autonomous regulatory domains that drive parasegment-specific expression of the Hox genes Ubx, abd-A, and Abd-B The Fab-7 boundary is located between the iab-6 and iab-7 domains and has two key functions: blocking cross-talk between these domains and at the same time promoting communication (boundary bypass) between iab-6 and the Abd-B promoter. Using a replacement strategy, we found that multimerized binding sites for the architectural proteins Pita, Su(Hw), and dCTCF function as conventional insulators and block cross-talk between the iab-6 and iab-7 domains; however, they lack bypass activity, and iab-6 is unable to regulate Abd-B Here we show that an ∼200-bp sequence of dHS1 from the Fab-7 boundary rescues the bypass defects of these multimerized binding sites. The dHS1 sequence is bound in embryos by a large multiprotein complex, Late Boundary Complex (LBC), that contains the zinc finger proteins CLAMP and GAF. Using deletions and mutations in critical GAGAG motifs, we show that bypass activity correlates with the efficiency of recruitment of LBC components CLAMP and GAF to the artificial boundary. These results indicate that LBC orchestrates long-distance communication between the iab-6 regulatory domain and the Abd-B gene, while the Pita, Su(Hw), and dCTCF proteins function to block local cross-talk between the neighboring regulatory domains iab-6 and iab-7.

The bithorax complex iab-7 Polycomb response element has a novel role in the functioning of the Fab-7 chromatin boundary.

Expression of the three bithorax complex homeotic genes is orchestrated by nine parasegment-specific regulatory domains. Autonomy of each domain is conferred by boundary elements (insulators). Here, we have used an in situ replacement strategy to reanalyze the sequences required for the functioning of one of the best-characterized fly boundaries, Fab-7. It was initially identified by a deletion, Fab-71, that transformed parasegment (PS) 11 into a duplicate copy of PS12. Fab-71 deleted four nuclease hypersensitive sites, HS*, HS1, HS2, and HS3, located between the iab-6 and iab-7 regulatory domains. Transgenic and P-element excision experiments mapped the boundary to HS*+HS1+HS2, while HS3 was shown to be the iab-7 Polycomb response element (PRE). Recent replacement experiments showed that HS1 is both necessary and sufficient for boundary activity when HS3 is also present in the replacement construct. Surprisingly, while HS1+HS3 combination has full boundary activity, we discovered that HS1 alone has only minimal function. Moreover, when combined with HS3, only the distal half of HS1, dHS1, is needed. A ~1,000 kD multiprotein complex containing the GAF protein, called the LBC, binds to the dHS1 sequence and we show that mutations in dHS1, that disrupt LBC binding in nuclear extracts, eliminate boundary activity and GAF binding in vivo. HS3 has binding sites for GAF and Pho proteins that are required for PRE silencing. In contrast, HS3 boundary activity only requires the GAF binding sites. LBC binding with HS3 in nuclear extracts, and GAF association in vivo, depend upon the HS3 GAF sites, but not the Pho sites. Consistent with a role for the LBC in HS3 boundary activity, the boundary function of the dHS1+HS3mPho combination is lost when the flies are heterozygous for a mutation in the GAF gene. Taken together, these results reveal a novel function for the iab-7 PREs in chromosome architecture.

In silico evolution of the hunchback gene indicates redundancy in cis-regulatory organization and spatial gene expression.

Biological development depends on the coordinated expression of genes in time and space. Developmental genes have extensive cis-regulatory regions which control their expression. These regions are organized in a modular manner, with different modules controlling expression at different times and locations. Both how modularity evolved and what function it serves are open questions. We present a computational model for the cis-regulation of the hunchback (hb) gene in the fruit fly (Drosophila). We simulate evolution (using an evolutionary computation approach from computer science) to find the optimal cis-regulatory arrangements for fitting experimental hb expression patterns. We find that the cis-regulatory region tends to readily evolve modularity. These cis-regulatory modules (CRMs) do not tend to control single spatial domains, but show a multi-CRM/multi-domain correspondence. We find that the CRM-domain correspondence seen in Drosophila evolves with a high probability in our model, supporting the biological relevance of the approach. The partial redundancy resulting from multi-CRM control may confer some biological robustness against corruption of regulatory sequences. The technique developed on hb could readily be applied to other multi-CRM developmental genes.

In silico evolution of gene cooption in pattern-forming gene networks.

Gene recruitment or cooption occurs when a gene, which may be part of an existing gene regulatory network (GRN), comes under the control of a new regulatory system. Such re-arrangement of pre-existing networks is likely more common for increasing genomic complexity than the creation of new genes. Using evolutionary computations (EC), we investigate how cooption affects the evolvability, outgrowth and robustness of GRNs. We use a data-driven model of insect segmentation, for the fruit fly Drosophila, and evaluate fitness by robustness to maternal variability-a major constraint in biological development. We compare two mechanisms of gene cooption: a simpler one with gene Introduction and Withdrawal operators; and one in which GRN elements can be altered by transposon infection. Starting from a minimal 2-gene network, insufficient for fitting the Drosophila gene expression patterns, we find a general trend of coopting available genes into the GRN, in order to better fit the data. With the transposon mechanism, we find co-evolutionary oscillations between genes and their transposons. These oscillations may offer a new technique in EC for overcoming premature convergence. Finally, we comment on how a differential equations (in contrast to Boolean) approach is necessary for addressing realistic continuous variation in biochemical parameters.