Cell-type identity of the avian utricle.

The avian utricle, a vestibular organ of the inner ear, displays turnover of sensory hair cells throughout life. This is in sharp contrast to the mammalian utricle, which shows limited regenerative capacity. Here, we use single-cell RNA sequencing to identify distinct marker genes for the different sensory hair cell subtypes of the chicken utricle, which we validated in situ. We provide markers for spatially distinct supporting cell populations and identify two transitional cell populations of dedifferentiating supporting cells and developing hair cells. Trajectory reconstruction resulted in an inventory of gene expression dynamics of natural hair cell generation in the avian utricle.

Spatiotemporal dynamics of inner ear sensory and non-sensory cells revealed by single-cell transcriptomics.

The utricle is a vestibular sensory organ that requires mechanosensitive hair cells to detect linear acceleration. In neonatal mice, new hair cells are derived from non-sensory supporting cells, yet cell type diversity and mechanisms of cell addition remain poorly characterized. Here, we perform computational analyses on single-cell transcriptomes to categorize cell types and resolve 14 individual sensory and non-sensory subtypes. Along the periphery of the sensory epithelium, we uncover distinct groups of transitional epithelial cells, marked by Islr, Cnmd, and Enpep expression. By reconstructing de novo trajectories and gene dynamics, we show that as the utricle expands, Islr+ transitional epithelial cells exhibit a dynamic and proliferative phase to generate new supporting cells, followed by coordinated differentiation into hair cells. Taken together, our study reveals a sequential and coordinated process by which non-sensory epithelial cells contribute to growth of the postnatal mouse sensory epithelium.

Cell-type identity of the avian cochlea.

In contrast to mammals, birds recover naturally from acquired hearing loss, which makes them an ideal model for inner ear regeneration research. Here, we present a validated single-cell RNA sequencing resource of the avian cochlea. We describe specific markers for three distinct types of sensory hair cells, including a previously unknown subgroup, which we call superior tall hair cells. We identify markers for the supporting cells associated with tall hair cells, which represent the facultative stem cells of the avian inner ear. Likewise, we present markers for supporting cells that are located below the short cochlear hair cells. We further infer spatial expression gradients of hair cell genes along the tonotopic axis of the cochlea. This resource advances neurobiology, comparative biology, and regenerative medicine by providing a basis for comparative studies with non-regenerating mammalian cochleae and for longitudinal studies of the regenerating avian cochlea.

Blinatumomab-induced T cell activation at single cell transcriptome resolution.

BACKGROUND: Bi-specific T-cell engager (BiTE) antibody is a class of bispecific antibodies designed for cancer immunotherapy. Blinatumomab is the first approved BiTE to treat acute B cell lymphoblastic leukemia (B-ALL). It brings killer T and target B cells into close proximity, activating patient's autologous T cells to kill malignant B cells via mechanisms such as cytolytic immune synapse formation and inflammatory cytokine production. However, the activated T-cell subtypes and the target cell-dependent T cell responses induced by blinatumomab, as well as the mechanisms of resistance to blinatumomab therapy are largely unknown. RESULTS: In this study, we performed single-cell sequencing analysis to identify transcriptional changes in T cells following blinatumomab-induced T cell activation using single cells from both, a human cell line model and a patient-derived model of blinatumomab-mediated cytotoxicity. In total, the transcriptome of 17,920 single T cells from the cell line model and 2271 single T cells from patient samples were analyzed. We found that CD8+ effector memory T cells, CD4+ central memory T cells, naïve T cells, and regulatory T cells were activated after blinatumomab treatment. Here, blinatumomab-induced transcriptional changes reflected the functional immune activity of the blinatumomab-activated T cells, including the upregulation of pathways such as the immune system, glycolysis, IFNA signaling, gap junctions, and IFNG signaling. Co-stimulatory (TNFRSF4 and TNFRSF18) and co-inhibitory (LAG3) receptors were similarly upregulated in blinatumomab-activated T cells, indicating ligand-dependent T cell functions. Particularly, B-ALL cell expression of TNFSF4, which encodes the ligand of T cell co-stimulatory receptor TNFRSF4, was found positively correlated with the response to blinatumomab treatment. Furthermore, recombinant human TNFSF4 protein enhanced the cytotoxic activity of blinatumomab against B-ALL cells. CONCLUSION: These results reveal a target cell-dependent mechanism of T-cell activation by blinatumomab and suggest that TNFSF4 may be responsible for the resistant mechanism and a potential target for combination therapy with blinatumomab, to treat B-ALL or other B-cell malignancies.

