Determining the effective dimensionality of the genetic variance-covariance matrix

Determining the dimensionality of G provides an important perspective on the genetic basis of a multivariate suite of traits. Since the introduction of Fisher's geometric model, the number of genetically independent traits underlying a set of functionally related phenotypic traits has been recognized as an important factor influencing the response to selection. Here, we show how the effective dimensionality of G can be established, using a method for the determination of the dimensionality of the effect space from a multivariate general linear model introduced by AMEMIYA (1985). We compare this approach with two other available methods, factor-analytic modeling and bootstrapping, using a half-sib experiment that estimated G for eight cuticular hydrocarbons of Drosophila serrata. In our example, eight pheromone traits were shown to be adequately represented by only two underlying genetic dimensions by Amemiya's approach and factor-analytic modeling of the covariance structure at the sire level. In, contrast, bootstrapping identified four dimensions with significant genetic variance. A simulation study indicated that while the performance of Amemiya's method was more sensitive to power constraints, it performed as well or better than factor-analytic modeling in correctly identifying the original genetic dimensions at moderate to high levels of heritability. The bootstrap approach consistently overestimated the number of dimensions in all cases and performed less well than Amemiya's method at subspace recovery.

Arthritic disease suppression and cartilage protection with glycosaminoglycan polypeptide complexes (Peptacans) derived from the cartilage extracellular matrix: A novel approach to therapy

Molecular fragments of cartilage are antigenic and can stimulate an autoimmune response. Oral administration of type II collagen prevents disease onset in animal models of arthritis but the effects of other matrix components have not been reported. We evaluated glycosaminoglycan polypeptides (GAG-P) and matrix proteins (CaP) from cartilage for a) mitigating disease activity in rats with collagen-induced arthritis (CIA) and adjuvant-induced arthritis (AIA) and b) stimulating proteoglycan (PG) synthesis by chondrocytes in-vitro. CIA and AIA were established in Wistar rats using standard methods. Agents were administered orally (10–200 mg/kg), either for seven days prior to disease induction (toleragenic protocol), or continuously for 15 days after injecting the arthritigen (prophylactic protocol). Joint swelling and arthritis scores were determined on day 15. Histological sections of joint tissues were assessed post-necropsy. In chondrocyte cultures, CaP + / − interleukin-1 stimulated PG biosynthesis. CaP was also active in preventing arthritis onset at 3.3, 10 or 20 mg/kg in the rat CIA model using the toleragenic protocol. It was only active at 20 and 200 mg/kg in the CIA prophylactic protocol. GAG-P was active in the CIA toleragenic protocol at 20 mg/kg but chondroitin sulfate and glucosamine hydrochloride or glucosamine sulfate were all inactive. The efficacy of CaP in the rat AIA model was less than in the CIA model. These findings lead us to suggest that oral CaP could be used as a disease-modifying anti-arthritic drug.

A novel extracellular role for tissue transglutaminase in matrix-bound VEGF-mediated angiogenesis

The importance of tissue transglutaminase (TG2) in angiogenesis is unclear and contradictory. Here we show that inhibition of extracellular TG2 protein crosslinking or downregulation of TG2 expression leads to inhibition of angiogenesis in cell culture, the aorta ring assay and in vivo models. In a human umbilical vein endothelial cell (HUVEC) co-culture model, inhibition of extracellular TG2 activity can halt the progression of angiogenesis, even when introduced after tubule formation has commenced and after addition of excess vascular endothelial growth factor (VEGF). In both cases, this leads to a significant reduction in tubule branching. Knockdown of TG2 by short hairpin (shRNA) results in inhibition of HUVEC migration and tubule formation, which can be restored by add back of wt TG2, but not by the transamidation-defective but GTP-binding mutant W241A. TG2 inhibition results in inhibition of fibronectin deposition in HUVEC monocultures with a parallel reduction in matrix-bound VEGFA, leading to a reduction in phosphorylated VEGF receptor 2 (VEGFR2) at Tyr1214 and its downstream effectors Akt and ERK1/2, and importantly its association with b1 integrin. We propose a mechanism for the involvement of matrix-bound VEGFA in angiogenesis that is dependent on extracellular TG2-related activity. © 2013 Macmillan Publishers Limited. All rights reserved.

