Oxygen transport in a three-dimensional microvascular network incorporated with early tumour growth and preexisting vessel cooption: Numerical simulation study

We propose a dynamic mathematical model of tissue oxygen transport by a preexisting three-dimensional microvascular network which provides nutrients for an in situ cancer at the very early stage of primary microtumour growth. The expanding tumour consumes oxygen during its invasion to the surrounding tissues and cooption of host vessels. The preexisting vessel cooption, remodelling and collapse are modelled by the changes of haemodynamic conditions due to the growing tumour. A detailed computational model of oxygen transport in tumour tissue is developed by considering (a) the time-varying oxygen advection diffusion equation within the microvessel segments, (b) the oxygen flux across the vessel walls, and (c) the oxygen diffusion and consumption with in the tumour and surrounding healthy tissue. The results show the oxygen concentration distribution at different time points of early tumour growth. In addition, the influence of preexisting vessel density on the oxygen transport has been discussed. The proposed model not only provides a quantitative approach for investigating the interactions between tumour growth and oxygen delivery, but also is extendable to model other molecules or chemotherapeutic drug transport in the future study.

Deep transverse metatarsal ligaments mechanical response during landing

The deep transverse metatarsal ligaments play an important role in stabilizing the metatarsal bones and manipulating foot transverse arch deformation. However, the biomechanical research about transverse metatarsal ligaments in the foot maneuver is quite few. Due to the difficulties and lack of better measurement technology for these ligaments experimental monitor, the load transfer mechanism and internal stress state also hadn't been well addressed. The purpose of this study was to develop a detailing foot finite element model including transverse metatarsal ligaments tissues, to investigate the mechanical response of transverse metatarsal ligaments during the landing condition. The transverse metatarsal ligaments were considered as hyperelastic material model was used to represent the nonlinear and nearly incompressible nature of the ligament tissue. From the simulation results, it is clearly to find that the peak maiximal principal stress of transverse metatarsal ligaments was between the third and fourth metatarsals. Meanwhile, it seems the transverse metatarsal ligaments in the middle position experienced higher tension than the sides transverse metatarsal ligaments.

A new method for the in vivo identification of mechanical properties in arteries from cine MRI images: Theoretical framework and validation

Quantifying the stiffness properties of soft tissues is essential for the diagnosis of many cardiovascular diseases such as atherosclerosis. In these pathologies it is widely agreed that the arterial wall stiffness is an indicator of vulnerability. The present paper focuses on the carotid artery and proposes a new inversion methodology for deriving the stiffness properties of the wall from cine-MRI (magnetic resonance imaging) data. We address this problem by setting-up a cost function defined as the distance between the modeled pixel signals and the measured ones. Minimizing this cost function yields the unknown stiffness properties of both the arterial wall and the surrounding tissues. The sensitivity of the identified properties to various sources of uncertainty is studied. Validation of the method is performed on a rubber phantom. The elastic modulus identified using the developed methodology lies within a mean error of 9.6%. It is then applied to two young healthy subjects as a proof of practical feasibility, with identified values of 625 kPa and 587 kPa for one of the carotid of each subject.

3D numerical simulation of avascular tumour growth: effect of hypoxic micro-environment in host tissue

A three-dimensional (3D) mathematical model of tumour growth at the avascular phase and vessel remodelling in host tissues is proposed with emphasis on the study of the interactions of tumour growth and hypoxic micro-environment in host tissues. The hybrid based model includes the continuum part, such as the distributions of oxygen and vascular endothelial growth factors (VEGFs), and the discrete part of tumour cells (TCs) and blood vessel networks. The simulation shows the dynamic process of avascular tumour growth from a few initial cells to an equilibrium state with varied vessel networks. After a phase of rapidly increasing numbers of the TCs, more and more host vessels collapse due to the stress caused by the growing tumour. In addition, the consumption of oxygen expands with the enlarged tumour region. The study also discusses the effects of certain factors on tumour growth, including the density and configuration of preexisting vessel networks and the blood oxygen content. The model enables us to examine the relationship between early tumour growth and hypoxic micro-environment in host tissues, which can be useful for further applications, such as tumour metastasis and the initialization of tumour angiogenesis.

Finite element analysis of deep transverse metatarsal ligaments mechanical response during landing

The deep transverse metatarsal ligaments (DTML) play an important role in stabilizing the metatarsal bones and manipulating foot transverse arch deformation. However, the biomechanical research about DTML in the foot maneuver is quite few. Due to the difficulties and lack of better measurement technology for these ligaments experimental monitor, the load transfer mechanism and internal stress state also hadn't been well addressed. The purpose of this study was to develop a detailing foot finite element model including DTML tissues, to investigate the mechanical response of DTML during the landing condition. The DTML was considered as hyperelastic material model was used to represent the nonlinear and nearly incompressible nature of the ligament tissue. From the simulation results, it is clearly to find that the peak maiximal principal stress of DTML was between the third and fourth metatarsals. Meanwhile, it seems the DTML in the middle position experienced higher tension than the sides DTML.