Prior activation state shapes the microglia response to antihuman TREM2 in a mouse model of Alzheimer's disease.

Triggering receptor expressed on myeloid cells 2 (TREM2) sustains microglia response to brain injury stimuli including apoptotic cells, myelin damage, and amyloid ß (Aß). Alzheimer's disease (AD) risk is associated with the TREM2R47H variant, which impairs ligand binding and consequently microglia responses to Aß pathology. Here, we show that TREM2 engagement by the mAb hT2AB as surrogate ligand activates microglia in 5XFAD transgenic mice that accumulate Aß and express either the common TREM2 variant (TREM2CV) or TREM2R47H scRNA-seq of microglia from TREM2CV-5XFAD mice treated once with control hIgG1 exposed four distinct trajectories of microglia activation leading to disease-associated (DAM), interferon-responsive (IFN-R), cycling (Cyc-M), and MHC-II expressing (MHC-II) microglia types. All of these were underrepresented in TREM2R47H-5XFAD mice, suggesting that TREM2 ligand engagement is required for microglia activation trajectories. Moreover, Cyc-M and IFN-R microglia were more abundant in female than male TREM2CV-5XFAD mice, likely due to greater Aß load in female 5XFAD mice. A single systemic injection of hT2AB replenished Cyc-M, IFN-R, and MHC-II pools in TREM2R47H-5XFAD mice. In TREM2CV-5XFAD mice, however, hT2AB brought the representation of male Cyc-M and IFN-R microglia closer to that of females, in which these trajectories had already reached maximum capacity. Moreover, hT2AB induced shifts in gene expression patterns in all microglial pools without affecting representation. Repeated treatment with a murinized hT2AB version over 10 d increased chemokines brain content in TREM2R47H-5XFAD mice, consistent with microglia expansion. Thus, the impact of hT2AB on microglia is shaped by the extent of TREM2 endogenous ligand engagement and basal microglia activation.

Systematic comparison of high-throughput single-cell RNA-seq methods for immune cell profiling.

BACKGROUND: Elucidation of immune populations with single-cell RNA-seq has greatly benefited the field of immunology by deepening the characterization of immune heterogeneity and leading to the discovery of new subtypes. However, single-cell methods inherently suffer from limitations in the recovery of complete transcriptomes due to the prevalence of cellular and transcriptional dropout events. This issue is often compounded by limited sample availability and limited prior knowledge of heterogeneity, which can confound data interpretation. RESULTS: Here, we systematically benchmarked seven high-throughput single-cell RNA-seq methods. We prepared 21 libraries under identical conditions of a defined mixture of two human and two murine lymphocyte cell lines, simulating heterogeneity across immune-cell types and cell sizes. We evaluated methods by their cell recovery rate, library efficiency, sensitivity, and ability to recover expression signatures for each cell type. We observed higher mRNA detection sensitivity with the 10x Genomics 5' v1 and 3' v3 methods. We demonstrate that these methods have fewer dropout events, which facilitates the identification of differentially-expressed genes and improves the concordance of single-cell profiles to immune bulk RNA-seq signatures. CONCLUSION: Overall, our characterization of immune cell mixtures provides useful metrics, which can guide selection of a high-throughput single-cell RNA-seq method for profiling more complex immune-cell heterogeneity usually found in vivo.

Fbxo2VHC mouse and embryonic stem cell reporter lines delineate in vitro-generated inner ear sensory epithelia cells and enable otic lineage selection and Cre-recombination.

