Tunable In Situ Synthesis of Ultrathin Extracellular Matrix-Derived Membranes in Organ-on-a-Chip Devices

Thin cell culture membranes in organ-on-a-chip (OOC) devices are used to model a wide range of thin tissues. While early and most current platforms use microporous or fibrous elastomeric or thermoplastic membranes, there is an emerging class of devices using extra-cellular matrix (ECM) protein-based membranes to improve their biological relevance. These ECM-based membranes present physiologically relevant properties, but they are difficult to integrate into OOC devices due to their relative fragility. Additionally, the specialized fabrication methods developed to date make comparison between methods difficult. This work presents the development and characterization of a method to produce ultrathin matrix-derived membranes (UMM) in OOC devices that requires only a preassembled thermoplastic device and a micropipette, decoupling the device and UMM fabrication processes. Control over the thickness and permeability of the UMM is demonstrated, along with integration of the UMM in a device enabling high-resolution on-chip microscopy. The reliability of the UMM fabrication method is leveraged to develop a medium-throughput well-plate format device with 32 independent UMM-integrated samples. Finally, proof-of-concept cell culture experiments are demonstrated. Due to its simplicity and controllability, the presented method has the potential to overcome technical barriers preventing wider adoption of physiologically relevant ECM-based membranes in OOC devices.

Quetiapine-Related Deaths: In Search of a Surrogate Endpoint

Quetiapine is a second-generation antipsychotic drug available for two and half decades. Due to increased misuse, prescription outside the approved indications, and availability on the black market, it is being encountered in medicolegal autopsies more frequently. For instance, it has been linked to increased mortality rates, most likely due to its adverse effects on the cardiovascular system. Its pharmacokinetic features and significant postmortem redistribution challenge traditional sampling in forensic toxicology. Therefore, a systematic literature review was performed, inclusive of PubMed, the Web of Science-core collection, and the Scopus databases; articles were screened for the terms "quetiapine", "death", and "autopsy" to reevaluate each matrix used as a surrogate endpoint in the forensic toxicology of quetiapine-related deaths. Ultimately, this review considers the results of five studies that were well presented (more than two matrices, data available for all analyses, for instance). The highest quetiapine concentrations were usually measured in the liver tissue. As interpreted by their authors, the results of the considered studies showed a strong correlation between some matrices, but, unfortunately, the studies presented models with poor goodness of fit. The distribution of quetiapine in distinct body compartments/tissues showed no statistically significant relationship with the length of the postmortem interval. Furthermore, this study did not confirm the anecdotal correlation of peripheral blood concentrations with skeletal muscle concentrations. Otherwise, there was no consistency regarding selecting an endpoint for analysis.

Piezoresistivity modeling of soft tissue electrical-mechanical properties: A validation study

In general, the electrical property of soft tissues is sensitive to the force applied to their surface. To further study the relationship between the force and the electrical property of soft tissues, this paper attempts to investigate the effect of static and higher-order stresses on electrical properties. Overall, a practical experimental platform is designed to acquire the force information and the electrical property of soft tissues during a contact procedure, which is featured different compression stimuli, such as constant pressing force, constant pressing speed, and step-force compression, etc. Furthermore, the piezoresistive characteristic is innovatively introduced to model the mechanical-electrical properties of soft tissue. Finite Element Modeling (FEM) is adopted to fit the static piezoresistivity of the soft tissue. Finally, experimental studies were performed to demonstrate the effect of stress on the electrical properties and the feasibility of the proposed piezoresistive model to describe soft tissues' mechanical and electrical properties.

Development and CT image-domain validation of a computational lung lesion model for use in virtual imaging trials

PURPOSE

Computational abnormalities (e.g., lesion models) for use in medical imaging simulation studies are frequently generated using data collected from clinical images. Although this approach allows for highly-customizable lesion detectability studies on clinical computed tomography (CT) data, the ground-truth lesion models produced with this method do not provide a sufficiently realistic lesion morphology for use with current anthropomorphic simulation studies. This work is intended to demonstrate that the new anatomically-informed lesion model presented here is not inferior to the previous lesion model under CT imaging, and can therefore provide a more biologically-informed model for use with simulated CT imaging studies.

METHODS

