A cellular and spatial map of the choroid plexus across brain ventricles and ages.

The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1ß (IL-1ß) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.

Structural and functional basis of mammalian microRNA biogenesis by Dicer.

MicroRNA (miRNA) and RNA interference (RNAi) pathways rely on small RNAs produced by Dicer endonucleases. Mammalian Dicer primarily supports the essential gene-regulating miRNA pathway, but how it is specifically adapted to miRNA biogenesis is unknown. We show that the adaptation entails a unique structural role of Dicer's DExD/H helicase domain. Although mice tolerate loss of its putative ATPase function, the complete absence of the domain is lethal because it assures high-fidelity miRNA biogenesis. Structures of murine Dicer•-miRNA precursor complexes revealed that the DExD/H domain has a helicase-unrelated structural function. It locks Dicer in a closed state, which facilitates miRNA precursor selection. Transition to a cleavage-competent open state is stimulated by Dicer-binding protein TARBP2. Absence of the DExD/H domain or its mutations unlocks the closed state, reduces substrate selectivity, and activates RNAi. Thus, the DExD/H domain structurally contributes to mammalian miRNA biogenesis and underlies mechanistical partitioning of miRNA and RNAi pathways.

Autoimmune amelogenesis imperfecta in patients with APS-1 and coeliac disease.

Ameloblasts are specialized epithelial cells in the jaw that have an indispensable role in tooth enamel formation-amelogenesis1. Amelogenesis depends on multiple ameloblast-derived proteins that function as a scaffold for hydroxyapatite crystals. The loss of function of ameloblast-derived proteins results in a group of rare congenital disorders called amelogenesis imperfecta2. Defects in enamel formation are also found in patients with autoimmune polyglandular syndrome type-1 (APS-1), caused by AIRE deficiency3,4, and in patients diagnosed with coeliac disease5-7. However, the underlying mechanisms remain unclear. Here we show that the vast majority of patients with APS-1 and coeliac disease develop autoantibodies (mostly of the IgA isotype) against ameloblast-specific proteins, the expression of which is induced by AIRE in the thymus. This in turn results in a breakdown of central tolerance, and subsequent generation of corresponding autoantibodies that interfere with enamel formation. However, in coeliac disease, the generation of such autoantibodies seems to be driven by a breakdown of peripheral tolerance to intestinal antigens that are also expressed in enamel tissue. Both conditions are examples of a previously unidentified type of IgA-dependent autoimmune disorder that we collectively name autoimmune amelogenesis imperfecta.

Bispecific IgG neutralizes SARS-CoV-2 variants and prevents escape in mice.

Neutralizing antibodies that target the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein are among the most promising approaches against COVID-191,2. A bispecific IgG1-like molecule (CoV-X2) has been developed on the basis of C121 and C135, two antibodies derived from donors who had recovered from COVID-193. Here we show that CoV-X2 simultaneously binds two independent sites on the RBD and, unlike its parental antibodies, prevents detectable spike binding to the cellular receptor of the virus, angiotensin-converting enzyme 2 (ACE2). Furthermore, CoV-X2 neutralizes wild-type SARS-CoV-2 and its variants of concern, as well as escape mutants generated by the parental monoclonal antibodies. We also found that in a mouse model of SARS-CoV-2 infection with lung inflammation, CoV-X2 protects mice from disease and suppresses viral escape. Thus, the simultaneous targeting of non-overlapping RBD epitopes by IgG-like bispecific antibodies is feasible and effective, and combines the advantages of antibody cocktails with those of single-molecule approaches.

Ovol2 promoter mutations in mice and human illuminate species-specific phenotypic divergence.

