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.

Characteristics and alteration origins of matrix minerals in volcaniclastic kimberlite of the Muskox pipe (Nunavut, Canada)

The matrix of volcaniclastic kimberlite (VK) from the Muskox pipe (Northern Slave Province, Nunavut, Canada) is interpreted to represent an overprint of an original clastic matrix. Muskox VK is subdivided into three different matrix mineral assemblages that reflect differences in the proportions of original primary matrix constituents, temperature of formation and nature of the altering fluids. Using whole rock X-ray fluorescence (XRF), whole rock X-ray diffraction (XRD), microprobe analyses, back-scatter electron (BSE) imaging, petrography and core logging, we find that most matrix minerals (serpentine, phlogopite, chlorite, saponite, monticellite, Fe-Ti oxides and calcite) lack either primary igneous or primary clastic textures. The mineralogy and textures are most consistent with formation through alteration overprinting of an original clastic matrix that form by retrograde reactions as the deposit cools, or, in the case of calcite, by precipitation from Ca-bearing fluids into a secondary porosity. The first mineral assemblage consists largely of serpentine, phlogopite, calcite, Fe-Ti oxides and monticellite and occurs in VK with relatively fresh framework clasts. Alteration reactions, driven by deuteric fluids derived from the juvenile constituents, promote the crystallisation of minerals that indicate relatively high temperatures of formation (> 400 °C). Lower-temperature minerals are not present because permeability was occluded before the deposit cooled to low temperatures, thus shielding the facies from further interaction with fluids. The other two matrix mineral assemblages consist largely of serpentine, phlogopite, calcite, +/- diopside, and +/- chlorite. They form in VK that contains more country rock, which may have caused the deposit to be cooler upon emplacement. Most framework components are completely altered, suggesting that larger volumes of fluids drove the alteration reactions. These fluids were likely of meteoric provenance and became heated by the volcaniclastic debris when they percolated into the VK infill. Most alteration reactions ceased at temperatures > 200 °C, as indicated by the absence or paucity of lower-temperature phases in most samples, such as saponite. Recognition that Muskox VK contains an original clastic matrix is a necessary first step for evaluating the textural configuration, which is important for reconstructing the physical processes responsible for the formation of the deposit.

Hierarchical clustering of the genetic connectivity matrix reveals the network topology of gene action on brain microstructure: An N=531 twin study

Genetic correlation (rg) analysis determines how much of the correlation between two measures is due to common genetic influences. In an analysis of 4 Tesla diffusion tensor images (DTI) from 531 healthy young adult twins and their siblings, we generalized the concept of genetic correlation to determine common genetic influences on white matter integrity, measured by fractional anisotropy (FA), at all points of the brain, yielding an NxN genetic correlation matrix rg(x,y) between FA values at all pairs of voxels in the brain. With hierarchical clustering, we identified brain regions with relatively homogeneous genetic determinants, to boost the power to identify causal single nucleotide polymorphisms (SNP). We applied genome-wide association (GWA) to assess associations between 529,497 SNPs and FA in clusters defined by hubs of the clustered genetic correlation matrix. We identified a network of genes, with a scale-free topology, that influences white matter integrity over multiple brain regions.

High threshold of β1 integrin inhibition required to block collagen I-induced membrane type-1 matrix metalloproteinase (MT1-MMP) activation of matrix metalloproteinase 2 (MMP-2)

