Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Prophylactic antibiotics or G(M)-CSF for the prevention of infections and improvement of survival in cancer patients receiving myelotoxic chemotherapy.

BACKGROUND Febrile neutropenia (FN) and other infectious complications are some of the most serious treatment-related toxicities of chemotherapy for cancer, with a mortality rate of 2% to 21%. The two main types of prophylactic regimens are granulocyte (macrophage) colony-stimulating factors (G(M)-CSF) and antibiotics, frequently quinolones or cotrimoxazole. Current guidelines recommend the use of colony-stimulating factors when the risk of febrile neutropenia is above 20%, but they do not mention the use of antibiotics. However, both regimens have been shown to reduce the incidence of infections. Since no systematic review has compared the two regimens, a systematic review was undertaken. OBJECTIVES To compare the efficacy and safety of G(M)-CSF compared to antibiotics in cancer patients receiving myelotoxic chemotherapy. SEARCH METHODS We searched The Cochrane Library, MEDLINE, EMBASE, databases of ongoing trials, and conference proceedings of the American Society of Clinical Oncology and the American Society of Hematology (1980 to December 2015). We planned to include both full-text and abstract publications. Two review authors independently screened search results. SELECTION CRITERIA We included randomised controlled trials (RCTs) comparing prophylaxis with G(M)-CSF versus antibiotics for the prevention of infection in cancer patients of all ages receiving chemotherapy. All study arms had to receive identical chemotherapy regimes and other supportive care. We included full-text, abstracts, and unpublished data if sufficient information on study design, participant characteristics, interventions and outcomes was available. We excluded cross-over trials, quasi-randomised trials and post-hoc retrospective trials. DATA COLLECTION AND ANALYSIS Two review authors independently screened the results of the search strategies, extracted data, assessed risk of bias, and analysed data according to standard Cochrane methods. We did final interpretation together with an experienced clinician. MAIN RESULTS In this updated review, we included no new randomised controlled trials. We included two trials in the review, one with 40 breast cancer patients receiving high-dose chemotherapy and G-CSF compared to antibiotics, a second one evaluating 155 patients with small-cell lung cancer receiving GM-CSF or antibiotics.We judge the overall risk of bias as high in the G-CSF trial, as neither patients nor physicians were blinded and not all included patients were analysed as randomised (7 out of 40 patients). We considered the overall risk of bias in the GM-CSF to be moderate, because of the risk of performance bias (neither patients nor personnel were blinded), but low risk of selection and attrition bias.For the trial comparing G-CSF to antibiotics, all cause mortality was not reported. There was no evidence of a difference for infection-related mortality, with zero events in each arm. Microbiologically or clinically documented infections, severe infections, quality of life, and adverse events were not reported. There was no evidence of a difference in frequency of febrile neutropenia (risk ratio (RR) 1.22; 95% confidence interval (CI) 0.53 to 2.84). The quality of the evidence for the two reported outcomes, infection-related mortality and frequency of febrile neutropenia, was very low, due to the low number of patients evaluated (high imprecision) and the high risk of bias.There was no evidence of a difference in terms of median survival time in the trial comparing GM-CSF and antibiotics. Two-year survival times were 6% (0 to 12%) in both arms (high imprecision, low quality of evidence). There were four toxic deaths in the GM-CSF arm and three in the antibiotics arm (3.8%), without evidence of a difference (RR 1.32; 95% CI 0.30 to 5.69; P = 0.71; low quality of evidence). There were 28% grade III or IV infections in the GM-CSF arm and 18% in the antibiotics arm, without any evidence of a difference (RR 1.55; 95% CI 0.86 to 2.80; P = 0.15, low quality of evidence). There were 5 episodes out of 360 cycles of grade IV infections in the GM-CSF arm and 3 episodes out of 334 cycles in the cotrimoxazole arm (0.8%), with no evidence of a difference (RR 1.55; 95% CI 0.37 to 6.42; P = 0.55; low quality of evidence). There was no significant difference between the two arms for non-haematological toxicities like diarrhoea, stomatitis, infections, neurologic, respiratory, or cardiac adverse events. Grade III and IV thrombopenia occurred significantly more frequently in the GM-CSF arm (60.8%) compared to the antibiotics arm (28.9%); (RR 2.10; 95% CI 1.41 to 3.12; P = 0.0002; low quality of evidence). Neither infection-related mortality, incidence of febrile neutropenia, nor quality of life were reported in this trial. AUTHORS' CONCLUSIONS As we only found two small trials with 195 patients altogether, no conclusion for clinical practice is possible. More trials are necessary to assess the benefits and harms of G(M)-CSF compared to antibiotics for infection prevention in cancer patients receiving chemotherapy.

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.

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.