Pathogenic variants in the highly conserved OVOL2 promoter region cause posterior polymorphous corneal dystrophy (PPCD) 1 by inducing an ectopic expression of the endothelial OVOL2 mRNA. Here we produced an allelic series of Ovol2 promoter mutations in the mouse model including the heterozygous c.-307T>C variant (RefSeq NM_021220.4) causing PPCD1 in humans. Despite the high evolutionary conservation of the Ovol2 promoter, only some alterations of its sequence had phenotypic consequences in mice. Four independent sequence variants in the distal part of the Ovol2 promoter had no significant effect on endothelial Ovol2 mRNA level or caused any ocular phenotype. In contrast, the mutation c.-307T>C resulted in increased Ovol2 expression in the corneal endothelium. However, only a small fraction of adult mice c.-307T>C heterozygotes developed ocular phenotypes such as irido-corneal adhesions, and corneal opacity. Interestingly, phenotypic penetrance was increased at embryonic stages. Notably, c.-307T>C mutation is located next to the Ovol1/Ovol2 transcription factor binding site. Mice carrying an allele with a deletion encompassing the Ovol2 binding site c.-307_-320del showed significant Ovol2 gene upregulation in the cornea endothelium and exhibited phenotypes similar to the c.-307T>C mutation. In conclusion, although the mutations c.-307T>C and -307_-320del lead to a comparably strong increase in endothelial Ovol2 expression as seen in PPCD1 patients, endothelial dystrophy was not observed in the mouse model, implicating species-specific differences in endothelial cell biology. Nonetheless, the emergence of dominant ocular phenotypes associated with Ovol2 promoter variants in mice implies a potential role of this gene in eye development and disease.

ABIN1 is a negative regulator of effector functions in cytotoxic T cells.

T cells are pivotal in the adaptive immune defense, necessitating a delicate balance between robust response against infections and self-tolerance. Their activation involves intricate cross-talk among signaling pathways triggered by the T-cell antigen receptors (TCR) and co-stimulatory or inhibitory receptors. The molecular regulation of these complex signaling networks is still incompletely understood. Here, we identify the adaptor protein ABIN1 as a component of the signaling complexes of GITR and OX40 co-stimulation receptors. T cells lacking ABIN1 are hyper-responsive ex vivo, exhibit enhanced responses to cognate infections, and superior ability to induce experimental autoimmune diabetes in mice. ABIN1 negatively regulates p38 kinase activation and late NF-κB target genes. P38 is at least partially responsible for the upregulation of the key effector proteins IFNG and GZMB in ABIN1-deficient T cells after TCR stimulation. Our findings reveal the intricate role of ABIN1 in T-cell regulation.

Functional canonical RNAi in mice expressing a truncated Dicer isoform and long dsRNA.

Canonical RNA interference (RNAi) is sequence-specific mRNA degradation guided by small interfering RNAs (siRNAs) made by RNase III Dicer from long double-stranded RNA (dsRNA). RNAi roles include gene regulation, antiviral immunity or defense against transposable elements. In mammals, RNAi is constrained by Dicer's adaptation to produce another small RNA class-microRNAs. However, a truncated Dicer isoform (ΔHEL1) supporting RNAi exists in mouse oocytes. A homozygous mutation to express only the truncated ΔHEL1 variant causes dysregulation of microRNAs and perinatal lethality in mice. Here, we report the phenotype and canonical RNAi activity in DicerΔHEL1/wt mice, which are viable, show minimal miRNome changes, but their endogenous siRNA levels are an order of magnitude higher. We show that siRNA production in vivo is limited by available dsRNA, but not by Protein kinase R, a dsRNA sensor of innate immunity. dsRNA expression from a transgene yields sufficient siRNA levels to induce efficient RNAi in heart and muscle. DicerΔHEL1/wt mice with enhanced canonical RNAi offer a platform for examining potential and limits of mammalian RNAi in vivo.

Characterization of mice carrying a neurodevelopmental disease-associated GluN2B(L825V) variant.