Background Matrix metalloproteinase-2 (MMP-2) is an endopeptidase that facilitates extracellular matrix remodeling and molecular regulation, and is implicated in tumor metastasis. Type I collagen (Col I) regulates the activation of MMP-2 through both transcriptional and post-transcriptional means; however gaps remain in our understanding of the involvement of collagen-binding ?1 integrins in collagen-stimulated MMP-2 activation. Methods Three ?1 integrin siRNAs were used to elucidate the involvement of ?1 integrins in the Col I-induced MMP-2 activation mechanism. ?1 integrin knockdown was analyzed by quantitative RT-PCR, Western Blot and FACS analysis. Adhesion assay and collagen gel contraction were used to test the biological effects of ?1 integrin abrogation. MMP-2 activation levels were monitored by gelatin zymography. Results All three ?1 integrin siRNAs were efficient at ?1 integrin knockdown and FACS analysis revealed commensurate reductions of integrins ?2 and ?3, which are heterodimeric partners of ?1, but not ?V, which is not. All three ?1 integrin siRNAs inhibited adhesion and collagen gel contraction, however only the siRNA showing the greatest magnitude of ?1 knockdown inhibited Col I-induced MMP-2 activation and reduced the accompanying upregulation of MT1-MMP, suggesting a dose response threshold effect. Re-transfection with codon-swapped ?1 integrin overcame the reduction in MMP-2 activation induced by Col-1, confirming the ?1 integrin target specificity. MMP-2 activation induced by TPA or Concanavalin A (Con A) was not inhibited by ?1 integrin siRNA knockdown. Conclusion Together, the data reveals that strong abrogation of ?1 integrin is required to block MMP-2 activation induced by Col I, which may have implications for the therapeutic targeting of ?1 integrin.

Effect of decellularization on the load-bearing characteristics of articular cartilage matrix

The application of decellularized extracellular matrices to aid tissue regeneration in reconstructive surgery and regenerative medicine has been promising. Several decellularization protocols for removing cellular materials from natural tissues such as heart valves are currently in use. This paper evaluates the feasibility of potential extension of this methodology relative to the desirable properties of load bearing joint tissues such as stiffness, porosity and ability to recover adequately after deformation to facilitate physiological function. Two decellularization protocols, namely: Trypsin and Triton X-100 were evaluated against their effects on bovine articular cartilage, using biomechanical, biochemical and microstructural techniques. These analyses revealed that decellularization with trypsin resulted in severe loss of mechanical stiffness including deleterious collapse of the collagen architecture which in turn significantly compromised the porosity of the construct. In contrast, triton X-100 detergent treatment yielded samples that retain mechanical stiffness relative to that of the normal intact cartilage sample, but the resulting construct contained ruminant cellular constituents. We conclude that both of these common decellularization protocols are inadequate for producing constructs that can serve as effective replacement and scaffolds to regenerate articular joint tissue.

Feasibility of autologous bone marrow mesenchymal stem cell-derived extracellular matrix scaffold for cartilage tissue engineering

Extracellular matrix (ECM) materials are widely used in cartilage tissue engineering. However, the current ECM materials are unsatisfactory for clinical practice as most of them are derived from allogenous or xenogenous tissue. This study was designed to develop a novel autologous ECM scaffold for cartilage tissue engineering. The autologous bone marrow mesenchymal stem cell-derived ECM (aBMSC-dECM) membrane was collected and fabricated into a three-dimensional porous scaffold via cross-linking and freeze-drying techniques. Articular chondrocytes were seeded into the aBMSC-dECM scaffold and atelocollagen scaffold, respectively. An in vitro culture and an in vivo implantation in nude mice model were performed to evaluate the influence on engineered cartilage. The current results showed that the aBMSC-dECM scaffold had a good microstructure and biocompatibility. After 4 weeks in vitro culture, the engineered cartilage in the aBMSC-dECM scaffold group formed thicker cartilage tissue with more homogeneous structure and higher expressions of cartilaginous gene and protein compared with the atelocollagen scaffold group. Furthermore, the engineered cartilage based on the aBMSC-dECM scaffold showed better cartilage formation in terms of volume and homogeneity, cartilage matrix content, and compressive modulus after 3 weeks in vivo implantation. These results indicated that the aBMSC-dECM scaffold could be a successful novel candidate scaffold for cartilage tissue engineering.

Discussion of “Generalized estimating equations: Notes on the choice of the working correlation matrix”

Objective To discuss generalized estimating equations as an extension of generalized linear models by commenting on the paper of Ziegler and Vens "Generalized Estimating Equations. Notes on the Choice of the Working Correlation Matrix". Methods Inviting an international group of experts to comment on this paper. Results Several perspectives have been taken by the discussants. Econometricians have established parallels to the generalized method of moments (GMM). Statisticians discussed model assumptions and the aspect of missing data Applied statisticians; commented on practical aspects in data analysis. Conclusions In general, careful modeling correlation is encouraged when considering estimation efficiency and other implications, and a comparison of choosing instruments in GMM and generalized estimating equations, (GEE) would be worthwhile. Some theoretical drawbacks of GEE need to be further addressed and require careful analysis of data This particularly applies to the situation when data are missing at random.