Assessment of different formulations of oral Mycobacterium bovis Bacille Calmette-Guerin (BCG) vaccine in rodent models for immunogenicity and protection against aerosol challenge with M-bovis

Bovine tuberculosis (bTB) caused by infection with Mycobacterium bovis is causing considerable economic loss to farmers and Government in the United Kingdom as its incidence is increasing. Efforts to control bTB in the UK are hampered by the infection in Eurasian badgers (Metes metes) that represent a wildlife reservoir and source of recurrent M. bovis exposure to cattle. Vaccination of badgers with the human TB vaccine, M. bovis Bacille Calmette-Guerin (BCG), in oral bait represents a possible disease control tool and holds the best prospect for reaching badger populations over a wide geographical area. Using mouse and guinea pig models, we evaluated the immunogenicity and protective efficacy, respectively, of candidate badger oral vaccines based on formulation of BCG in lipid matrix, alginate beads, or a novel microcapsular hybrid of both lipid and alginate. Two different oral doses of BCG were evaluated in each formulation for their protective efficacy in guinea pigs, while a single dose was evaluated in mice. In mice, significant immune responses (based on lymphocyte proliferation and expression of IFN-gamma) were only seen with the lipid matrix and the lipid in alginate microcapsular formulation, corresponding to the isolation of viable BCG from alimentary tract lymph nodes. In guinea pigs, only BCG formulated in lipid matrix conferred protection to the spleen and lungs following aerosol route challenge with M. bovis. Protection was seen with delivery doses in the range 10(6)-10(7) CFU, although this was more consistent in the spleen at the higher dose. No protection in terms of organ CFU was seen with BCG administered in alginate beads or in lipid in alginate microcapsules, although 10(7) in the latter formulation conferred protection in terms of increasing body weight after challenge and a smaller lung to body weight ratio at necropsy. These results highlight the potential for lipid, rather than alginate, -based vaccine formulations as suitable delivery vehicles for an oral BCG vaccine in badgers.

Some Properties of Mittag-Leffler Functions and Matrix-Variate Analogues: A Statistical Perspective

Mathematical Subject Classification 2010:26A33, 33E99, 15A52, 62E15.

Matrix-Variate Statistical Distributions and Fractional Calculus

MSC 2010: 15A15, 15A52, 33C60, 33E12, 44A20, 62E15 Dedicated to Professor R. Gorenflo on the occasion of his 80th birthday

An improved framework for dynamic origin-destination (O-D) matrix estimation

This dissertation aims to improve the performance of existing assignment-based dynamic origin-destination (O-D) matrix estimation models to successfully apply Intelligent Transportation Systems (ITS) strategies for the purposes of traffic congestion relief and dynamic traffic assignment (DTA) in transportation network modeling. The methodology framework has two advantages over the existing assignment-based dynamic O-D matrix estimation models. First, it combines an initial O-D estimation model into the estimation process to provide a high confidence level of initial input for the dynamic O-D estimation model, which has the potential to improve the final estimation results and reduce the associated computation time. Second, the proposed methodology framework can automatically convert traffic volume deviation to traffic density deviation in the objective function under congested traffic conditions. Traffic density is a better indicator for traffic demand than traffic volume under congested traffic condition, thus the conversion can contribute to improving the estimation performance. The proposed method indicates a better performance than a typical assignment-based estimation model (Zhou et al., 2003) in several case studies. In the case study for I-95 in Miami-Dade County, Florida, the proposed method produces a good result in seven iterations, with a root mean square percentage error (RMSPE) of 0.010 for traffic volume and a RMSPE of 0.283 for speed. In contrast, Zhou's model requires 50 iterations to obtain a RMSPE of 0.023 for volume and a RMSPE of 0.285 for speed. In the case study for Jacksonville, Florida, the proposed method reaches a convergent solution in 16 iterations with a RMSPE of 0.045 for volume and a RMSPE of 0.110 for speed, while Zhou's model needs 10 iterations to obtain the best solution, with a RMSPE of 0.168 for volume and a RMSPE of 0.179 for speed. The successful application of the proposed methodology framework to real road networks demonstrates its ability to provide results both with satisfactory accuracy and within a reasonable time, thus establishing its potential usefulness to support dynamic traffic assignment modeling, ITS systems, and other strategies.