N-methyl-D-aspartate receptors (NMDARs), encoded by GRIN genes, are ionotropic glutamate receptors playing a critical role in synaptic transmission, plasticity, and synapse development. Genome sequence analyses have identified variants in GRIN genes in patients with neurodevelopmental disorders, but the underlying disease mechanisms are not well understood. Here, we have created and evaluated a transgenic mouse line carrying a missense variant Grin2bL825V , corresponding to a de-novo GRIN2B variant encoding GluN2B(L825V) found in a patient with intellectual disability (ID) and autism spectrum disorder (ASD). We used HEK293T cells expressing recombinant receptors and primary hippocampal neurons prepared from heterozygous Grin2bL825V/+ (L825V/+) and wild-type Grin2b+/+ (+/+) male and female mice to assess the functional impact of the variant. Whole-cell NMDAR currents were reduced in neurons prepared from L825V/+ compared to +/+ mice. Peak amplitude of NMDAR-mediated evoked excitatory postsynaptic currents (NMDAR-eEPSC) was not changed, but NMDAR-eEPSCs in L825V/+ neurons had faster deactivation compared to +/+ neurons and were less sensitive to a GluN2B-selective antagonist ifenprodil. Together, these results suggest a decreased functional contribution of GluN2B subunits to synaptic NMDAR currents in hippocampal neurons from L825V/+ mice. The analysis of the GluN2B(L825V) subunit surface expression and synaptic localization revealed no differences compared to wild-type GluN2B. Behavioral testing of mice of both sexes demonstrated hypoactivity, anxiety, and impaired sensorimotor gating in the L825V/+ strain, particularly affecting males, as well as cognitive symptoms. The heterozygous L825V/+ mouse offers a clinically relevant model of GRIN2B-related ID/ASD and our results suggest synaptic-level functional changes that may contribute to neurodevelopmental pathology.Significance statement Variants in genes for subunits of N-methyl-D-aspartate receptors (NMDARs), a subtype of ionotropic glutamate receptors, are associated with neurodevelopmental disorders. Here we have generated a transgenic mouse model of a de-novo missense GRIN2B gene variant, identified in a patient with intellectual disability and autism, that introduces a single amino acid substitution (L825V) in the NMDAR GluN2B subunit. Di- and triheteromeric NMDARs containing the GluN2B(L825V) subunit have a reduced channel open probability. Synaptic NMDAR currents in neurons from heterozygous L825V/+ mice have accelerated deactivation and reduced ifenprodil sensitivity, suggesting synaptic loss of GluN2B function. L825V/+ mice show increased anxiety, impaired sensorimotor gating, and cognitive deficits, consistent with patient symptoms. Our study describes a clinically relevant mouse model of GRIN2B-related neurodevelopmental pathology.

Expression of human-specific ARHGAP11B in mice leads to neocortex expansion and increased memory flexibility.

Neocortex expansion during human evolution provides a basis for our enhanced cognitive abilities. Yet, which genes implicated in neocortex expansion are actually responsible for higher cognitive abilities is unknown. The expression of human-specific ARHGAP11B in embryonic/foetal mouse, ferret and marmoset neocortex was previously found to promote basal progenitor proliferation, upper-layer neuron generation and neocortex expansion during development, features commonly thought to contribute to increased cognitive abilities. However, a key question is whether this phenotype persists into adulthood and if so, whether cognitive abilities are indeed increased. Here, we generated a transgenic mouse line with physiological ARHGAP11B expression that exhibits increased neocortical size and upper-layer neuron numbers persisting into adulthood. Adult ARHGAP11B-transgenic mice showed altered neurobehaviour, notably increased memory flexibility and a reduced anxiety level. Our data are consistent with the notion that neocortex expansion by ARHGAP11B, a gene implicated in human evolution, underlies some of the altered neurobehavioural features observed in the transgenic mice, such as the increased memory flexibility, a neocortex-associated trait, with implications for the increase in cognitive abilities during human evolution.

Early perturbation of Wnt signaling reveals patterning and invagination-evagination control points in molar tooth development.

Tooth formation requires complex signaling interactions both within the oral epithelium and between the epithelium and the underlying mesenchyme. Previous studies of the Wnt/ß-catenin pathway have shown that tooth formation is partly inhibited in loss-of-function mutants, and gain-of-function mutants have perturbed tooth morphology. However, the stage at which Wnt signaling is first important in tooth formation remains unclear. Here, using an Fgf8-promoter-driven, and therefore early, deletion of ß-catenin in mouse molar epithelium, we found that loss of Wnt/ß-catenin signaling completely deletes the molar tooth, demonstrating that this pathway is central to the earliest stages of tooth formation. Early expression of a dominant-active ß-catenin protein also perturbs tooth formation, producing a large domed evagination at early stages and supernumerary teeth later on. The early evaginations are associated with premature mesenchymal condensation marker, and are reduced by inhibition of condensation-associated collagen synthesis. We propose that invagination versus evagination morphogenesis is regulated by the relative timing of epithelial versus mesenchymal cell convergence regulated by canonical Wnt signaling. Together, these studies reveal new aspects of Wnt/ß-catenin signaling in tooth formation and in epithelial morphogenesis more broadly.

MEIS-WNT5A axis regulates development of fourth ventricle choroid plexus.

