Skin single-cell transcriptomics reveals a core of sebaceous gland-relevant genes shared by mice and humans.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) has been widely applied to dissect cellular heterogeneity in normal and diseased skin. Sebaceous glands, essential skin components with established functions in maintaining skin integrity and emerging roles in systemic energy metabolism, have been largely neglected in scRNA-seq studies. METHODS: Departing from mouse and human skin scRNA-seq datasets, we identified gene sets expressed especially in sebaceous glands with the open-source R-package oposSOM. RESULTS: The identified gene sets included sebaceous gland-typical genes as Scd3, Mgst1, Cidea, Awat2 and KRT7. Surprisingly, however, there was not a single overlap among the 100 highest, exclusively in sebaceous glands expressed transcripts in mouse and human samples. Notably, both species share a common core of only 25 transcripts, including mitochondrial and peroxisomal genes involved in fatty acid, amino acid, and glucose processing, thus highlighting the intense metabolic rate of this gland. CONCLUSIONS: This study highlights intrinsic differences in sebaceous lipid synthesis between mice and humans, and indicates an important role for peroxisomal processes in this context. Our data also provides attractive starting points for experimentally addressing novel candidates regulating sebaceous gland homeostasis.

Correction: Winkler et al. The Adhesion G-Protein-Coupled Receptor GPR115/ADGRF4 Regulates Epidermal Differentiation and Associates with Cytoskeletal KRT1. Cells 2022, 11, 3151.

In the original publication [...].

Indirect Virus Transmission via Fomites Can Counteract Lock-Down Effectiveness.

The spread of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has raised major health policy questions. Direct transmission via respiratory droplets seems to be the dominant route of its transmission. However, indirect transmission via shared contact of contaminated objects may also occur. The contribution of each transmission route to epidemic spread might change during lock-down scenarios. Here, we simulate viral spread of an abstract epidemic considering both routes of transmission by use of a stochastic, agent-based SEIR model. We show that efficient contact tracing (CT) at a high level of incidence can stabilize daily cases independently of the transmission route long before effects of herd immunity become relevant. CT efficacy depends on the fraction of cases that do not show symptoms. Combining CT with lock-down scenarios that reduce agent mobility lowers the incidence for exclusive direct transmission scenarios and can even eradicate the epidemic. However, even for small fractions of indirect transmission, such lockdowns can impede CT efficacy and increase case numbers. These counterproductive effects can be reduced by applying measures that favor distancing over reduced mobility. In summary, we show that the efficacy of lock-downs depends on the transmission route. Our results point to the particular importance of hygiene measures during mobility lock-downs.

The Adhesion G-Protein-Coupled Receptor GPR115/ADGRF4 Regulates Epidermal Differentiation and Associates with Cytoskeletal KRT1.

Among the 33 human adhesion G-protein-coupled receptors (aGPCRs), a unique subfamily of GPCRs, only ADGRF4, encoding GPR115, shows an obvious skin-dominated transcriptomic profile, but its expression and function in skin is largely unknown. Here, we report that GPR115 is present in a small subset of basal and in most suprabasal, noncornified keratinocytes of the stratified epidermis, supporting epidermal transcriptomic data. In psoriatic skin, characterized by hyperproliferation and delayed differentiation, the expression of GPR115 and KRT1/10, the fundamental suprabasal keratin dimer, is delayed. The deletion of ADGRF4 in HaCaT keratinocytes grown in an organotypic mode abrogates KRT1 and reduces keratinocyte stratification, indicating a role of GPR115 in epidermal differentiation. Unexpectedly, endogenous GPR115, which is not glycosylated and is likely not proteolytically processed, localizes intracellularly along KRT1/10-positive keratin filaments in a regular pattern. Our data demonstrate a hitherto unknown function of GPR115 in the regulation of epidermal differentiation and KRT1.

Reduced adhesion of aged intestinal stem cells contributes to an accelerated clonal drift.

