Whole-brain imaging with single-cell resolution using chemical cocktails and computational analysis

Systems-level identification and analysis of cellular circuits in the brain will require the development of whole-brain imaging with single-cell resolution. To this end, we performed comprehensive chemical screening to develop a whole-brain clearing and imaging method, termed CUBIC (clear, unobstructed brain imaging cocktails and computational analysis). CUBIC is a simple and efficient method involving the immersion of brain samples in chemical mixtures containing aminoalcohols, which enables rapid whole-brain imaging with single-photon excitation microscopy. CUBIC is applicable to multicolor imaging of fluorescent proteins or immunostained samples in adult brains and is scalable from a primate brain to subcellular structures. We also developed a whole-brain cell-nuclear counterstaining protocol and a computational image analysis pipeline that, together with CUBIC reagents, enable the visualization and quantification of neural activities induced by environmental stimulation. CUBIC enables time-course expression profiling of whole adult brains with single-cell resolution.

Whole-Body Imaging with Single-Cell Resolution by Tissue Decolorization

The development of whole-body imaging at single-cell resolution enables system-level approaches to studying cellular circuits in organisms. Previous clearing methods focused on homogenizing mismatched refractive indices of individual tissues, enabling reductions in opacity but falling short of achieving transparency. Here, we show that an aminoalcohol decolorizes blood by efficiently eluting the heme chromophore from hemoglobin. Direct transcardial perfusion of an aminoalcohol-containing cocktail that we previously termed CUBIC coupled with a 10 day to 2 week clearing protocol decolorized and rendered nearly transparent almost all organs of adult mice as well as the entire body of infant and adult mice. This CUBIC-perfusion protocol enables rapid whole-body and whole-organ imaging at single-cell resolution by using light-sheet fluorescent microscopy. The CUBIC protocol is also applicable to 3D pathology, anatomy, and immunohistochemistry of various organs. These results suggest that whole-body imaging of colorless tissues at high resolution will contribute to organism-level systems biology.

Stability and convergence of a new finite volume method for a two-sided space-fractional diffusion equation

In this paper, we consider a two-sided space-fractional diffusion equation with variable coefficients on a finite domain. Firstly, based on the nodal basis functions, we present a new fractional finite volume method for the two-sided space-fractional diffusion equation and derive the implicit scheme and solve it in matrix form. Secondly, we prove the stability and convergence of the implicit fractional finite volume method and conclude that the method is unconditionally stable and convergent. Finally, some numerical examples are given to show the effectiveness of the new numerical method, and the results are in excellent agreement with theoretical analysis.

Reduced intensity conditioning for allogeneic hematopoietic stem cell transplant determines the kinetics of acute Graft-Versus-Host Disease

Background Preparative myeloablative conditioning regimens for allogeneic hematopoietic stem-cell transplantation (HSCT) may control malignancy and facilitate engraftment but also contribute to transplant related mortality, cytokine release, and acute graft-versus-host disease (GVHD). Reduced intensity conditioning (RIC) regimens have decreased transplant related mortality but the incidence of acute GVHD, while delayed, remains unchanged. There are currently no in vivo allogeneic models of RIC HSCT, limiting studies into the mechanism behind RIC-associated GVHD. Methods We developed two RIC HSCT models that result in delayed onset GVHD (major histocompatibility complex mismatched (UBI-GFP/BL6 [H-2b]→BALB/c [H-2d]) and major histocompatibility complex matched, minor histocompatibility mismatched (UBI-GFP/BL6 [H-2b]→BALB.B [H-2b])) enabling the effect of RIC on chimerism, dendritic cell (DC) chimerism, and GVHD to be investigated. Results In contrast with myeloablative conditioning, we observed that RIC-associated delayed-onset GVHD is characterized by low production of tumor necrosis factor-α, maintenance of host DC, phenotypic DC activation, increased T-regulatory cell numbers, and a delayed emergence of activated donor DC. Furthermore, changes to the peritransplant milieu in the recipient after RIC lead to the altered activation of DC and the induction of T-regulatory responses. Reduced intensity conditioning recipients suffer less early damage to GVHD target organs. However, as donor cells engraft, activated donor DC and rising levels of tumor necrosis factor-α are associated with a later onset of severe GVHD. Conclusions Delineating the mechanisms underlying delayed onset GVHD in RIC HSCT recipients is vital to improve the prediction of disease onset and allow more targeted interventions for acute GVHD.

Developing a therapeutic anti-dendritic cell antibody to prevent graft versus host disease

Host and donor dendritic cells (DC) stimulate alloreactive donor T lymphocytes, and initiate GVHD. We have shown that polyclonal antibody to the DC surface activation marker human CD83 (anti hCD83), which depletes activated DC, can prevent human DC and T cell induced lethal xenogeneic GVHD in SCID mice without impairing T cell mediated anti-leukaemic and anti-viral (CMV and influenza) immunity (J Exp Med 2009; 206: 387). Therefore, we made and tested a polyclonal anti mouse CD83 (RAM83) antibody in murine HSCT models and developed a human mAb against hCD83 as a potential new therapeutic immunosuppressive agent.

