Covariance structure in the skull of Catarrhini: a case of pattern stasis and magnitude evolution

The study of the genetic variance/covariance matrix (G-matrix) is a recent and fruitful approach in evolutionary biology, providing a window of investigating for the evolution of complex characters. Although G-matrix studies were originally conducted for microevolutionary timescales, they could be extrapolated to macroevolution as long as the G-matrix remains relatively constant, or proportional, along the period of interest. A promising approach to investigating the constancy of G-matrices is to compare their phenotypic counterparts (P-matrices) in a large group of related species; if significant similarity is found among several taxa, it is very likely that the underlying G-matrices are also equivalent. Here we study the similarity of covariance and correlation structure in a broad sample of Old World monkeys and apes (Catarrhini). We made phylogenetically structured comparisons of correlation and covariance matrices derived from 39 skull traits, ranging from between species to the superfamily level. We also compared the overall magnitude of integration between skull traits (r(2)) for all Catarrhim genera. Our results show that P-matrices were not strictly constant among catarrhines, but the amount of divergence observed among taxa was generally low. There was significant and positive correlation between the amount of divergence in correlation and covariance patterns among the 30 genera and their phylogenetic distances derived from a recently proposed phylogenetic hypothesis. Our data demonstrate that the P-matrices remained relatively similar along the evolutionary history of catarrhines, and comparisons with the G-matrix available for a New World monkey genus (Saguinus) suggests that the same holds for all anthropoids. The magnitude of integration, in contrast, varied considerably among genera, indicating that evolution of the magnitude, rather than the pattern of inter-trait correlations, might have played an important role in the diversification of the catarrhine skull. (C) 2009 Elsevier Ltd. All rights reserved.

The Evolution of Modularity in the Mammalian Skull I: Morphological Integration Patterns and Magnitudes

Morphological integration refers to the modular structuring of inter-trait relationships in an organism, which could bias the direction and rate of morphological change, either constraining or facilitating evolution along certain dimensions of the morphospace. Therefore, the description of patterns and magnitudes of morphological integration and the analysis of their evolutionary consequences are central to understand the evolution of complex traits. Here we analyze morphological integration in the skull of several mammalian orders, addressing the following questions: are there common patterns of inter-trait relationships? Are these patterns compatible with hypotheses based on shared development and function? Do morphological integration patterns and magnitudes vary in the same way across groups? We digitized more than 3,500 specimens spanning 15 mammalian orders, estimated the correspondent pooled within-group correlation and variance/covariance matrices for 35 skull traits and compared those matrices among the orders. We also compared observed patterns of integration to theoretical expectations based on common development and function. Our results point to a largely shared pattern of inter-trait correlations, implying that mammalian skull diversity has been produced upon a common covariance structure that remained similar for at least 65 million years. Comparisons with a rodent genetic variance/covariance matrix suggest that this broad similarity extends also to the genetic factors underlying phenotypic variation. In contrast to the relative constancy of inter-trait correlation/covariance patterns, magnitudes varied markedly across groups. Several morphological modules hypothesized from shared development and function were detected in the mammalian taxa studied. Our data provide evidence that mammalian skull evolution can be viewed as a history of inter-module parcellation, with the modules themselves being more clearly marked in those lineages with lower overall magnitude of integration. The implication of these findings is that the main evolutionary trend in the mammalian skull was one of decreasing the constraints to evolution by promoting a more modular architecture.

Genetic Polymorphisms of Matrix Metalloproteinases: Susceptibility and Prognostic Implications for Prostate Cancer

Purpose: Prostate cancer is the most common tumor in males in Brazil. Single nucleotide polymorphisms have been demonstrated to exist in the promoter regions of matrix metalloproteinase genes and they are associated with the development and progression of some cancers. We investigated the correlation between MMP1, 2, 7 and 9 polymorphisms with susceptibility to prostate cancer, and classic prognostic parameters of prostate cancer. Materials and Methods: Genomic DNA was extracted using conventional protocols. The DNA sequence containing the polymorphic site was amplified by realtime polymerase chain reaction using TaqMan (R) fluorescent probes. Results: For the MMP1 gene the polymorphic allele was more common in the control group than in the prostate cancer group (p <0.001). For the MMP9 gene the incidence of the polymorphic homozygote genotype was higher in the prostate cancer group (p <0.001). For higher stage tumors (pT3) a polymorphic allele in the MMP2 gene was more common (p = 0.026). When considering Gleason score, the polymorphic homozygote genotype of MMP9 was more common in Gleason 6 or less tumors (p = 0.003), while a polymorphic allele in the MMP2 gene was more common in Gleason 7 or greater tumors (p = 0.042). Conclusions: MMP1 and MMP2 may protect against prostate cancer development and MMP9 may be related to higher risk. In contrast, MMP9 polymorphism was associated with a lower Gleason score and MMP2 polymorphism was associated with nonorgan confined disease.