The choroid plexus (ChP) produces cerebrospinal fluid and forms an essential brain barrier. ChP tissues form in each brain ventricle, each one adopting a distinct shape, but remarkably little is known about the mechanisms underlying ChP development. Here, we show that epithelial WNT5A is crucial for determining fourth ventricle (4V) ChP morphogenesis and size in mouse. Systemic Wnt5a knockout, or forced Wnt5a overexpression beginning at embryonic day 10.5, profoundly reduced ChP size and development. However, Wnt5a expression was enriched in Foxj1-positive epithelial cells of 4V ChP plexus, and its conditional deletion in these cells affected the branched, villous morphology of the 4V ChP. We found that WNT5A was enriched in epithelial cells localized to the distal tips of 4V ChP villi, where WNT5A acted locally to activate non-canonical WNT signaling via ROR1 and ROR2 receptors. During 4V ChP development, MEIS1 bound to the proximal Wnt5a promoter, and gain- and loss-of-function approaches demonstrated that MEIS1 regulated Wnt5a expression. Collectively, our findings demonstrate a dual function of WNT5A in ChP development and identify MEIS transcription factors as upstream regulators of Wnt5a in the 4V ChP epithelium.

Mitochondrial HER2 stimulates respiration and promotes tumorigenicity.

BACKGROUND: Amplification of HER2, a receptor tyrosine kinase and a breast cancer-linked oncogene, is associated with aggressive disease. HER2 protein is localised mostly at the cell membrane, but a fraction translocates to mitochondria. Whether and how mitochondrial HER2 contributes to tumorigenicity is currently unknown. METHODS: We enriched the mitochondrial (mt-)HER2 fraction in breast cancer cells using an N-terminal mitochondrial targeting sequence and analysed how this manipulation impacts bioenergetics and tumorigenic properties. The role of the tyrosine kinase activity of mt-HER2 was assessed in wild type, kinase-dead (K753M) and kinase-enhanced (V659E) mtHER2 constructs. RESULTS: We document that mt-HER2 associates with the oxidative phosphorylation system, stimulates bioenergetics and promotes larger respiratory supercomplexes. mt-HER2 enhances proliferation and invasiveness in vitro and tumour growth and metastatic potential in vivo, in a kinase activity-dependent manner. On the other hand, constitutively active mt-HER2 provokes excessive mitochondria ROS generation, sensitises to cell death, and restricts growth of primary tumours, suggesting that regulation of HER2 activity in mitochondria is required for the maximal pro-tumorigenic effect. CONCLUSIONS: mt-HER2 promotes tumorigenicity by supporting bioenergetics and optimal redox balance.

A review of standardized high-throughput cardiovascular phenotyping with a link to metabolism in mice.

Cardiovascular diseases cause a high mortality rate worldwide and represent a major burden for health care systems. Experimental rodent models play a central role in cardiovascular disease research by effectively simulating human cardiovascular diseases. Using mice, the International Mouse Phenotyping Consortium (IMPC) aims to target each protein-coding gene and phenotype multiple organ systems in single-gene knockout models by a global network of mouse clinics. In this review, we summarize the current advances of the IMPC in cardiac research and describe in detail the diagnostic requirements of high-throughput electrocardiography and transthoracic echocardiography capable of detecting cardiac arrhythmias and cardiomyopathies in mice. Beyond that, we are linking metabolism to the heart and describing phenotypes that emerge in a set of known genes, when knocked out in mice, such as the leptin receptor (Lepr), leptin (Lep), and Bardet-Biedl syndrome 5 (Bbs5). Furthermore, we are presenting not yet associated loss-of-function genes affecting both, metabolism and the cardiovascular system, such as the RING finger protein 10 (Rfn10), F-box protein 38 (Fbxo38), and Dipeptidyl peptidase 8 (Dpp8). These extensive high-throughput data from IMPC mice provide a promising opportunity to explore genetics causing metabolic heart disease with an important translational approach.

Comprehensive ECG reference intervals in C57BL/6N substrains provide a generalizable guide for cardiac electrophysiology studies in mice.