The lesion model was simulated initially from a seed cell with 10 μm diameter placed in an anatomical location within segmented lung CT and was allowed to reproduce locally within the available solid angle in discrete time-intervals (corresponding to synchronous cell cycles) up to a size of ∼200 μm in diameter. Daughter cells of generation G were allowed also to reproduce on the next available time-step given sufficient space. At lesion sizes beyond 200 μm in diameter, the health of subregions of cells were tracked with a Markov chain technique, indicating which regions were likely to continue growing, which were likely stable, and which were likely to develop necrosis given their proximity to anatomical features and other lesion cells. For lesion sizes beyond 500 μm, the lesion was represented with three nested, triangulated surfaces (corresponding to proliferating, dormant, and necrotic regions), indicating how discrete volumes of the lesion were behaving at a particular time. Lesions were then assigned smoothly-varying material properties based on their cellular level health in each region, resulting in a multi-material lesion model. The lesions produced with this model were then voxelized and placed into lung CT images for comparison with both prior work and clinical data. This model was subject to an observer study in which cardiothoracic imaging radiologists assessed the realism of both clinical and synthetic lesions in CT images.

RESULTS

The useable outputs of this work were voxel- or surface-based, validated, computational lesions, at a scale clearly visible on clinical CT (3-4 cm). Analysis of the observer study results indicated that the computationally-generated lesions were indistinguishable from clinical lesions (AUC = 0.49, 95% CI = [0.36, 0.61]) and non-inferior to an earlier image-based lesion model-indicating the advantage of the model for use in both hybrid CT images and in simulated CT imaging of the lungs.

CONCLUSIONS

Results indicated the non-inferiority of this model as compared to previous methods, indicating the utility of the model for use in both hybrid CT images and in simulated CT imaging.

Customizable Microfluidic Origami Liver-on-a-Chip (oLOC)

The design and manufacture of an origami-based liver-on-a-chip device are presented, together with demonstrations of the chip's effectiveness at recapitulating some of the liver's key in vivo architecture, physical microenvironment, and functions. Laser-cut layers of polyimide tape are folded together with polycarbonate nanoporous membranes to create a stack of three adjacent flow chambers separated by the membranes. Endothelial cells are seeded in the upper and lower flow chambers to simulate sinusoids, and hepatocytes are seeded in the middle flow chamber. Nutrients and metabolites flow through the simulated sinusoids and diffuse between the vascular pathways and the hepatocyte layers, mimicking physiological microcirculation. Studies of cell viability, metabolic functions, and hepatotoxicity of pharmaceutical compounds show that the endothelialized liver-on-a-chip model is conducive to maintaining hepatocyte functions and evaluation of the hepatotoxicity of drugs. Our unique origami approach speeds chip development and optimization, effectively simplifying the laboratory-scale fabrication of on-chip models of human tissues without necessarily reducing their structural and functional sophistication.

Anatomically- and physiologically-informed computational model of hepatic contrast perfusion for virtual imaging trials

PURPOSE

Virtual (in silico) imaging trials (VITs), involving computerized phantoms and models of the imaging process, provide a modern alternative to clinical imaging trials. VITs are faster, safer, and enable otherwise-impossible investigations. Current phantoms used in VITs are limited in their ability to model functional behavior such as contrast perfusion which is an important determinant of dose and image quality in CT imaging. In our prior work with the XCAT computational phantoms, we determined and modeled inter-organ (organ to organ) intravenous contrast concentration as a function of time from injection. However, intra-organ concentration, heterogeneous distribution within a given organ, was not pursued. We extend our methods in this work to model intra-organ concentration within the XCAT phantom with a specific focus on the liver.

METHODS

Intra-organ contrast perfusion depends on the organ's vessel network. We modeled the intricate vascular structures of the liver, informed by empirical and theoretical observations of anatomy and physiology. The developed vessel generation algorithm modeled a dual-input-single-output vascular network as a series of bifurcating surfaces to optimally deliver flow within the bounding surface of a given XCAT liver. Using this network, contrast perfusion was simulated within voxelized versions of the phantom by using knowledge of the blood velocities in each vascular structure, vessel diameters and length, and the time since the contrast entered the hepatic artery. The utility of the enhanced phantom was demonstrated through a simulation study with the phantom voxelized prior to CT simulation with the relevant liver vasculature prepared to represent blood and iodinated contrast media. The spatial extent of the blood-contrast mixture was compared to clinical data.