A genetic polymorphism of matrix metalloproteinase 9 (MMP-9) affects the changes in circulating MMP-9 levels induced by highly active antiretroviral therapy in HIV patients

We examined whether two functional polymorphisms (g.-1562C>T and g.-90(CA)14-24) in the matrix metalloproteinase (MMP)-9 gene or MMP-9 haplotypes affect the circulating levels of pro-MMP-9 and pro-MMP-9/TIMP-1 (tissue inhibitor of metalloproteinase-1) ratios in AIDS patients, and modulate alterations in these biomarkers after highly active antiretroviral therapy (HAART). We studied 82 patients commencing HAART. Higher pro-MMP-9 concentrations and pro-MMP-9/TIMP-1 ratios were found in CT/TT patients compared with CC patients. HAART decreased pro-MMP-9 levels and pro-MMP-9/TIMP-1 ratios in CT/TT patients, it did not modify pro-MMP-9 levels and it increased pro-MMP-9/TIMP-1 ratios in CC patients. The g.-90(CA)14-24 polymorphism, however, produced no significant effects. Moreover, we found no significant differences in HAART-induced changes in plasma pro-MMP-9, TIMP-1 and pro-MMP-9/TIMP-1 ratios when different MMP-9 haplotypes were compared. These findings suggest that the g.-1562C>T polymorphism affects pro-MMP-9 levels in patients with AIDS and modulates the alterations in pro-MMP-9 levels caused by HAART, thus possibly affecting the risk of cardiovascular complications. The Pharmacogenomics Journal (2009) 9, 265-273; doi: 10.1038/tpj.2009.13; published online 21 April 2009

Covariance matrix formula for Birnbaum-Saunders regression models

The Birnbaum-Saunders regression model is commonly used in reliability studies. We derive a simple matrix formula for second-order covariances of maximum-likelihood estimators in this class of models. The formula is quite suitable for computer implementation, since it involves only simple operations on matrices and vectors. Some simulation results show that the second-order covariances can be quite pronounced in small to moderate sample sizes. We also present empirical applications.

Comparison of different models to estimate genetic parameters for carcass traits in a commercial broiler line

The effect of genetic and non-genetic factors for carcass, breast meat and leg weights, and yields of a commercial broiler line were investigated using the restricted maximum likelihood method, considering four different animal models, including or excluding maternal genetic effect with covariance between direct and maternal genetic effects, and maternal permanent environmental effect. The likelihood ratio test was used to determine the most adequate model for each trait. For carcass, breast, and leg weight, and for carcass and breast yield, maternal genetic and permanent environmental effects as well as the covariance between direct and maternal genetic effects were significant. The estimates of direct and maternal heritability were 0.17 and 0.04 for carcass weight, 0.26 and 0.06 for breast weight, 0.22 and 0.02 for leg weight, 0.32 and 0.02 for carcass yield, and 0.52 and 0.04 for breast yield, respectively. For leg yield, maternal permanent environmental effect was important, in addition to direct genetic effects. For that trait, direct heritability and maternal permanent environmental variance as a proportion of the phenotypic variance were 0.43 and 0.02, respectively. The results indicate that ignoring maternal effects in the models, even though they were of small magnitude (0.02 to 0.06), tended to overestimate direct genetic variance and heritability for all traits.

Additive genetic relationships between scrotal circumference, heifer pregnancy, and stayability in Nellore cattle