While the mouse has been a productive model for inner ear studies, a lack of highly specific genes and tools has presented challenges. The absence of definitive otic lineage markers and tools is limiting in vitro studies of otic development, where innate cellular heterogeneity and disorganization increase the reliance on lineage-specific markers. To address this challenge in mice and embryonic stem (ES) cells, we targeted the lineage-specific otic gene Fbxo2 with a multicistronic reporter cassette (Venus/Hygro/CreER = VHC). In otic organoids derived from ES cells, Fbxo2VHC specifically delineates otic progenitors and inner ear sensory epithelia. In mice, Venus expression and CreER activity reveal a cochlear developmental gradient, label the prosensory lineage, show enrichment in a subset of type I vestibular hair cells, and expose strong expression in adult cerebellar granule cells. We provide a toolbox of multiple spectrally distinct reporter combinations for studies that require use of fluorescent reporters, hygromycin selection, and conditional Cre-mediated recombination.

Correction to: Aminoglycoside Damage and Hair Cell Regeneration in the Chicken Utricle.

This article was updated to correct a formatting error in Table 1.

Transcriptional Dynamics of Hair-Bundle Morphogenesis Revealed with CellTrails.

Protruding from the apical surface of inner ear sensory cells, hair bundles carry out mechanotransduction. Bundle growth involves sequential and overlapping cellular processes, which are concealed within gene expression profiles of individual cells. To dissect such processes, we developed CellTrails, a tool for uncovering, analyzing, and visualizing single-cell gene-expression dynamics. Utilizing quantitative gene-expression data for key bundle proteins from single cells of the developing chick utricle, we reconstructed de novo a bifurcating trajectory that spanned from progenitor cells to mature striolar and extrastriolar hair cells. Extraction and alignment of developmental trails and association of pseudotime with bundle length measurements linked expression dynamics of individual genes with bundle growth stages. Differential trail analysis revealed high-resolution dynamics of transcripts that control striolar and extrastriolar bundle development, including those that encode proteins that regulate [Ca2+]i or mediate crosslinking and lengthening of actin filaments.

Correction to: Aminoglycoside Damage and Hair Cell Regeneration in the Chicken Utricle.

This article was updated to correct a formatting error in Table 1.

Aminoglycoside Damage and Hair Cell Regeneration in the Chicken Utricle.

In this study, we present a systematic characterization of hair cell loss and regeneration in the chicken utricle in vivo. A single unilateral surgical delivery of streptomycin caused robust decline of hair cell numbers in striolar as well as extrastriolar regions, which in the striola was detected very early, 6 h post-insult. During the initial 12 h of damage response, we observed global repression of DNA replication, in contrast to the natural, mitotic hair cell production in undamaged control utricles. Regeneration of hair cells in striolar and extrastriolar regions occurred via high rates of asymmetric supporting cell divisions, accompanied by delayed replenishment by symmetric division. While asymmetric division of supporting cells is the main regenerative response to aminoglycoside damage, the detection of symmetric divisions supports the concept of direct transdifferentiation where supporting cells need to be replenished after their phenotypic conversion into new hair cells. Supporting cell divisions appear to be well coordinated because total supporting cell numbers throughout the regenerative process were invariant, despite the initial large-scale loss of hair cells. We conclude that a single ototoxic drug application provides an experimental framework to study the precise onset and timing of utricle hair cell regeneration in vivo. Our findings indicate that initial triggers and signaling events occur already within a few hours after aminoglycoside exposure. Direct transdifferentiation and asymmetric division of supporting cells to generate new hair cells subsequently happen largely in parallel and persist for several days.

Single-cell analysis delineates a trajectory toward the human early otic lineage.

Efficient pluripotent stem cell guidance protocols for the production of human posterior cranial placodes such as the otic placode that gives rise to the inner ear do not exist. Here we use a systematic approach including defined monolayer culture, signaling modulation, and single-cell gene expression analysis to delineate a developmental trajectory for human otic lineage specification in vitro. We found that modulation of bone morphogenetic protein (BMP) and WNT signaling combined with FGF and retinoic acid treatments over the course of 18 days generates cell populations that develop chronological expression of marker genes of non-neural ectoderm, preplacodal ectoderm, and early otic lineage. Gene expression along this differentiation path is distinct from other lineages such as endoderm, mesendoderm, and neural ectoderm. Single-cell analysis exposed the heterogeneity of differentiating cells and allowed discrimination of non-neural ectoderm and otic lineage cells from off-target populations. Pseudotemporal ordering of human embryonic stem cell and induced pluripotent stem cell-derived single-cell gene expression profiles revealed an initially synchronous guidance toward non-neural ectoderm, followed by comparatively asynchronous occurrences of preplacodal and otic marker genes. Positive correlation of marker gene expression between both cell lines and resemblance to mouse embryonic day 10.5 otocyst cells implied reasonable robustness of the guidance protocol. Single-cell trajectory analysis further revealed that otic progenitor cell types are induced in monolayer cultures, but further development appears impeded, likely because of lack of a lineage-stabilizing microenvironment. Our results provide a framework for future exploration of stabilizing microenvironments for efficient differentiation of stem cell-generated human otic cell types.