Reference ranges provide a powerful tool for diagnostic decision-making in clinical medicine and are enormously valuable for understanding normality in pre-clinical scientific research that uses in vivo models. As yet, there are no published reference ranges for electrocardiography (ECG) in the laboratory mouse. The first mouse-specific reference ranges for the assessment of electrical conduction are reported herein generated from an ECG dataset of unprecedented scale. International Mouse Phenotyping Consortium data from over 26,000 conscious or anesthetized C57BL/6N wildtype control mice were stratified by sex and age to develop robust ECG reference ranges. Interesting findings include that heart rate and key elements from the ECG waveform (RR-, PR-, ST-, QT-interval, QT corrected, and QRS complex) demonstrate minimal sexual dimorphism. As expected, anesthesia induces a decrease in heart rate and was shown for both inhalation (isoflurane) and injectable (tribromoethanol) anesthesia. In the absence of pharmacological, environmental, or genetic challenges, we did not observe major age-related ECG changes in C57BL/6N-inbred mice as the differences in the reference ranges of 12-week-old compared to 62-week-old mice were negligible. The generalizability of the C57BL/6N substrain reference ranges was demonstrated by comparison with ECG data from a wide range of non-IMPC studies. The close overlap in data from a wide range of mouse strains suggests that the C57BL/6N-based reference ranges can be used as a robust and comprehensive indicator of normality. We report a unique ECG reference resource of fundamental importance for any experimental study of cardiac function in mice.

Batch alignment via retention orders for preprocessing large-scale multi-batch LC-MS experiments.

MOTIVATION: Meticulous selection of chromatographic peak detection parameters and algorithms is a crucial step in preprocessing liquid chromatography-mass spectrometry (LC-MS) data. However, as mass-to-charge ratio and retention time shifts are larger between batches than within batches, finding apt parameters for all samples of a large-scale multi-batch experiment with the aim of minimizing information loss becomes a challenging task. Preprocessing independent batches individually can curtail said problems but requires a method for aligning and combining them for further downstream analysis. RESULTS: We present two methods for aligning and combining individually preprocessed batches in multi-batch LC-MS experiments. Our developed methods were tested on six sets of simulated and six sets of real datasets. Furthermore, by estimating the probabilities of peak insertion, deletion and swap between batches in authentic datasets, we demonstrate that retention order swaps are not rare in untargeted LC-MS data. AVAILABILITY AND IMPLEMENTATION: kmersAlignment and rtcorrectedAlignment algorithms are made available as an R package with raw data at https://metabocombiner.img.cas.cz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Bardet-Biedl Syndrome ciliopathy is linked to altered hematopoiesis and dysregulated self-tolerance.

Bardet-Biedl Syndrome (BBS) is a pleiotropic genetic disease caused by the dysfunction of primary cilia. The immune system of patients with ciliopathies has not been investigated. However, there are multiple indications that the impairment of the processes typically associated with cilia may have influence on the hematopoietic compartment and immunity. In this study, we analyze clinical data of BBS patients and corresponding mouse models carrying mutations in Bbs4 or Bbs18. We find that BBS patients have a higher prevalence of certain autoimmune diseases. Both BBS patients and animal models have altered red blood cell and platelet compartments, as well as elevated white blood cell levels. Some of the hematopoietic system alterations are associated with BBS-induced obesity. Moreover, we observe that the development and homeostasis of B cells in mice is regulated by the transport complex BBSome, whose dysfunction is a common cause of BBS. The BBSome limits canonical WNT signaling and increases CXCL12 levels in bone marrow stromal cells. Taken together, our study reveals a connection between a ciliopathy and dysregulated immune and hematopoietic systems.

ATF2 loss promotes tumor invasion in colorectal cancer cells via upregulation of cancer driver TROP2.

In cancer, the activating transcription factor 2 (ATF2) has pleiotropic functions in cellular responses to growth stimuli, damage, or inflammation. Due to only limited studies, the significance of ATF2 in colorectal cancer (CRC) is not well understood. We report that low ATF2 levels correlated with worse prognosis and tumor aggressiveness in CRC patients. NanoString gene expression and ChIP analysis confirmed trophoblast cell surface antigen 2 (TROP2) as a novel inhibitory ATF2 target gene. This inverse correlation was further observed in primary human tumor tissues. Immunostainings revealed that high intratumoral heterogeneity for ATF2 and TROP2 expression was sustained also in liver metastasis. Mechanistically, our in vitro data of CRISPR/Cas9-generated ATF2 knockout (KO) clones revealed that high TROP2 levels were critical for cell de-adhesion and increased cell migration without triggering EMT. TROP2 was enriched in filopodia and displaced Paxillin from adherens junctions. In vivo imaging, micro-computer tomography, and immunostainings verified that an ATF2KO/TROP2high status triggered tumor invasiveness in in vivo mouse and chicken xenograft models. In silico analysis provided direct support that ATF2low/TROP2high expression status defined high-risk CRC patients. Finally, our data demonstrate that ATF2 acts as a tumor suppressor by inhibiting the cancer driver TROP2. Therapeutic TROP2 targeting might prevent particularly the first steps in metastasis, i.e., the de-adhesion and invasion of colon cancer cells.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density.

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