Upon aging, the function of the intestinal epithelium declines with a concomitant increase in aging-related diseases. ISCs play an important role in this process. It is known that ISC clonal dynamics follow a neutral drift model. However, it is not clear whether the drift model is still valid in aged ISCs. Tracking of clonal dynamics by clonal tracing revealed that aged crypts drift into monoclonality substantially faster than young ones. However, ISC tracing experiments, in vivo and ex vivo, implied a similar clonal expansion ability of both young and aged ISCs. Single-cell RNA sequencing for 1,920 high Lgr5 ISCs from young and aged mice revealed increased heterogeneity among subgroups of aged ISCs. Genes associated with cell adhesion were down-regulated in aged ISCs. ISCs of aged mice indeed show weaker adhesion to the matrix. Simulations applying a single cell-based model of the small intestinal crypt demonstrated an accelerated clonal drift at reduced adhesion strength, implying a central role for reduced adhesion for affecting clonal dynamics upon aging.

Organoid Cultures In Silico: Tools or Toys?

The implementation of stem-cell-based organoid culture more than ten years ago started a development that created new avenues for diagnostic analyses and regenerative medicine. In parallel, computational modelling groups realized the potential of this culture system to support their theoretical approaches to study tissues in silico. These groups developed computational organoid models (COMs) that enabled testing consistency between cell biological data and developing theories of tissue self-organization. The models supported a mechanistic understanding of organoid growth and maturation and helped linking cell mechanics and tissue shape in general. What comes next? Can we use COMs as tools to complement the equipment of our biological and medical research? While these models already support experimental design, can they also quantitatively predict tissue behavior? Here, we review the current state of the art of COMs and discuss perspectives for their application.

Epigenetic Drifts during Long-Term Intestinal Organoid Culture.

Organoids retain the morphological and molecular patterns of their tissue of origin, are self-organizing, relatively simple to handle and accessible to genetic engineering. Thus, they represent an optimal tool for studying the mechanisms of tissue maintenance and aging. Long-term expansion under standard growth conditions, however, is accompanied by changes in the growth pattern and kinetics. As a potential explanation of these alterations, epigenetic drifts in organoid culture have been suggested. Here, we studied histone tri-methylation at lysine 4 (H3K4me3) and 27 (H3K27me3) and transcriptome profiles of intestinal organoids derived from mismatch repair (MMR)-deficient and control mice and cultured for 3 and 20 weeks and compared them with data on their tissue of origin. We found that, besides the expected changes in short-term culture, the organoids showed profound changes in their epigenomes also during the long-term culture. The most prominent were epigenetic gene activation by H3K4me3 recruitment to previously unmodified genes and by H3K27me3 loss from originally bivalent genes. We showed that a long-term culture is linked to broad transcriptional changes that indicate an ongoing maturation and metabolic adaptation process. This process was disturbed in MMR-deficient mice, resulting in endoplasmic reticulum (ER) stress and Wnt activation. Our results can be explained in terms of a mathematical model assuming that epigenetic changes during a long-term culture involve DNA demethylation that ceases if the metabolic adaptation is disturbed.

How to Obtain a Mega-Intestine with Normal Morphology: In Silico Modelling of Postnatal Intestinal Growth in a Cd97-Transgenic Mouse.

Intestinal cylindrical growth peaks in mice a few weeks after birth, simultaneously with crypt fission activity. It nearly stops after weaning and cannot be reactivated later. Transgenic mice expressing Cd97/Adgre5 in the intestinal epithelium develop a mega-intestine with normal microscopic morphology in adult mice. Here, we demonstrate premature intestinal differentiation in Cd97/Adgre5 transgenic mice at both the cellular and molecular levels until postnatal day 14. Subsequently, the growth of the intestinal epithelium becomes activated and its maturation suppressed. These changes are paralleled by postnatal regulation of growth factors and by an increased expression of secretory cell markers, suggesting growth activation of non-epithelial tissue layers as the origin of enforced tissue growth. To understand postnatal intestinal growth mechanistically, we study epithelial fate decisions during this period with the use of a 3D individual cell-based computer model. In the model, the expansion of the intestinal stem cell (SC) population, a prerequisite for crypt fission, is largely independent of the tissue growth rate and is therefore not spontaneously adaptive. Accordingly, the model suggests that, besides the growth activation of non-epithelial tissue layers, the formation of a mega-intestine requires a released growth control in the epithelium, enabling accelerated SC expansion. The similar intestinal morphology in Cd97/Adgre5 transgenic and wild type mice indicates a synchronization of tissue growth and SC expansion, likely by a crypt density-controlled contact inhibition of growth of intestinal SC proliferation. The formation of a mega-intestine with normal microscopic morphology turns out to originate in changes of autonomous and conditional specification of the intestinal cell fate induced by the activation of Cd97/Adgre5.