RESULTS

The vascular structures of the liver were generated with size and orientation which resulted in minimal energy expenditure required to maintain blood flow. Intravenous contrast was simulated as having known concentration and known total volume in the liver as calibrated from time-concentration curves (TCC). Measurements of simulated CT ROIs were found to agree with clinically-observed values of early arterial phase contrast enhancement of the parenchyma (∼5 HU). Similarly, early enhancement in the hepatic artery was found to agree with average clinical enhancement (180 HU).

CONCLUSIONS

The computational methods presented here furthered the development of the XCAT phantoms allowing for multi-timepoint contrast perfusion simulations, enabling more anthropomorphic virtual clinical trials intended for optimization of current clinical imaging technologies and applications. This article is protected by copyright. All rights reserved.

Quantification of Myocyte Disarray in Human Cardiac Tissue

Proper three-dimensional (3D)-cardiomyocyte orientation is important for an effective tension production in cardiac muscle. Cardiac diseases can cause severe remodeling processes in the heart, such as cellular misalignment, that can affect both the electrical and mechanical functions of the organ. To date, a proven methodology to map and quantify myocytes disarray in massive samples is missing. In this study, we present an experimental pipeline to reconstruct and analyze the 3D cardiomyocyte architecture in massive samples. We employed tissue clearing, staining, and advanced microscopy techniques to detect sarcomeres in relatively large human myocardial strips with micrometric resolution. Z-bands periodicity was exploited in a frequency analysis approach to extract the 3D myofilament orientation, providing an orientation map used to characterize the tissue organization at different spatial scales. As a proof-of-principle, we applied the proposed method to healthy and pathologically remodeled human cardiac tissue strips. Preliminary results suggest the reliability of the method: strips from a healthy donor are characterized by a well-organized tissue, where the local disarray is log-normally distributed and slightly depends on the spatial scale of analysis; on the contrary, pathological strips show pronounced tissue disorganization, characterized by local disarray significantly dependent on the spatial scale of analysis. A virtual sample generator is developed to link this multi-scale disarray analysis with the underlying cellular architecture. This approach allowed us to quantitatively assess tissue organization in terms of 3D myocyte angular dispersion and may pave the way for developing novel predictive models based on structural data at cellular resolution.

Technical Note: Histological validation of anatomical imaging for breast modeling using a novel cryo-microtome

PURPOSE

Precise correlation between three-dimensional (3D) imaging and histology can aid biomechanical modeling of the breast. We develop a framework to register ex vivo images to histology using a novel cryo-fluorescence tomography (CFT) device.

METHODS

A formalin-fixed cadaveric breast specimen, including chest wall, was subjected to high-resolution magnetic resonance (MR) imaging. The specimen was then frozen and embedded in an optimal cutting temperature (OCT) compound. The OCT block was placed in a CFT device with an overhead camera and 50 μm thick slices were successively shaved off the block. After each shaving, the block-face was photographed. At select locations including connective/adipose tissue, muscle, skin, and fibroglandular tissue, 20 μm sections were transferred onto cryogenic tape for manual hematoxylin and eosin staining, histological assessment, and image capture. A 3D white-light image was automatically reconstructed from the photographs by aligning fiducial markers embedded in the OCT block. The 3D MR image, 3D white-light image, and photomicrographs were rigidly registered. Target registration errors (TREs) were computed based on 10 pairs of points marked at fibroglandular intersections. The overall MR-histology registration was used to compare the MR intensities at tissue extraction sites with a one-way analysis of variance.

RESULTS

The MR image to CFT-captured white-light image registration achieved a mean TRE of 0.73 ± 0.25 mm (less than the 1 mm MR slice resolution). The block-face white-light image and block-face photomicrograph registration showed visually indistinguishable alignment of anatomical structures and tissue boundaries. The MR intensities at the four tissue sites identified from histology differed significantly (p < 0.01). Each tissue pair, except the skin-connective/adipose tissue pair, also had significantly different MR intensities (p < 0.01).

CONCLUSIONS

