Search

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2 + TMPRSS2 + cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Although single-cell RNA sequencing studies have begun to provide compendia of cell expression profiles 1-9 , it has been difficult to systematically identify and localize all molecular cell types in individual organs to create a full molecular cell atlas. Here, using droplet- and plate-based single-cell RNA sequencing of approximately 75,000 human cells across all lung tissue compartments and circulating blood, combined with a multi-pronged cell annotation approach, we create an extensive cell atlas of the human lung. We define the gene expression profiles and anatomical locations of 58 cell populations in the human lung, including 41 out of 45 previously known cell types and 14 previously unknown ones. This comprehensive molecular atlas identifies the biochemical functions of lung cells and the transcription factors and markers for making and monitoring them; defines the cell targets of circulating hormones and predicts local signalling interactions and immune cell homing; and identifies cell types that are directly affected by lung disease genes and respiratory viruses. By comparing human and mouse data, we identified 17 molecular cell types that have been gained or lost during lung evolution and others with substantially altered expression profiles, revealing extensive plasticity of cell types and cell-type-specific gene expression during organ evolution including expression switches between cell types. This atlas provides the molecular foundation for investigating how lung cell identities, functions and interactions are achieved in development and tissue engineering and altered in disease and evolution.

Aging is a pleiotropic process affecting many aspects of mammalian physiology. Mammals are composed of distinct cell type identities and tissue environments, but the influence of these cell identities and environments on the trajectory of aging in individual cells remains unclear. Here, we performed single-cell RNA-seq on >50,000 individual cells across three tissues in young and old mice to allow for direct comparison of aging phenotypes across cell types. We found transcriptional features of aging common across many cell types, as well as features of aging unique to each type. Leveraging matrix factorization and optimal transport methods, we found that both cell identities and tissue environments exert influence on the trajectory and magnitude of aging, with cell identity influence predominating. These results suggest that aging manifests with unique directionality and magnitude across the diverse cell identities in mammals., (© 2019 Kimmel et al.; Published by Cold Spring Harbor Laboratory Press.)