Fighting Against Promoter DNA Hyper-Methylation: Protective Histone Modification Profiles of Stress-Resistant Intestinal Stem Cells.

Aberrant DNA methylation in stem cells is a hallmark of aging and tumor development. Recently, we have suggested that promoter DNA hyper-methylation originates in DNA repair and that even successful DNA repair might confer this kind of epigenetic long-term change. Here, we ask for interrelations between promoter DNA methylation and histone modification changes observed in the intestine weeks after irradiation and/or following Msh2 loss. We focus on H3K4me3 recruitment to the promoter of H3K27me3 target genes. By RNA- and histone ChIP-sequencing, we demonstrate that this recruitment occurs without changes of the average gene transcription and does not involve H3K9me3. Applying a mathematical model of epigenetic regulation of transcription, we show that the recruitment can be explained by stronger DNA binding of H3K4me3 and H3K27me3 histone methyl-transferases as a consequence of lower DNA methylation. This scenario implicates stable transcription despite of H3K4me3 recruitment, in agreement with our RNA-seq data. Following several kinds of stress, including moderate irradiation, stress-sensitive intestinal stem cell (ISCs) are known to become replaced by more resistant populations. Our simulation results suggest that the stress-resistant ISCs are largely protected against promoter hyper-methylation of H3K27me3 target genes.

Loss of Msh2 and a single-radiation hit induce common, genome-wide, and persistent epigenetic changes in the intestine.

BACKGROUND: Mismatch repair (MMR)-deficiency increases the risk of colorectal tumorigenesis. To determine whether the tumors develop on a normal or disturbed epigenetic background and how radiation affects this, we quantified genome-wide histone H3 methylation profiles in macroscopic normal intestinal tissue of young radiated and untreated MMR-deficient VCMsh2LoxP/LoxP (Msh2-/-) mice months before tumor onset. RESULTS: Histone H3 methylation increases in Msh2-/- compared to control Msh2+/+ mice. Activating H3K4me3 and H3K36me3 histone marks frequently accumulate at genes that are H3K27me3 or H3K4me3 modified in Msh2+/+ mice, respectively. The genes recruiting H3K36me3 enrich in gene sets associated with DNA repair, RNA processing, and ribosome biogenesis that become transcriptionally upregulated in the developing tumors. A similar epigenetic effect is present in Msh2+/+ mice 4 weeks after a single-radiation hit, whereas radiation of Msh2-/- mice left their histone methylation profiles almost unchanged. CONCLUSIONS: MMR deficiency results in genome-wide changes in histone H3 methylation profiles preceding tumor development. Similar changes constitute a persistent epigenetic signature of radiation-induced DNA damage.

Linking DNA Damage and Age-Related Promoter DNA Hyper-Methylation in the Intestine.

Aberrant DNA methylation in stem cells is a hallmark of aging and tumor development. Here, we explore whether and how DNA damage repair might impact on these time-dependent changes, in particular in proliferative intestinal stem cells. We introduce a 3D multiscale computer model of intestinal crypts enabling simulation of aberrant DNA and histone methylation of gene promoters during aging. We assume histone state-dependent activity of de novo DNA methyltransferases (DNMTs) and methylation-dependent binding of maintenance DNMTs to CpGs. We simulate aging with and without repeated DNA repair. Motivated by recent findings on the histone demethylase KDM2b, we consider that DNA repair is associated with chromatin opening and improved recruitment of de novo DNMTs. Our results suggest that methylation-dependent binding of maintenance DNMTs to CpGs, establishing bistable DNA methylation states, is a prerequisite to promoter hyper-methylation following DNA repair. With this, the transient increase in de novo DNMT activity during repair can induce switches from low to high methylation states. These states remain stable after repair, leading to an epigenetic drift. The switches are most frequent in genes with H3K27me3 modified promoters. Our model provides a mechanistic explanation on how even successful DNA repair might confer long term changes of the epigenome.