Fine sectioning in a highly controlled imaging/sectioning environment enables accurate registration between the MR image and histology. Statistically significant differences in MR signal intensities between histological tissues are indicators for the specificity of correlation between MRI and histology.

When Stiffness Matters: Mechanosensing in Heart Development and Disease

During embryonic morphogenesis, the heart undergoes a complex series of cellular phenotypic maturations (e.g., transition of myocytes from proliferative to quiescent or maturation of the contractile apparatus), and this involves stiffening of the extracellular matrix (ECM) acting in concert with morphogenetic signals. The maladaptive remodeling of the myocardium, one of the processes involved in determination of heart failure, also involves mechanical cues, with a progressive stiffening of the tissue that produces cellular mechanical damage, inflammation, and ultimately myocardial fibrosis. The assessment of the biomechanical dependence of the molecular machinery (in myocardial and non-myocardial cells) is therefore essential to contextualize the maturation of the cardiac tissue at early stages and understand its pathologic evolution in aging. Because systems to perform multiscale modeling of cellular and tissue mechanics have been developed, it appears particularly novel to design integrated mechano-molecular models of heart development and disease to be tested in ex vivo reconstituted cells/tissue-mimicking conditions. In the present contribution, we will discuss the latest implication of mechanosensing in heart development and pathology, describe the most recent models of cell/tissue mechanics, and delineate novel strategies to target the consequences of heart failure with personalized approaches based on tissue engineering and induced pluripotent stem cell (iPSC) technologies.

Modeling Emergent Properties in the Brain Using Tissue Models to Investigate Neurodegenerative Disease

Here we describe emergent properties of the brain and the key challenges associated with modelling them in vitro. Modeling emergent properties of the brain will provide insights into brain function, development, and disease.

Estimating Vocal Fold Contact Pressure from Raw Laryngeal High-Speed Videoendoscopy Using a Hertz Contact Model

The development of trauma-induced lesions of the vocal folds (VFs) has been linked to a high collision pressure on the VF surface. However, there are no direct methods for the clinical assessment of VF collision, thus limiting the objective assessment of these disorders. In this study, we develop a video processing technique to directly quantify the mechanical impact of the VFs using solely laryngeal kinematic data. The technique is based on an edge tracking framework that estimates the kinematic sequence of each VF edge with a Kalman filter approach and a Hertzian impact model to predict the contact force during the collision. The proposed formulation overcomes several limitations of prior efforts since it uses a more relevant VF contact geometry, it does not require calibrated physical dimensions, it is normalized by the tissue properties, and it applies a correction factor for using a superior view only. The proposed approach is validated against numerical models, silicone vocal fold models, and prior studies. A case study with high-speed videoendoscopy recordings provides initial insights between the sound pressure level and contact pressure. Thus, the proposed method has a high potential in clinical practice and could also be adapted to operate with laryngeal stroboscopic systems.

A method to directly estimate the contact pressure of the vocal folds using uncalibrated laryngeal kinematic data is presented. The approach is promising in enhancing the objective assessment of vocal function in clinical settings, especially for studying same-subject variations.

Modeling Lung Architecture in the XCAT Series of Phantoms: Physiologically Based Airways, Arteries and Veins

The purpose of this paper was to extend the extended cardiac-torso (XCAT) series of computational phantoms to include a detailed lung architecture including airways and pulmonary vasculature. Eleven XCAT phantoms of varying anatomy were used in this paper. The lung lobes and initial branches of the airways, pulmonary arteries, and veins were previously defined in each XCAT model. These models were extended from the initial branches of the airways and vessels to the level of terminal branches using an anatomically-based volume-filling branching algorithm. This algorithm grew the airway and vasculature branches separately and iteratively without intersecting each other using cylindrical models with diameters estimated by order-based anatomical measurements. Geometrical features of the extended branches were compared with the literature anatomy values to quantitatively evaluate the models. These features include branching angle, length to diameter ratio, daughter to parent diameter ratio, asymmetrical branching pattern, diameter, and length ratios. The XCAT phantoms were then used to simulate CT images to qualitatively compare them with the original phantom images. The proposed growth model produced 46369 ± 12521 airways, 44737 ± 11773 arteries, and 39819 ± 9988 veins to the XCAT phantoms. Furthermore, the growth model was shown to produce asymmetrical airway, artery, and vein networks with geometrical attributes close to morphometry and model based studies. The simulated CT images of the phantoms were judged to be more realistic, including more airways and pulmonary vessels compared with the original phantoms. Future work will seek to add a heterogeneous parenchymal background into the XCAT lungs to make the phantoms even more representative of human anatomy, paving the way towards the use of XCAT models as a tool to virtually evaluate the current and emerging medical imaging technologies.