Scrotal circumference data from 47,605 Nellore young bulls, measured at around 18 mo of age (SC18), were analyzed simultaneously with 27,924 heifer pregnancy (HP) and 80,831 stayability (STAY) records to estimate their additive genetic relationships. Additionally, the possibility that economically relevant traits measured directly in females could replace SC18 as a selection criterion was verified. Heifer pregnancy was defined as the observation that a heifer conceived and remained pregnant, which was assessed by rectal palpation at 60 d. Females were exposed to sires for the first time at about 14 mo of age (between 11 and 16 mo). Stayability was defined as whether or not a cow calved every year up to 5 yr of age, when the opportunity to breed was provided. A Bayesian linear-threshold-threshold analysis via Gibbs sampler was used to estimate the variance and covariance components of the multitrait model. Heritability estimates were 0.42 +/- 0.01, 0.53 +/- 0.03, and 0.10 +/- 0.01, for SC18, HP, and STAY, respectively. The genetic correlation estimates were 0.29 +/- 0.05, 0.19 +/- 0.05, and 0.64 +/- 0.07 between SC18 and HP, SC18 and STAY, and HP and STAY, respectively. The residual correlation estimate between HP and STAY was -0.08 +/- 0.03. The heritability values indicate the existence of considerable genetic variance for SC18 and HP traits. However, genetic correlations between SC18 and the female reproductive traits analyzed in the present study can only be considered moderate. The small residual correlation between HP and STAY suggests that environmental effects common to both traits are not major. The large heritability estimate for HP and the high genetic correlation between HP and STAY obtained in the present study confirm that EPD for HP can be used to select bulls for the production of precocious, fertile, and long-lived daughters. Moreover, SC18 could be incorporated in multitrait analysis to improve the prediction accuracy for HP genetic merit of young bulls.

SIZE AS A LINE OF LEAST RESISTANCE II: DIRECT SELECTION ON SIZE OR CORRELATED RESPONSE DUE TO CONSTRAINTS?

Evolutionary change in New World Monkey (NWM) skulls occurred primarily along the line of least resistance defined by size (including allometric) variation (g(max)). Although the direction of evolution was aligned with this axis, it was not clear whether this macroevolutionary pattern results from the conservation of within population genetic covariance patterns (long-term constraint) or long-term selection along a size dimension, or whether both, constraints and selection, were inextricably involved. Furthermore, G-matrix stability can also be a consequence of selection, which implies that both, constraints embodied in g(max) and evolutionary changes observed on the trait averages, would be influenced by selection Here, we describe a combination of approaches that allows one to test whether any particular instance of size evolution is a correlated by-product due to constraints (g(max)) or is due to direct selection on size and apply it to NWM lineages as a case study. The approach is based on comparing the direction and amount of evolutionary change produced by two different simulated sets of net-selection gradients (beta), a size (isometric and allometric size) and a nonsize set. Using this approach it is possible to distinguish between the two hypotheses (indirect size evolution due to constraints or direct selection on size), because although both may produce an evolutionary response aligned with g(max), the amount of change produced by random selection operating through the variance/covariance patterns (constraints hypothesis) will be much smaller than that produced by selection on size (selection hypothesis). Furthermore, the alignment of simulated evolutionary changes with g(max) when selection is not on size is not as tight as when selection is actually on size, allowing a statistical test of whether a particular observed case of evolution along the line of least resistance is the result of selection along it or not. Also, with matrix diagonalization (principal components [PC]) it is possible to calculate directly the net-selection gradient on size alone (first PC [PC1]) by dividing the amount of phenotypic difference between any two populations by the amount of variation in PC1, which allows one to benchmark whether selection was on size or not

The Poisson-exponential lifetime distribution

In this paper we proposed a new two-parameters lifetime distribution with increasing failure rate. The new distribution arises on a latent complementary risk problem base. The properties of the proposed distribution are discussed, including a formal proof of its probability density function and explicit algebraic formulae for its reliability and failure rate functions, quantiles and moments, including the mean and variance. A simple EM-type algorithm for iteratively computing maximum likelihood estimates is presented. The Fisher information matrix is derived analytically in order to obtaining the asymptotic covariance matrix. The methodology is illustrated on a real data set. (C) 2010 Elsevier B.V. All rights reserved.

Genetic variability in S(1) families from different maize populations.

Genetic variability in S(1) families from different maize populations. The objectives of the present work were directed towards the study of genetic: variablilty In seven maize populations with a broad genetic base, as a guide for population improvement. The field evaluation was conducted in completely randomized blocks, at one location (Anhembi, Sao Paulo state) with different groups, of S(1) families Obtained from seven populations (GO-D: dent type, GO-F: flint type, GO-L: long car, GO-G: thick Car; and composites G3, G4 and GO-S). Estimates were obtained for genetic variance (progeny mean basis), phenotypic variance of families means, and coefficient of heritability (broad sense) for progeny means. Estimates of heritability were high for Car weight (0,89 to 0.94), car length (0.77 to 0.88) and car diameter (0.77 to 0.92); and lower for plant height (0.58 to 0.80) and Car height (0.54 to 0.84), thus showing the high Potential of the populations for recurrent selection based oil S, families. Ear yield in the base populations used as controls varied front 11,200 kg ha(-1) (GO-D) to 12,800 kg ha(-1) (G3). The means of S(1) families varied from 6,070 kg ha(-1) (GO-F) to 7,380 kg ha(-1) (G4); the Inbreeding depression in S(1) Families varied front 37.5% (G4) to 48.0% (G3) relative to the non-inbred population.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Skull Modularity in Neotropical Marsupials and Monkeys: Size Variation and Evolutionary Constraint and Flexibility