Normalization procedure for relaxation studies in NMR quantum information processing

NMR quantum information processing studies rely on the reconstruction of the density matrix representing the so-called pseudo-pure states (PPS). An initially pure part of a PPS state undergoes unitary and non-unitary (relaxation) transformations during a computation process, causing a ""loss of purity"" until the equilibrium is reached. Besides, upon relaxation, the nuclear polarization varies in time, a fact which must be taken into account when comparing density matrices at different instants. Attempting to use time-fixed normalization procedures when relaxation is present, leads to various anomalies on matrices populations. On this paper we propose a method which takes into account the time-dependence of the normalization factor. From a generic form for the deviation density matrix an expression for the relaxing initial pure state is deduced. The method is exemplified with an experiment of relaxation of the concurrence of a pseudo-entangled state, which exhibits the phenomenon of sudden death, and the relaxation of the Wigner function of a pseudo-cat state.

Permutation-based stochastic ordering using pairwise comparisons

The topic of this work concerns nonparametric permutation-based methods aiming to find a ranking (stochastic ordering) of a given set of groups (populations), gathering together information from multiple variables under more than one experimental designs. The problem of ranking populations arises in several fields of science from the need of comparing G>2 given groups or treatments when the main goal is to find an order while taking into account several aspects. As it can be imagined, this problem is not only of theoretical interest but it also has a recognised relevance in several fields, such as industrial experiments or behavioural sciences, and this is reflected by the vast literature on the topic, although sometimes the problem is associated with different keywords such as: "stochastic ordering", "ranking", "construction of composite indices" etc., or even "ranking probabilities" outside of the strictly-speaking statistical literature. The properties of the proposed method are empirically evaluated by means of an extensive simulation study, where several aspects of interest are let to vary within a reasonable practical range. These aspects comprise: sample size, number of variables, number of groups, and distribution of noise/error. The flexibility of the approach lies mainly in the several available choices for the test-statistic and in the different types of experimental design that can be analysed. This render the method able to be tailored to the specific problem and the to nature of the data at hand. To perform the analyses an R package called SOUP (Stochastic Ordering Using Permutations) has been written and it is available on CRAN.

Skeletal matrices, muci, and the origin of invertebrate calcification.

The sudden appearance of calcified skeletons among many different invertebrate taxa at the Precambrian-Cambrian transition may have required minor reorganization of preexisting secretory functions. In particular, features of the skeletal organic matrix responsible for regulating crystal growth by inhibition may be derived from mucous epithelial excretions. The latter would have prevented spontaneous calcium carbonate overcrusting of soft tissues exposed to the highly supersaturated Late Proterozoic ocean [Knoll, A. H., Fairchild, I. J. & Swett, K. (1993) Palaios 8, 512-525], a putative function for which we propose the term "anticalcification." We tested this hypothesis by comparing the serological properties of skeletal water-soluble matrices and mucous excretions of three invertebrates--the scleractinian coral Galaxea fascicularis and the bivalve molluscs Mytilus edulis and Mercenaria mercenaria. Crossreactivities recorded between muci and skeletal water-soluble matrices suggest that these different secretory products have a high degree of homology. Furthermore, freshly extracted muci of Mytilus were found to inhibit calcium carbonate precipitation in solution.