Applicability of pebble matrix filtration for the pre-treatment of surface waters containing high turbidity and NOM

Purification of drinking water is routinely achieved by use of conventional coagulants and disinfection procedures. However, there are instances such as flood events when the level of turbidity reaches extreme levels while NOM may be an issue throughout the year. Consequently, there is a need to develop technologies which can effectively treat water of high turbidity during flood events and natural organic matter (NOM) content year round. It was our hypothesis that pebble matrix filtration potentially offered a relatively cheap, simple and reliable means to clarify such challenging water samples. Therefore, a laboratory scale pebble matrix filter (PMF) column was used to evaluate the turbidity and natural organic matter (NOM) pre-treatment performance in relation to 2013 Brisbane River flood water. Since the high turbidity was only a seasonal and short term problem, the general applicability of pebble matrix filters for NOM removal was also investigated. A 1.0 m deep bed of pebbles (the matrix) partly in-filled with either sand or crushed glass was tested, upon which was situated a layer of granular activated carbon (GAC). Turbidity was measured as a surrogate for suspended solids (SS), whereas, total organic carbon (TOC) and UV Absorbance at 254 nm were measured as surrogate parameters for NOM. Experiments using natural flood water showed that without the addition of any chemical coagulants, PMF columns achieved at least 50% turbidity reduction when the source water contained moderate hardness levels. For harder water samples, above 85% turbidity reduction was obtained. The ability to remove 50% turbidity without chemical coagulants may represent significant cost savings to water treatment plants and added environmental benefits accrue due to less sludge formation. A TOC reduction of 35-47% and UV-254 nm reduction of 24-38% was also observed. In addition to turbidity removal during flood periods, the ability to remove NOM using the pebble matrix filter throughout the year may have the benefit of reducing disinfection by-products (DBP) formation potential and coagulant demand at water treatment plants. Final head losses were remarkably low, reaching only 11 cm at a filtration velocity of 0.70 m/h.

The politics of foreign aid: A median voter perspective

The paper presents a model where the median voter in the donor country determines the support of foreign aid. It is first established that an individual in the donor country is affected by the direct benefits (due to altruism) and costs (due to taxes) of giving aid, and by the indirect benefits or costs of a change in the terms of trade. Then it is shown that the latter effect works through changing both the donor country's average income and its distribution of income. Given the stylized facts of a capital-abundant donor country and relatively capital-poor median voter, it is shown how redistribution-of-income effects soften the impact of terms-of-trade changes on the political support for foreign aid.

3D extracellular matrix interactions modulate tumour cell growth, invasion and angiogenesis in engineered tumour microenvironments

Interactions between tumour cells and extracellular matrix proteins of the tumour microenvironment play crucial roles in cancer progression. So far, however, there are only a few experimental platforms available that allow us to study these interactions systematically in a mechanically defined three-dimensional (3D) context. Here, we have studied the effect of integrin binding motifs found within common extracellular matrix (ECM) proteins on 3D breast (MCF-7) and prostate (PC-3, LNCaP) cancer cell cultures, and co-cultures with endothelial and mesenchymal stromal cells. For this purpose, matrix metalloproteinase-degradable biohybrid poly(ethylene) glycol-heparin hydrogels were decorated with the peptide motifs RGD, GFOGER (collagen I), or IKVAV (laminin-111). Over 14 days, cancer spheroids of 100-200µm formed. While the morphology of poorly invasive MCF-7 and LNCaP cells was not modulated by any of the peptide motifs, the aggressive PC-3 cells exhibited an invasive morphology when cultured in hydrogels comprising IKVAV and GFOGER motifs compared to RGD motifs or nonfunctionalised controls. PC-3 (but not MCF-7 and LNCaP) cell growth and endothelial cell infiltration were also significantly enhanced in IKVAV and GFOGER presenting gels. Taken together, we have established a 3D culture model that allows for dissecting the effect of biochemical cues on processes relevant to early cancer progression. These findings provide a basis for more mechanistic studies that may further advance our understanding of how ECM modulates cancer cell invasion and how to ultimately interfere with this process.