Considerations for Using the Vasculature as a Coordinate System to Map All the Cells in the Human Body

Introduction

When examining coordinate systems that could be used for different biological entities, specific characteristics should be considered. For example, chromosome number and base pair position form a linear coordinate system that is sufficient to define location in the reference genome (1). However, a coordinate system for the entire human body, which can serve the needs of HuBMAP, HCA, HPA, and other mapping efforts, is far more complex for several reasons (2–5). First, cells live in three dimensions; and, additional information is needed to describe the size and rotational orientation of any tissue sample used to generate the data. Second, while genome base pair position does not depend on how DNA is folded and packed into the nucleus, the positions of individual elements (e.g., cells, organs) within the body are dynamic, as they may change with skeletal and muscle movements, gravity, breathing, heart beating, and other functions and forces that distort tissue. Third, although the variability of DNA between people is <1%, human body geometry spans a wide range of heights, weights, and shapes, which vary by sex and race and change significantly over a person's lifespan.

One option is to use a 3D Cartesian coordinate system to describe position within the body. Cartesian coordinate systems have the benefit of being familiar and easily understood. With the body in standard anatomical position, left-right, caudal-cranial, and posterior-anterior form three perpendicular (x-y-z) axes. For smaller scale analyses, centering the origin on an anatomic landmark might be most useful and potentially less variable between different people. For example, it would be easier to measure the distance of a tissue sample from the kidney to the superior pole of the kidney than to the top of a person's head. The disadvantage of Cartesian coordinates is that its axes do not follow the natural shape of the body. As a result, creating a reference body map by digitally “stitching” together tissue specimens collected from people of different sizes would require a complex process of warping, rotating, and aligning the data. Doing this in a way that scales to the entire body while maintaining cellular-level precision would be extremely challenging.

Other 3D coordinate systems have been developed for different applications. For example, spherical coordinates, with axes representing latitude, longitude, and elevation, conveniently describe location on the surface of the Earth. A simple path in spherical coordinates, such as “1 km east,” is much more difficult to communicate using Cartesian coordinates. Similarly, coordinate systems for the human body should ideally have axes that follow anatomical structures, genetic expression patterns, chemical gradients, and/or other biologically relevant pathways. Several organ-specific coordinate systems have been developed, such as Talairach coordinates for functional brain imaging, but these do not extend to the entire body (6).

Vascular Based Coordinate System

In 2017, a Common Coordinate Framework (CCF) Meeting was organized by the National Institutes of Health, the Broad Institute, the Sanger Institute, and the Chan Zuckerberg Initiative. Meeting participants, which included anatomists, pathologists, clinicians, and technology experts from around the world, suggested several approaches for mapping the cells of the human body, including anatomically based coordinate systems (7). One of these, which uses the vasculature, was further discussed at a 2019 CCF Workshop, specifically for localizing kidney cells (8). In this manuscript, we present a conceptual overview of this vascular based coordinate system and describe its benefits and limitations (Figure 1).

FIGURE 1

