Effects of Interspecific Gene Flow on the Phenotypic Variance-Covariance Matrix in Lake Victoria Cichlids

Quantitative genetics theory predicts adaptive evolution to be constrained along evolutionary lines of least resistance. In theory, hybridization and subsequent interspecific gene flow may however rapidly change the evolutionary constraints of a population and eventually change its evolutionary potential, but empirical evidence is still scarce. Using closely related species pairs of Lake Victoria cichlids sampled from four different islands with different levels of interspecific gene flow, we tested for potential effects of introgressive hybridization on phenotypic evolution in wild populations. We found that these effects differed among our study species. Constraints measured as the eccentricity of phenotypic variance-covariance matrices declined significantly with increasing gene flow in the less abundant species for matrices that have a diverged line of least resistance. In contrast we find no such decline for the more abundant species. Overall our results suggest that hybridization can change the underlying phenotypic variance-covariance matrix, potentially increasing the adaptive potential of such populations.

DEVCON: Stata module to apply the deviation contrast transform to estimation results

devcon transforms the coefficients of 0/1 dummy variables so that they reflect deviations from the "grand mean" rather than deviations from the reference category (the transformed coefficients are equivalent to those obtained by the so called "effects coding") and adds the coefficient for the reference category. The variance-covariance matrix of the estimates is transformed accordingly. The transformed estimated can be used with post estimation procedures. In particular, devcon can be used to solve the identification problem for dummy variable effects in the so-called Blinder-Oaxaca decomposition (see the oaxaca package).

Brain areas involved in medial temporal lobe seizures: a principal component analysis of ictal SPECT data

The study describes brain areas involved in medial temporal lobe (mTL) seizures of 12 patients. All patients showed so-called oro-alimentary behavior within the first 20 s of clinical seizure manifestation characteristic of mTL seizures. Single photon emission computed tomography (SPECT) images of regional cerebral blood flow (rCBF) were acquired from the patients in ictal and interictal phases and from normal volunteers. Image analysis employed categorical comparisons with statistical parametric mapping and principal component analysis (PCA) to assess functional connectivity. PCA supplemented the findings of the categorical analysis by decomposing the covariance matrix containing images of patients and healthy subjects into distinct component images of independent variance, including areas not identified by the categorical analysis. Two principal components (PCs) discriminated the subject groups: patients with right or left mTL seizures and normal volunteers, indicating distinct neuronal networks implicated by the seizure. Both PCs were correlated with seizure duration, one positively and the other negatively, confirming their physiological significance. The independence of the two PCs yielded a clear clustering of subject groups. The local pattern within the temporal lobe describes critical relay nodes which are the counterpart of oro-alimentary behavior: (1) right mesial temporal zone and ipsilateral anterior insula in right mTL seizures, and (2) temporal poles on both sides that are densely interconnected by the anterior commissure. Regions remote from the temporal lobe may be related to seizure propagation and include positively and negatively loaded areas. These patterns, the covarying areas of the temporal pole and occipito-basal visual association cortices, for example, are related to known anatomic paths.

Morphological Diversity and the Roles of Contingency, Chance and Determinism in African Cichlid Radiations

Background: Deterministic evolution, phylogenetic contingency and evolutionary chance each can influence patterns of morphological diversification during adaptive radiation. In comparative studies of replicate radiations, convergence in a common morphospace implicates determinism, whereas non-convergence suggests the importance of contingency or chance. Methodology/Principal Findings: The endemic cichlid fish assemblages of the three African great lakes have evolved similar sets of ecomorphs but show evidence of non-convergence when compared in a common morphospace, suggesting the importance of contingency and/or chance. We then analyzed the morphological diversity of each assemblage independently and compared their axes of diversification in the unconstrained global morphospace. We find that despite differences in phylogenetic composition, invasion history, and ecological setting, the three assemblages are diversifying along parallel axes through morphospace and have nearly identical variance-covariance structures among morphological elements. Conclusions/Significance: By demonstrating that replicate adaptive radiations are diverging along parallel axes, we have shown that non-convergence in the common morphospace is associated with convergence in the global morphospace. Applying these complimentary analyses to future comparative studies will improve our understanding of the relationship between morphological convergence and non-convergence, and the roles of contingency, chance and determinism in driving morphological diversification.

Estimating Particle Shape and Orientation Using Volume Tensors

We develop statistical procedures for estimating shape and orientation of arbitrary three-dimensional particles. We focus on the case where particles cannot be observed directly, but only via sections. Volume tensors are used for describing particle shape and orientation, and we derive stereological estimators of the tensors. These estimators are combined to provide consistent estimators of the moments of the so-called particle cover density. The covariance structure associated with the particle cover density depends on the orientation and shape of the particles. For instance, if the distribution of the typical particle is invariant under rotations, then the covariance matrix is proportional to the identity matrix. We develop a non-parametric test for such isotropy. A flexible Lévy-based particle model is proposed, which may be analysed using a generalized method of moments in which the volume tensors enter. The developed methods are used to study the cell organization in the human brain cortex.

Longitudinal evaluation of cartilage composition of matrix-associated autologous chondrocyte transplants with 3-T delayed gadolinium-enhanced MRI of cartilage