Diatom-based Models for Inferring Hydrology and Periphyton Abundance in a Subtropical Karstic Wetland: Implications for Ecosystem-Scale Bioassessment

We developed diatom-based prediction models of hydrology and periphyton abundance to inform assessment tools for a hydrologically managed wetland. Because hydrology is an important driver of ecosystem change, hydrologic alterations by restoration efforts could modify biological responses, such as periphyton characteristics. In karstic wetlands, diatoms are particularly important components of mat-forming calcareous periphyton assemblages that both respond and contribute to the structural organization and function of the periphyton matrix. We examined the distribution of diatoms across the Florida Everglades landscape and found hydroperiod and periphyton biovolume were strongly correlated with assemblage composition. We present species optima and tolerances for hydroperiod and periphyton biovolume, for use in interpreting the directionality of change in these important variables. Predictions of these variables were mapped to visualize landscape-scale spatial patterns in a dominant driver of change in this ecosystem (hydroperiod) and an ecosystem-level response metric of hydrologic change (periphyton biovolume). Specific diatom assemblages inhabiting periphyton mats of differing abundance can be used to infer past conditions and inform management decisions based on how assemblages are changing. This study captures diatom responses to wide gradients of hydrology and periphyton characteristics to inform ecosystem-scale bioassessment efforts in a large wetland.

Realized Matrix-Exponential Stochastic Volatility with Asymmetry, Long Memory and Spillovers

The paper develops a novel realized matrix-exponential stochastic volatility model of multivariate returns and realized covariances that incorporates asymmetry and long memory (hereafter the RMESV-ALM model). The matrix exponential transformation guarantees the positivedefiniteness of the dynamic covariance matrix. The contribution of the paper ties in with Robert Basmann’s seminal work in terms of the estimation of highly non-linear model specifications (“Causality tests and observationally equivalent representations of econometric models”, Journal of Econometrics, 1988, 39(1-2), 69–104), especially for developing tests for leverage and spillover effects in the covariance dynamics. Efficient importance sampling is used to maximize the likelihood function of RMESV-ALM, and the finite sample properties of the quasi-maximum likelihood estimator of the parameters are analysed. Using high frequency data for three US financial assets, the new model is estimated and evaluated. The forecasting performance of the new model is compared with a novel dynamic realized matrix-exponential conditional covariance model. The volatility and co-volatility spillovers are examined via the news impact curves and the impulse response functions from returns to volatility and co-volatility.

Linear Hyperspectral Unmixing Using L0-norm Approximations and Nonnegative Matrix Factorization

Spectral unmixing (SU) is a technique to characterize mixed pixels of the hyperspectral images measured by remote sensors. Most of the existing spectral unmixing algorithms are developed using the linear mixing models. Since the number of endmembers/materials present at each mixed pixel is normally scanty compared with the number of total endmembers (the dimension of spectral library), the problem becomes sparse. This thesis introduces sparse hyperspectral unmixing methods for the linear mixing model through two different scenarios. In the first scenario, the library of spectral signatures is assumed to be known and the main problem is to find the minimum number of endmembers under a reasonable small approximation error. Mathematically, the corresponding problem is called the $\ell_0$-norm problem which is NP-hard problem. Our main study for the first part of thesis is to find more accurate and reliable approximations of $\ell_0$-norm term and propose sparse unmixing methods via such approximations. The resulting methods are shown considerable improvements to reconstruct the fractional abundances of endmembers in comparison with state-of-the-art methods such as having lower reconstruction errors. In the second part of the thesis, the first scenario (i.e., dictionary-aided semiblind unmixing scheme) will be generalized as the blind unmixing scenario that the library of spectral signatures is also estimated. We apply the nonnegative matrix factorization (NMF) method for proposing new unmixing methods due to its noticeable supports such as considering the nonnegativity constraints of two decomposed matrices. Furthermore, we introduce new cost functions through some statistical and physical features of spectral signatures of materials (SSoM) and hyperspectral pixels such as the collaborative property of hyperspectral pixels and the mathematical representation of the concentrated energy of SSoM for the first few subbands. Finally, we introduce sparse unmixing methods for the blind scenario and evaluate the efficiency of the proposed methods via simulations over synthetic and real hyperspectral data sets. The results illustrate considerable enhancements to estimate the spectral library of materials and their fractional abundances such as smaller values of spectral angle distance (SAD) and abundance angle distance (AAD) as well.