In vivo velocity vector imaging and time-resolved strain rate measurements in the wall of blood vessels using MRI

In this paper, we present a new approach for velocity vector imaging and time-resolved measurements of strain rates in the wall of human arteries using MRI and we prove its feasibility on two examples: in vitro on a phantom and in vivo on the carotid artery of a human subject. Results point out the promising potential of this approach for investigating the mechanics of arterial tissues in vivo.

Study on tracheal collapsibility, compliance, and stress by considering nonlinear mechanical property of cartilage

Tracheal cartilage has been widely regarded as a linear elastic material either in experimental studies or in analytic and numerical models. However, it has been recently demonstrated that, like other fiber-oriented biological tissues, tracheal cartilage is a nonlinear material, which displays higher strength in compression than in extension. Considering the nonlinearity requires a more complex theoretical frame work and costs more to simulate. This study aims to quantify the deviation due to the simplified treatment of the tracheal cartilage as a linear material. It also evaluates the improved accuracy gained by considering the nonlinearity. Pig tracheal rings were used to exam the mechanical properties of cartilage and muscular membrane. By taking into account the asymmetric shape of tracheal cartilage, the collapse behavior of complete rings was simulated, and the compliance of airway and stress in the muscular membrane were discussed. The results obtained were compared with those assuming linear mechanical properties. The following results were found: (1) Models based on both types of material properties give a small difference in representing collapse behavior; (2) regarding compliance, the relative difference is big, ranging from 10 to 40% under negative pressure conditions; and (3) the difference in determining stress in the muscular membrane is small too: <5%. In conclusion, treating tracheal cartilage as a linear material will not cause big deviations in representing the collapse behavior, and mechanical stress in the muscular part, but it will induce a big deviation in predicting the compliance, particularly when the transmural pressure is lower than -0.5 kPa. The results obtained in this study may be useful in both understanding the collapse behavior of trachea and in evaluating the error induced by the simplification of treating the tracheal cartilage as a linear elastic material.

Exploring methods of preparing functional cartilage-bone xenografts for joint repair

This thesis explores the feasibility of donor-receiver concept for joint replacement where cartilage-bone tissues can be taken from either human or other mammals and prepared scientifically for repairing focal joint defects in knees, hips and shoulders. The manufactured construct is immunologically inert and is capable of acting as a scaffold for engineering new cartilage-bone laminates when placed in the joint. Innovative manufacturing procedures and assessment techniques were developed for appraising this tissue-based scaffold. This research has demonstrated that tissue replacement technology can be applied in situations where blood vessels are absent such as in articular cartilage.

Use of electronic portal imaging devices for electron treatment verification

This study aims to help broaden the use of electronic portal imaging devices (EPIDs) for pre-treatment patient positioning verification, from photon-beam radiotherapy to photon- and electron-beam radiotherapy, by proposing and testing a method for acquiring clinicallyuseful EPID images of patient anatomy using electron beams, with a view to enabling and encouraging further research in this area. EPID images used in this study were acquired using all available beams from a linac configured to deliver electron beams with nominal energies of 6, 9, 12, 16 and 20 MeV, as well as photon beams with nominal energies of 6 and 10 MV. A widely-available heterogeneous, approximately-humanoid, thorax phantom was used, to provide an indication of the contrast and noise produced when imaging different types of tissue with comparatively realistic thicknesses. The acquired images were automatically calibrated, corrected for the effects of variations in the sensitivity of individual photodiodes, using a flood field image. For electron beam imaging, flood field EPID calibration images were acquired with and without the placement of blocks of water-equivalent plastic (with thicknesses approximately equal to the practical range of electrons in the plastic) placed upstream of the EPID, to filter out the primary electron beam, leaving only the bremsstrahlung photon signal. While the electron beam images acquired using a standard (unfiltered) flood field calibration were observed to be noisy and difficult to interpret, the electron beam images acquired using the filtered flood field calibration showed tissues and bony anatomy with levels of contrast and noise that were similar to the contrast and noise levels seen in the clinically acceptable photon beam EPID images. The best electron beam imaging results (highest contrast, signal-to-noise and contrast-to-noise ratios) were achieved when the images were acquired using the higher energy electron beams (16 and 20 MeV) when the EPID was calibrated using an intermediate (12 MeV) electron beam energy. These results demonstrate the feasibility of acquiring clinically-useful EPID images of patient anatomy using electron beams and suggest important avenues for future investigation, thus enabling and encouraging further research in this area. There is manifest potential for the EPID imaging method proposed in this work to lead to the clinical use of electron beam imaging for geometric verification of electron treatments in the future.

