SARS‐CoV‐2 receptor ACE 2 and TMPRSS 2 are primarily expressed in bronchial transient secretory cells

The SARS-CoV-2 pandemic affecting the human respiratory system severely challenges public health and urgently demands for increasing our understanding of COVID-19 pathogenesis, especially host factors facilitating virus infection and replication. SARS-CoV-2 was reported to enter cells via binding to ACE2, followed by its priming by TMPRSS2. Here, we investigate ACE2 and TMPRSS2 expression levels and their distribution across cell types in lung tissue (twelve donors, 39,778 cells) and in cells derived from subsegmental bronchial branches (four donors, 17,521 cells) by single nuclei and single cell RNA sequencing, respectively. While TMPRSS2 is strongly expressed in both tissues, in the subsegmental bronchial branches ACE2 is predominantly expressed in a transient secretory cell type. Interestingly, these transiently differentiating cells show an enrichment for pathways related to RHO GTPase function and viral processes suggesting increased vulnerability for SARS-CoV-2 infection. Our data provide a rich resource for future investigations of COVID-19 infection and pathogenesis

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno's original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways

Comparison of single cell transcriptomics technologies and their application to investigate cellular heterogeneity in healthy and diseased lung

Multicellular organisms rely on the concerted interaction of a multitude of cells, which are often highly specialised and give rise to complex tissues. As vastly different cellular phenotypes emerge from the same genotype that is shared across all cells of an organism, the transcriptome represents a key mediator driving different cell types and cell states that give rise to functional tissues. These are also subject to environmental factors or intrinsic changes that may disrupt homeostasis and lead to disease. In the human lung, the effects of tobacco smoke exposure, still the greatest risk factor for lung cancer, have not been fully resolved at the cellular level. Moreover, cellular heterogeneity may be significant for the emergence of lung cancer in never smokers, a growing proportion of global cases. A focused investigation of cellular heterogeneity in the healthy lung and lung cancers is therefore highly warranted. During the last decade, technological advancements have made it possible to interrogate the transcriptome of single cells by novel next generation sequencing approaches. While previous studies were limited to averaging transcriptome information over many cells, single cell RNA sequencing (scRNA-seq) technologies are now enabling the investigation of cellular phenotypes in healthy and diseased tissues at unprecedented resolution. In this thesis, I adapt different scRNA-seq technologies to process fresh or biobanked samples from different tissues and species, thus enabling comparisons across diverse origins. We identify specific advantages, limitations and experimental challenges associated with each technology. I then perform a comprehensive single-cell transcriptomics study of healthy lung and lung adenocarcinoma (LADC). Based on twelve healthy lung samples, we generate a reference cell atlas that provides a rich resource for investigating cellular diversity in the human alveolar lung. Its utility is demonstrated by probing the expression of genes that are implicated in host cell entry of SARS-CoV-2 virus, thereby vi contributing to our understanding of coronavirus infections. By comparing single cell profiles from smokers and never smokers, we also resolve the involvement of distinct cell types in the maintenance of an inflammatory state in smoker lungs, and we identify key mediators of inflammatory processes induced by tobacco smoke exposure in fibroblasts and endothelial cells. To investigate cell type diversity and microenvironment interactions in LADC, I analyse 26 tumour tissue samples and resolve functional malignant cell subpopulations linked by a differentiation hierarchy in both smokers and never smokers. They comprise proliferating and intermediate undifferentiated cells as well as two differentiated tumour cell states implicated in cancer progression and invasiveness. Distinct macrophage and fibroblast subpopulations which contribute to a tumourigenic environment are also detected. A subset of proliferating tumour cells show differential immune modulating activity dependent on smoking status, with implications for future treatment approaches. Taken together, these results provide a comparison of rapidly developing scRNA-seq technologies for use in further studies and demonstrate their utility to dissect cellular heterogeneity and identify transcriptional programmes in the healthy and diseased lung. By applying these technologies, I add to our understanding of SARS-CoV-2 entry into human lung cells, define the alveolar lung cell types affected by tobacco smoke exposure, and provide deeper insight into cellular heterogeneity of LADC and the tumour microenvironment. These findings represent a valuable reference for future translational studies

Major changes of cell function and toxicant sensitivity in cultured cells undergoing mild, quasi-natural genetic drift

Genomic drift affects the functional properties of cell lines, and the reproducibility of data from in vitro studies. While chromosomal aberrations and mutations in single pivotal genes are well explored, little is known about effects of minor, possibly pleiotropic, genome changes. We addressed this question for the human dopaminergic neuronal precursor cell line LUHMES by comparing two subpopulations (SP) maintained either at the American-Type-Culture-Collection (ATCC) or by the original provider (UKN). Drastic differences in susceptibility towards the specific dopaminergic toxicant 1-methyl-4-phenylpyridinium (MPP+) were observed. Whole-genome sequencing was performed to identify underlying genetic differences. While both SP had normal chromosome structures, they displayed about 70 differences on the level of amino acid changing events. Some of these differences were confirmed biochemically, but none offered a direct explanation for the altered toxicant sensitivity pattern. As second approach, markers known to be relevant for the intended use of the cells were specifically tested. The “ATCC” cells rapidly down-regulated the dopamine-transporter and tyrosine-hydroxylase after differentiation, while “UKN” cells maintained functional levels. As the respective genes were not altered themselves, we conclude that polygenic complex upstream changes can have drastic effects on biochemical features and toxicological responses of relatively similar SP of cells