Linking stem cell function and growth pattern of intestinal organoids.

Intestinal stem cells (ISCs) require well-defined signals from their environment in order to carry out their specific functions. Most of these signals are provided by neighboring cells that form a stem cell niche, whose shape and cellular composition self-organize. Major features of this self-organization can be studied in ISC-derived organoid culture. In this system, manipulation of essential pathways of stem cell maintenance and differentiation results in well-described growth phenotypes. We here provide an individual cell-based model of intestinal organoids that enables a mechanistic explanation of the observed growth phenotypes. In simulation studies of the 3D structure of expanding organoids, we investigate interdependences between Wnt- and Notch-signaling which control the shape of the stem cell niche and, thus, the growth pattern of the organoids. Similar to in vitro experiments, changes of pathway activities alter the cellular composition of the organoids and, thereby, affect their shape. Exogenous Wnt enforces transitions from branched into a cyst-like growth pattern; known to occur spontaneously during long term organoid expansion. Based on our simulation results, we predict that the cyst-like pattern is associated with biomechanical changes of the cells which assign them a growth advantage. The results suggest ongoing stem cell adaptation to in vitro conditions during long term expansion by stabilizing Wnt-activity. Our study exemplifies the potential of individual cell-based modeling in unraveling links between molecular stem cell regulation and 3D growth of tissues. This kind of modeling combines experimental results in the fields of stem cell biology and cell biomechanics constituting a prerequisite for a better understanding of tissue regeneration as well as developmental processes.

Adhesion forces and cortical tension couple cell proliferation and differentiation to drive epidermal stratification.

To establish and maintain organ structure and function, tissues need to balance stem cell proliferation and differentiation rates and coordinate cell fate with position. By quantifying and modelling tissue stress and deformation in the mammalian epidermis, we find that this balance is coordinated through local mechanical forces generated by cell division and delamination. Proliferation within the basal stem/progenitor layer, which displays features of a jammed, solid-like state, leads to crowding, thereby locally distorting cell shape and stress distribution. The resulting decrease in cortical tension and increased cell-cell adhesion trigger differentiation and subsequent delamination, reinstating basal cell layer density. After delamination, cells establish a high-tension state as they increase myosin II activity and convert to E-cadherin-dominated adhesion, thereby reinforcing the boundary between basal and suprabasal layers. Our results uncover how biomechanical signalling integrates single-cell behaviours to couple proliferation, cell fate and positioning to generate a multilayered tissue.

The Regulatory Capacity of Bivalent Genes-A Theoretical Approach.

Bivalent genes are frequently associated with developmental and lineage specification processes. Resolving their bivalency enables fast changes in their expression, which potentially can trigger cell fate decisions. Here, we provide a theoretical model of bivalency that allows for predictions on the occurrence, stability and regulatory capacity of this prominent modification state. We suggest that bivalency enables balanced gene expression heterogeneity that constitutes a prerequisite of robust lineage priming in somatic stem cells. Moreover, we demonstrate that interactions between the histone and DNA methylation machineries together with the proliferation activity control the stability of the bivalent state and can turn it into an unmodified state. We suggest that deregulation of these interactions underlies cell transformation processes as associated with acute myeloid leukemia (AML) and provide a model of AML blast formation following deregulation of the Ten-eleven Translocation (TET) pathway.

Stem cell competition in the gut: insights from multi-scale computational modelling.

