Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution.

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4+, CD8+, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-ß signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

A single nucleotide polymorphism of cyclin-dependent kinase inhibitor 1B (p27Kip1) associated with human vein graft failure affects growth of human venous adventitial cells but not smooth muscle cells.

BACKGROUND: Cyclin-dependent kinase inhibitor 1B (p27Kip1) is a cell-cycle inhibitor whose -838C>A single nucleotide polymorphism (rs36228499; hereafter called p27 SNP) has been associated with the clinical failure of peripheral vein grafts, but the functional effects of this SNP have not been demonstrated. METHODS: Human saphenous vein adventitial cells and intimal/medial smooth muscle cells (SMCs) were derived from explants obtained at the time of lower extremity bypass operations. We determined the following in adventitial cells and SMCs as a function of the p27 SNP genotype: (1) p27 promoter activity, (2) p27 messenger (m)RNA and protein levels, and (3) growth and collagen gel contraction. Deoxyribonuclease I footprinting was also performed in adventitial cells and SMCs. RESULTS: p27 promoter activity, deoxyribonuclease I footprinting, p27 mRNA levels, and p27 protein levels demonstrated that the p27 SNP is functional in adventitial cells and SMCs. Both cell types with the graft failure protective AA genotype had more p27 mRNA and protein. As predicted because of higher levels of p27 protein, adventitial cells with the AA genotype grew slower than those of the CC genotype. Unexpectedly, SMCs did not show this genotype-dependent growth response. CONCLUSIONS: These results support the functionality of the p27 SNP in venous SMCs and adventitial cells, but an effect of the SNP on cell proliferation is limited to only adventitial cells. These data point to a potential role for adventitial cells in human vein graft failure and also suggest that SMCs express factors that interfere with the activity of p27.

Selecting Optimal Peptides for Targeted Proteomic Experiments in Human Plasma Using In Vitro Synthesized Proteins as Analytical Standards.

In targeted proteomics, the development of robust methodologies is dependent upon the selection of a set of optimal peptides for each protein-of-interest. Unfortunately, predicting which peptides and respective product ion transitions provide the greatest signal-to-noise ratio in a particular assay matrix is complicated. Using in vitro synthesized proteins as analytical standards, we report here an empirically driven method for the selection of said peptides in a human plasma assay matrix.

Using Data Independent Acquisition (DIA) to Model High-responding Peptides for Targeted Proteomics Experiments.

Targeted mass spectrometry is an essential tool for detecting quantitative changes in low abundant proteins throughout the proteome. Although selected reaction monitoring (SRM) is the preferred method for quantifying peptides in complex samples, the process of designing SRM assays is laborious. Peptides have widely varying signal responses dictated by sequence-specific physiochemical properties; one major challenge is in selecting representative peptides to target as a proxy for protein abundance. Here we present PREGO, a software tool that predicts high-responding peptides for SRM experiments. PREGO predicts peptide responses with an artificial neural network trained using 11 minimally redundant, maximally relevant properties. Crucial to its success, PREGO is trained using fragment ion intensities of equimolar synthetic peptides extracted from data independent acquisition experiments. Because of similarities in instrumentation and the nature of data collection, relative peptide responses from data independent acquisition experiments are a suitable substitute for SRM experiments because they both make quantitative measurements from integrated fragment ion chromatograms. Using an SRM experiment containing 12,973 peptides from 724 synthetic proteins, PREGO exhibits a 40-85% improvement over previously published approaches at selecting high-responding peptides. These results also represent a dramatic improvement over the rules-based peptide selection approaches commonly used in the literature.

Developmental fate and cellular maturity encoded in human regulatory DNA landscapes.

Cellular-state information between generations of developing cells may be propagated via regulatory regions. We report consistent patterns of gain and loss of DNase I-hypersensitive sites (DHSs) as cells progress from embryonic stem cells (ESCs) to terminal fates. DHS patterns alone convey rich information about cell fate and lineage relationships distinct from information conveyed by gene expression. Developing cells share a proportion of their DHS landscapes with ESCs; that proportion decreases continuously in each cell type as differentiation progresses, providing a quantitative benchmark of developmental maturity. Developmentally stable DHSs densely encode binding sites for transcription factors involved in autoregulatory feedback circuits. In contrast to normal cells, cancer cells extensively reactivate silenced ESC DHSs and those from developmental programs external to the cell lineage from which the malignancy derives. Our results point to changes in regulatory DNA landscapes as quantitative indicators of cell-fate transitions, lineage relationships, and dysfunction.

Personal and population genomics of human regulatory variation.

The characteristics and evolutionary forces acting on regulatory variation in humans remains elusive because of the difficulty in defining functionally important noncoding DNA. Here, we combine genome-scale maps of regulatory DNA marked by DNase I hypersensitive sites (DHSs) from 138 cell and tissue types with whole-genome sequences of 53 geographically diverse individuals in order to better delimit the patterns of regulatory variation in humans. We estimate that individuals likely harbor many more functionally important variants in regulatory DNA compared with protein-coding regions, although they are likely to have, on average, smaller effect sizes. Moreover, we demonstrate that there is significant heterogeneity in the level of functional constraint in regulatory DNA among different cell types. We also find marked variability in functional constraint among transcription factor motifs in regulatory DNA, with sequence motifs for major developmental regulators, such as HOX proteins, exhibiting levels of constraint comparable to protein-coding regions. Finally, we perform a genome-wide scan of recent positive selection and identify hundreds of novel substrates of adaptive regulatory evolution that are enriched for biologically interesting pathways such as melanogenesis and adipocytokine signaling. These data and results provide new insights into patterns of regulatory variation in individuals and populations and demonstrate that a large proportion of functionally important variation lies beyond the exome.

SRC Homology 2 Domain Binding Sites in Insulin, IGF-1 and FGF receptor mediated signaling networks reveal an extensive potential interactome.

Specific peptide ligand recognition by modular interaction domains is essential for the fidelity of information flow through the signal transduction networks that control cell behavior in response to extrinsic and intrinsic stimuli. Src homology 2 (SH2) domains recognize distinct phosphotyrosine peptide motifs, but the specific sites that are phosphorylated and the complement of available SH2 domains varies considerably in individual cell types. Such differences are the basis for a wide range of available protein interaction microstates from which signaling can evolve in highly divergent ways. This underlying complexity suggests the need to broadly map the signaling potential of systems as a prerequisite for understanding signaling in specific cell types as well as various pathologies that involve signal transduction such as cancer, developmental defects and metabolic disorders. This report describes interactions between SH2 domains and potential binding partners that comprise initial signaling downstream of activated fibroblast growth factor (FGF), insulin (Ins), and insulin-like growth factor-1 (IGF-1) receptors. A panel of 50 SH2 domains screened against a set of 192 phosphotyrosine peptides defines an extensive potential interactome while demonstrating the selectivity of individual SH2 domains. The interactions described confirm virtually all previously reported associations while describing a large set of potential novel interactions that imply additional complexity in the signaling networks initiated from activated receptors. This study of pTyr ligand binding by SH2 domains provides valuable insight into the selectivity that underpins complex signaling networks that are assembled using modular protein interaction domains.

Rapid empirical discovery of optimal peptides for targeted proteomics.

We report a method for high-throughput, cost-efficient empirical discovery of optimal proteotypic peptides and fragment ions for targeted proteomics applications using in vitro-synthesized proteins. We demonstrate the approach using human transcription factors, which are typically difficult, low-abundance targets and empirically derived proteotypic peptides for 98% of the target proteins. We show that targeted proteomic assays developed using our approach facilitate robust in vivo quantification of human transcription factors.