An organism is built through a series of contingent factors, yet it is determined by historical, physical, and developmental constraints. A constraint should not be understood as an absolute obstacle to evolution, as it may also generate new possibilities for evolutionary change. Modularity is, in this context, an important way of organizing biological information and has been recognized as a central concept in evolutionary biology bridging on developmental, genetics, morphological, biochemical, and physiological studies. In this article, we explore how modularity affects the evolution of a complex system in two mammalian lineages by analyzing correlation, variance/covariance, and residual matrices (without size variation). We use the multivariate response to selection equation to simulate the behavior of Eutheria and Metharia skulls in terms of their evolutionary flexibility and constraints. We relate these results to classical approaches based on morphological integration tests based on functional/developmental hypotheses. Eutherians (Neotropical primates) showed smaller magnitudes of integration compared with Metatheria (didelphids) and also skull modules more clearly delimited. Didelphids showed higher magnitudes of integration and their modularity is strongly influenced by within-groups size variation to a degree that evolutionary responses are basically aligned with size variation. Primates still have a good portion of the total variation based on size; however, their enhanced modularization allows a broader spectrum of responses, more similar to the selection gradients applied (enhanced flexibility). Without size variation, both groups become much more similar in terms of modularity patterns and magnitudes and, consequently, in their evolutionary flexibility. J. Exp. Zool. (Mol. Dev. Evol.) 314B:663-683, 2010. (C) 2010 Wiley-Liss, Inc.

Evaluating gene by sex and age interactions on cardiovascular risk factors in Brazilian families

Background: In family studies, it is important to evaluate the impact of genes and environmental factors on traits of interest. In particular, the relative influences of both genes and the environment may vary in different strata of the population of interest, such as young and old individuals, or males and females. Methods: In this paper, extensions of the variance components model are used to evaluate heterogeneity in the genetic and environmental variance components due to the effects of sex and age (the cutoff between young and old was 43 yrs). The data analyzed were from 81 Brazilian families (1,675 individuals) of the Baependi Family Heart Study. Results: The models allowing for heterogeneity of variance components by sex suggest that genetic and environmental variances are not different in males and females for diastolic blood pressure, LDL-cholesterol, and HDL-cholesterol, independent of the covariates included in the models. However, for systolic blood pressure, fasting glucose and triglycerides, the evidence for heterogeneity was dependent on the covariates in the model. For instance, in the presence of sex and age covariates, heterogeneity in the genetic variance component was suggested for fasting glucose. But, for systolic blood pressure, there was no evidence of heterogeneity in any of the two variance components. Except for the LDL-cholesterol, models allowing for heterogeneity by age provide evidence of heterogeneity in genetic variance for triglycerides and systolic and diastolic blood pressure. There was evidence of heterogeneity in environmental variance in fasting glucose and HDL-cholesterol. Conclusions: Our results suggest that heterogeneity in trait variances should not be ignored in the design and analyses of gene-finding studies involving these traits, as it may generate additional information about gene effects, and allow the investigation of more sophisticated models such as the model including sex-specific oligogenic variance components.

PARTIAL STABILITY FOR A CLASS OF NONLINEAR SYSTEMS

This paper studies a nonlinear, discrete-time matrix system arising in the stability analysis of Kalman filters. These systems present an internal coupling between the state components that gives rise to complex dynamic behavior. The problem of partial stability, which requires that a specific component of the state of the system converge exponentially, is studied and solved. The convergent state component is strongly linked with the behavior of Kalman filters, since it can be used to provide bounds for the error covariance matrix under uncertainties in the noise measurements. We exploit the special features of the system-mainly the connections with linear systems-to obtain an algebraic test for partial stability. Finally, motivated by applications in which polynomial divergence of the estimates is acceptable, we study and solve a partial semistability problem.

Multifractality in the random parameter model for multivariate time series

The Random Parameter model was proposed to explain the structure of the covariance matrix in problems where most, but not all, of the eigenvalues of the covariance matrix can be explained by Random Matrix Theory. In this article, we explore the scaling properties of the model, as observed in the multifractal structure of the simulated time series. We use the Wavelet Transform Modulus Maxima technique to obtain the multifractal spectrum dependence with the parameters of the model. The model shows a scaling structure compatible with the stylized facts for a reasonable choice of the parameter values. (C) 2009 Elsevier B.V. All rights reserved.