Limitations on the number of unique protein and DNA molecules that can be characterized microscopically in a single tissue specimen impede advances in understanding the biological basis of health and disease. Here we present a multiplexed fluorescence microscopy method (MxIF) for quantitative, single-cell, and subcellular characterization of multiple analytes in formalin-fixed paraffinembedded tissue. Chemical inactivation of fluorescent dyes after each image acquisition round allows reuse of common dyes in iterative staining and imaging cycles. The mild inactivation chemistry is compatible with total and phosphoprotein detection, as well as DNA FISH. Accurate computational registration of sequential images is achieved by aligning nuclear counterstain-derived fiducial points. Individual cells, plasma membrane, cytoplasm, nucleus, tumor, and stromal regions are segmented to achieve cellular and subcellular quantification of multiplexed targets. In a comparison of pathologist scoring of diaminobenzidine staining of serial sections and automated MxIF scoring of a single section, human epidermal growth factor receptor 2, estrogen receptor, p53, and androgen receptor staining by diaminobenzidine and MxIF methods yielded similar results. Single-cell staining patterns of 61 protein antigens by MxIF in 747 colorectal cancer subjects reveals extensive tumor heterogeneity, and cluster analysis of divergent signaling through ERK1/2, S6 kinase 1, and 4E binding protein 1 provides insights into the spatial organization of mechanistic target of rapamycin and MAPK signal transduction. Our results suggest MxIF should be broadly applicable to problems in the fields of basic biological research, drug discovery and development, and clinical diagnostics.cancer diagnostics | high-content cellular analysis | image analysis | mTOR | multiplexing

No abstract

SummaryEpithelial organ morphogenesis involves reciprocal interactions between epithelial and mesenchymal cell types to balance progenitor cell retention and expansion with cell differentiation for evolution of tissue architecture. Underlying submandibular salivary gland branching morphogenesis is the regulated proliferation and differentiation of perhaps several progenitor cell populations, which have not been characterized throughout development, and yet are critical for understanding organ development, regeneration, and disease. Here we applied a serial multiplexed fluorescent immunohistochemistry technology to map the progressive refinement of the epithelial and mesenchymal cell populations throughout development from embryonic day 14 through postnatal day 20. Using computational single cell analysis methods, we simultaneously mapped the evolving temporal and spatial location of epithelial cells expressing subsets of differentiation and progenitor markers throughout salivary gland development. We mapped epithelial cell differentiation markers, including aquaporin 5, PSP, SABPA, and mucin 10 (acinar cells); cytokeratin 7 (ductal cells); and smooth muscle α-actin (myoepithelial cells) and epithelial progenitor cell markers, cytokeratin 5 and c-kit. We used pairwise correlation and visual mapping of the cells in multiplexed images to quantify the number of single- and double-positive cells expressing these differentiation and progenitor markers at each developmental stage. We identified smooth muscle α-actin as a putative early myoepithelial progenitor marker that is expressed in cytokeratin 5-negative cells. Additionally, our results reveal dynamic expansion and redistributions of c-kit- and K5-positive progenitor cell populations throughout development and in postnatal glands. The data suggest that there are temporally and spatially discreet progenitor populations that contribute to salivary gland development and homeostasis.

A method for flexible fitting of molecular models into three-dimensional electron microscopy (3D-EM) reconstructions at a resolution range of 8-12 A is proposed. The approach uses the evolutionarily related structural variability existing among the protein domains of a given superfamily, according to structural databases such as CATH. A structural alignment of domains belonging to the superfamily, followed by a principal components analysis, is performed, and the first three principal components of the decomposition are explored. Using rigid body transformations for the secondary structure elements (SSEs) plus the cyclic coordinate descent algorithm to close the loops, stereochemically correct models are built for the structure to fit. All of the models are fitted into the 3D-EM map, and the best one is selected based on crosscorrelation measures. This work applies the method to both simulated and experimental data and shows that the flexible fitting was able to produce better results than rigid body fitting.

The challenges faced in analyzing optical imaging data from neurons include a low signal-to-noise ratio of the acquired images and the multiscale nature of the tubular structures that range in size from hundreds of microns to hundreds of nanometers. In this paper, we address these challenges and present a computational framework for an automatic, three-dimensional (3D) morphological reconstruction of live nerve cells. The key aspects of this approach are: (i) detection of neuronal dendrites through learning 3D tubular models, and (ii) skeletonization by a new algorithm using a morphology-guided deformable model for extracting the dendritic centerline. To represent the neuron morphology, we introduce a novel representation, the Minimum Shape-Cost (MSC) Tree that approximates the dendrite centerline with sub-voxel accuracy and demonstrate the uniqueness of such a shape representation as well as its computational efficiency. We present extensive quantitative and qualitative results that demonstrate the accuracy and robustness of our method.