Scuphr: A probabilistic framework for cell lineage tree reconstruction.

Cell lineage tree reconstruction methods are developed for various tasks, such as investigating the development, differentiation, and cancer progression. Single-cell sequencing technologies enable more thorough analysis with higher resolution. We present Scuphr, a distance-based cell lineage tree reconstruction method using bulk and single-cell DNA sequencing data from healthy tissues. Common challenges of single-cell DNA sequencing, such as allelic dropouts and amplification errors, are included in Scuphr. Scuphr computes the distance between cell pairs and reconstructs the lineage tree using the neighbor-joining algorithm. With its embarrassingly parallel design, Scuphr can do faster analysis than the state-of-the-art methods while obtaining better accuracy. The method's robustness is investigated using various synthetic datasets and a biological dataset of 18 cells.

Clonally heritable gene expression imparts a layer of diversity within cell types.

Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information.

Long-read whole-genome analysis of human single cells.

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

Divergent clonal differentiation trajectories establish CD8+ memory T cell heterogeneity during acute viral infections in humans.

The CD8+ T cell response to an antigen is composed of many T cell clones with unique T cell receptors, together forming a heterogeneous repertoire of effector and memory cells. How individual T cell clones contribute to this heterogeneity throughout immune responses remains largely unknown. In this study, we longitudinally track human CD8+ T cell clones expanding in response to yellow fever virus (YFV) vaccination at the single-cell level. We observed a drop in clonal diversity in blood from the acute to memory phase, suggesting that clonal selection shapes the circulating memory repertoire. Clones in the memory phase display biased differentiation trajectories along a gradient from stem cell to terminally differentiated effector memory fates. In secondary responses, YFV- and influenza-specific CD8+ T cell clones are poised to recapitulate skewed differentiation trajectories. Collectively, we show that the sum of distinct clonal phenotypes results in the multifaceted human T cell response to acute viral infections.

Expansions of adaptive-like NK cells with a tissue-resident phenotype in human lung and blood.

Human adaptive-like "memory" CD56dimCD16+ natural killer (NK) cells in peripheral blood from cytomegalovirus-seropositive individuals have been extensively investigated in recent years and are currently explored as a treatment strategy for hematological cancers. However, treatment of solid tumors remains limited due to insufficient NK cell tumor infiltration, and it is unknown whether large expansions of adaptive-like NK cells that are equipped for tissue residency and tumor homing exist in peripheral tissues. Here, we show that human lung and blood contains adaptive-like CD56brightCD16- NK cells with hallmarks of tissue residency, including expression of CD49a. Expansions of adaptive-like lung tissue-resident NK (trNK) cells were found to be present independently of adaptive-like CD56dimCD16+ NK cells and to be hyperresponsive toward target cells. Together, our data demonstrate that phenotypically, functionally, and developmentally distinct subsets of adaptive-like NK cells exist in human lung and blood. Given their tissue-related character and hyperresponsiveness, human lung adaptive-like trNK cells might represent a suitable alternative for therapies targeting solid tumors.

Unique transcriptional and protein-expression signature in human lung tissue-resident NK cells.

Human lung tissue-resident NK cells (trNK cells) are likely to play an important role in host responses towards viral infections, inflammatory conditions and cancer. However, detailed insights into these cells are still largely lacking. Here we show, using RNA sequencing and flow cytometry-based analyses, that subsets of human lung CD69+CD16- NK cells display hallmarks of tissue-residency, including high expression of CD49a, CD103, and ZNF683, and reduced expression of SELL, S1PR5, and KLF2/3. CD49a+CD16- NK cells are functionally competent, and produce IFN-γ, TNF, MIP-1ß, and GM-CSF. After stimulation with IL-15, they upregulate perforin, granzyme B, and Ki67 to a similar degree as CD49a-CD16- NK cells. Comparing datasets from trNK cells in human lung and bone marrow with tissue-resident memory CD8+ T cells identifies core genes co-regulated either by tissue-residency, cell-type or location. Together, our data indicate that human lung trNK cells have distinct features, likely regulating their function in barrier immunity.

Conbase: a software for unsupervised discovery of clonal somatic mutations in single cells through read phasing.

Accurate variant calling and genotyping represent major limiting factors for downstream applications of single-cell genomics. Here, we report Conbase for the identification of somatic mutations in single-cell DNA sequencing data. Conbase leverages phased read data from multiple samples in a dataset to achieve increased confidence in somatic variant calls and genotype predictions. Comparing the performance of Conbase to three other methods, we find that Conbase performs best in terms of false discovery rate and specificity and provides superior robustness on simulated data, in vitro expanded fibroblasts and clonal lymphocyte populations isolated directly from a healthy human donor.

Thyroid hormone drives the expression of mouse carbonic anhydrase Car4 in kidney, lung and brain.

Thyroid hormone is a well-known regulator of brain, lung and kidney development and function. However, the molecular mechanisms by which the hormone exerts its function have remained largely enigmatic, and only a limited set of target genes have been identified in these tissues. Using a mouse model with a mutation in thyroid hormone receptor α1 (TRα1), we here demonstrate that the expression of carbonic anhydrase 4 in lung and brain of the adult animal depends on intact TRα1 signaling. In the kidney, carbonic anhydrase 4 mRNA and protein are not affected by the mutant TRα1, but are acutely repressed by thyroid hormone. However, neither lung function--as measured by respiration rate and oxygen saturation--nor urine pH levels were affected by altered carbonic anhydrase 4 levels, suggesting that other carbonic anhydrases are likely to compensate. Taken together, our findings identify a previously unknown marker of TRα1 action in brain and lung, and provide a novel negatively regulated target gene to assess renal thyroid hormone status.

Transplanted Bone Marrow-Derived Cells Contribute to Human Adipogenesis.

Because human white adipocytes display a high turnover throughout adulthood, a continuous supply of precursor cells is required to maintain adipogenesis. Bone marrow (BM)-derived progenitor cells may contribute to mammalian adipogenesis; however, results in animal models are conflicting. Here we demonstrate in 65 subjects who underwent allogeneic BM or peripheral blood stem cell (PBSC) transplantation that, over the entire lifespan, BM/PBSC-derived progenitor cells contribute ∼10% to the subcutaneous adipocyte population. While this is independent of gender, age, and different transplantation-related parameters, body fat mass exerts a strong influence, with up to 2.5-fold increased donor cell contribution in obese individuals. Exome and whole-genome sequencing of single adipocytes suggests that BM/PBSC-derived progenitors contribute to adipose tissue via both differentiation and cell fusion. Thus, at least in the setting of transplantation, BM serves as a reservoir for adipocyte progenitors, particularly in obese subjects.