The receptor-type protein tyrosine phosphatase CD45 promotes onset and severity of IL-1ß-mediated autoinflammatory osteomyelitis.

A number of human autoinflammatory diseases manifest with severe inflammatory bone destruction. Mouse models of these diseases represent valuable tools that help us to understand molecular mechanisms triggering this bone autoinflammation. The Pstpip2cmo mouse strain is among the best characterized of these; it harbors a mutation resulting in the loss of adaptor protein PSTPIP2 and development of autoinflammatory osteomyelitis. In Pstpip2cmo mice, overproduction of interleukin-1ß (IL-1ß) and reactive oxygen species by neutrophil granulocytes leads to spontaneous inflammation of the bones and surrounding soft tissues. However, the upstream signaling events leading to this overproduction are poorly characterized. Here, we show that Pstpip2cmo mice deficient in major regulator of Src-family kinases (SFKs) receptor-type protein tyrosine phosphatase CD45 display delayed onset and lower severity of the disease, while the development of autoinflammation is not affected by deficiencies in Toll-like receptor signaling. Our data also show deregulation of pro-IL-1ß production by Pstpip2cmo neutrophils that are attenuated by CD45 deficiency. These data suggest a role for SFKs in autoinflammation. Together with previously published work on the involvement of protein tyrosine kinase spleen tyrosine kinase, they point to the role of receptors containing immunoreceptor tyrosine-based activation motifs, which after phosphorylation by SFKs recruit spleen tyrosine kinase for further signal propagation. We propose that this class of receptors triggers the events resulting in increased pro-IL-1ß synthesis and disease initiation and/or progression.

Semantic clustering analysis of E3-ubiquitin ligases in gastrointestinal tract defines genes ontology clusters with tissue expression patterns.

BACKGROUND: Ubiquitin ligases (Ub-ligases) are essential intracellular enzymes responsible for the regulation of proteome homeostasis, signaling pathway crosstalk, cell differentiation and stress responses. Individual Ub-ligases exhibit their unique functions based on the nature of their substrates. They create a complex regulatory network with alternative and feedback pathways to maintain cell homeostasis, being thus important players in many physiological and pathological conditions. However, the functional classification of Ub-ligases needs to be revised and extended. METHODS: In the current study, we used a novel semantic biclustering technique for expression profiling of Ub-ligases and ubiquitination-related genes in the murine gastrointestinal tract (GIT). We accommodated a general framework of the algorithm for finding tissue-specific gene expression clusters in GIT. In order to test identified clusters in a biological system, we used a model of epithelial regeneration. For this purpose, a dextran sulfate sodium (DSS) mouse model, following with in situ hybridization, was used to expose genes with possible compensatory features. To determine cell-type specific distribution of Ub-ligases and ubiquitination-related genes, principal component analysis (PCA) and Uniform Manifold Approximation and Projection technique (UMAP) were used to analyze the Tabula Muris scRNA-seq data of murine colon followed by comparison with our clustering results. RESULTS: Our established clustering protocol, that incorporates the semantic biclustering algorithm, demonstrated the potential to reveal interesting expression patterns. In this manner, we statistically defined gene clusters consisting of the same genes involved in distinct regulatory pathways vs distinct genes playing roles in functionally similar signaling pathways. This allowed us to uncover the potentially redundant features of GIT-specific Ub-ligases and ubiquitination-related genes. Testing the statistically obtained results on the mouse model showed that genes clustered to the same ontology group simultaneously alter their expression pattern after induced epithelial damage, illustrating their complementary role during tissue regeneration. CONCLUSIONS: An optimized semantic clustering protocol demonstrates the potential to reveal a readable and unique pattern in the expression profiling of GIT-specific Ub-ligases, exposing ontologically relevant gene clusters with potentially redundant features. This extends our knowledge of ontological relationships among Ub-ligases and ubiquitination-related genes, providing an alternative and more functional gene classification. In a similar way, semantic cluster analysis could be used for studding of other enzyme families, tissues and systems.