OBJECTIVE: The purposes of this study were to use delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) to evaluate the zonal distribution of glycosaminoglycans (GAGs) in normal cartilage and repair tissue and to use 3-T MRI to monitor the GAG content in matrix-associated autologous chondrocyte transplants. SUBJECTS AND METHODS: Fifteen patients who underwent matrix-associated autologous chondrocyte transplantation in the knee joint underwent MRI at baseline and 3-T follow-up MRI 1 year later. Total and zonal changes in longitudinal relaxivity (deltaR1) and relative deltaR1 were calculated for repair tissue and normal hyaline cartilage and compared by use of analysis of variance. RESULTS: There was a significant difference between the mean deltaR1 of repair tissue and that of reference cartilage at baseline and follow-up (p < 0.001). There was a significant increase in deltaR1 value and a decrease in GAG content from the deep layer to the superficial layer in the reference cartilage and almost no variation and significantly higher values for the repair tissue at both examinations. At 1-year follow-up imaging, there was a 22.7% decrease in deltaR1 value in the deep zone of the transplant. CONCLUSION: T1 mapping with dGEMRIC at 3 T shows the zonal structure of normal hyaline cartilage, highly reduced zonal variations in repair tissue, and a tendency toward an increase in global and zonal GAG content 1 year after transplantation.

Quantitative T2 mapping during follow-up after matrix-associated autologous chondrocyte transplantation (MACT): full-thickness and zonal evaluation to visualize the maturation of cartilage repair tissue

The purpose of this article was to evaluate the potential of in vivo zonal T2-mapping as a noninvasive tool in the longitudinal visualization of cartilage repair tissue maturation after matrix-associated autologous chondrocyte transplantation (MACT). Fifteen patients were treated with MACT and evaluated cross-sectionally, with a baseline MRI at a follow-up of 19.7 +/- 12.1 months after cartilage transplantation surgery of the knee. In the same 15 patients, 12 months later (31.7 +/- 12.0 months after surgery), a longitudinal 1-year follow-up MRI was obtained. MRI was performed on a 3 Tesla MR scanner; morphological evaluation was performed using a double-echo steady-state sequence; T2 maps were calculated from a multiecho, spin-echo sequence. Quantitative mean (full-thickness) and zonal (deep and superficial) T2 values were calculated in the cartilage repair area and in control cartilage sites. A statistical analysis of variance was performed. Full-tickness T2 values showed no significant difference between sites of healthy cartilage and cartilage repair tissue (p < 0.05). Using zonal T2 evaluation, healthy cartilage showed a significant increase from the deep to superficial cartilage layers (p < 0.05). Cartilage repair tissue after MACT showed no significant zonal increase from deep to superficial cartilage areas during baseline MRI (p > 0.05); however, during the 1-year follow-up, a significant zonal stratification could be observed (p < 0.05). Morphological evaluation showed no significant difference between the baseline and the 1-year follow-up MRI. T2 mapping seems to be more sensitive in revealing changes in the repair tissue compared to morphological MRI. In vivo zonal T2 assessment may be sensitive enough to characterize the maturation of cartilage repair tissue.

Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

Cell and matrix morphology in articular cartilage from adult human knee and ankle joints suggests depth-associated adaptations to biomechanical and anatomical roles

OBJECTIVE Marked differences exist between human knee and ankle joints regarding risks and progression of osteoarthritis (OA). Pathomechanisms of degenerative joint disease may therefore differ in these joints, due to differences in tissue structure and function. Focussing on structural issues which are design goals for tissue engineering, we compared cell and matrix morphologies in different anatomical sites of adult human knee and ankle joints. METHODS Osteochondral explants were acquired from knee and ankle joints of deceased persons aged 20 to 40 years and analyzed for cell, matrix and tissue morphology using confocal and electron microscopy and unbiased stereological methods. Variations associated with joint (knee versus ankle) and biomechanical role (convex versus concave articular surfaces) were identified by 2-way analysis of variance and post-hoc analysis. RESULTS Knee cartilage exhibited higher cell densities in the superficial zone than ankle cartilage. In the transitional zone, higher cell densities were observed in association with convex versus concave articular surfaces, without significant differences between knee and ankle cartilage. Highly uniform cell and matrix morphologies were evident throughout the radial zone in the knee and ankle, regardless of tissue biomechanical role. Throughout the knee and ankle cartilage sampled, chondron density was remarkably constant at approximately 4.2×10(6) chondrons/cm(3). CONCLUSION Variation of cartilage cell and matrix morphologies with changing joint and biomechanical environments suggests that tissue structural adaptations are performed primarily by the superficial and transitional zones. Data may aid the development of site-specific cartilage tissue engineering, and help identify conditions where OA is likely to occur.

A Thalamic-Fronto-Parietal Structural Covariance Network Emerging in the Course of Recovery from Hand Paresis after Ischemic Stroke.

AIM To describe structural covariance networks of gray matter volume (GMV) change in 28 patients with first-ever stroke to the primary sensorimotor cortices, and to investigate their relationship to hand function recovery and local GMV change. METHODS Tensor-based morphometry maps derived from high-resolution structural images were subject to principal component analyses to identify the networks. We calculated correlations between network expression and local GMV change, sensorimotor hand function and lesion volume. To verify which of the structural covariance networks of GMV change have a significant relationship to hand function, we performed an additional multivariate regression approach. RESULTS Expression of the second network, explaining 9.1% of variance, correlated with GMV increase in the medio-dorsal (md) thalamus and hand motor skill. Patients with positive expression coefficients were distinguished by significantly higher GMV increase of this structure during stroke recovery. Significant nodes of this network were located in md thalamus, dorsolateral prefrontal cortex, and higher order sensorimotor cortices. Parameter of hand function had a unique relationship to the network and depended on an interaction between network expression and lesion volume. Inversely, network expression is limited in patients with large lesion volumes. CONCLUSION Chronic phase of sensorimotor cortical stroke has been characterized by a large scale co-varying structural network in the ipsilesional hemisphere associated specifically with sensorimotor hand skill. Its expression is related to GMV increase of md thalamus, one constituent of the network, and correlated with the cortico-striato-thalamic loop involved in control of motor execution and higher order sensorimotor cortices. A close relation between expression of this network with degree of recovery might indicate reduced compensatory resources in the impaired subgroup.