Dynamic morphometric characterization of local connective tissue network structure in humans using ultrasound
© Langevin et al; licensee BioMed Central Ltd. 2007
Received: 21 November 2006
Accepted: 05 June 2007
Published: 05 June 2007
In humans, connective tissue forms a complex, interconnected network throughout the body that may have mechanosensory, regulatory and signaling functions. Understanding these potentially important phenomena requires non-invasive measurements of collagen network structure that can be performed in live animals or humans. The goal of this study was to show that ultrasound can be used to quantify dynamic changes in local connective tissue structure in vivo. We first performed combined ultrasound and histology examinations of the same tissue in two subjects undergoing surgery: in one subject, we examined the relationship of ultrasound to histological images in three dimensions; in the other, we examined the effect of a localized tissue perturbation using a previously developed robotic acupuncture needling technique. In ten additional non-surgical subjects, we quantified changes in tissue spatial organization over time during needle rotation vs. no rotation using ultrasound and semi-variogram analyses.
3-D renditions of ultrasound images showed longitudinal echogenic sheets that matched with collagenous sheets seen in histological preparations. Rank correlations between serial 2-D ultrasound and corresponding histology images resulted in high positive correlations for semi-variogram ranges computed parallel (r = 0.79, p < 0.001) and perpendicular (r = 0.63, p < 0.001) to the surface of the skin, indicating concordance in spatial structure between the two data sets. Needle rotation caused tissue displacement in the area surrounding the needle that was mapped spatially with ultrasound elastography and corresponded to collagen bundles winding around the needle on histological sections. In semi-variograms computed for each ultrasound frame, there was a greater change in the area under the semi-variogram curve across successive frames during needle rotation compared with no rotation. The direction of this change was heterogeneous across subjects. The frame-to-frame variability was 10-fold (p < 0.001) greater with rotation than with no rotation indicating changes in tissue structure during rotation.
The combination of ultrasound and semi-variogram analyses allows quantitative assessment of dynamic changes in the structure of human connective tissue in vivo.
Connective tissue, including the fasciae surrounding all muscles and organs, forms a continuous network throughout the body. This network also includes superficial subcutaneous connective tissue (pannicular fascia) which forms interconnecting longitudinal fibrous sheets intermixed with adipose tissue and linked to the skin by string-like "skin ligaments" [1–9]. Recent evidence suggests that this body-wide connective tissue network is not only composed of a web-like collagenous matrix, but also includes a network of fibroblasts that are dynamically responsive to mechanical stimulation [10–13]. These fibroblasts are linked to the collagen matrix by specialized proteins, as well as to each other by abundant cell-to-cell contacts [10, 14]. Modern understanding of complex systems suggests that, given these two types of interconnections, the dynamic responsiveness of fibroblasts within connective tissue may be affected by the structural characteristics of the collagenous network both locally and at a distance [15, 16]. Understanding these potentially important phenomena requires non-invasive measurements of both collagen network structure and cell function that can be performed in live animals or humans. Such measurements could be used to quantify changes over time and in response to specific experimental stimuli. Recently developed techniques using injectable fluorescent tagged proteins allow functional assessment of cells in live animals . Non-invasive in vivo measurement of collagen network microstructure has been previously achieved using polarization sensitive optical coherence tomography (PS-OCT) [18, 19]. However, this technique can only visualize structures within 1–2 millimeters of the skin surface, and is not well suited to the examination of deeper subcutaneous and perimuscular connective tissues. A technique allowing quantitative evaluation of both superficial and deep connective tissue network structure therefore may open up a new area of investigation in systems biology.
For many organ systems (e.g. cardiovascular), ultrasound is increasingly used to visualize human anatomy non-invasively and obtain in vivo tissue measurements [20–24]. In addition, elasticity imaging techniques use signal processing methods on consecutive ultrasound frames to detect tissue motion and strain, providing biomechanical as well as morphological information from ultrasound data . Elastography also has had an important impact on several clinical and research applications, such as the detection of differences in tissue stiffness associated with pathology in the breast [26, 27], prostate [28, 29], blood vessels , heart , skin and cartilage . During ultrasound imaging of biological materials, echoes generated by homogeneous material (e.g. fat) produce diffusely scattered signals, while echoes generated by interfaces of organized tissues with different acoustic impedances (product of density and speed of sound) produce more correlated "specular" signals [34–37]. It is generally assumed that the echogenic bands seen in human subcutaneous tissue correspond to connective tissue. This notion is supported by studies of connective tissue structures using both ultrasound and magnetic resonance imaging (MRI) [21, 38], as well as the successful use of echogenic bands as landmarks during ultrasound-guided percutaneous procedures such as regional anesthesia [39, 40]. However, to our knowledge, the correspondence between connective tissue structures revealed by ultrasound and "gold-standard" histological tissue identification has not been demonstrated quantitatively. Since human connective tissue contains a combination of adipose and fibrous (collagenous) tissues, the majority of specular echoes forming ultrasound images of connective tissue are expected to originate from interfaces between regions with predominantly fat (density = 0.917–0.940 g/cm3, speed of sound = 1,410–1,480 m/s), regions with predominantly collagen (density = 1.1–1.2 g/cm3, speed of sound = 1,540–1,570 m/s), and regions with predominantly glycosaminoglycans and water (density 0.998–1.05 g/cm3, speed of sound = 1480–1,500 m/s). Several parameters, including the ratio of specular-to-diffuse backscatter coefficients and the ultrasound intensity signal-to-noise ratio, have been used to assess tissue organization by quantifying the presence of organizational scatterers within tissues [35, 41–45]. These measures, however, do not reflect structural continuity in a given plane or direction, which is one of principal characteristics of connective tissue. Semi-variogram (also termed "variogram") analysis is a technique so far little applied to biology but commonly used in the hydrogeological and geophysical communities to evaluate the organization of laminar structures beneath the earth's surface [46–50]. Statistical analysis of ultrasound images using semi-variograms therefore offers a new way to quantify the organizational pattern of connective tissue in vivo, specifically regarding the structural continuity of longitudinal collagenous sheets, as well as changes in this structural organization over time.
The goals of this study were 1) examine the degree of correspondence between the three dimensional pattern of echoes generated by ultrasound imaging of connective tissue and the three dimensional pattern formed by fibrous connective tissue in histological serial sections of the same piece of tissue and 2) show that ultrasound can be used to quantify dynamic changes in connective tissue structure in vivo induced using a previously developed robotic acupuncture needling technique [51, 52]. We have previously demonstrated in animal experiments that rotation of an inserted ultra-fine acupuncture needle will cause winding and gathering of collagen fibers from the periphery toward the needle, and therefore can be used to achieve highly specific mechanical perturbations of fibrous connective tissue [52–54]. In this study, we have combined robotic needling with ultrasound elasticity imaging techniques, ultrasound semi-variogram analyses and histology to quantify structural changes occurring in connective tissue as a result of needle rotation. We used 1) elasticity imaging techniques to map areas of tissue motion, 2) histology to identify tissue structures in relation to ultrasound images and 3) semi-variograms to quantify changes in tissue spatial organization over time.
Results and discussion
Comparison of three dimensional ultrasound and histology image visualizations
Dynamic changes in connective tissue structure during needle rotation
Together, our results show that 1) needle rotation produced changes in connective tissue structure and 2) individuals responded differently from one another in response to the same mechanical input (i.e. needle rotation). This heterogeneous response of tissues observed among subjects was interesting and may be due to differences in initial tissue structure. Increased or decreased structural continuity over time may correspond to collagenous bundles either becoming more straight and tightly organized (i.e. less wavy) or "breaking up" due to the mechanical pull resulting from the tissue winding around the needle. Whether collagen bundles straighten or break in response to a given mechanical input may depend on a number of factors including collagen bundle size, cross-linking, tissue glycosaminoglycan concentration and water content [55–57]. Further studies will be needed to examine dynamic tissue behavior in relation to initial tissue factors and may lead to important new insights on how connective tissue behaves in the presence of pathology including injury, inflammation, scarring and fibrosis.
A potential limitation of structural tissue analysis based on ultrasound imaging is that the strength of reflected ultrasound signals is dependent to some degree on the orientation of a given structure relative to the direction of the incident wave. In this study, because the ultrasound transducer was placed perpendicular to the skin, organized structures oriented parallel to the skin would be expected to generate stronger echoes than structures perpendicular to the skin. Figure 2 indeed suggests that the small transverse septae within subcutaneous tissue are less well seen in the ultrasound images compared with the large longitudinal sheets. However, the close similarity between ultrasound and histology semi-variograms for both directions (parallel and perpendicular to the skin) indicates that ultrasound, nevertheless, captured the dominant structural characteristics of the tissue. Further studies examining other types of connective tissue should consider this limitation and control for ultrasound transducer position when comparing different experimental conditions.
The complementary techniques used in this study show that the laminar collagenous structures seen with ultrasound reorganized during rotation of an inserted needle. These changes were localized to the area of the needle and were accompanied by a "whorl" of collagen in corresponding biopsy sections. The combination of ultrasound and semi-variogram analysis used in this study is novel in biology. These imaging methods may be used along with dynamic visualization of molecular markers and computational modeling to investigate connective tissue as a functional network system capable of responding to local perturbations.
Twelve human subjects (11 females, 1 male; age range 24–74) were recruited and signed informed consent for this study. The Institutional Review Board at the University of Vermont approved the study methods and procedures. Two subjects were patients undergoing spinal fusion surgery (one 74 year old female and one 47 year old female). In the first surgical subject, we examined the degree of correspondence between serial ultrasound and histology images of the same tissue sample. In the second surgical subject, we examined structural changes occurring during robotic needle rotation using both ultrasound and histology, (ultrasound data was acquired with the transducer in a fixed position while an inserted needle was rotated, followed by tissue excision with the needle still in place). In the ten non-surgical subjects, we performed robotic acupuncture needling at the same anatomical location on the thigh bilaterally to compare the effect of needle rotation vs. needle insertion alone without rotation.
Comparison of ultrasound and histology three dimensional image visualizations (N = 1 subject)
With the patient anesthetized, a 4 cm × 4 cm skin surface area located at the bone graft site, (centered on an area 4 cm lateral to the posterior superior iliac spine and 4 cm caudal to the iliac crest) was imaged by ultrasound (Figure 1A, B). An ultrasound cine-recording was acquired using a GE System Five (Vingmed) scanner equipped with a linear array ultrasound transducer (frequency 10 MHz). The transducer was translated across the defined area at a rate of 2.5 mm/second while acquiring ultrasound images at 10 frames/second with an ultrasound image depth of 4 cm, yielding serial 3-d data on a 64 cm3 cube of tissue with pixel intensity ranging from 0–255. After ultrasound imaging, a 4 cm × 2 cm × 4 cm tissue sample was excised at the site of the surgical incision, centered on the middle of the area scanned by ultrasound (Figure 1A, C). The tissue was cauterized at the rostral end for orientation purposes, and then fixed in 10% formalin. After fixation the tissue sample was cut transversely into seven 4 mm thick blocks (Figure 1C, D). The seven tissue blocks were paraffin-embedded and serial-sectioned at 10 μm thickness. Every 25th section was stained with hematoxylin/eosin (H&E). Sections were photographed using an Olympus Bx51 reflected light microscope (Melville, NY). The H&E images were digitized and converted to gray scale. A smaller number of additional serial sections also were stained with Masson Trichrome, and compared with the corresponding H&E images to verify the accuracy of H&E staining in defining the outline of the collagenous subcutaneous connective tissue structures (Figure 1E, F). In order to assess changes in tissue dimensions during processing, we measured the dimensions of 1) the freshly excised whole tissue biopsy, 2) the whole fixed biopsy, 3) each trimmed tissue block before paraffin embedding 4) paraffin embedded tissue blocks and 5) tissue sections on histology slides. We found that minimal tissue shrinkage occurred during this overall process (on average, tissue measurements increased by 6.1% after paraformaldehyde fixation and decreased by 6.8% after paraffin embedding). We therefore did not include a shrinkage factor in the 3-d rendition of histological images.
The ultrasound cine-recording included 160 individual 2-D frames. Each ultrasound frame was cropped to match the area of the histology image. Histology slides were stacked sequentially using major landmarks of the tissue (e.g. blood vessels) for orientation. However, tissue landmarks were not used to match ultrasound and histology images with each other. Rather, ultrasound and histology data sets were matched as follows: the histology section corresponding to the middle of the biopsy (middle histology section) was first identified based on serial section numbers; this histology section was then matched with the ultrasound image corresponding to the middle of the ultrasound recording (80th frame); next, x coordinates were matched by lining up the medial-lateral midpoint of the middle histology section with the medial-lateral midpoint of the ultrasound image; finally, y coordinates were matched by lining up the deep border of the dermis in the middle histology section to the deep border of the dermis in the middle ultrasound image. Thus, histology sections were stacked relative to each other before matching the whole histology data set with the serial ultrasound images. This avoided the potential bias associated with using connective tissue layers themselves as landmarks. Three-dimensional visualizations of ultrasound and histology image stacks were performed using Environmental Visualization EVS-Pro V (8.0) software.
Effect of needle rotation using ultrasound and histology (N = 1 subject)
With the patient anesthetized, an ultrasound cine-recording was acquired at the same location as above, with the transducer in a fixed position while an inserted needle was rotated, followed by tissue excision with the needle still in place. Needle insertion and rotation were performed using a computer-controlled needling instrument as previously described [51, 52] and disposable stainless steel acupuncture needles (Seirin, Japan, 0.25 mm in diameter and 50 mm in length). The needling instrument and ultrasound transducer were placed in a clamp system such that the transducer was perpendicular to the skin and the needling instrument was at a 20° angle with respect to the transducer. Needle insertion depth was 40 mm. Needle rotation and axial translational speeds were 6 revolutions/sec and 10 mm/s respectively. The needle was rotated for 4 seconds for a total of 24 continuous uni-directional revolutions.
An ultrasound recording consisting of seventy successive frames was acquired continuously during the needling procedure at 13.2 frames/s. Ultrasound frequency was 6.9 MHz. Raw ultrasonic radio frequency (RF) data were collected and stored following digitization by the ultrasound system at a sampling frequency of 20 MHz.
After ultrasound data acquisition, the needling instrument was disconnected from the needle, leaving the needle in place within the tissue. A 4 cm × 2 cm × 4 cm tissue sample was then excised as above, centered on the needle. After fixation, the position of the needle in the tissue was marked by dipping the tip of the needle (protruding from the end of the tissue) into green ink before pulling it out of the tissue, thus dragging the ink through the tissue as the needle was being pulled back. The tissue sample was then cut into five blocks perpendicular to the needle as shown in Figure 4A, D, E, F and processed for histology.
Comparison of needle rotation vs. needle insertion without rotation (N = 10 subjects)
Ultrasound cine-recordings were acquired at 13.2 frames/sec as described above on the anterior thigh (18 cm superior to the middle of the superior edge of the patella) bilaterally. The two testing locations were randomized to needle rotation (8 clockwise uni-directional revolutions) vs. no rotation. Needle insertion depth was calculated for each individual so that the needle tip extended 5 mm past the peri-muscular fascia, and the needle insertion angle was 45° to the transducer. Both raw radio-frequency (RF) ultrasound data and B-scan image were saved for quantification of tissue displacement (elasticity imaging) and changes in tissue structure (semi-variogram analysis).
Ultrasound data post-processing for measurement of tissue displacement during needle rotation (using elasticity imaging techniques)
The displacement occurring between each of the seventy successive frames acquired during needle rotation was calculated using cross-correlation techniques with a 3 mm window and a window overlap of 90% as previously described [25, 52]. The term 'displacement' refers to the incremental axial (i.e. along the propagation of the ultrasound beam, or top to bottom, on all images shown) motion that the tissue undergoes between two successively acquired ultrasound frames (i.e. after 76 ms have elapsed).
where N(h) is the number of pixel pairs separated by the distance, h; and u(a) and u(a+h) are the respective pixel intensities at locations (a) and some distance away (a+h). The difference in pixel intensity between each data point and every other data point in a given direction was calculated and these differences were binned according to separation distance (h) (see horizontal axes in Figure 3). For example, all pairs of points separated by a distance ranging between 0 and 10 pixels were included in the first bin and represented by the first data point of Figure 3; the y axis represents the average semi-variance in pixel intensity associated with these pairs of data. Semi-variance values were calculated and averaged over each bin using MATLAB version 6.1 (MathWorks, Inc, Natick, MA).
In this study, we used semi-variograms to 1) compare B-scan ultrasound images with histology based on quantitative measures representing tissue structure and 2) quantify changes in tissue structure over time during needle rotation. For the ultrasound/histology comparison, we used a standard semi-variogram parameter (the range) derived using non-linear regression models fit using JMP statistical software (SAS Institute Inc., Cary, NC) [58–60]. The range is the most commonly used measure of spatial structure for spatially auto-correlated data sets and is defined as the distance at which γ(h) (the variability in the difference in pixel intensity between pairs of points) reaches a plateau. Larger ranges correspond to less spatial variability and therefore greater spatial correlation (i.e. points located large distances apart are still correlated) while smaller ranges imply more spatial variation and therefore weaker spatial correlation. Spearmans's rank correlation was used to examine the concordance between estimates obtained from ultrasound and histology.
The area under the semi-variogram curve (AUC) was used to examine structural changes in tissue during needle rotation. The AUC was calculated using the trapezoidal rule from the raw semi-variogram data corresponding to each B-scan ultrasound image acquired during needle rotation (or no rotation). We used the AUC as a surrogate measure for the range as it provides a quantitative measure of spatial structure without necessitating non-linear curve fitting for each of the 300 ultrasound frames (15 frames × 2conditions × 10 subjects) in parallel and perpendicular directions. The AUC is sensitive to change in data structure, and was found to be highly correlated with the range in the ultrasound/histology comparison (r = 0.88, P < 0.001 and r = 0.90, P < 0.001 for directions parallel and perpendicular to the skin respectively).
Repeated measures analysis of variance was used to compare mean AUC between rotation and no-rotation conditions. The mean AUC for each condition was based on the average of the 15 frames across the 10 subjects. Analysis of variance was also used to compare rotation and no-rotation conditions with respect to the mean absolute difference in AUC between successive frames. This outcome measure examines structural changes while accounting for the temporal ordering of the frames. Variance component analysis (SAS, PROC VARCOMP) was used the obtain estimates of the inter-frame variance in AUC (VAUC) for rotation and no-rotation conditions. Variance estimates for the two conditions were compared based on an F-test. Statistical significance was determined based on α = 0.05.
The authors thank Brendan H. Kennedy and Alexander T. Danco for technical assistance. This study was funded by the National Institutes of Health Center for Complementary and Alternative Medicine Research Grant RO1 AT01121. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health.
- Abu-Hijleh MF, Roshier AL, Al-Shboul Q, Dharap AS, Harris PF: The membranous layer of superficial fascia: evidence for its widespread distribution in the body. Surg Radiol Anat. 2006, 28 (6): 606-619. 10.1007/s00276-006-0142-8.PubMedView ArticleGoogle Scholar
- Singer E: Fasciae of the human body and their relations to the organs they envelop. Fasciae of the human body and their relations to the organs they envelop. 1935, : 1-34. Baltimore, The Williams & Wilkins CompanyGoogle Scholar
- Tobin CE, Benjamin JA: Anatomic and clinical re-evaluation of Camper's, Scarpa's, and Colles' fasciae. Surgery Today. 1949, 88: 545-559.Google Scholar
- Congdon ED, Edson J, Yanitelli S: Gross structure of the subcutaneous layer of the anterior and lateral trunk in the male. Am J Anat. 1946, 79: 399-429. 10.1002/aja.1000790304.PubMedView ArticleGoogle Scholar
- Nash LG, Phillips MN, Nicholson H, Barnett R, Zhang M: Skin ligaments: regional distribution and variation in morphology. Clin Anat. 2004, 17: 287-293. 10.1002/ca.10203.PubMedView ArticleGoogle Scholar
- Ghassemi A, Prescher A, Riediger D, Axer H: Anatomy of the SMAS revisited. Aesthetic Plast Surg. 2003, 27: 258-264. 10.1007/s00266-003-3065-3.PubMedView ArticleGoogle Scholar
- Lockwood TE: Superficial fascial system (SFS) of the trunk and extremities: a new concept. Plast Reconstr Surg. 1991, 87 (6): 1009-1018.PubMedView ArticleGoogle Scholar
- Markman B: Anatomy and physiology of adipose tissue. Clin Plast Surg. 1989, 16: 235-244.PubMedGoogle Scholar
- Markman B, Barton FE: Anatomy of the subcutaneous tissue of the trunk and lower extremity. Plast Reconstr Surg. 1987, 80: 248-254. 10.1097/00006534-198708000-00015.PubMedView ArticleGoogle Scholar
- Langevin HM, Cornbrooks CJ, Taatjes DJ: Fibroblasts form a body-wide cellular network. Histochem Cell Biol. 2004, 122: 7-15. 10.1007/s00418-004-0667-z.PubMedView ArticleGoogle Scholar
- Langevin HM, Bouffard NA, Badger GJ, Iatridis JC, Howe AK: Dynamic fibroblast cytoskeletal response to subcutaneous tissue stretch ex vivo and in vivo. Am J Physiol Cell Physiol. 2005, 288: C747-56. 10.1152/ajpcell.00420.2004.PubMedView ArticleGoogle Scholar
- Langevin HM, Bouffard NA, Badger GJ, Churchill DL, Howe AK: Subcutaneous tissue fibroblast cytoskeletal remodeling induced by acupuncture: Evidence for a mechanotransduction-based mechanism. J Cell Physiol. 2006, 207: 767-774. 10.1002/jcp.20623.PubMedView ArticleGoogle Scholar
- Langevin HM, Storch KN, Cipolla MJ, White SL, Buttolph TR, Taatjes DJ: Fibroblast spreading induced by connective tissue stretch involves intracellular redistribution of alpha- and beta-actin. Histochem Cell Biol. 2006, 125: 487-495. 10.1007/s00418-005-0138-1.PubMedView ArticleGoogle Scholar
- Giancotti FG, Ruoslahti E: Integrin signaling. Science. 1999, 285: 1028-1032. 10.1126/science.285.5430.1028.PubMedView ArticleGoogle Scholar
- Langevin HM: Connective tissue: a body-wide signaling network?. Med Hypotheses. 2006, 66: 1074-1077. 10.1016/j.mehy.2005.12.032.PubMedView ArticleGoogle Scholar
- Bar-Yam Y: Dynamics of complex systems. Studies in nonlinearity. 2003, xvi, 848 p.-Boulder, CO, Westview PressGoogle Scholar
- Prescher JA, Bertozzi CR: Chemistry in living systems. Nat Chem Biol. 2005, 1: 13-21. 10.1038/nchembio0605-13.PubMedView ArticleGoogle Scholar
- Liu B, Harman M, Giattina S, Stamper DL, Demakis C, Chilek M, Raby S, Brezinski ME: Characterizing of tissue microstructure with single-detector polarization-sensitive optical coherence tomography. Appl Opt. 2006, 45: 4464-4479. 10.1364/AO.45.004464.PubMedView ArticleGoogle Scholar
- Pierce MC, Strasswimmer J, Hyle Park B, Cense B, De Boer JF: Birefringence measurements in human skin using polarization-sensitive optical coherence tomography. J Biomed Opt. 2004, 9: 287-291. 10.1117/1.1645797.PubMedView ArticleGoogle Scholar
- Martin DC, Medri MK, Chow RS, Oxorn V, Leekam RN, Agur AM, McKee NH: Comparing human skeletal muscle architectural parameters of cadavers with in vivo ultrasonographic measurements. J Anat. 2001, 199: 429-434. 10.1046/j.1469-7580.2001.19940429.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Theobald P, Bydder G, Dent C, Nokes L, Pugh N, Benjamin M: The functional anatomy of Kager's fat pad in relation to retrocalcaneal problems and other hindfoot disorders. J Anat. 2006, 208: 91-97. 10.1111/j.1469-7580.2006.00510.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Maganaris CN, Kawakami Y, Fukunaga T: Changes in aponeurotic dimensions upon muscle shortening: in vivo observations in man. J Anat. 2001, 199: 449-456. 10.1046/j.1469-7580.2001.19940449.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Mellor RH, Bush NL, Stanton AW, Bamber JC, Levick JR, Mortimer PS: Dual-frequency ultrasound examination of skin and subcutis thickness in breast cancer-related lymphedema. Breast J. 2004, 10: 496-503. 10.1111/j.1075-122X.2004.21458.x.PubMedView ArticleGoogle Scholar
- Turnbull DH, Starkoski BG, Harasiewicz KA, Semple JL, From L, Gupta AK, Sauder DN, Foster FS: A 40-100 MHz B-scan ultrasound backscatter microscope for skin imaging. Ultrasound Med Biol. 1995, 21: 79-88. 10.1016/0301-5629(94)00083-2.PubMedView ArticleGoogle Scholar
- Ophir J, Cespedes I, Ponnekanti H, Yazdi Y, Li X: Elastography: a quantitative method for imaging the elasticity of biological tissues. Ultrason Imaging. 1991, 13: 111-134. 10.1016/0161-7346(91)90079-W.PubMedView ArticleGoogle Scholar
- Garra BS, Cespedes EI, Ophir J, Spratt SR, Zuurbier RA, Magnant CM, Pennanen MF: Elastography of breast lesions: initial clinical results. Radiology. 1997, 202: 79-86.PubMedView ArticleGoogle Scholar
- Konofagou EE, Harrigan T, Ophir J: Shear strain estimation and lesion mobility assessment in elastography. Ultrasonics. 2000, 38: 400-404. 10.1016/S0041-624X(99)00125-0.PubMedView ArticleGoogle Scholar
- Kallel F, Price RE, Konofagou E, Ophir J: Elastographic imaging of the normal canine prostate in vitro. Ultrason Imaging. 1999, 21: 201-215.PubMedView ArticleGoogle Scholar
- Souchon R, Rouviere O, Gelet A, Detti V, Srinivasan S, Ophir J, Chapelon JY: Visualisation of HIFU lesions using elastography of the human prostate in vivo: preliminary results. Ultrasound Med Biol. 2003, 29: 1007-1015. 10.1016/S0301-5629(03)00065-6.PubMedView ArticleGoogle Scholar
- de Korte CL, Ignacio Cespedes EI, van der Steen AF, Lancee CT: Intravascular elasticity imaging using ultrasound: feasibility studies in phantoms. Ultrasound Med Biol. 1997, 23: 735-746. 10.1016/S0301-5629(97)00004-5.PubMedView ArticleGoogle Scholar
- Konofagou EE, D'Hooge J, Ophir J: Myocardial elastography--a feasibility study in vivo. Ultrasound Med Biol. 2002, 28: 475-482. 10.1016/S0301-5629(02)00488-X.PubMedView ArticleGoogle Scholar
- Gennisson JL, Baldeweck T, Tanter M, Catheline S, Fink M, Sandrin L, Cornillon C, Querleux B: Assessment of elastic parameters of human skin using dynamic elastography. IEEE Trans Ultrason Ferroelectr Freq Control. 2004, 51: 980-989. 10.1109/TUFFC.2004.1324402.PubMedView ArticleGoogle Scholar
- Konofagou EE, Harrigan TP, Ophir J, Krouskop TA: Poroelastography: imaging the poroelastic properties of tissues. Ultrasound Med Biol. 2001, 27: 1387-1397. 10.1016/S0301-5629(01)00433-1.PubMedView ArticleGoogle Scholar
- Kremkau FW: Diagnostic ultrasound : principles and instruments. 1998, xi, 590 p.-Philadelphia, Pa., W.B. Saunders, 5thGoogle Scholar
- Garra BS: In vivo liver and splenic tissue characterization by scattering. Ultrasonic scattering in biological tissues. Edited by: Shung KK and Thieme GA. 1993, 347-391. Boca Raton, CRC PressGoogle Scholar
- Insana MF, Wagner RF, Garra BS, Brown DG, Shawker TH: Analysis of ultrasound image texture via generalized Rician statistics. Proceedings of the International Society for Optical Engineering. 1985, 556: 153-159.Google Scholar
- Lizzi FL, Alam SK, Mikaelian S, Lee P, Feleppa EJ: On the statistics of ultrasonic spectral parameters. Ultrasound Med Biol. 2006, 32: 1671-1685. 10.1016/j.ultrasmedbio.2006.09.002.PubMedView ArticleGoogle Scholar
- Gregg JM, Silberstein M, Schneider T, Kerr JB, Marks P: Sonography of plantar plates in cadavers: correlation with MRI and histology. AJR Am J Roentgenol. 2006, 186: 948-955. 10.2214/AJR.04.1481.PubMedView ArticleGoogle Scholar
- Popovic J, Morimoto M, Wambold D, Blanck TJ, Rosenberg AD: Current practice of ultrasound-assisted regional anesthesia. Pain Pract. 2006, 6: 127-134.PubMedView ArticleGoogle Scholar
- Marhofer P, Greher M, Kapral S: Ultrasound guidance in regional anaesthesia. Br J Anaesth. 2005, 94: 7-17. 10.1093/bja/aei002.PubMedView ArticleGoogle Scholar
- Wagner RF, Insana MF, Brown DG: Unified approach to the detection and classification of speckle texture in diagnostic ultrasound. Proceedings of SPIE-The International Society for Optical Engineering. 1985, 556: 146-152.Google Scholar
- Hall CS, Scott MJ, Lanza GM, Miller JG, Wickline SA: The extracellular matrix is an important source of ultrasound backscatter from myocardium. J Acoust Soc Am. 2000, 107: 612-619. 10.1121/1.428327.PubMedView ArticleGoogle Scholar
- Wilson T, Chen Q, Zagzebski JA, Varghese T, VanMiddlesworth L: Initial clinical experience imaging scatterer size and strain in thyroid nodules. J Ultrasound Med. 2006, 25: 1021-1029.PubMedGoogle Scholar
- Bridal SL, Fournier C, Coron A, Leguerney I, Laugier P: Ultrasonic backscatter and attenuation (11-27 MHz) variation with collagen fiber distribution in ex vivo human dermis. Ultrason Imaging. 2006, 28: 23-40.PubMedView ArticleGoogle Scholar
- Gaitini D, Lederman M, Baruch Y, Ghersin E, Veitsman E, Kerner H, Shalem B, Yaniv G, Sarfaty C, Azhari H: Computerised analysis of liver texture with correlation to needle biopsy. Ultraschall Med. 2005, 26: 197-202. 10.1055/s-2005-858267.PubMedView ArticleGoogle Scholar
- Journel AG, Huijbregts CJ: Mining geostatistics. 1978, x, 600 p.-London ; New York, Academic PressGoogle Scholar
- Diaz G, Zucca A, Setzu MD, Cappai C: Chromatin pattern by variogram analysis. Microsc Res Tech. 1997, 39: 305-311. 10.1002/(SICI)1097-0029(19971101)39:3<305::AID-JEMT9>3.0.CO;2-I.PubMedView ArticleGoogle Scholar
- Gough JJ, Kent JT, O'Higgins P, Ellison LT: Variogram methods for the analysis of bony trabecular shadows in plain radiographs. Int J Biomed Comput. 1994, 35: 141-153. 10.1016/0020-7101(94)90063-9.PubMedView ArticleGoogle Scholar
- McShane LM, Meier KL, Wassermann EM: A comparison of spatial prediction techniques for an exploratory analysis of human cortical motor representations. Stat Med. 1997, 16: 1337-1355. 10.1002/(SICI)1097-0258(19970630)16:12<1337::AID-SIM564>3.0.CO;2-W.PubMedView ArticleGoogle Scholar
- Rossi RE, Mulla DJ, Journel AG, Franz EH: Geostatistical tools for modeling and interpreting ecological spatial dependence. Ecological Monographs. 1992, 62: 277-314. 10.2307/2937096.View ArticleGoogle Scholar
- Langevin HM, Churchill DL, Fox JR, Badger GJ, Garra BS, Krag MH: Biomechanical response to acupuncture needling in humans. J Appl Physiol. 2001, 91: 2471-2478.PubMedGoogle Scholar
- Langevin HM, Konofagou EE, Badger GJ, Churchill DL, Fox JR, Ophir J, Garra BS: Tissue displacements during acupuncture using ultrasound elastography techniques. Ultrasound Med Biol. 2004, 30: 1173-1183. 10.1016/j.ultrasmedbio.2004.07.010.PubMedView ArticleGoogle Scholar
- Langevin HM, Churchill DL, Cipolla MJ: Mechanical signaling through connective tissue: a mechanism for the therapeutic effect of acupuncture. FASEB J. 2001, 15: 2275-2282. 10.1096/fj.01-0015hyp.PubMedView ArticleGoogle Scholar
- Langevin HM, Churchill DL, Wu J, Badger GJ, Yandow JA, Fox JR, Krag MH: Evidence of connective tissue involvement in acupuncture. Faseb J. 2002, 16: 872-874.PubMedGoogle Scholar
- Armstrong CG, Mow VC: Variations in the intrinsic mechanical properties of human articular cartilage with age, degeneration, and water content. J Bone Joint Surg Am. 1982, 64: 88-94.PubMedGoogle Scholar
- Lyons G, Eisenstein SM, Sweet MB: Biochemical changes in intervertebral disc degeneration. Biochim Biophys Acta. 1981, 673: 443-453.PubMedView ArticleGoogle Scholar
- Woo SLY, Livesay GA, Runco TJ, Young EP: Structure and function of tendons and ligaments. Basic orthopaedic biomechanics. Edited by: Mow VC and Hayes WC. 1997, 209-251. Philadelphia, Lippincott-Raven, 2ndGoogle Scholar
- Marsily G: Quantitative hydrogeology : groundwater hydrology for engineers. 1986, xix, 440 p.-Orlando, FL, Academic PressGoogle Scholar
- Goovaerts P: Geostatistical tools for characterizing the spatial variability of microbiological and physico-chemical soil properties. Biology and Fertility of Soils. 1998, 27: 315-334. 10.1007/s003740050439.View ArticleGoogle Scholar
- Isaaks EH, Srivastava RM: Applied geostatistics. 1989, xix, 561 p.-New York, Oxford University PressGoogle Scholar