Fusion transcript loci share many genomic features with non-fusion loci

Background Fusion transcripts are found in many tissues and have the potential to create novel functional products. Here, we investigate the genomic sequences around fusion junctions to better understand the transcriptional mechanisms mediating fusion transcription/splicing. We analyzed data from prostate (cancer) cells as previous studies have shown extensively that these cells readily undergo fusion transcription. Results We used the FusionMap program to identify high-confidence fusion transcripts from RNAseq data. The RNAseq datasets were from our (N = 8) and other (N = 14) clinical prostate tumors with adjacent non-cancer cells, and from the LNCaP prostate cancer cell line that were mock-, androgen- (DHT), and anti-androgen- (bicalutamide, enzalutamide) treated. In total, 185 fusion transcripts were identified from all RNAseq datasets. The majority (76 %) of these fusion transcripts were ‘read-through chimeras’ derived from adjacent genes in the genome. Characterization of sequences at fusion loci were carried out using a combination of the FusionMap program, custom Perl scripts, and the RNAfold program. Our computational analysis indicated that most fusion junctions (76 %) use the consensus GT-AG intron donor-acceptor splice site, and most fusion transcripts (85 %) maintained the open reading frame. We assessed whether parental genes of fusion transcripts have the potential to form complementary base pairing between parental genes which might bring them into physical proximity. Our computational analysis of sequences flanking fusion junctions at parental loci indicate that these loci have a similar propensity as non-fusion loci to hybridize. The abundance of repetitive sequences at fusion and non-fusion loci was also investigated given that SINE repeats are involved in aberrant gene transcription. We found few instances of repetitive sequences at both fusion and non-fusion junctions. Finally, RT-qPCR was performed on RNA from both clinical prostate tumors and adjacent non-cancer cells (N = 7), and LNCaP cells treated as above to validate the expression of seven fusion transcripts and their respective parental genes. We reveal that fusion transcript expression is similar to the expression of parental genes. Conclusions Fusion transcripts maintain the open reading frame, and likely use the same transcriptional machinery as non-fusion transcripts as they share many genomic features at splice/fusion junctions.

Epithelial-mesenchymal axis in head and neck cancer cell lines

Objectives: Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous tumour type which necessitates multiple invitro models to attain an appreciation of its multiple subtypes. The phenomenon of epithelial-mesenchymal transition (EMT) isimportant to the development of a metastatic cancer cell phenotype being relevant to the ability of cancer cells to intravasate intovasculature and to invade tissues. The role of EMT in human papilloma virus (HPV) positive HNSCC is not well understood. Thispaper aims to characterize seven HNSCC cell lines (FaDu, SCC-25, SCC-15, CAL27, RPMI2650) including two new HPV-16positive HNSCC cell lines (UD-SCC2, 93-VU-147T) for their epithelial and mesenchymal properties. Materials and methods: A panel of HNSCC cell lines from multiple head and neck anatomical sites were profiled for basalexpression of epithelial and mesenchymal characteristics at mRNA, protein and functional levels (proliferative, migratory andinvasive properties). Furthermore, 3D spheroid forming capabilities were investigated. Results: We found that the HPV-16 positive cell line, in particular UD-SCC2 demonstrated a more invasive and mesenchymalphenotype at the molecular and functional levels suggesting HPV infection may mediate some of these cellular properties.Moreover, HPV-negative cell lines were not strictly epithelial presenting with a dynamic range of expression. Conclusions: This study presents the molecular and phenotypic diversity of HNSCC cell lines. It highlights the need formore studies in this field and a scoring system where HNSCC cell lines are ranked according to their respective epithelial andmesenchymal nature. This data will be useful to anyone modelling HNSCC behaviour, providing a molecular context which willenable them to decipher cell phenotypes and to develop therapies which block EMT progression.

Prolonged hyperinsulinemia affects metabolic signal transduction markers in a tissue specific manner