On Queueinginventory Models – product form solution; reservation, Cancellation and common life time

Queueing theory is the mathematical study of ‘queue’ or ‘waiting lines’ where an item from inventory is provided to the customer on completion of service. A typical queueing system consists of a queue and a server. Customers arrive in the system from outside and join the queue in a certain way. The server picks up customers and serves them according to certain service discipline. Customers leave the system immediately after their service is completed. For queueing systems, queue length, waiting time and busy period are of primary interest to applications. The theory permits the derivation and calculation of several performance measures including the average waiting time in the queue or the system, mean queue length, traffic intensity, the expected number waiting or receiving service, mean busy period, distribution of queue length, and the probability of encountering the system in certain states, such as empty, full, having an available server or having to wait a certain time to be served.

The third dimension: new developments in cell culture models for colorectal research

Cellular models are important tools in various research areas related to colorectal biology and associated diseases. Herein, we review the most widely used cell lines and the different techniques to grow them, either as cell monolayer, polarized two-dimensional epithelia on membrane filters, or as three-dimensional spheres in scaffoldfree or matrix-supported culture conditions. Moreover, recent developments, such as gut-on-chip devices or the ex vivo growth of biopsy-derived organoids, are also discussed. We provide an overview on the potential applications but also on the limitations for each of these techniques, while evaluating their contribution to provide more reliable cellular models for research, diagnostic testing, or pharmacological validation related to colon physiology and pathophysiology.

Animal Models of Friedreich Ataxia

Friedreich ataxia (FRDA) is the most common form of autosomal-recessive ataxia. Common nonmotor features include cardiomyopathy and diabetes mellitus. At present, no effective treatments are available to prevent disease progression. Age of onset varies from infancy to adulthood. In the majority of patients, FRDA is caused by intronic GAA expansions in FXN, which encodes a highly-conserved small mitochondrial matrix protein, frataxin. A mouse model of FRDA has been difficult to generate because complete loss of frataxin causes early embryonic lethality. Although there are some controversies about the function of frataxin, recent biochemical and structural studies have confirmed that it is a component of the multiprotein complex that assembles iron-sulfur clusters in the mitochondrial matrix. The main consequences of frataxin deficiency are energy deficit, altered iron metabolism, and oxidative damage.

A new method for conditioning stochastic groundwater flow models in fractured media

Many geological formations consist of crystalline rocks that have very low matrix permeability but allow flow through an interconnected network of fractures. Understanding the flow of groundwater through such rocks is important in considering disposal of radioactive waste in underground repositories. A specific area of interest is the conditioning of fracture transmissivities on measured values of pressure in these formations. This is the process where the values of fracture transmissivities in a model are adjusted to obtain a good fit of the calculated pressures to measured pressure values. While there are existing methods to condition transmissivity fields on transmissivity, pressure and flow measurements for a continuous porous medium there is little literature on conditioning fracture networks. Conditioning fracture transmissivities on pressure or flow values is a complex problem because the measurements are not linearly related to the fracture transmissivities and they are also dependent on all the fracture transmissivities in the network. We present a new method for conditioning fracture transmissivities on measured pressure values based on the calculation of certain basis vectors; each basis vector represents the change to the log transmissivity of the fractures in the network that results in a unit increase in the pressure at one measurement point whilst keeping the pressure at the remaining measurement points constant. The fracture transmissivities are updated by adding a linear combination of basis vectors and coefficients, where the coefficients are obtained by minimizing an error function. A mathematical summary of the method is given. This algorithm is implemented in the existing finite element code ConnectFlow developed and marketed by Serco Technical Services, which models groundwater flow in a fracture network. Results of the conditioning are shown for a number of simple test problems as well as for a realistic large scale test case.