Porcine skin as a source of biodegradable matrices: alkaline treatment and glutaraldehyde crosslinking

In this work, the modifications promoted by alkaline hydrolysis and glutaraldehyde (GA) crosslinking on type I collagen found in porcine skin have been studied. Collagen matrices were obtained from the alkaline hydrolysis of porcine skin, with subsequent GA crosslinking in different concentrations and reaction times. The elastin content determination showed that independent of the treatment, elastin was present in the matrices. Results obtained from in vitro trypsin degradation indicated that with the increase of GA concentration and reaction time, the degradation rate decreased. From thermogravimetry and differential scanning calorimetry analysis it can be observed that the collagen in the matrices becomes more resistant to thermal degradation as a consequence of the increasing crosslink degree. Scanning electron microscopy analysis indicated that after the GA crosslinking, collagen fibers become more organized and well-defined. Therefore, the preparations of porcine skin matrices with different degradation rates, which can be used in soft tissue reconstruction, are viable.

Protein-polysaccharide viscoelastic matrices: synergic effects of amylose on lysozyme physical gelation in aqueous dimethylsulfoxide

A synergic effect of amylose on rheological characteristics of lysozyme physical gels evolved out of dimethylsulfoxide-water was verified and analyzed. The dynamics of the gels were experimentally approached by oscillatory rheology. The synergic effect was characterized by a decrease in the threshold DMSO volume fraction required for lysozyme gelation, and by a significant strengthening of the gel structure at over-critical solvent and protein concentrations. Drastic changes in the relaxation and creep curve patterns for systems in the presence of amylose were verified. Complex viscosity dependence on temperature was found to conform to an Arrhenius-like equation, allowing the determination of an activation energy term (Ea, apparent) for discrimination of gel rigidity. A dilatant effect was found to be induced by temperature on the flow behavior of lysozyme dispersions in DMSO-H(2)O in sub-critical conditions for gelation, which was greatly intensified by the presence of amylose in the samples. That transition was interpreted as reflecting a change from a predominant colloidal flow regime, where globular components are the prevailing structural elements, to a mainly fibrillar, polymeric flow behavior.

Morphological and mechanical properties of hybrid matrices of polysiloxane-polyvinyl alcohol prepared by sol-gel technique and their potential for immobilizing enzyme

Hybrid matrices of polysiloxane-polyvinyl alcohol (POS-PVA) were prepared by sol-gel technique using different concentrations of the organic component (polyvinyl alcohol, PVA) in the synthesis medium. The goal was to prepare carriers for immobilizing enzyme by taking into consideration properties as hardness, mean pore diameter, specific surface area and pore size distribution. The matrices were activated with sodium metaperiodate to render functional groups for binding the lipase from Candida rugosa, used here as a study model. Results showed that low proportion of PVA gave POS-PVA with low surface area and pore volume, although with higher hardness. The chemical activation decreased the pore volume and increased the pore size with a decrease on the surface area of about 60-75%. The matrices for enzyme immobilization were chosen considering the best combination of high surface area and hardness. Thus, the POS-PVA prepared with 5.56 x 10(-5) M of PVA with a surface area of 123 m(2)/g and hardness of 71 HV (50 gf 30 s) was shown to be suitable to immobilize the lipase, with an immobilization yield of about 40%. (c) 2008 Elsevier B.V. All rights reserved.

Sources of Variations When Comparing Concrete Carbonation Results

Carbonation is one of the main concerns for concrete service life in tropical countries. The mechanism and materials that produce it have been widely studied as well as natural and accelerated methods to report and analyze it. In spite of reported investigations, there is a need for information that could allow an adequate interpretation of the results of the standardization process. This lack of information can produce variations not only in the interpretation but also in the predictions of service life. The purpose of this paper is to analyze and discuss variables that could be sources of error, especially when performing accelerated tests. As a result, a methodologies to minimize variations when interpreting and comparing results is proposed, such as specimen geometry and preconditioning, spacing, relative humidity, and CO(2) concentration.