Drug accumulation into single drug-sensitive and drug-resistant prostate cancer cells conducted on the single cell bioanalyzer

Multidrug resistance (MDR) occurs in prostate cancer, and this happens when the cancer cells resist chemotherapeutic drugs by pumping them out of the cells. MDR inhibitors such as cyclosporin A (CsA) can stop the pumping and enhance the drugs accumulated in the cells. The cellular drug accumulation is monitored using a microfluidic chip mounted on a single cell bioanalyzer. This equipment has been developed to measure accumulation of drugs such as doxorubicin (DOX) and fluorescently labeled paclitaxel (PTX) in single prostate cancer cells. The inhibition of drug efflux on the same prostate cell was examined in drug-sensitive and drug-resistant cells. Accumulation of these drug molecules was not found in the MDR cells, PC-3 RX-DT2R cells. Enhanced drug accumulation was observed only after treating the MDR cell in the presence of 5 μM of CsA as the MDR inhibitor. We envision this monitoring of the accumulation of fluorescent molecules (drug or fluorescent molecules), if conducted on single patient cancer cells, can provide information for clinical monitoring of patients undergoing chemotherapy in the future.

Numerical investigation of motion and deformation of a single red blood cell in a stenosed capillary

It is generally assumed that influence of the red blood cells (RBCs) is predominant in blood rheology. The healthy RBCs are highly deformable and can thus easily squeeze through the smallest capillaries having internal diameter less than their characteristic size. On the other hand, RBCs infected by malaria or other diseases are stiffer and so less deformable. Thus it is harder for them to flow through the smallest capillaries. Therefore, it is very important to critically and realistically investigate the mechanical behavior of both healthy and infected RBCs which is a current gap in knowledge. The motion and the steady state deformed shape of the RBCs depend on many factors, such as the geometrical parameters of the capillary through which blood flows, the membrane bending stiffness and the mean velocity of the blood flow. In this study, motion and deformation of a single two-dimensional RBC in a stenosed capillary is explored by using smoothed particle hydrodynamics (SPH) method. An elastic spring network is used to model the RBC membrane, while the RBC's inside fluid and outside fluid are treated as SPH particles. The effect of RBC's membrane stiffness (kb), inlet pressure (P) and geometrical parameters of the capillary on the motion and deformation of the RBC is studied. The deformation index, RBC's mean velocity and the cell membrane energy are analyzed when the cell passes through the stenosed capillary. The simulation results demonstrate that the kb, P and the geometrical parameters of the capillary have a significant impact on the RBCs' motion and deformation in the stenosed section.

Melanoma cell colony expansion parameters revealed by approximate Bayesian computation

In vitro studies and mathematical models are now being widely used to study the underlying mechanisms driving the expansion of cell colonies. This can improve our understanding of cancer formation and progression. Although much progress has been made in terms of developing and analysing mathematical models, far less progress has been made in terms of understanding how to estimate model parameters using experimental in vitro image-based data. To address this issue, a new approximate Bayesian computation (ABC) algorithm is proposed to estimate key parameters governing the expansion of melanoma cell (MM127) colonies, including cell diffusivity, D, cell proliferation rate, λ, and cell-to-cell adhesion, q, in two experimental scenarios, namely with and without a chemical treatment to suppress cell proliferation. Even when little prior biological knowledge about the parameters is assumed, all parameters are precisely inferred with a small posterior coefficient of variation, approximately 2–12%. The ABC analyses reveal that the posterior distributions of D and q depend on the experimental elapsed time, whereas the posterior distribution of λ does not. The posterior mean values of D and q are in the ranges 226–268 µm2h−1, 311–351 µm2h−1 and 0.23–0.39, 0.32–0.61 for the experimental periods of 0–24 h and 24–48 h, respectively. Furthermore, we found that the posterior distribution of q also depends on the initial cell density, whereas the posterior distributions of D and λ do not. The ABC approach also enables information from the two experiments to be combined, resulting in greater precision for all estimates of D and λ.

A cellular automata model to investigate immune cell–tumor cell interactions in growing tumors in two spatial dimensions

We develop a hybrid cellular automata model to describe the effect of the immune system and chemokines on a growing tumor. The hybrid cellular automata model consists of partial differential equations to model chemokine concentrations, and discrete cellular automata to model cell–cell interactions and changes. The computational implementation overlays these two components on the same spatial region. We present representative simulations of the model and show that increasing the number of immature dendritic cells (DCs) in the domain causes a decrease in the number of tumor cells. This result strongly supports the hypothesis that DCs can be used as a cancer treatment. Furthermore, we also use the hybrid cellular automata model to investigate the growth of a tumor in a number of computational “cancer patients.” Using these virtual patients, the model can explain that increasing the number of DCs in the domain causes longer “survival.” Not surprisingly, the model also reflects the fact that the parameter related to tumor division rate plays an important role in tumor metastasis.