B-mode and acoustic radiation force impulse (ARFI) imaging of prostate zonal anatomy: comparison with 3T T2-weighted MR imaging

Prostate cancer (PCa) is the most common non-cutaneous malignancy among men in the United States and the second leading cause of cancer-related death. Multi-parametric magnetic resonance imaging (mpMRI) has gained recent popularity to characterize PCa. Acoustic Radiation Force Impulse (ARFI) imaging has the potential to aid PCa diagnosis and management by using tissue stiffness to evaluate prostate zonal anatomy and lesions. MR and B-mode/ARFI in vivo imaging datasets were compared with one another and with gross pathology measurements made immediately after radical prostatectomy. Images were manually segmented in 3D Slicer to delineate the central gland (CG) and prostate capsule, and 3D models were rendered to evaluate zonal anatomy dimensions and volumes. Both imaging modalities showed good correlation between estimated organ volume and gross pathologic weights. Ultrasound and MR total prostate volumes were well correlated (R(2) = 0.77), but B-mode images yielded prostate volumes that were larger (16.82% ± 22.45%) than MR images, due to overestimation of the lateral dimension (18.4% ± 13.9%), with less significant differences in the other dimensions (7.4% ± 17.6%, anterior-to-posterior, and -10.8% ± 13.9%, apex-to-base). ARFI and MR CG volumes were also well correlated (R(2) = 0.85). CG volume differences were attributed to ARFI underestimation of the apex-to-base axis (-28.8% ± 9.4%) and ARFI overestimation of the lateral dimension (21.5% ± 14.3%). B-mode/ARFI imaging yielded prostate volumes and dimensions that were well correlated with MR T2-weighted image (T2WI) estimates, with biases in the lateral dimension due to poor contrast caused by extraprostatic fat. B-mode combined with ARFI imaging is a promising low-cost, portable, real-time modality that can complement mpMRI for PCa diagnosis, treatment planning, and management.

Identification of reduced-order thermal therapy models using thermal MR images: theory and validation

In this paper, we develop and validate a method to identify computationally efficient site- and patient-specific models of ultrasound thermal therapies from MR thermal images. The models of the specific absorption rate of the transduced energy and the temperature response of the therapy target are identified in the reduced basis of proper orthogonal decomposition of thermal images, acquired in response to a mild thermal test excitation. The method permits dynamic reidentification of the treatment models during the therapy by recursively utilizing newly acquired images. Such adaptation is particularly important during high-temperature therapies, which are known to substantially and rapidly change tissue properties and blood perfusion. The developed theory was validated for the case of focused ultrasound heating of a tissue phantom. The experimental and computational results indicate that the developed approach produces accurate low-dimensional treatment models despite temporal and spatial noises in MR images and slow image acquisition rate.

Analysis of human fibroadenomas using three-dimensional impedance maps

Three-dimensional impedance maps (3DZMs) are virtual volumes of acoustic impedance values constructed from histology to represent tissue microstructure acoustically. From the 3DZM, the ultrasonic backscattered power spectrum can be predicted and model based scatterer properties, such as effective scatterer diameter (ESD), can be estimated. Additionally, the 3DZM can be exploited to visualize and identify possible scattering sites, which may aid in the development of more effective scattering models to better represent the ultrasonic interaction with underlying tissue microstructure. In this study, 3DZMs were created from a set of human fibroadenoma samples. ESD estimates were made assuming a fluid-filled sphere form factor model from 3DZMs of volume 300×300×300 μm. For a collection of 33 independent human fibroadenoma tissue samples, the ESD was estimated to be 111±40.7 μm. The 3DZMs were then investigated visually to identify possible scattering sources which conformed to the estimated model scatterer dimensions. This estimation technique allowed a better understanding of the spatial distribution and variability of the estimates throughout the volume.