An unbridgeable gap? Comparing the HIV/AIDS epidemics in Australia and sub-Saharan Africa

OBJECTIVE: To compare the HIV/AIDS epidemics in Australia and sub-Saharan Africa, to outline reasons for differences, and to consider implications for the Asia and Pacific region. METHODS: Comparison of key indicators of the epidemic in Australia, and Africa viewed largely through the experience of the Hlabisa health district, South Africa. RESULTS: To the end of 1997, for all Australia, the estimated cumulative number of HIV infections was approximately 19,000, whereas in Hlabisa 31,000 infections are estimated to have occurred. Compared with the low and declining incidence of HIV in Australia (<1%), estimated incidence in Hlabisa rose to 10% in 1997. In all, 94% of Australian infections have been amongst men; in Hlabisa equal numbers of males and females are infected. Consequently, whereas 3000 children were perinatally exposed to HIV in Hlabisa in 1998 alone, 160 Australian children have been exposed this way. In Australia, HIV-related disease is characterised by opportunistic infection whereas in Hlabisa tuberculosis and wasting dominate. Surveys among gay men in Sydney and Melbourne indicate >80% of HIV infected people receive antiretroviral therapy whereas in Hlabisa these drugs are not available. IMPLICATIONS: It seems possible that Asia and the Pacific will experience a similar HIV/AIDS epidemic to that in Africa. Levels of HIV are already high in parts of Asia, and social conditions in parts of the region might be considered ripe for the spread of HIV. As Australia strengthens economic and political ties within the region, so should more be done to help Pacific and Asian neighbours to prevent and respond to the HIV epidemic.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.

A randomized trial comparing hormone replacement therapy (HRT) and HRT plus calcitriol in the treatment of postmenopausal osteoporosis with vertebral fractures: Benefit of the combination on total body and hip density

We report a prospective, randomized, multi-center, open-label 2-year trial of 81 postmenopausal women aged 53-79 years with at least one minimal-trauma vertebral fracture (VF) and low (T-score below 2) lumbar bone mineral density (BMD). Group HRT received piperazine estrone sulfate (PES) 0.625 - 1.25 mg/d +/- medroxyprogesterone acetate (MPA) 2.5 - 5 mg/d,- group HRT/D received HRT plus calcitriol 0.25 mug bd. All with a baseline dietary calcium (Ca) of < I g/d received Ca carbonate 0.6 g nocte. Final data were on 66 - 70 patients. On HRT/D, significant (P < 0.001) BNID increases from baseline by DXA were at total body - head, trochanter, Ward's, total hip, inter-trochanter and femoral shaft (% group mean Delta 4.2, 6.1, 9.3. 3.7. 3.3 and 3.3%, respectively). On HRT, at these significant Deltas were restricted to the trochanter and sites. si Wards. Significant advantages of HRT/D over HRT were in BMD of total body (- head), total hip and trochanter (all P = 0.01). The differences in mean Delta at these sites were 1.3, 2.6 and 3.9%. At the following, both groups Improved significantly -lumbar spine (AP and lateral), forearm shaft and ultradistal tibia/fibula. The weightbearing, site - specific benefits of the combination associated with significant suppression of parathyroid hormone-suggest a beneficial effect on cortical bone. Suppression of bone turnover was significantly greater on HRT/D (serum osteocalcin P = 0.024 and urinary hydroxyproline/creatinine ratio P = 0.035). There was no significant difference in the number of patients who developed fresh VFs during the trial (HRT 8/36, 22%; HRT/D 4/34, 12% - intention to treat); likewise in the number who developed incident nonvertebral fractures. This Is the first study comparing the 2 treatments in a fracture population. The results indicate a significant benefit of calcitriol combined with HRT on total body BMD and on BNID at the hip, the major site of osteoporotic fracture.