Cell type discovery and representation in the era of high-content single cell phenotyping

Background

A fundamental characteristic of multicellular organisms is the specialization of functional cell types through the process of differentiation. These specialized cell types not only characterize the normal functioning of different organs and tissues, they can also be used as cellular biomarkers of a variety of different disease states and therapeutic/vaccine responses. In order to serve as a reference for cell type representation, the Cell Ontology has been developed to provide a standard nomenclature of defined cell types for comparative analysis and biomarker discovery. Historically, these cell types have been defined based on unique cellular shapes and structures, anatomic locations, and marker protein expression. However, we are now experiencing a revolution in cellular characterization resulting from the application of new high-throughput, high-content cytometry and sequencing technologies. The resulting explosion in the number of distinct cell types being identified is challenging the current paradigm for cell type definition in the Cell Ontology.

Results

In this paper, we provide examples of state-of-the-art cellular biomarker characterization using high-content cytometry and single cell RNA sequencing, and present strategies for standardized cell type representations based on the data outputs from these cutting-edge technologies, including “context annotations” in the form of standardized experiment metadata about the specimen source analyzed and marker genes that serve as the most useful features in machine learning-based cell type classification models. We also propose a statistical strategy for comparing new experiment data to these standardized cell type representations.

Conclusion

The advent of high-throughput/high-content single cell technologies is leading to an explosion in the number of distinct cell types being identified. It will be critical for the bioinformatics community to develop and adopt data standard conventions that will be compatible with these new technologies and support the data representation needs of the research community. The proposals enumerated here will serve as a useful starting point to address these challenges.

     

Interaction network of the ribosome assembly machinery from a eukaryotic thermophile

Ribosome biogenesis in eukaryotic cells is a highly dynamic and complex process innately linked to cell proliferation. The assembly of ribosomes is driven by a myriad of biogenesis factors that shape pre‐ribosomal particles by processing and folding the ribosomal RNA and incorporating ribosomal proteins. Biochemical approaches allowed the isolation and characterization of pre‐ribosomal particles from Saccharomyces cerevisiae, which lead to a spatiotemporal map of biogenesis intermediates along the path from the nucleolus to the cytoplasm. Here, we cloned almost the entire set (～180) of ribosome biogenesis factors from the thermophilic fungus Chaetomium thermophilum in order to perform an in‐depth analysis of their protein–protein interaction network as well as exploring the suitability of these thermostable proteins for structural studies. First, we performed a systematic screen, testing about 80 factors for crystallization and structure determination. Next, we performed a yeast 2‐hybrid analysis and tested about 32,000 binary combinations, which identified more than 1000 protein–protein contacts between the thermophilic ribosome assembly factors. To exemplary verify several of these interactions, we performed biochemical reconstitution with the focus on the interaction network between 90S pre‐ribosome factors forming the ctUTP‐A and ctUTP‐B modules, and the Brix‐domain containing assembly factors of the pre‐60S subunit. Our work provides a rich resource for biochemical reconstitution and structural analyses of the conserved ribosome assembly machinery from a eukaryotic thermophile.     

Leucine‐rich repeat receptor‐like kinase II phylogenetics reveals five main clades throughout the plant kingdom

Receptor‐like kinases (RLKs) represent the largest group of cell surface receptors in plants. The monophyletic leucine‐rich repeat (LRR)‐RLK subfamily II is considered to contain the somatic embryogenesis receptor kinases (SERKs) and NSP‐interacting kinases known to be involved in developmental processes and cellular immunity in plants. There are only a few published studies on the phylogenetics of LRR‐RLKII; unfortunately these suffer from poor taxon/gene sampling. Hence, it is not clear how many and what main clades this family contains, let alone what structure–function relationships exist. We used 1342 protein sequences annotated as ‘SERK’ and ‘SERK‐like’ plus related sequences in order to estimate phylogeny within the LRR‐RLKII clade, using the nematode protein kinase Pelle as an outgroup. We reconstruct five main clades (LRR‐RLKII 1–5), in each of which the main pattern of land plant relationships re‐occurs, confirming previous hypotheses that duplication events happened in this gene subfamily prior to divergence among land plant lineages. We show that domain structures and intron–exon boundaries within the five clades are well conserved in evolution. Furthermore, phylogenetic patterns based on the separate LRR and kinase parts of LRR‐RLKs are incongruent: whereas the LRR part supports a LRR‐RLKII 2/3 sister group relationship, the kinase part supports clades 1/2. We infer that the kinase part includes few ‘radical’ amino acid changes compared with the LRR part. Finally, our results confirm that amino acids involved in each LRR‐RLKII–receptor complex interaction are located at N‐capping residues, and that the short amino acid motifs of this interaction domain are highly conserved throughout evolution within the five LRR‐RLKII clades.     

Potential motion related bias in the worn dosimeter measurements of two childhood leukemia studies