Major changes of cell function and toxicant sensitivity in cultured cells undergoing mild, quasi-natural genetic drift

Genomic drift affects the functional properties of cell lines, and the reproducibility of data from in vitro studies. While chromosomal aberrations and mutations in single pivotal genes are well explored, little is known about effects of minor, possibly pleiotropic, genome changes. We addressed this question for the human dopaminergic neuronal precursor cell line LUHMES by comparing two subpopulations (SP) maintained either at the American-Type-Culture-Collection (ATCC) or by the original provider (UKN). Drastic differences in susceptibility towards the specific dopaminergic toxicant 1-methyl-4-phenylpyridinium (MPP+) were observed. Whole-genome sequencing was performed to identify underlying genetic differences. While both SP had normal chromosome structures, they displayed about 70 differences on the level of amino acid changing events. Some of these differences were confirmed biochemically, but none offered a direct explanation for the altered toxicant sensitivity pattern. As second approach, markers known to be relevant for the intended use of the cells were specifically tested. The "ATCC" cells rapidly down-regulated the dopamine-transporter and tyrosine-hydroxylase after differentiation, while "UKN" cells maintained functional levels. As the respective genes were not altered themselves, we conclude that polygenic complex upstream changes can have drastic effects on biochemical features and toxicological responses of relatively similar SP of cells.publishe

Development of a neural rosette formation assay (RoFA) to identify neurodevelopmental toxicants and to characterize their transcriptome disturbances

The first in vitro tests for developmental toxicity made use of rodent cells. Newer teratology tests, e.g. developed during the ESNATS project, use human cells and measure mechanistic endpoints (such as transcriptome changes). However, the toxicological implications of mechanistic parameters are hard to judge, without functional/morphological endpoints. To address this issue, we developed a new version of the human stem cell-based test STOP-tox(UKN). For this purpose, the capacity of the cells to self-organize to neural rosettes was assessed as functional endpoint: pluripotent stem cells were allowed to differentiate into neuroepithelial cells for 6 days in the presence or absence of toxicants. Then, both transcriptome changes were measured (standard STOP-tox(UKN)) and cells were allowed to form rosettes. After optimization of staining methods, an imaging algorithm for rosette quantification was implemented and used for an automated rosette formation assay (RoFA). Neural tube toxicants (like valproic acid), which are known to disturb human development at stages when rosette-forming cells are present, were used as positive controls. Established toxicants led to distinctly different tissue organization and differentiation stages. RoFA outcome and transcript changes largely correlated concerning (1) the concentration-dependence, (2) the time dependence, and (3) the set of positive hits identified amongst 24 potential toxicants. Using such comparative data, a prediction model for the RoFA was developed. The comparative analysis was also used to identify gene dysregulations that are particularly predictive for disturbed rosette formation. This ‘RoFA predictor gene set’ may be used for a simplified and less costly setup of the STOP-tox(UKN) assay.publishe

Development of a neural rosette formation assay (RoFA) to identify neurodevelopmental toxicants and to characterize their transcriptome disturbances

The first in vitro tests for developmental toxicity made use of rodent cells. Newer teratology tests, e.g. developed during the ESNATS project, use human cells and measure mechanistic endpoints (such as transcriptome changes). However, the toxicological implications of mechanistic parameters are hard to judge, without functional/morphological endpoints. To address this issue, we developed a new version of the human stem cell-based test STOP-tox((UKN)). For this purpose, the capacity of the cells to self-organize to neural rosettes was assessed as functional endpoint: pluripotent stem cells were allowed to differentiate into neuroepithelial cells for 6 days in the presence or absence of toxicants. Then, both transcriptome changes were measured (standard STOP-tox((UKN))) and cells were allowed to form rosettes. After optimization of staining methods, an imaging algorithm for rosette quantification was implemented and used for an automated rosette formation assay (RoFA). Neural tube toxicants (like valproic acid), which are known to disturb human development at stages when rosette-forming cells are present, were used as positive controls. Established toxicants led to distinctly different tissue organization and differentiation stages. RoFA outcome and transcript changes largely correlated concerning (1) the concentration-dependence, (2) the time dependence, and (3) the set of positive hits identified amongst 24 potential toxicants. Using such comparative data, a prediction model for the RoFA was developed. The comparative analysis was also used to identify gene dysregulations that are particularly predictive for disturbed rosette formation. This 'RoFA predictor gene set' may be used for a simplified and less costly setup of the STOP-tox((UKN)) assay

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals1,2,3,4,5,6, probably owing to the evolutionary pressure on males to be reproductively successful7. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals1-6, probably owing to the evolutionary pressure on males to be reproductively successful7. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals