Abstract

Current advances in staining and imaging of biological tissues provide a wealth of data, but there are few metrics to quantitatively describe these findings. In particular there is a need for a metric that would characterize the correlation and consistency of orientation of different biological constructs within a tissue. We aimed to create such a metric and to demonstrate its use with images of cardiac tissues. The co-orientational order parameter (COOP) was based on the mathematical framework of a classical parameter, the orientational order parameter (OOP). Theorems were proven to illustrate the properties and boundaries of the COOP, which was then applied to both synthetic and experimental data. We showed the COOP to be useful for quantifying the correlation of orientation of constructs such as actin filaments and sarco-meric Z-Iines. As expected, cardiac tissues showed perfect correlation between actin filaments and Z-Iines. We also demonstrated the use of COOP to quantify the consistency of construct orientation within cells of the same shape. The COOP provides a quantitative tool to characterize tissues beyond co-Iocalization or single construct orientation distribution. In the future, this new parameter could be used to represent the quantitative changes during maturation of cardiac tissue, pathological malformation, and other processes.

Author Summary

In tissue engineering, recreating the in vivo architecture is an important aspect of in vitro experiments. An ability to quantify organization of cellular constructs, specifically the correlation in their orientations, would greatly enhance both our understanding of the function of each construct and provide a tool to evaluate engineered tissues. In this work, we have developed a parameter to evaluate the correlation and consistency of construct orientation. We have extensively characterized this co-orientational order parameter analytically, validated it using synthetic data, and demonstrated its use with experimental data of Z-lines and actin fibrils architectures in engineered cardiac tissues. As long as the orientation angles and location of constructs are known, the parameter can quantify both orientational correlation and consistency. Thus, the co-orientational order parameter has a Wide scope of application both in biology and potentially outside of it.

Introduction

The study of architecture is therefore key to understanding how the cellular microenvironment potentiates function, and may provide new insights in the study of physiological mechanisms. Furthermore, for proper function, different components of the cytoskeleton, cell, or tissue need to co-localize and orient properly with respect to each other [5, 6]. Quantifying the degree of orientation of cells and subcellular components, both relative to themselves and to other components, is thus crucial for evaluating the quality of engineered tissues [7].

The quantification of orientation from biological images is in principle straightforward and can be either done manually [8] or with a variety of computational algorithms [6, 9—11]. As far as the second problem is concerned, summarizing the overall organization after image analysis involves selecting a metric, which is more controversial. As a result, a wide variety of metrics are utilized in the bio-imag-ing field. For example, some assume that the parameter can be described as the standard deviation of a truncated Gaussian, or normal, distribution [12]. Others use the von Mises circular distribution [9, 10], which is a wrapped normal distribution. However, cellular and cytoskele-ton distributions are often non-Gaussian, and their being non-Gaussian may be of crucial importance [13].

The OOP is a mathematical construct developed to quantify the degree of order in anisotropic medias [14]. Mathematically, the OOP is equivalent to resultant vector length from the circular distribution with a period of 71’ [15]. In biology, the OOP has been successfully employed to characterize organization of bacteria [16] , fibroblasts [17] , vascular smooth muscle [18], actin fibrils alignment in valve endothelial cells [19], and Z-lines in cardiac muscle [20]. However, there is a lack of a robust correlation metric that has been characterized for use with biological images. The suite of correlation parameters provided by circular statistics are either too limited to be used with cytoskeleton organization or so complex the results are hard to connect back to biological phenotypes [21]. Other correlation metrics are also not ideal for correlating orientations of the cytoskeleton components, and to date no metric has been developed or specifically characterized for this purpose [22—24].

As an example, we apply this Co-Orientational Order Parameter (COOP) to compare orientation of Z-lines and actin filaments in a neonatal rat ventricular myocyte (NRVM) monolayer. Lastly, we show how the COOP can be used to measure the consistency of building a cardiomyocyte on a triangular island of extracellular matrix (ECM).

Results

Theoretical Results

In designing the new metric, we aimed to overcome the challenge of analyzing orientation of multiple pseudo vectors, i.e. the metric needed to be symmetric with the period of 71’ for both vectors. The OOP was designed to analyze the organization of pseudo vectors, and has become a standard parameter for use in liquid crystals [20]. The OOP ranges from zero, for isotropic, to one, for aligned mediums (81 Fig), and it has been applied to various biological systems [14, 16, 20]. However, the OOP was not designed to evaluate the correlation of orientation of coupled constructs such as actin fibrils and Z-lines. The first step in creating the COOP was to formally define the problem. Let the first construct be P, a set of pseudo vectors 15;, and the second construct be Q, a set of pseudo vectors (Fig. 1A). The order tensor and OOP of each field is: defined a new field F (a set of pseudo vectors ): Physically, field F represent the angle (0) between the two biological constructs, and The metric was then calculated similarly to the COP: The analytical solution of the COOP is:

The director represents the mean angle between the two constructs. Alternatively, the COOP can be written in the expended form: The COOP was designed to range between zero and one. Here we present a series of theorems that illustrate the various properties of the COOP. Table 1. Symmetry permutations.

Symmetry is an important characteristic of both the GOP and the COOP because it alleviates calculation errors that may arise when there is a random choice of signs for the pseudo vectors (Fig. 1B). Symmetry can be easily shown for GOP (Sl Supplemental Text). As can be seen from the anyaltical solution of the onstrate the pseudo-symmetry of the COOP.

All of these can be reduced to the same equation with no difference in sign, which proves pseudo-vector symmetry of the COOP. We also demonstrated that the COOP is symmetric to the switch of P and Q:

There are four valid results for an angles between two pseudo-vectors: 90, —90, 71—90, 90—71 (Fig. 1B). For any symmetry permutation the director will correspond to one of these four angles. However, it is essential that the translation from f1 to 90 is handled with this symmetry in mind.

To verify that rotation of Q with respect to P does not affect COOP (Fig. 1C), let If each was rotated by angle 1/, the rotated field th would be defined by: Then the rotated field Frat is: As the angle 1/ is constant:

As a result COOPmt = COOP. Thus we have proven that field rotation does not affect COOP, and without loosing generality we can assume 1% 2 fig 2 [1,0]. This proves that the mean p angle between fibers cannot be used to evaluate the correlation of orientations.

We next proved that the same field compared to itself would obtain a COOP of one (Fig. 1D). Imagine two sets of pseudo vectors distributed in a 2D space, and Let Then, the field of this system can be written as: The mean tensor: Therefore, the COOP of constructs P and Q: We obtained a COOP of 1 thus constructs P and Q are a perfectly co-oriented system Which is expected as P = Q.

For any given pair of fields With OOPp and OOPQ there is a range of possible COOP values, With a maximum range of zero to one. To aid in interpreting the meaning of the COOP, we derived the value for the COOP (COOPu) as a function of the OOPs for the maximally uncorrelated system. For the COOP to be uncorrelated the two fields needed to have no correlation between their orientations (Fig. 1E). If and only if and are assumed independent and f1 p 2 f1 q (Theorem 2), the analytical solution Equation (7) can Using the classical probability distribution property: Based on the definition of the COP the uncorrelated COOP is: COOPu is the lowest limit of the COOP if and only if the two constructs are independent. However, if P and Q are correlated in multiple ways it is possible to achieve a smaller value of

Theorem 5: Anti-correlated COOP. To prove that COOPu is not necessarily the minimum COOP, we constructed an example where P was composed of p: and p3 , and construct Q was composed of q: and q}: The mean tensor and COOP of this system are:

. . }, COOPu>O. Thus, the COOP can be lower than COOPu. Such a case can be graphically imagined if you have constructs correlated in two different ways in the same cell (Fig. 1F), similar to a single cardiomyocyte with punctate alpha-actinin at the ends and well defined Z-lines in the middle. This demonstrates that, similarly to the OOP, the COOP will not fully capture second order correlations (81 Supplemental Text). Theorem 6: Correlated COOP limit. For the upper limit, we derived the value of the COOP as a function of the OOPs for the maXimally correlated constructs. To determine the correlated COOP (COOPC), we assumed the two fields, P and Q, were almost identical, except that Q was rotated by a random noise angle 9 (Fig. 1G). We can assume, without loosing gener-As 9 was assumed to be random noise generated, 9 and a are independent, and therefore:

OOPPCOOPC.

Solving for COOPC and rewriting it in a more general form we obtained: This is the upper limit of COOP if the two constructs are correlated but are subject to random biological variance (noise). This would not be the limit in a system where the variance is not random.

We also showed that the correlated COOP is not necessarily the maximum. To prove this we defined P and Q as: If the expression in Equation (36) is positive, then COOP 2 COOPC. An example of such a case is When 6 is constant and Equation (36) becomes:

Based on Theorems 4—7 we can diVide the range of COOP into three regions (Fig. 11): anti-correlated, normal, and ultra-correlated. The boundaries of these regions

Experimentally, it might be more relevant to know how close the COOP is to the uncorrelated and correlated boundaries as these carry biological implications. We therefore defined a Normalized COOP, Which is a measure of how close a parameter is to COOPu or COOPC:

11). Any parameter has limitations. For example, the OOP of the red line-segments in the anti-correlat-ed schematic (Fig. 1F) would be equivalent to an isotropic organization (OOP ~ 0) even though the red line-segments are well organized. Similarly, the COOP cannot be used to identify correlations in tissues that have very nontrivial spatially dependent correlations (more complex than the ultra-correlated case (Fig. 1H)). We believe that most biological constructs for which the COOP has been developed will never exhibit this behavior. However, if the COOP is ever found to be statistically significantly greater than COOPC, it will be essential to reevaluate the applicability of the parameter.

To interpret the information provided by the COOP, we would need to know which region our tissue falls under. Our ability to do so will be limited by the error inherent in any measurement and the width of the normal COOP range. COOPu approaches the maximum (COOPu = 1) as both OOPp and OOPq approach one (Fig. 1](i)). COOPC, however, approaches its maximum when the two order parameters are close to being equal (Fig. 1](ii)). This shows, for example, that if OOPp = 0 and OOPq = 0, then the normal region of the COOP ranges [0, 1]. However, if OOPp>O and OOPq = O, the normal range does not exist and a COOP>O would indicate an ultra-correlated system. If the error in the system is so large that there is no statistically significant difference between the boundaries (COOPu and COOPC), it would not be possible to differentiate between the regions. Therefore, for the parameter to be useful, the maximum allowable error and minimum sample size have to be experimentally realistic. To estimate the error and sample size, we calculated the propagation of error in COOPu and COOPC, and used them in the student t-test to calculate statistical significance. Assuming OOPp and OOPQ are normally distributed with the standard deviation UOOPP and O'OOPQ, respectively, and OOPP and OOPQ are independent, then the variance: and the standard deviation: Assuming OOPP>OOPQ, the variance for the correlated COOP: formulated the estimate for the the standard deviation of COOPC:

We calculated the t-value and degrees of freedom to complete the t-test. We assumed significance for a p-value less than 0.05 for the two sample t-test. There are two variables that impact significance as a function of OOPs: the error (000p) and the sample size (N). For visualization, we calculated the maximum allowable experimental error (max (00013)) at sample size of four (Fig. 1](iii)). The maximal allowable error at N = 4 was at its highest, max(aoop) = 0.18, for OOPP = OOPQ = 0.60, which is greater than experimental error reported for GOP organization [13]. While the maximal allowable error approaches zero as COP—>0, it rapidly increases for OOP>0.

The minimum sample size ranges from two to infinity (Fig. 1](iv)). Naturally, when COOPC—>COOPu, the min(N) —> 00 as it is not possible to find them significantly different. For 000p = 0.04, the minimum sample size is between two and five for most of the COP values. For convenience, we also calculated the minimum sample size for a range of higher experimental errors (82 Fig.). Obviously higher error would necessitate more samples to maintain significance. However, errors reported for GOP experimentally correspond to normal requirements in sample size. The errors and sample sizes were confirmed to be experimentally realistic, thus we next moved to testing the parameter with synthetic and experimental data.

Synthetic Results

The first step to validate the new parameter was to construct four limiting cases that should lead to specific COOP values. We constructed a custom MatLab code that could be interfaced with experimental or synthetic (computer generated) data. The code was used to confirm the COOP for four synthetic limiting cases. The first test was to compare two perfectly aligned constructs (Fig. 2A(i)). It was clear from the image that these two constructs were perfectly correlated, and we therefore expected COOP = 1. We first confirmed that both constructs were perfectly aligned with OOPP = OOPQ = 1 (Fig. 2A(ii)). Analytically, it is clear that COOP = 1, by first applying Theorem 2 and then Theorem 3. When the synthetic data for this condition was analyzed by the code, the COOP was confirmed to be one (Fig. 2A(iii)). It is worthwhile to note that the director generated by the COOP code corresponds to the angle, 90, between the constructs.

2B(i)). In this case, there was no correlation, and we expected COOP = O. Analytically, even if we supposed that the two constructs were maximally correlated (COOP = COOPC) and applied Theorem 6 based on OOPQ = 1 and OOPp = 0 (Fig. 2B(ii)), the COOP = O. This was confirmed by the result of the synthetic case (Fig. 2B(iii)). Synthetic Image

A Perfect

No

C

No

D Perfect

Synthetic data for four limiting cases. A) Two perfectly organized samples are always perfectly correlated; (B) A perfectly organized and an isotropic construct cannot be correlated; Two isotropic cases can be (C) completely uncorrelated and (D) perfectly correlated. For (AD) (i) Schematic of a small section of the synthetic data; (ii) GOP and normalized occurrence for both constructs; (iii) COOP results with 00 designating the mean angle between constructs.

2C (i)-(ii)). We expected the COOP for uncorrelated constructs to be zero, and analytically this was confirmed by Theorem 4. These findings were also confirmed by the results of the synthetic data (Fig. 2C(iii)).

2D(i)-(ii)). We eXpected the COOP to be one as there was perfect correlation, and this was proven by Theorem 6. Synthetically we showed the COOP = 1, and again the director gave the average angle, 90, between constructs. The limiting conditions validated the parameter for four simple cases. However, the interpretation of the COOP gets more complex When neither construct is perfectly aligned or isotropic.

To understand the limits of the COOP parameter, we constructed a series of cases With different organizations by truncation of Gaussian distributions With specified standard deViations (Fig. 3A and 3B). We created synthetic data for the uncorrelated case (dark blue Fig. 3C) by generating two separate random number sets using the appropriate truncated Gaussian distribution for each. For the correlated case (brown Fig. 3C), we generated the first, more organized, data set by the same method. Then the noise was generated such that When it was added to the first data set, the new set would have the target distribution. Both methods of creating P and Q data sets lead to the same desired indiVidual distributions (Fig. 3A and 3B). To Visualize the results we constructed a slider, sketched in Fig. 11, for each case (Fig. 3C) color-coded to indicate the boundaries, COOPu (dark blue) and COOPC (brown), as well as the three regions: anti-correlated (light blue), normal (ranging from dark blue to brown), and ultra-correlated (bright green). We eXpected that When OOPP = 1 (Fig. 3A(i)), COOPu = COOPC = OOPQ for all OOPQ values, Which was confirmed by the results of the synthetic data (Fig. 3C(i)-(iii)). Also, if one OOP = 0 (Fig. 3B(i)) and the other OOP 7S 0, the COOPu = COOPC = O. This was also confirmed With synthetic results (Fig. 3C(i) and 3C(iv)). For a case Where OOP 7S {0 or 1} the normal range of the COOP was greater if the OOPs were smaller and closer to each other, Which can be seen by comparing Fig. 3C(V) and Fig. 3C(Vi). These synthetic results confirm the results of Fig. 1, Which showed that the COOP parameter is best applied in situations Where the correlation is not obviously dictated by the organization of the indiVidual constructs.

Experimental Results

We used actin fibers and sarcomeric Z-lines in NRVMs to validate the parameter and code. The tissues were stained for a-actinin and phalloidin to identify Z-line and actin fibril directions, respectively (Fig. 4A). The program we used to identify the direction of the construct was based on a fingerprint identification code [6, 25, 26] , and it assigned a direction to every nonempty pixel in the image. However, computing correlation of constructs based on individual pixels introduced too many errors. Indeed, most images of cytoskeleton constructs are obtained via immunostaining and imaging. The accuracy of the COOP will be a direct consequence of image quality. If the images are of poor quality (poor contrast, dead cells, etc. ), it will not be possible to accurately extract the construct direction data, and thus, the COOP will not be accurate. However, if it is possible to accurately extract directionality data, the COOP can be used. The images we have collected for the proof-of-principle data set are representative of the images normally used to study tissue architecture [6, 20]. The images may be resolvable to different degrees, and image acquisition procedures can introduce inaccuracies at smaller scales. For example, during the collection of this data, it is customary to ensure that each channel results in sharp images of the corresponding construct, such as actin and Z-lines. This sometimes requires focusing on slightly different planes, and as a result, the images of the same field-of-view could be slightly offset from each other. This along with other imaging inaccuracies lead us to develop a procedure to average out these small errors by calculating the direction of each construct within grid-squares. While any consistent small grid can be used for comparing results across multiple tissues, we recommend picking a grid size that would correspond to a natural biological unit. In this case, the image was partitioned into a grid (Fig. 4B), which was chosen such that at minimum, two Z-lines could be expected to fit within each grid-square, ensuring that at least one “sarcomere” complex is within each square (S3 Fig.). The spacing between Z-lines in NRVM tissues is 1.9—2.1 um [27], thus we chose a grid size close to 4.2 pm (~ 30 pixels). For each grid (1'), we calculated the average direction for both constructs (i.e. Z-lines or actin fibrils) (Fig. 4C). To account for varying densities and partial grid-squares at the edges of the images, the area density (p) was calculated for each grid-square using the number of non zero angles in the grid-square divided by its area. Each grid-square (i) was assigned a weight (W1) based on the OOP and density of constructs (p):

Additionally, grid-squares with better alignment have higher weights than grid-squares with isotropic organization, which prevents loss of consistency.

However, histological samples may not be perfect, with some eXample imperfections identified in the zoomed-in image of Fig. 4C, which correspond to the non-perpendicular arrows identified in Fig. 4D. Note that the lengths of the arrows correspond to the weights assigned to each square, so these imperfections can significantly alter the resulting COOP. To test the parameter, we took four coverslips with 10 fields of view imaged for each and identified noiseless regions with minimal imperfections (Fig. 4E-F). Before any angle detection was done, Image] was used to merge fields of view containing Z-lines and actin f1-brils, and regions with minimal imperfections were chosen by an experienced user. Only merged images that contained four or more regions of minimal imperfections were used (Fig. 4E). This was done prior to the organizational analysis to eliminate bias. The implementation of the code is summarized in a flow chart (S4 Fig. ).

We measured no significant difference in the values of COOPu or COOPC between the full and noiseless data sets (Fig. 4G). There was also no significant difference between COOP and COOPC for the masked, noiseless images. In contrast there was a significant difference between COOP and COOPC (p<0.001) for the raw images. This is reflected in the Normalized COOP which is 1 for the noiseless images and less than 0.5 for the full data (Fig. 4H). This illustrates that while the parameter is capable of capturing the eXpected correlation in a Table 2. Pairwise comparison of sarcomeric Z-lines consistency. Cell 1 Cell 2 Cell 3 Cell 4 Cell 5 Cell 1 Cell 2 Cell 3 Cell 4 Cell 5 sarcomere between actin fibril and Z-line orientation, cardiac tissues may have imperfections that result in a lower value of COOP. The mean angle for the noiseless images also shows the expected perpendicular correlation between orientations of sarcomere Z-lines and actin fibrils with minimal error (Fig. 4H). While, they were also approximately perpendicular in the raw images, the error was greater than for the noiseless images.

As an example we seeded NRVMs on triangular islands, and stained the samples for nucleus, Z-lines, actin fibrils, and fibronectin (Fig. 41). For analysis we ensured that the fibronectin islands lined up for all five cells (SS Fig). Then the Z-line (solid bars) or actin fibril (dashed bars) images were compared in a pairwise manner (Tables 2 and 3). This showed an additional experimental confirmation of Theorem 3: COOP = 1 for each cell when it was compared with itself (Fig. 41). The COOP was calculated for the same grid as in the isotropic images. The results showed that although there is an overall consistency between cells, i.e., myofibrils were bundled along the edges of the triangle, the orientation was not fully consistent at a smaller length scale (grid size). Indeed, the average COOP for ten pairwise comparisons

4K). This demonstrated another potential use of this parameter.

Discussion

For example, second-order stereology has been used to analyze spatial arrangements of constructs in images. Noorafshan et al. used second-order stereology to examine the correlation between the spatial arrangements of cardiomyocytes and microvessels [22]. Their method involved pair correlation and cross-correlation functions to determine positive or negative correlation at different distances. However, neither simple colocalization nor second-order stere-ology do not analyze the relative orientations of biological constructs. Table 3. Pairwise comparison of actin fibrils consistency. COOP actin Cell 1 Cell 2 Cell 3 Cell 4 Cell 5 Cell 1 Cell 2 Cell 3 Cell 4 Cell 5

The COOP would allow for investigation of mechanisms or functions that rely on not only spatial proximity, but also specific organizational schemes. To properly interpret the meaning of the parameter values, we characterized it through a series of analytical theorems. As a result, we defined three regimes—normal, ultra-correlated, and anti-correlated—that have biological implications. After validating the parameter with synthetic data, we demonstrated its use with experimental images by showing that perfect portions of cardiac tissues have the expected correlation of the orientation between sarcomeric z-discs and actin fibrils. The reduction in the COOP for unmasked (raw) data suggests that the defects in the architecture will be distinguished by our new method (Fig. 4G-H). The code we have developed can also be used to calculate the mean angle between constructs thus allowing for tracking of mean angle changes as a function of experimental conditions. Furthermore, the parameter can be used to calculate the consistency of orientational organization to help evaluate the importance of orientational order.

The best metrics place the least number of constraints on the distribution of orientations. For example, the standard deviation is not an appropriate metric for quantifying orientation distribution of Z-lines as they are not distributed normally. The OOP works with all types of distributions, and it has an additional benefit of being symmetric to pseudo-vectors [14]. As the COOP was designed with similar math, it shares the same benefits as the OOP such as pseudo-vector symmetry.

For example, Feng et al. use normalized cross-correlation to compare two images with a possible rotation or change in scale [23]. Their method involves identifying a relatively small number of points of interest and matching the comparison based on them. The Feng et al. method is insensitive to the rotation of the whole image (i.e. rotation plus translation), while the COOP method is insensitive to the rotation of all vectors without translation. The normalized cross-correlation is a powerful tool, but not appropriate whole cell architectural metric, the COOP is therefore useful for comparing consistency in similarly shaped cells with matching ECM islands, but cannot be used to identify the same cell that has been re-scaled and rotated. Another example of mathematical tools for image analysis is a set of a nonparametric circular statistics tests such as Watson’s U3 test, which is designed to evaluate the probability that a sample comes from a specific distribution or that two samples come from the same distribution. For instance, nonparametric circular statistics has been utilized to evaluate if a pattern of migration of different objects is the same [30, 31]. However, these tests do not consider the location of each sample pair, thus while the results can correspond to the COOP in very special situations these parameters address fundamentally different questions. Thus, these circular statistics tests are not a good tool to evaluate orientational correlation of co-localized pseudo-vectors.

For example, the OOP has been used to quantify the organization of the bacterial population in a channel with respect to the channel direction [16]. Indeed, this is equivalent to a rudimentary case of the COOP where one of the constructs, the channel, is perfectly organized (Fig. 3C(ii)-(iii)). The COOP is more general in that it can be used when neither construct is perfectly organized. Circular statistics tool-sets include some correlation metrics [24], such as the circular correlation coefficient [21] which corresponds to the COOP in the same case. Specificically the circular correlation coefficient can only be used for uniform distributions (i.e. isotropic tissues). In that special case, the circular correlation coefficient and the COOP converge to the same equation (SI Supplemental Text). However, the more general vector formulation of the circular correlation coefficient, while not constrained to a uniform distribution, is very complex, and thus cannot be easily characterized the way we have done for the COOP. This circular correlation coefficient would not be a convenient metric for cytoskeleton or cellular orientation quantification. The COOP can be calculated so long as two sets of angular distributions and their locations are known, and it has been extensively characterized. Thus, this new parameter can be used with a multitude of biological systems.

In disease, loss of this organization leads to reduction in function, such as the myofibril organization changes in dilated cardiomyopathy [32—34]. However, proper organization of one element in a cell or tissue is not sufficient. The multiple constructs have to be properly organized with respect to each other, and that organization can have biological implications. This has been shown to occur during maturation of myocytes where the a-actinin is initially punctate and parallel to actin fibrils, but, in mature cardiomyocytes, becomes part of the newly formed sarcomeric Z-lines, which are perpendicular to actin fibrils [35]. In this case the relative orientation of a-actinin staining and actin fibrils indicates the progressive maturation of myofibrils. An additional example where orientation of different constructs affects each other is when the organization of the extracellular matrix can be used to control the architecture of cells [2, 13, 36]. Conversely, cells have been shown to change the orientation of the extracellular matrix fibrils [37]. Another instance of organization correlation can be found in the different cell types and collagen fibrils within heart valves [38]. The common use of such metrics as the OOP and COOP for biological and medical sciences will allow for a quantitative evaluation of tissue engineered substrates from a variety of cell sources. Combining such metrics with histology will create a universal evaluation metric between in vitro and in vivo systems favorably impacting our ability to design replacement tissues, to create in vitro drug testing platforms, and to evaluate pathological reports in the clinic.

Materials and Methods

Ethics Statement

This protocol met the guidelines for the use of vertebrate animals in research and teaching of the Faculty of Arts and Sciences of UCI. It also followed recommendations of the NIH Guide for the Care and Use of Laboratory Animals and was in accordance With existing federal (9 CFR Parts 1, 2 & 3), state, and city laws and regulations governing the use of animals in research and teaching.

Implementing COOP Calculation

To facilitate the calculation of the COOP we created a custom MATLAB code. The code was designed to have an input of angles for P and Q organized such that the information of which

Synthetic Data

Synthetic data of isotropic constructs for limiting conditions (Fig. 2) was generated using a random number generator (rand) that provides a uniform distribution of at least 106 random values in MATLAB. Each construct used in testing COOPu and COOPC (Fig. 3) contained 108 random numbers (MATLAB function normrnd) that were normally distributed With the specified mean and standard deviation. We have included the codes to create synthetic data as supporting codes (81 Code and 82 Code).

Experimental Data

To make the substrates 25 mm glass coverslips were coated with PDMS (Ellsworth Adhesives, Germantown, WI) and cured for 12 hours in a 60°C oven. To create triangular myocytes we utilized a microcontact printing procedure similar to that described by Tan et al [39]. A mask with the desired pattern was designed using Adobe Illustrator (Adobe Systems Incorporated, San Jose, CA) and made by Front Range Photomask (Palmer Lake, CO). The mask was used to make a silicone wafer (Integrated Nanosystems Research Facility, Irvine, CA). A polydimethylsiloxane (PDMS) stamp, cast from a silicon master, was used to contact transfer the extracellular matrix (ECM) protein fibronectin (FN) (Fisher Scientific Company, Hanover Park, IL) onto a UV-sterilized (UVO, Ielight Company, Inc. Irvine, CA) PDMS-coated coverslip. Fabricated substrates underwent one 10 minute pluronics (250g of Pluronics F- 127, Sigma-Aldrich, Inc., Saint Louis, MO) wash and three rinses of phosphate buffer-saline (PBS) (Life Technologies, Carlsbad, CA). To make isotropic substrates, UV-sterilized PDMS-coated coverslip were coated with FN for 10 minutes and underwent three PBS washes. The substrates were stored at 4°C prior to NRVM seeding.

Cell cultures of NRVMs were prepared from two-day old Sprague-Dawley rats (Charles River Laboratories, Wilmington, MA). A mid-sternal incision was made in order to expose the heart of the neonatal rat for dissection. Ventricular tissue was removed and rinsed in a Hanks balanced salt solution buffer (Life Technologies, Carlsbad, CA) and placed in 1 mg/mL trypsin solution (Sigma-Aldrich, Inc., Saint Louis, M0) to be shaken overnight (12 hour incubation) at 4°C. The next day, isolated tissue was dissociated into individual cells by treatment with four separate washes of 1 mg/mL collagenase type II (Worthington Biochemical, Lakewood, NI) for two minutes at 37°C. Isolated cardiomyocytes were resuspended in M199 culture medium (Invitrogen, Carlsbad, CA) supplemented with 10% heat-inactivated Fetal Bovine Serum, 10 mM HEPES, 20 mM glucose, 2 mM L-glutamine, 1.5 uM vitamin B-12 and 50 U/ml penicillin. The cell solution was filtered with a 40 um filter (Thermo Fisher Scientific, Waltham, MA), and the remaining cells were pre-plated multiple times to eliminate fibroblast contamination. Immediately after purification, myocytes were plated on substrates (prepared as detailed above) at a density of 106 or 105 cells per well in a standard six-well plate for confluent or sparse cultures, respectively. These were incubated at 37°C with a 5% C02 atmosphere. Seeded cultures underwent a wash with PBS 24 hours after plating to remove unattached and dead myocytes. They were then cultured in 10% serum media for another 24 hours at which point the media was changed to 2% serum media. After a total of 72 hours in culture, the samples were fixed and immunostained.

After each incubation period, cover-slips were rinsed three times with PBS for 5—10 min. Each coverslip was then mounted onto a microscope slide preserved with prolong gold antifade reagent (Life Technologies, Carlsbad, CA). The images were collected using an IX-83 inverted motorized microscope (Olympus America, Center Valley, PA) with an UPLFLN 40x oil immersion objective (Olympus America, Center Valley, PA) and a digital CCD camera ORCA-R2 C10600-10B (Hamamatsu Photonics, Shizuoka Prefecture, Japan). For isotropic monolayers, at least ten fields of view were collected for every sample.

A macro was created in Image] that allows the user to select regions without imperfections in an image. The regions that were not selected became masked, resulting in a series of images with only regions of interest displayed in the new masked images for every channel imaged (for example: DAPI, m-cherry, GFP)(Fig. 4). The masked images could then be analyzed using the same codes used for raw images.

To determine construct angles, we adapted a previous MATLAB code that detects ridges of a fingerprint [20, 25, 26]. This code was used to detect Z-lines and actin fibers in the images. In the code, a binary mask applied to the image determined the constructs and a filter was applied to clean up the constructs that were identified in the images. The code took pixel information from the images and for every nonempty pixel in the image, a pseudo-vector was calculated and used to determine the GDP for Z-lines and actin f1-brils, as well as a new set of pseudo-vectors for each square in the grid (Fig. 4B, D) These new pseudo-vectors were then utilized to calculate the COOP between two constructs (i.e. Z-lines and actin fibrils) or two cells (Fig. 4).

To calculate the average angle between the constructs ((60)) and the standard deviation of those angles (060) across multiple conver-slips, it is essential to keep in mind that the angle period is 71. The simplest way, but not the only way, to generate (60) and 060 is to calculate the director of the director pseudo-vectors fiALL of each coverslip. Meaning that in Equation (1) 2 f1 i where i is the coverslip, and the fiALL is the the eigenvector of the tensor from Equation (1) that corresponds to the eigenvalue from Equation (2). The angle for each coverslip is then determined as follows

Obviously if the COOP —>O then this procedure is useless as the angles will be inconsistent between the constructs. However, this procedure is a convenient way to determine the range in which it is most convenient to report the angle (i.e. 0—HT or — gag). Additionally, there are two possible angles that can be calculated (Fig. 1B). We chose to provide the clockwise angle from 13 to q, but it is also possible to calculate only the acute angle instead. To compare the COOP, COOPC, and COOPu in the analysis of the experimental data, the one way ANOVA with the Student-New-man-Keuls test was used.

Supporting Information

The supporting text provides the reader With more information about standard methods, such as the COP, circular statistics, and circular correlation. This information is provided in the same notation used for the manuscript. (PDF)

For (AD) schematic of the construct is on the left, and the orientation distribution With the GOP and standard deviation is on the right. (A) Perfect alignment; (B) almost perfect organization; (C) somewhat anisotropic; (D)

Statistical significance at p< 0.05, With OOP error of O'OOP =

For (AD) schematic of square grids (dashed black line outlines) on Z-lines(red). (AB) Grid size equivalent to one sarcomere complex length. (CD) Grid size equivalent to two sarcomere complex length. (A, C) Grid by chance aligns With Z-lines. (B, D) Grid does not align With Z-lines. we choose the grids shown in CD because we cannot control the alignment of grid to sarcomere complex.

Flow chart sketching the implementation of the neW method for experimental data. The additional steps for generating “noiseless” images With minimal imperfections is highlighted in blue. (TIF)

Image of the fibronectin island for each cell in Fig. 41 were cropped such that the triangular islands aligned With each other. (TIF)

Variables and parameters. A table of definitions for all variables and parameters. (PDF)

Code for synthetic data for four limiting cases. Code that creates synthetic data for four limiting cases.

Code for boundaries of normal region of the COOP. Code that creates synthetic data for boundaries of normal region of the COOP.

Author Contributions

Performed the experiments: NKD MAE DB. Analyzed the data: NKD AG. Contributed reagents/materials/analysis tools: NKD MAE AG. Wrote the paper: NKD MAE DB AG KKP.