This work presents the workflow for generating single cell transcriptomes derived from primary human uterine and cervical tissue obtained during planned cesarean hysterectomies. In total, a catalogue of 310 single cell transcriptomes are obtained, cell types present in these biopsies are inferred, and specific genes defining each of the cellular types present in the tissue are identified. Further validation of the inferred cell identity is also demonstrated via meta-analysis of independent repositories in literature generated by bulk sequenced data of fluorescence-activated cell sorting sorted cells., (© 2019 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Background: High-fidelity preservation strategies for primary tissues are in great demand in the single cell RNAseq community. A reliable method would greatly expand the scope of feasible multi-site collaborations and maximize the utilization of technical expertise. When choosing a method, standardizability and fidelity are important factors to consider due to the susceptibility of single-cell RNAseq analysis to technical noise. Existing approaches such as cryopreservation and chemical fixation are less than ideal for failing to satisfy either or both of these standards., Results: Here we propose a new strategy that leverages preservation schemes developed for organ transplantation. We evaluated the strategy by storing intact mouse kidneys in organ transplant preservative solution at hypothermic temperature for up to 4 days (6 h, 1, 2, 3, and 4 days), and comparing the quality of preserved and fresh samples using FACS and single cell RNAseq. We demonstrate that the strategy effectively maintained cell viability, transcriptome integrity, cell population heterogeneity, and transcriptome landscape stability for samples after up to 3 days of preservation. The strategy also facilitated the definition of the diverse spectrum of kidney resident immune cells, to our knowledge the first time at single cell resolution., Conclusions: Hypothermic storage of intact primary tissues in organ transplant preservative maintains the quality and stability of the transcriptome of cells for single cell RNAseq analysis. The strategy is readily generalizable to primary specimens from other tissue types for single cell RNAseq analysis.

Background: It remains difficult to predict and to measure the efficacy of pharmacological immunosuppression. We hypothesized that measuring the B-cell repertoire would enable assessment of the overall level of immunosuppression after heart transplantation., Methods and Findings: In this proof-of-concept study, we implemented a molecular-barcode-based immune repertoire sequencing assay that sensitively and accurately measures the isotype and clonal composition of the circulating B cell repertoire. We used this assay to measure the temporal response of the B cell repertoire to immunosuppression after heart transplantation. We selected a subset of 12 participants from a larger prospective cohort study (ClinicalTrials.gov NCT01985412) that is ongoing at Stanford Medical Center and for which enrollment started in March 2010. This subset of 12 participants was selected to represent post-heart-transplant events, with and without acute rejection (six participants with moderate-to-severe rejection and six without). We analyzed 130 samples from these patients, with an average follow-up period of 15 mo. Immune repertoire sequencing enables the measurement of a patient's net state of immunosuppression (correlation with tacrolimus level, r = -0.867, 95% CI -0.968 to -0.523, p = 0.0014), as well as the diagnosis of acute allograft rejection, which is preceded by increased immune activity with a sensitivity of 71.4% (95% CI 30.3% to 94.9%) and a specificity of 82.0% (95% CI 72.1% to 89.1%) (cell-free donor-derived DNA as noninvasive gold standard). To illustrate the potential of immune repertoire sequencing to monitor atypical post-transplant trajectories, we analyzed two more patients, one with chronic infections and one with amyloidosis. A larger, prospective study will be needed to validate the power of immune repertoire sequencing to predict rejection events, as this proof-of-concept study is limited to a small number of patients who were selected based on several criteria including the availability of a large number of samples and the absence or presence of rejection events., Conclusions: If confirmed in larger, prospective studies, the method described here has potential applications in the tailored management of post-transplant immunosuppression and, more broadly, as a method for assessing the overall activity of the immune system.

Antibody heavy chains contain a variable and a constant region. The constant region of the antibody heavy chain is encoded by multiple groups of exons which define the isotype and therefore many functional characteristics of the antibody. We performed both single B cell RNAseq and long read single molecule sequencing of antibody heavy chain transcripts and were able to identify novel exons for IGHA1 and IGHA2 as well as novel isoforms for IGHM antibody heavy chain.

Histocompatibility is the basis by which multicellular organisms of the same species distinguish self from nonself. Relatively little is known about the mechanisms underlying histocompatibility reactions in lower organisms. Botryllus schlosseri is a colonial urochordate, a sister group of vertebrates, that exhibits a genetically determined natural transplantation reaction, whereby self-recognition between colonies leads to formation of parabionts with a common vasculature, whereas rejection occurs between incompatible colonies. Using genetically defined lines, whole-transcriptome sequencing, and genomics, we identified a single gene that encodes self-nonself and determines "graft" outcomes in this organism. This gene is significantly up-regulated in colonies poised to undergo fusion and/or rejection, is highly expressed in the vasculature, and is functionally linked to histocompatibility outcomes. These findings establish a platform for advancing the science of allorecognition.

Botryllus schlosseri is a colonial urochordate that follows the chordate plan of development following sexual reproduction, but invokes a stem cell-mediated budding program during subsequent rounds of asexual reproduction. As urochordates are considered to be the closest living invertebrate relatives of vertebrates, they are ideal subjects for whole genome sequence analyses. Using a novel method for high-throughput sequencing of eukaryotic genomes, we sequenced and assembled 580 Mbp of the B. schlosseri genome. The genome assembly is comprised of nearly 14,000 intron-containing predicted genes, and 13,500 intron-less predicted genes, 40% of which could be confidently parceled into 13 (of 16 haploid) chromosomes. A comparison of homologous genes between B. schlosseri and other diverse taxonomic groups revealed genomic events underlying the evolution of vertebrates and lymphoid-mediated immunity. The B. schlosseri genome is a community resource for studying alternative modes of reproduction, natural transplantation reactions, and stem cell-mediated regeneration. DOI:http://dx.doi.org/10.7554/eLife.00569.001.

The human antibody repertoire is one of the most important defenses against infectious disease, and the development of vaccines has enabled the conferral of targeted protection to specific pathogens. However, there are many challenges to measuring and analyzing the immunoglobulin sequence repertoire, including that each B cell's genome encodes a distinct antibody sequence, that the antibody repertoire changes over time, and the high similarity between antibody sequences. We have addressed these challenges by using high-throughput long read sequencing to perform immunogenomic characterization of expressed human antibody repertoires in the context of influenza vaccination. Informatic analysis of 5 million antibody heavy chain sequences from healthy individuals allowed us to perform global characterizations of isotype distributions, determine the lineage structure of the repertoire, and measure age- and antigen-related mutational activity. Our analysis of the clonal structure and mutational distribution of individuals' repertoires shows that elderly subjects have a decreased number of lineages but an increased prevaccination mutation load in their repertoire and that some of these subjects have an oligoclonal character to their repertoire in which the diversity of the lineages is greatly reduced relative to younger subjects. We have thus shown that global analysis of the immune system's clonal structure provides direct insight into the effects of vaccination and provides a detailed molecular portrait of age-related effects.

It is thought that the adaptive immune system of immature organisms follows a more deterministic program of antibody creation than is found in adults. We used high-throughput sequencing to characterize the diversifying antibody repertoire in zebrafish over five developmental time points. We found that the immune system begins in a highly stereotyped state with preferential use of a small number of V (variable) D (diverse) J (joining) gene segment combinations, but that this stereotypy decreases dramatically as the zebrafish mature, with many of the top VDJ combinations observed in 2-wk-old zebrafish virtually disappearing by 1 mo. However, we discovered that, in the primary repertoire, there are strong correlations in VDJ use that increase with zebrafish maturity, suggesting that VDJ recombination involves a level of deterministic programming that is unexpected. This stereotypy is masked by the complex diversification processes of antibody maturation; the variation and lack of correlation in full repertoires between individuals appears to be derived from randomness in clonal expansion during the affinity maturation process. These data provide a window into the mechanisms of VDJ recombination and diversity creation and allow us to better understand how the adaptive immune system achieves diversity.

Mapping genes that underlie complex genetic traits, including genes that determine susceptibility to common diseases, requires an efficient method for high-resolution genotyping. Single-nucleotide differences between pairs of allelic sequences from unrelated individuals occur approximately once in every kilobase. Genomic mismatch scanning (GMS), by analyzing numerous single-nucleotide polymorphisms in a single genome-wide step, offers a potentially powerful and efficient approach to linkage analysis. GMS, originally developed in a yeast system, is shown here to be applicable to the more complex mouse and human genomes.

The development and progression of cancer and the experimental reversal of tumorigenicity are accompanied by complex changes in patterns of gene expression. Microarrays of cDNA provide a powerful tool for studying these complex phenomena. The tumorigenic properties of a human melanoma cell line, UACC-903, can be suppressed by introduction of a normal human chromosome 6, resulting in a reduction of growth rate, restoration of contact inhibition, and suppression of both soft agar clonogenicity and tumorigenicity in nude mice. We used a high density microarray of 1,161 DNA elements to search for differences in gene expression associated with tumour suppression in this system. Fluorescent probes for hybridization were derived from two sources of cellular mRNA [UACC-903 and UACC-903(+6)] which were labelled with different fluors to provide a direct and internally controlled comparison of the mRNA levels corresponding to each arrayed gene. The fluorescence signals representing hybridization to each arrayed gene were analysed to determine the relative abundance in the two samples of mRNAs corresponding to each gene. Previously unrecognized alterations in the expression of specific genes provide leads for further investigation of the genetic basis of the tumorigenic phenotype of these cells.

A disease which has blinded 50,000 Americans and threatens to blind an additional estimated one million victims who are unaware that they have the disease warrants serious consideration in the health curriculum. That disease, glaucoma, can be neither prevented nor cured. It results when abnormal fluid dynamics within the eye cause increased pressure. The increased pressure can usually be controlled to prevent blindness from occurring. However, prevention of blindness from glaucoma requires early diagnosis and treatment. Since glaucoma is rarely accompanied by pain or warning signs, the best hope for early diagnosis is with regular eye examinations. Health education is the way to teach people the nature and consequences of glaucoma and the role each individual can play in minimizing his or her chances of becoming blind from the disease.