miR-92a regulates TGF-ß1-induced WISP1 expression in pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is the most common and fatal form of idiopathic interstitial pneumonia. MicroRNAs (miRNAs), short, single-stranded RNAs that regulate protein expression in a post-transcriptional manner, have recently been demonstrated to contribute to IPF pathogenesis. We have previously identified WNT1-inducible signaling pathway protein 1 (WISP1) as a highly expressed pro-fibrotic mediator in IPF, but the underlying mechanisms resulting in increased WISP1 expression, remain elusive. Here, we investigated whether WISP1 is a target of miRNA regulation. We applied a novel supervised machine learning approach, which predicted miR-30a/d and miR-92a target sites in regions of the human WISP1 3'UTR preferentially bound by the miRNA ribonucleoprotein complex. Both miRNAs were decreased in IPF samples, whereas WISP1 protein was increased. We demonstrated further that transforming growth factor (TGF)-ß1-induced WISP1 expression in primary lung fibroblasts in vitro and lung homogenates in vivo. Notably, miR-30a and miR-92a reversed TGF-ß1-induced WISP1 mRNA expression in lung fibroblasts. Moreover, miR-92a inhibition increased WISP1 protein expression in lung fibroblasts. An inverse relationship for WISP1 and miR-92a was found in a TGF-ß1 dependent lung fibrosis model in vivo. Finally, we found significantly increased WISP1 expression in primary IPF fibroblasts, which negatively correlated with miR-92a level ex vivo. Altogether, our findings indicate a regulatory role of miR-92a for WISP1 expression in pulmonary fibrosis.

Rare variants in LRRK1 and Parkinson's disease.

Approximately 20 % of individuals with Parkinson's disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset PD to identify 15 potentially causal variants. Segregation analysis and frequency assessment in 862 PD cases and 1,014 ethnically matched controls highlighted variants in EEF1D and LRRK1 as the best candidates. Mutation screening of the coding regions of these genes in 862 cases and 1,014 controls revealed several novel non-synonymous variants in both genes in cases and controls. An in silico multi-model bioinformatics analysis was used to prioritize identified variants in LRRK1 for functional follow-up. However, protein expression, subcellular localization, and cell viability were not affected by the identified variants. Although it has yet to be proven conclusively that variants in LRRK1 are indeed causative of PD, our data strengthen a possible role for LRRK1 in addition to LRRK2 in the genetic underpinnings of PD but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.

Rare variants in PLXNA4 and Parkinson's disease.

Approximately 20% of individuals with Parkinson's disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset familial PD followed by frequency assessment in 975 PD cases and 1014 ethnically-matched controls and linkage analysis to identify potentially causal variants. Based on the predicted penetrance and the frequencies, a variant in PLXNA4 proved to be the best candidate and PLXNA4 was screened for additional variants in 862 PD cases and 940 controls, revealing an excess of rare non-synonymous coding variants in PLXNA4 in individuals with PD. Although we cannot conclude that the variant in PLXNA4 is indeed the causative variant, these findings are interesting in the light of a surfacing role of axonal guidance mechanisms in neurodegenerative disorders but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.

Profibrotic role of miR-154 in pulmonary fibrosis.