Insulin dysregulation is common in horses although the mechanisms of metabolic dysfunction are poorly understood. We hypothesized that insulin signaling in striated (cardiac and skeletal) muscle and lamellae may be mediated through different receptors as a result of receptor content, and that transcriptional regulation of downstream signal transduction and glucose transport may also differ between tissues sites during hyperinsulinemia. Archived samples from horses treated with a prolonged insulin infusion or a balanced electrolyte solution were used. All treated horses developed marked hyperinsulinemia and clinical laminitis. Protein expression was compared across tissues for the insulin receptor and insulin-like growth factor 1 receptor (IGF-1R) by immunoblotting. Gene expression of metabolic insulin-signaling markers (insulin receptor substrate 1, Akt2, and glycogen synthase kinase 3 beta [GSK-3β]) and glucose transport (basal glucose transporter 1 and insulin-sensitive glucose transporter 4) was evaluated using real-time reverse transcription polymerase chain reaction. Lamellar tissue contained significantly more IGF-1R protein than skeletal muscle, indicating the potential significance of IGF-1R signaling for this tissue. Gene expression of the selected markers of insulin signaling and glucose transport in skeletal muscle and lamellar tissues was unaffected by prolonged hyperinsulinemia. In contrast, the significant upregulation of Akt2, GSK-3β, GLUT1, and GLUT4 gene expression in cardiac tissue suggested that the prolonged hyperinsulinemia induced an increase in insulin sensitivity and a transcriptional activation of glucose transport. Responses to insulin are tissue-specific, and extrapolation of data across tissue sites is inappropriate.

Diet-induced hypercholesterolemia promotes androgen-independent prostate cancer metastasis via IQGAP1 and caveolin-1

Obesity and metabolic syndrome are associated with several cancers, however, the molecular mechanisms remain to be fully elucidated. Recent studies suggest that hypercholesterolemia increases intratumoral androgen signaling in prostate cancer, but it is unclear whether androgenindependent mechanisms also exist. Since hypercholesterolemia is associated with advanced, castrate-resistant prostate cancer, in this study, we aimed to determine whether and how hypercholesterolemia affects prostate cancer progression in the absence of androgen signaling. We demonstrate that diet-induced hypercholesterolemia promotes orthotopic xenograft PC-3 cell metastasis, concomitant with elevated expression of caveolin-1 and IQGAP1 in xenograft tumor tissues. In vitro cholesterol treatment of PC-3 cells stimulated migration and increased IQGAP1 and caveolin-1 protein level and localization to a detergent-resistant fraction. Down-regulation of caveolin-1 or IQGAP1 in PC-3 cells reduced migration and invasion in vitro, and hypercholesterolemia-induced metastasis in vivo. Double knock-down of caveolin-1 and IQGAP1 showed no additive effect, suggesting that caveolin-1 and IQGAP1 act via the same pathway. Taken together, our data show that hypercholesterolemia promotes prostate cancer metastasis independent of the androgen pathway, in part by increasing IQGAP1 and caveolin-1. These results have broader implications for managing metastasis of cancers in general as IQGAP1 and hypercholesterolemia are implicated in the progression of several cancers.

Miniscleral lens wear influences corneal curvature and optics

Purpose Transient changes in corneal topography associated with soft and conventional or reverse geometry rigid contact lens wear have been well documented; however, only a few studies have examined the influence of scleral contact lens wear upon the cornea. Therefore, in this study, we examined the influence of modern miniscleral contact lenses, which land entirely on the sclera and overlying tissues, upon anterior corneal curvature and optics. Methods Anterior corneal topography and elevation data were acquired using Scheimpflug imaging (Pentacam HR, Oculus) immediately prior to and following 8 hours of miniscleral contact lens wear in 15 young healthy adults (mean age 22 ± 3 years, 8 East Asian, 7 Caucasian) with normal corneae. Corneal diurnal variations were accounted for using data collected on a dedicated measurement day without contact lens wear. Corneal clearance was quantified using an optical coherence tomographer (RS-3000, Nidek) following lens insertion and after 8 hours of lens wear. Results Although corneal clearance was maintained throughout the 8 hour lens wear period, significant corneal flattening (up to 0.08 ± 0.04 mm) was observed, primarily in the superior mid-peripheral cornea, which resulted in a slight increase in against-the-rule corneal astigmatism (mean +0.02/-0.15 x 94 for an 8 mm diameter). Higher order aberration terms of horizontal coma, vertical coma and spherical aberration all underwent significant changes for an 8 mm corneal diameter (p ≤ 0.01), which typically resulted in a decrease in RMS error values (mean change in total higher order RMS -0.035 ± 0.046 µm for an 8 mm diameter). There was no association between the magnitude of change in central or mid-peripheral corneal clearance during lens wear and the observed changes in corneal curvature (p > 0.05). However, Asian participants displayed a significantly greater reduction in corneal clearance (p = 0.04) and greater superior-nasal corneal flattening compared to Caucasians (p = 0.048). Conclusions Miniscleral contact lenses that vault the cornea induce significant changes in anterior corneal surface topography and higher order aberrations following 8 hours of lens wear. The region of greatest corneal flattening was observed in the superior-nasal mid-periphery, more so in Asian participants. Practitioners should be aware that corneal measurements obtained following miniscleral lens removal may mask underlying corneal steepening.

Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.