Mesenchymal stem cells as mediators of the neuronal cell niche

This study examined the role of heparan sulfate proteoglycans (HSPGs) in neural lineage differentiation of human mesenchymal stem cells (hMSCs). Several HSPGs were identified as potential new targets controlling neural fate specification and may be applied to the development of improved models to examine and repair brain damage. hMSCs were characterised throughout extended in vitro expansion for neural lineage potential (neurons, astrocytes, oligodendrocytes) and differentiated using terminal differentiation and intermediate sphere formation. Brain damage and neurological disorders caused by injury or disease affect a large number of people often resulting in lifelong disabilities. Multipotent mesenchymal stem cells have a large capacity for self-renewal and provide an excellent model to examine the regulation and contribution of both stem cells and their surrounding microenvironment to the repair of neural tissue damage.

Treatment of silk fibroin with poly(ethylene glycol) for the enhancement of corneal epithelial cell growth

A silk protein, fibroin, was isolated from the cocoons of the domesticated silkworm (Bombyx mori) and cast into membranes to serve as freestanding templates for tissue-engineered corneal cell constructs to be used in ocular surface reconstruction. In this study, we sought to enhance the attachment and proliferation of corneal epithelial cells by increasing the permeability of the fibroin membranes and the topographic roughness of their surface. By mixing the fibroin solution with poly(ethylene glycol) (PEG) of molecular weight 300 Da, membranes were produced with increased permeability and with topographic patterns generated on their surface. In order to enhance their mechanical stability, some PEG-treated membranes were also crosslinked with genipin. The resulting membranes were thoroughly characterized and compared to the non-treated membranes. The PEG-treated membranes were similar in tensile strength to the non-treated ones, but their elastic modulus was higher and elongation lower, indicating enhanced rigidity. The crosslinking with genipin did not induce a significant improvement in mechanical properties. In cultures of a human-derived corneal epithelial cell line (HCE-T), the PEG treatment of the substratum did not improve the attachment of cells and it enhanced only slightly the cell proliferation in the longer term. Likewise, primary cultures of human limbal epithelial cells grew equally well on both non-treated and PEG-treated membranes, and the stratification of cultures was consistently improved in the presence of an underlying culture of irradiated 3T3 feeder cells, irrespectively of PEG-treatment. Nevertheless, the cultures grown on the PEG-treated membranes in the presence of feeder cells did display a higher nuclear-to-cytoplasmic ratio suggesting a more proliferative phenotype. We concluded that while the treatment with PEG had a significant effect on some structural properties of the B. mori silk fibroin (BMSF) membranes, there were minimal gains in the performance of these materials as a substratum for corneal epithelial cell growth. The reduced mechanical stability of freestanding PEG-treated membranes makes them a less viable choice than the non-treated membranes.

Numerical treatment of a two-dimensional variable-order fractional nonlinear reaction-diffusion model

A two-dimensional variable-order fractional nonlinear reaction-diffusion model is considered. A second-order spatial accurate semi-implicit alternating direction method for a two-dimensional variable-order fractional nonlinear reaction-diffusion model is proposed. Stability and convergence of the semi-implicit alternating direct method are established. Finally, some numerical examples are given to support our theoretical analysis. These numerical techniques can be used to simulate a two-dimensional variable order fractional FitzHugh-Nagumo model in a rectangular domain. This type of model can be used to describe how electrical currents flow through the heart, controlling its contractions, and are used to ascertain the effects of certain drugs designed to treat arrhythmia.

GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

Pressure placed on paediatric haematopoietic stem cell donors: Views from health professionals

Aim Paediatric haematopoietic stem cell donors undergo non-therapeutic procedures and endure known and unknown physical and psychosocial risks for the benefit of a family member. One ethical concern is the risk they may be pressured by parents or health professionals to act as a donor. This paper adds to what is known about this topic by presenting the views of health professionals. Methods This qualitative study involved semi-structured interviews with 14 health professionals in Australasia experienced in dealing with paediatric donors. Transcripts were analysed using established qualitative methodologies. Results Health professionals considered that some paediatric donors experience pressure to donate. Situations were identified that were likely to increase the risk of pressure being placed on donors and views were expressed about the ethical ‘appropriateness’ of these practices within the family setting. Conclusions Children may be subject to pressure from family and health professionals to be tested and act as donors, Therefore, our ethical obligation to these children extends to implementing donor focused processes – including independent health professionals and the appointment of a donor advocate – to assist in detecting and addressing instances of inappropriate pressure being placed on a child.