In this study, we explored the regulation and the role of up-regulated microRNAs in idiopathic pulmonary fibrosis (IPF), a progressive interstitial lung disease of unknown origin. We analyzed the expression of microRNAs in IPF lungs and identified 43 significantly up-regulated microRNAs. Twenty-four of the 43 increased microRNAs were localized to the chromosome 14q32 microRNA cluster. We validated the increased expression of miR-154, miR-134, miR-299-5p, miR-410, miR-382, miR-409-3p, miR-487b, miR-31, and miR-127 by quantitative RT-PCR and determined that they were similarly expressed in embryonic lungs. We did not find evidence for differential methylation in this region, but analysis of transcription factor binding sites identified multiple SMAD3-binding elements in the 14q32 microRNA cluster. TGF-ß1 stimulation of normal human lung fibroblasts (NHLF) caused up-regulation of microRNAs on chr14q32 that were also increased in IPF lungs. Chromatin immunoprecipitation confirmed binding of SMAD3 to the putative promoter of miR-154. Mir-154 was increased in IPF fibroblasts, and transfection of NHLF with miR-154 caused significant increases in cell proliferation and migration. The increase in proliferation induced by TGF-ß was not observed when NHLF or IPF fibroblasts were transfected with a mir-154 inhibitor. Transfection with miR-154 caused activation of the WNT pathway in NHLF. ICG-001 and XAV939, inhibitors of the WNT/ß-catenin pathway, reduced the proliferative effect of miR-154. The potential role of miR-154, one of multiple chr14q32 microRNA cluster members up-regulated in IPF and a regulator of fibroblast migration and proliferation, should be further explored in IPF.

Cis-acting polymorphisms affect complex traits through modifications of microRNA regulation pathways.

Genome-wide association studies (GWAS) have become an effective tool to map genes and regions contributing to multifactorial human diseases and traits. A comparably small number of variants identified by GWAS are known to have a direct effect on protein structure whereas the majority of variants is thought to exert their moderate influences on the phenotype through regulatory changes in mRNA expression. MicroRNAs (miRNAs) have been identified as powerful posttranscriptional regulators of mRNAs. Binding to their target sites, which are mostly located within the 3'-untranslated region (3'-UTR) of mRNA transcripts, they modulate mRNA expression and stability. Until today almost all human mRNA transcripts are known to harbor at least one miRNA target site with an average of over 20 miRNA target sites per transcript. Among 5,101 GWAS-identified sentinel single nucleotide polymorphisms (SNPs) that correspond to 18,884 SNPs in linkage disequilibrium (LD) with the sentinels (r2 ≥ 0.8) we identified a significant overrepresentation of SNPs that affect the 3'-UTR of genes (OR = 2.33, 95% CI = 2.12-2.57, P < 10(-52)). This effect was even stronger considering all SNPs in one LD bin a single signal (OR = 4.27, 95% CI = 3.84-4.74, P < 10(-114)). Based on crosslinking immunoprecipitation data we identified four mechanisms affecting miRNA regulation by 3'-UTR mutations: (i) deletion or (ii) creation of miRNA recognition elements within validated RNA-induced silencing complex binding sites, (iii) alteration of 3'-UTR splicing leading to a loss of binding sites, and (iv) change of binding affinity due to modifications of 3'-UTR folding. We annotated 53 SNPs of a total of 288 trait-associated 3'-UTR SNPs as mediating at least one of these mechanisms. Using a qualitative systems biology approach, we demonstrate how our findings can be used to support biological interpretation of GWAS results as well as to provide new experimentally testable hypotheses.

The sufficient minimal set of miRNA seed types.

MOTIVATION: Pairing between the target sequence and the 6-8 nt long seed sequence of the miRNA presents the most important feature for miRNA target site prediction. Novel high-throughput technologies such as Argonaute HITS-CLIP afford meanwhile a detailed study of miRNA:mRNA duplices. These interaction maps enable a first discrimination between functional and non-functional target sites in a bulky fashion. Prediction algorithms apply different seed paradigms to identify miRNA target sites. Therefore, a quantitative assessment of miRNA target site prediction is of major interest. RESULTS: We identified a set of canonical seed types based on a transcriptome wide analysis of experimentally verified functional target sites. We confirmed the specificity of long seeds but we found that the majority of functional target sites are formed by less specific seeds of only 6 nt indicating a crucial role of this type. A substantial fraction of genuine target sites arenon-conserved. Moreover, the majority of functional sites remain uncovered by common prediction methods.