Figure 1. The hub-and-spoke structure of a vascular coordinate system. The vascular coordinate system consists of representative vessel loops that begin and end in the chambers of the heart and extend to the functional units of each organ (counterclockwise from upper-left: alveoli, glomeruli, and hepatic lobules). The “address” of a cell includes of the name of the loop (e.g., “R” for renal and “H” for hepatic) and the branching level of the nearest vessel (e.g., “RA1” for afferent arterioles in the renal loop). Note that the walls of the heart lie along the coronary (“C”) vascular loop, and the cells of some structures internal to the heart are actually positioned on the periphery of the coordinate system, near the capillaries that provide them with blood. There are multiple ways to pinpoint an even more precise “GPS-like” position, including, going clockwise from upper-right, the “hypoxia fingerprint” that results from decreasing oxygen levels farther from the vessel, the unique “histologic fingerprint” and “gene expression fingerprint” in different types of vascular endothelial cells, and the distinctive “vascular architecture fingerprint” found in different types of tissue. The vasculature extends to all parts of the body (bottom) and frames all organs at all scales. Source Files: Images of the hepatic lobules, glomeruli, alveoli, gene expression heat map, and whole body vasculature are adapted from Wikimedia Commons: https://commons.wikimedia.org/wiki/File:2423_Microscopic_Anatomy_of_Liver.jpg (Public Domain); https://commons.wikimedia.org/wiki/File:Juxtaglomerular_Apparatus_and_Glomerulus.jpg (CC BY 3.0 license); https://commons.wikimedia.org/wiki/File:Alveolus_diagram.svg (Public Domain); https://commons.wikimedia.org/wiki/File:Iris_dendrogram.png (CC BY-SA 4.0 license); and, https://commons.wikimedia.org/wiki/File:Circulatory_System_en.svg (Public Domain). Gene expression heat map recreated from Guo et al. (9) using simulated data. Vascular architecture based on Bosetti et al. (10) (CC BY 4.0 license).

The goal of a vascular coordinate system is not to build a map of the vasculature. Rather, the idea is to use known vascular pathways through the body as an axis in a coordinate system that can describe the position of anatomical elements, such as cells, in the tissue surrounding the vasculature. In this proposed framework, the chambers of the heart may be viewed as the origin in the whole-body map. The vascular axis contains vascular loops (comparable to “roads” in the Google Map metaphor) that extend out through the aorta or pulmonary artery, gradually taper off to single cell size as they approach every other cell in every organ in the body, and return to the heart. Along these pathways, the vasculature follows the unique biology of the various organs and tissues; and, likewise, the organs and tissues of the body could not exist without the unique properties of their closely integrated vasculature.

Just as the genome map unfolds the DNA into a linear sequence, we can imagine unfolding the complex 3D twists and turns of the vasculature axis into a simpler 2D “hub-and-spoke” shape with the chambers of the heart in the center. This schematic representation of the vasculature makes it easier to (1) describe location in the body; (2) align vascular pathways with varying 3D shapes but equivalent function from different people; (3) identify patterns such as changes in cell type and gene expression as one transitions from larger to smaller vessels within an organ (along a spoke); and, (4) compare tissue at the same level (e.g., capillaries) in different organs (across spokes).

An organ, such as the kidney, has tens of thousands of vasculature pathways if every capillary is considered distinctly. However, in order to make the construction of a CCF feasible, we will make an additional simplifying assumption that many of these pathways are indistinguishable from each other. Notably, the smallest vessels (capillaries) are at the center of the small specialized functional units in each organ (e.g., hepatic lobules in the liver, glomeruli in the kidney, and alveoli in the lung). A whole organ requires many of these vascular-centered functional units that are anatomically and physiologically similar (11). This similarity can be used to collapse numerous vascular pathways into a few representative ones (Figure 2). In the future, as we learn more about the biomolecular profiles of these pathways, we can refine the CCF by splitting the pathways into any subtypes that are discovered.

FIGURE 2

Figure 2. Unfolding the renal vascular pathway down to the single cell level. A representative vascular pathway enters the kidney through the renal artery, passes through the glomerulus, and returns out the renal vein. At the macroscale (left) and mesoscale (center), labels in yellow boxes are anatomical structures that correspond to different vessels along this loop. At the microscale (right), labels in yellow boxes are different types of cells that are within a short distance of a nearby vessel in the glomerulus. The position of individual cells can be described by the nearest vessel's address (e.g., “RA0” for glomerulus capillary), the distance in micrometers or number of cells from the vessel endothelium and an angle perpendicular to the length of the vessel (e.g., “RA0-10 μm-135°” or “RA0-1c-135°”). On the far right, each vessel and its surrounding tissue has been extracted from the original image and aligned vertically to make it easy to see how cell types and distribution change along the vascular pathway. In this example, the aligned images have also been rotated so that the center of the glomerulus is always on the right. Source Files: Images of the kidney and glomerulus are adapted from Wikimedia Commons: https://commons.wikimedia.org/wiki/File:KidneyStructures_PioM.svg (CC BY 3.0 license); and, https://commons.wikimedia.org/wiki/File:Renal_corpuscle-en.svg (CC BY-SA 4.0 license).

An “address” of an individual cell in a vascular coordinate system has four components. (1) The first is the nearest anatomically named vessel that most precisely identifies a vascular loop (the vascular territory or “watershed”), such as the right renal artery. (2) The next component represents the branching level of the vessel. In the kidney, examples from proximal to distal include the interlobar arteries, cortical radiate arteries, afferent arterioles, and glomerular capillaries. The branching level can also be defined numerically, with capillaries at the zero level, and higher levels with larger numbers. (3) The third component is the perpendicular distance in micrometers or number of cells from the endothelial layer of the vessel. (4) An optional fourth component could be the rotational angle perpendicular to the length of the vessel. Zero degrees can be defined many ways depending on the application. For example, it can be relative to a marker on the tissue specimen, or it can point toward an anatomical structure, such as the center of the nearest glomerulus in the kidney. The utility and need for these components would need to be confirmed by experimental testing.

Another approach is to use endothelial cells as “anchoring cells” to identify the relative positions of the other cells in a GPS-like fashion. Endothelium lines the entire vasculature. Distinct endothelial cell (EC) types appear in different parts of the vasculature, which can be detected histologically through unique cell characteristics and neighboring cells (the “histologic fingerprint”), as well as through analysis of their specific gene profiles (the “gene expression fingerprint,” Figure 1 heat map) (9, 12, 13). Vascular organizational pattern, branching complexity, diameter, etc. provide additional recognizable features that form a “vascular architecture fingerprint” for each tissue (Figure 1 false-color fluorescent vasculature images) (10, 14). Each tissue has specific oxygen gradients, which can be used to determine distance from the vessel. When PO2 cannot be measured directly, one can infer the cell distance from the induction/expression gradients of hypoxia-inducible genes in various types of cells (the “hypoxia fingerprint,” Figure 1 transverse view of a vessel) (15–17).

Benefits of a Vascular CCF

A vasculature coordinate system makes sense biologically. Every living cell must be within a small radius of the closest blood vessel (100 μm to 1 mm, depending on the tissue) in order to receive oxygen (18); and, every vessel is lined by the same continuous layer of endothelium (19). Thus, the vasculature forms an unbroken pathway that reaches all parts of the body, seamlessly tapering down across scales, from macroscale (whole body/clinical) to mesoscale to microscale (single cell-size capillary vessels). As part of this, it handles the gradual transition from large arteries and veins, which are conserved in most people, down to the millions of microscopic vessels that are recognizable by category (e.g., glomerulus capillary) but not by individual names. The vascular system adapts to individual body size and shape, accounting for inter-individual variations (20, 21). It defines the shape of functional units in different organs, such as liver lobules, kidney glomeruli, and lung alveoli. Normal development of new tissue begins with vessel formation before additional cell types can grow around it (22).

The vascular system also reaches all the organs, tissues, and cells in a contiguous fashion. Other anatomical structures do not seem to have similar versatility. For example, the skeletal system extends to the entire body; however, the distance from a cell within an organ to the nearest bone can be several centimeters, and that bone is unlikely to be part of the organ, and even less likely to be part of a tissue specimen. The nervous system follows many of the same pathways as larger vessels invested with contractile properties, but it is not known whether they extend to follow the smaller vessels that reach all tissue cells. Thus, unlike cells' dependence on oxygen, there isn't a guarantee that every cell is in close proximity to a nerve. Organ-specific reference systems are based on unique structural characteristics that exist only in that individual organ (e.g., specific regions of the brain) and are unlikely to be useful for representing cell locations elsewhere in the body.

There are many practical benefits of a vascular coordinate system. Vascular pathways in all organs have been studied extensively and are described in the literature in detail. They are useful, well-known, and have standardized names across many specialty domains. Tissue vascularization patterns are used clinically to diagnose disease; and, vessels are used in surgery and biopsies as the main anatomical landmarks and to define “vascular territories” (23, 24). Although the exact position in 3D Cartesian space might be difficult to determine for a resected tissue specimen, the surgeon will typically be able to indicate which vascular watershed included the specimen.

Limitations of a Vascular CCF

Only knowing the base pairs in a genome sequence does not provide enough information to determine how the DNA is folded in a particular cell. Similarly, although the unfolded vascular coordinate system can describe the 3D spatial relationship of nearby cells at the microscale, the winding paths that vessels take and the relative distance between cells at larger scales is lost. Information about larger scale structures could be obtained by integrating clinical imaging and histological and molecular patterning.

Some additional positional information is lost by collapsing vascular pathways. For example, a vascular coordinate system might initially treat all glomerulus capillaries in the kidney as the same structure, potentially masking differences in glomeruli at the superior and inferior renal poles. Though, as we learn more about the biomolecular profiles of these pathways, we can iteratively refine the vascular CCF over time by splitting pathways into any subtypes that are discovered.

Despite these limitations, a vascular coordinate system retains enough positional information to answer many types of research questions, by allowing researchers to precisely localize individual cells within functional units, tissues and organs, or compare the effects of context on various cell types across the whole body at the macro- and meso- scales.

Next Steps and Future Directions

In Figure 1 we present a highly simplified and stylized representation of the vascular CCF in order to introduce the general concept. The actual vascular pathways are far more complex, with the full vascular CCF eventually consisting of hundreds or thousands of structures as well as branching patterns that do not follow a strict hub-and-spoke layout. For example, the veins in the hepatic portal system, which carry blood from multiple organs to the liver, and circulatory anastomoses, such as the circle of Willis in the brain, would run perpendicular to and cross over the “spokes” in Figure 1.

The vascular CCF will also need to incorporate normal anatomical variants of the vasculature, such as the thyroidea ima and double cystic arteries, in the same way that single nucleotide polymorphisms (SNPs) describe normal variations in the genome. As we expand and fill in the details of the vascular CCF, it might take on a more snowflake or web-like appearance, with probabilities assigned to different sections to indicate the likelihood of finding them in a given person.

The actual vascular CCF will require a more sophisticated coding scheme for naming vessels than what is shown in Figure 1, in order to represent its complexity more accurately and make it “computable” by software algorithms. One approach is to build on existing anatomical ontologies, such as UBERON, which already contains unique Uniform Resource Identifiers (URIs) for many vessels and defines various relationship types, such as “branching_part_of,” to describe how vessels connect (25).

A cell can be nearby more than one vessel, in which case it is possible to assign it a separate address relative to each vessel. Future research is needed to determine what the appropriate “cutoff” distance should be. It might be desirable to store multiple addresses for cells to help triangulate their position in 3D space.

Although our initial motivation for a vascular CCF is to define a coordinate system for specifying the location of cells, the novel radial layout of the human body around the vascular, as depicted in Figure 1, could have other scientific applications. For example, positional information within this framework has analogies to embryologic development and might provide insight into mechanisms that cause disease or dysregulations that arise from a group of cells. For instance, in the vascular CCF, macrovascular diseases associated with type 2 diabetes (e.g., ischemic heart disease, peripheral vascular disease, and cerebrovascular disease) are clustered near the center, while the microvascular diseases (e.g., retinopathy, nephropathy, and neuropathy) form a ring around the outside (26).

Summary

The vascular system has several properties that respond to the characteristics previously deemed as desirable for a CCF to enable mapping all the cells in the human body (7). (1) It works across several scales. Through a continuous endothelium it makes a seamless transition from large anatomical structures down to the cellular level. (2) It is applicable to all body tissues. Because organs develop around vessels, the vasculature frames organ architecture at all scales. Vascular patterns are so distinctive that images of vessels alone, with all other cells removed, can easily be used to identify the organ. (3) It accounts for donor differences. Vascularization naturally adapts to individual variations in tissue, body size and shape. Instead of having to map the 3D coordinates of each glomerulus in one person to those in another person, the vascular CCF would position each cell within a representative glomerulus capillary from each person. (4) The vasculature is well-known to practitioners of many different clinical specialties, helping both diagnose disease and guide surgeries, and used as landmarks by pathologists, radiologists, and other clinicians. This would facilitate positioning tissue samples in the body and registering the location of individual cells. While no single CCF is ideal for all use cases, a vascular CCF defines a natural coordinate system that would make it easy to combine biomolecular data from multiple people and ask biologically relevant research questions.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.proteinatlas.org.

Author Contributions

GW conceived this study based on reports from the 2017 CCF Meeting and 2019 CCF Workshop. KB and YJ provided critical feedback. GW wrote the initial draft, with additional text transcribed from the CCFWS-01 slide presentation. All authors contributed to editing the manuscript.

Funding

This research was supported by National Institutes of Health (NIH) award OT2OD026671. This work does not represent the opinion of the NIH or the National Heart, Lung, and Blood Institute (NHLBI).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Zorina S. Galis from the National Institutes of Health for helping develop the idea of a vascular-based coordinate system at the 2017 CCF Meeting and 2019 CCF Workshop, providing us with presentations from these events, and participating in numerous discussions that formed the basis of this manuscript. Expert comments by the reviewers helped clarify and detail the proposed vascular based common coordinate system.

References

1. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, et al. Initial sequencing and analysis of the human genome. Nature. (2001) 409:860–921. doi: 10.1038/35057062

PubMed Abstract | CrossRef Full Text | Google Scholar

2. HuBMAP Consortium. The human body at cellular resolution: the NIH Human Biomolecular Atlas Program. Nature. (2019) 574:187–92. doi: 10.1038/s41586-019-1629-x

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Uhlén M, Björling E, Agaton C, Szigyarto CA, Amini B, Andersen E, et al. A human protein atlas for normal and cancer tissues based on antibody proteomics. Mol Cell Proteomics. (2005) 4:1920–32. doi: 10.1074/mcp.M500279-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Rozenblatt-Rosen O, Stubbington MJT, Regev A, Teichmann SA. The Human Cell Atlas: from vision to reality. Nature. (2017) 550:451–3. doi: 10.1038/550451a

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Regev A, Teichmann SA, Lander ES, Amit I, Benoist C, Birney E, et al. The human cell atlas. Elife. (2017) 6:e27041. doi: 10.7554/eLife.27041

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Talairach J, Szikla G. Application of stereotactic concepts to the surgery of epilepsy. Acta Neurochir Suppl. (1980) 30:35–54. doi: 10.1007/978-3-7091-8592-6_5

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Common Coordinate Framework (CCF) Meeting. (2017). Bethesda, M. D. Meeting Summary. Available online at: https://commonfund.nih.gov/sites/default/files/CCF%20summary%20final%20.pdf (accessed October 29, 2019).

8. Common Coordinate Framework Workshop (CCFWS-01). (2019). Bloomington, IN: Indiana University. Available online at: https://ccfws.cns.iu.edu/ (accessed October 29, 2019).

9. Guo S, Zhou Y, Xing C, Lok J, Som AT, Ning M, et al. The vasculome of the mouse brain. PLoS ONE. (2012) 7:e52665. doi: 10.1371/journal.pone.0052665

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Bosetti F, Galis ZS, Bynoe MS, Charette M, Cipolla MJ, del Zoppa GJ, et al. (2016). “Small blood vessels: big health problems?”: scientific recommendations of the National Institutes of Health Workshop. J Am Heart Assoc. 5:e004389. doi: 10.1161/JAHA.116.004389

PubMed Abstract | CrossRef Full Text | Google Scholar

11. de Bono B, Grenon P, Baldock R, Hunter P. Functional tissue units and their primary tissue motifs in multi-scale physiology. J Biomed Semantics. (2013) 4:22. doi: 10.1186/2041-1480-4-22

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Vanlandewijck M, He L, Mäe MA, Andrae J, Ando K, Del Gaudio F, et al. A molecular atlas of cell types and zonation in the brain vasculature. Nature. (2018) 554:475–80. doi: 10.1038/nature25739

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Jourde-Chiche N, Fakhouri F, Dou L, Bellien J, Burtey S, Frimat M, et al. Endothelium structure and function in kidney health and disease. Nat Rev Nephrol. (2019) 15:87–108. doi: 10.1038/s41581-018-0098-z

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Johnson C, Sung HJ, Lessner SM, Fini ME, Galis ZS. Matrix metalloproteinase-9 is required for adequate angiogenic revascularization of ischemic tissues: potential role in capillary branching. Circ Res. (2004) 94:262–8. doi: 10.1161/01.RES.0000111527.42357.62

PubMed Abstract | CrossRef Full Text | Google Scholar

15. De Santis V, Singer M. Tissue oxygen tension monitoring of organ perfusion: rationale, methodologies, and literature review. Br J Anaesth. (2015) 115:357–65. doi: 10.1093/bja/aev162

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Patnaik MM, Tefferi A. The complete evaluation of erythrocytosis: congenital and acquired. Leukemia. (2009) 23:834–44. doi: 10.1038/leu.2009.54

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Koch CJ. Measurement of absolute oxygen levels in cells and tissues using oxygen sensors and 2-nitroimidazole EF5. Methods Enzymol. (2002) 352:3–31. doi: 10.1016/S0076-6879(02)52003-6

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tsai AG, Johnson PC, Intaglietta M. Oxygen gradients in the microcirculation. Physiol Rev. (2003) 83:933–63. doi: 10.1152/physrev.00034.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hacking C, Gaillard F. (2019). Aorta. Radiopaedia.org. Available online at: https://radiopaedia.org/articles/aorta (accessed Sep 29, 2019).

20. Paruchuri V, Salhab KF, Kuzmik G, Gubernikoff G, Fang H, Rizzo JA, et al. Aortic size distribution in the general population: explaining the size paradox in aortic dissection. Cardiology. (2015) 131:265–72. doi: 10.1159/000381281

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Abadi E, Segars WP, Sturgeon GM, Roos JE, Ravin CE, Samei E. Modeling lung architecture in the XCAT series of phantoms: physiologically based airways, arteries and veins. IEEE Trans Med Imaging. (2018) 37:693–702. doi: 10.1109/TMI.2017.2769640

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Si-Tayeb K, Lemaigre FP, Duncan SA. Organogenesis and development of the liver. Dev Cell. (2010) 18:175–89. doi: 10.1016/j.devcel.2010.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Taylor GI, Palmer JH. The vascular territories (angiosomes) of the body: experimental study and clinical applications. Br J Plast Surg. (1987) 40:113–41. doi: 10.1016/0007-1226(87)90185-8

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Majno P, Mentha G, Toso C, Morel P, Peitgen HO, Fasel JH. Anatomy of the liver: an outline with three levels of complexity–a further step towards tailored territorial liver resections. J Hepatol. (2014) 60:654–62. doi: 10.1016/j.jhep.2013.10.026

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Mungall CJ, Torniai C, Gkoutos GV, Lewis SE, Haendel MA. Uberon, an integrative multi-species anatomy ontology. Genome Biol. (2012) 13:R5. doi: 10.1186/gb-2012-13-1-r5

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Cade WT. Diabetes-related microvascular and macrovascular diseases in the physical therapy setting. Phys Ther. (2008) 88:1322–35. doi: 10.2522/ptj.20080008

PubMed Abstract | CrossRef Full Text | Google Scholar