Three-dimensional (3D) computational tissue models can provide a comprehensive description of tissue dynamics at the molecular, cellular and tissue level. Moreover, they can support the development of hypotheses about cellular interactions and about synergies between major signalling pathways. We exemplify these capabilities by simulation of a 3D single-cell-based model of mouse small intestinal crypts. We analyse the impact of lineage specification, distribution and cellular lifespan on clonal competition and study effects of Notch- and Wnt activation on fixation of mutations within the tissue. Based on these results, we predict that experimentally observed synergistic effects between autonomous Notch- and Wnt signalling in triggering intestinal tumourigenesis originate in the suppression of Wnt-dependent secretory lineage specification by Notch, giving rise to an increased fixation probability of Wnt-activating mutations. Our study demonstrates that 3D computational tissue models can support a mechanistic understanding of long-term tissue dynamics under homeostasis and during transformation.

A Branch-and-Bound Approach for Tautomer Enumeration.

Knowledge about tautomer forms of a structure is important since, e.g., a property prediction for a molecule can yield to different results which depend on the individual tautomer. Tautomers are isomers that can be transformed to each other through chemical equilibrium reactions. In this paper the first exact Branch-and-Bound (B&B) algorithm to calculate tautomer structures is proposed. The algorithm is complete in the sense of tautomerism and generates all possible tautomers of a structure according to the tautomer definition, it is initialized with. To be efficient, the algorithm takes advantage of symmetric and formation properties. Some restrictions are used to enable an early pruning of some branches of the B&B tree. This is important, since a simple enumeration strategy would lead to number of candidate tautomers that is exponentially increasing with the number of hydrogen atoms and their attachment sites. The proposed implementation of the B&B algorithm covers the majority of the prototropic tautomer cases, but can be adapted to other kinds of tautomerism too. Furthermore, a computer processable definition of tautomerism is given in the form of the moving hydrogen atom problem.

AgeFactDB--the JenAge Ageing Factor Database--towards data integration in ageing research.

AgeFactDB (http://agefactdb.jenage.de) is a database aimed at the collection and integration of ageing phenotype data including lifespan information. Ageing factors are considered to be genes, chemical compounds or other factors such as dietary restriction, whose action results in a changed lifespan or another ageing phenotype. Any information related to the effects of ageing factors is called an observation and is presented on observation pages. To provide concise access to the complete information for a particular ageing factor, corresponding observations are also summarized on ageing factor pages. In a first step, ageing-related data were primarily taken from existing databases such as the Ageing Gene Database--GenAge, the Lifespan Observations Database and the Dietary Restriction Gene Database--GenDR. In addition, we have started to include new ageing-related information. Based on homology data taken from the HomoloGene Database, AgeFactDB also provides observation and ageing factor pages of genes that are homologous to known ageing-related genes. These homologues are considered as candidate or putative ageing-related genes. AgeFactDB offers a variety of search and browse options, and also allows the download of ageing factor or observation lists in TSV, CSV and XML formats.

Tautomer identification and tautomer structure generation based on the InChI code.

An algorithm is introduced that enables a fast generation of all possible prototropic tautomers resulting from the mobile H atoms and associated heteroatoms as defined in the InChI code. The InChI-derived set of possible tautomers comprises (1,3)-shifts for open-chain molecules and (1,n)-shifts (with n being an odd number >3) for ring systems. In addition, our algorithm includes also, as extension to the InChI scope, those larger (1,n)-shifts that can be constructed from joining separate but conjugated InChI sequences of tautomer-active heteroatoms. The developed algorithm is described in detail, with all major steps illustrated through explicit examples. Application to approximately 72,500 organic compounds taken from EINECS (European Inventory of Existing Commercial Chemical Substances) shows that around 11% of the substances occur in different heteroatom-prototropic tautomeric forms. Additional QSAR (quantitative structure-activity relationship) predictions of their soil sorption coefficient and water solubility reveal variations across tautomers up to more than two and 4 orders of magnitude, respectively. For a small subset of nine compounds, analysis of quantum chemically predicted tautomer energies supports the view that among all tautomers of a given compound, those restricted to H atom exchanges between heteroatoms usually include the thermodynamically most stable structures.