Time averaged field measurements produced by a Positron dosimeter worn by study subjects was the primary method of exposure evaluation in two Canadian studies of childhood leukemia and AC magnetic field exposure. Statistically significant but mutually contradictory results obtained in the two studies, done in different locales but under similar study conditions, have not been explained. This report examines operational features of the Positron meter, including an unanticipated sensitivity to wearer motion. If the convalescent cases studied were less active than their healthy controls, as one might expect, then the meter's characteristic responses to motion, particularly as they would affect case-control distributions above and below the different referent group cutpoints used in the two studies, could help to explain both the unexpected inverse risks reported in the larger study and the unusually high risks reported in the smaller study.     

Sustained effects of biofeedback-assisted relaxation therapy in essential hypertension

The usefulness of biofeedback-assisted relaxation as an adjunct or substitute for pharmacotherapy in essential hypertension can be enhanced if the effects are shown to persist after formal treatment has ended. Patients with essential hypertension successfully treated with biofeedback-assisted relaxation were recalled for follow-up yearly after the termination of treatment. Twenty-six of 40 patients met the BP criterion for success. At one-, two-, and three-year follow-up, 31%, 38%, and 27% of the successful completers continued to meet the criterion for success. The pretreatment-posttreatment decreases in BP were accompanied by decreases in forehead muscle tension and urinary cortisol. Forehead muscle tension, urinary cortisol, and anxiety levels were significantly lower than pretreatment one year after the end of treatment. Self-report data were used to assess continued relaxation practice. No relationship was found between practice and any other dependent measure. It appears that some patients trained in biofeedback-assisted relaxation can maintain lowered blood pressure, muscle tension, anxiety, and cortisol levels over the long term; however, the role of relaxation practice in maintaining these lowered levels remains unclear.     

Restoration of the CD4 T cell compartment after long-term highly active antiretroviral therapy without phenotypical signs of accelerated immunological aging

It remains uncertain whether full T cell reconstitution can be established in HIV-infected children and adults with long-term sustained virological control by highly active antiretroviral therapy (HAART). In this study, we comprehensively analyzed various phenotypical markers of CD4 T cell recovery. In addition to measuring T cell activation and proliferation markers, CD4 T cell generation and aging of the CD4 T cell compartment were assessed by measuring TCR excision circles and the fraction of CD31-expressing naive CD4 T cells. In all children and in adults with relatively high CD4 T cell counts at start of therapy (>200 cells/microl), total CD4 T cell numbers normalized within 1 year of therapy. After long-term HAART (4.4-9.6 years), naive CD4 T cell counts had normalized in both groups. Although in adults with low baseline CD4 T cell counts (<200 cells/microl) total CD4 T cell numbers normalized eventually after at least 7 years of HAART, naive CD4 T cell counts had still not recovered. TCR excision circle data showed that thymic T cell production contributed to naive T cell recovery at all ages. The fraction of CD31-expressing naive CD4 T cells was found to be normal, suggesting that the CD4 T cell repertoire was diverse after long-term HAART. Hence, under sustained viral suppression during long-term HAART, the T cell compartment has the potential to fully recover by generating new naive T cells both in children and in adults with high baseline CD4 T cells counts. Irrespective of baseline CD4 T cell counts, reconstitution occurred without a significant effect on T cell aging as reflected by markers for replicative history.     

Immune-complex-trapping cells in the spleen of the oriental fire-bellied toad,Bombina orientalis

Summary The spleen of the oriental fire-bellied toad, Bombina orientalis, consists of well-developed white pulp, separated from the lymphocytic marginal zone by the connective tissue boundary layer. Injection of peroxidase-conjugated rabbit anti-peroxidase revealed that these immune complexes were localized on the surface of acid-phosphatase-positive and non-specific-esterasepositive cells in the white pulp. The majority of immunecomplex-trapping cells were present around the blood vessels. Cell processes of some of these cells penetrated into the wall of blood vessels. The significance of the present findings is discussed with respect to the evolution of immune-complex-trapping cells in the spleen.     

Beitrag zur Flora von Bosnien und der Hercegovina

Ohne Zusammenfassung     

Quantification of hypoxic regions distant from occlusions in cerebral penetrating arteriole trees

The microvasculature plays a key role in oxygen transport in the mammalian brain. Despite the close coupling between cerebral vascular geometry and local oxygen demand, recent experiments have reported that microvascular occlusions can lead to unexpected distant tissue hypoxia and infarction. To better understand the spatial correlation between the hypoxic regions and the occlusion sites, we used both in vivo experiments and in silico simulations to investigate the effects of occlusions in cerebral penetrating arteriole trees on tissue hypoxia. In a rat model of microembolisation, 25 μm microspheres were injected through the carotid artery to occlude penetrating arterioles. In representative models of human cortical columns, the penetrating arterioles were occluded by simulating the transport of microspheres of the same size and the oxygen transport was simulated using a Green’s function method. The locations of microspheres and hypoxic regions were segmented, and two novel distance analyses were implemented to study their spatial correlation. The distant hypoxic regions were found to be present in both experiments and simulations, and mainly due to the hypoperfusion in the region downstream of the occlusion site. Furthermore, a reasonable agreement for the spatial correlation between hypoxic regions and occlusion sites is shown between experiments and simulations, which indicates the good applicability of in silico models in understanding the response of cerebral blood flow and oxygen transport to microemboli.