Novel method to estimate the phenotypic variation explained by genome-wide association studies reveals large fraction of the missing heritability.

Genome-wide association studies (GWAS) are conducted with the promise to discover novel genetic variants associated with diverse traits. For most traits, associated markers individually explain just a modest fraction of the phenotypic variation, but their number can well be in the hundreds. We developed a maximum likelihood method that allows us to infer the distribution of associated variants even when many of them were missed by chance. Compared to previous approaches, the novelty of our method is that it (a) does not require having an independent (unbiased) estimate of the effect sizes; (b) makes use of the complete distribution of P-values while allowing for the false discovery rate; (c) takes into account allelic heterogeneity and the SNP pruning strategy. We applied our method to the latest GWAS meta-analysis results of the GIANT consortium. It revealed that while the explained variance of genome-wide (GW) significant SNPs is around 1% for waist-hip ratio (WHR), the observed P-values provide evidence for the existence of variants explaining 10% (CI=[8.5-11.5%]) of the phenotypic variance in total. Similarly, the total explained variance likely to exist for height is estimated to be 29% (CI=[28-30%]), three times higher than what the observed GW significant SNPs give rise to. This methodology also enables us to predict the benefit of future GWA studies that aim to reveal more associated genetic markers via increased sample size.

Environmental factors, spatial variation, and specific requirements of Chironomidae in Mediterranean reference streams

Chironomidae spatial distribution was investigated at 63 near-pristine sites in 22 catchments of the Iberian Mediterranean coast. We used partial redundancy analysis to study Chironomidae community responses to a number of environmental factors acting at several spatial scales. The percentage of variation explained by local factors (23.3%) was higher than that explained by geographical (8.5%) or regional factors(8%). Catchment area, longitude, pH, % siliceous rocks in the catchment, and altitude were the best predictors of Chironomidae assemblages. We used a k-means cluster analysis to classified sites into 3 major groups based on Chironomidae assemblages. These groups were explained mainly by longitudinal zonation and geographical position, and were defined as 1) siliceous headwater streams, 2) mid-altitude streams with small catchment areas, and 3) medium-sized calcareous streams. Distinct species assemblages with associated indicator taxa were established for each stream category using IndVal analysis. Species responses to previously identified key environmental variables were determined, and optima and tolerances were established by weighted average regression. Distinct ecological requirements were observed among genera and among species of the same genus. Some genera were restricted to headwater systems (e.g., Diamesa), whereas others (e.g., Eukiefferiella) had wider ecological preferences but with distinct distributions among congenerics. In the present period of climate change, optima and tolerances of species might be a useful tool to predict responses of different species to changes in significant environmental variables, such as temperature and hydrology.

Estimates of patient costs related with population morbidity: Can indirect costs affect the results?

A number of health economics works require patient cost estimates as a basic information input.However the accuracy of cost estimates remains in general unspecified. We propose to investigate howthe allocation of indirect costs or overheads can affect the estimation of patient costs in order to allow forimprovements in the analysis of patient costs estimates. Instead of focusing on the costing method, thispaper proposes to highlight changes in variance explained observed when a methodology is chosen. Wecompare three overhead allocation methods for a specific Spanish population adjusted using the ClinicalRisk Groups (CRG), and we obtain different series of full-cost group estimates. As a result, there aresignificant gains in the proportion of the variance explained, depending upon the methodology used.Furthermore, we find that the global amount of variation explained by risk adjustment models dependsmainly on direct costs and is independent of the level of aggregation used in the classification system.

Étude pharmacogénomique sur l’utilisation de la warfarine en pratique clinique réelle

Contexte: Bien que plusieurs algorithmes pharmacogénétiques de prédiction de doses de warfarine aient été publiés, peu d’études ont comparé la validité de ces algorithmes en pratique clinique réelle. Objectif: Évaluer trois algorithmes pharmacogénomiques dans une population de patients qui initient un traitement à la warfarine et qui souffrent de fibrillation auriculaire ou de problèmes de valves cardiaques. Analyser la performance des algorithmes de Gage et al., de Michaud et al. ainsi que de l’IWPC quant à la prédiction de la dose de warfarine permettant d’atteindre l’INR thérapeutique. Méthodes: Un devis de cohorte rétrospectif fut utilisé afin d’évaluer la validité des algorithmes chez 605 patients ayant débuté une thérapie de warfarine à l’Institut de Cardiologie de Montréal. Le coefficient de corrélation de Pearson ainsi que l’erreur absolue moyenne ont été utilisés pour évaluer la précision des algorithmes. L’exactitude clinique des prédictions de doses fut évaluée en calculant le nombre de patients pour qui la dose prédite était sous-estimée, idéalement estimée ou surestimée. Enfin, la régression linéaire multiple a été utilisée pour évaluer la validité d’un modèle de prédiction de doses de warfarine obtenu en ajoutant de nouvelles covariables. Résultats : L’algorithme de Gage a obtenu la proportion de variation expliquée la plus élevée (R2 ajusté = 44 %) ainsi que la plus faible erreur absolue moyenne (MAE = 1.41 ± 0.06). De plus, la comparaison des proportions de patients ayant une dose prédite à moins de 20 % de la dose observée a confirmé que l’algorithme de Gage était également le plus performant. Conclusion : Le modèle publié par Gage en 2008 est l’algorithme pharmacogénétique le plus exact dans notre population pour prédire des doses thérapeutiques de warfarine.

Evaluation of milk urea nitrogen as a diagnostic of protein feeding

An evaluation of milk urea nitrogen (MUN) as a diagnostic of protein feeding in dairy cows was performed using mean treatment data (n = 306) from 50 production trials conducted in Finland (n = 48) and Sweden (n = 2). Data were used to assess the effects of diet composition and certain animal characteristics on MUN and to derive relationships between MUN and the efficiency of N utilization for milk production and urinary N excretion. Relationships were developed using regression analysis based on either models of fixed factors or using mixed models that account for between-experiment variations. Dietary crude protein (CP) content was the best single predictor of MUN and accounted for proportionately 0.778 of total variance [ MUN (mg/dL) = -14.2 + 0.17 x dietary CP content (g/kg dry matter)]. The proportion of variation explained by this relationship increased to 0.952 when a mixed model including the random effects of study was used, but both the intercept and slope remained unchanged. Use of rumen degradable CP concentration in excess of predicted requirements, or the ratio of dietary CP to metabolizable energy as single predictors, did not explain more of the variation in MUN (R-2 = 0.767 or 0.778, respectively) than dietary CP content. Inclusion of other dietary factors with dietary CP content in bivariate models resulted in only marginally better predictions of MUN (R-2 = 0.785 to 0.804). Closer relationships existed between MUN and dietary factors when nutrients (CP to metabolizable energy) were expressed as concentrations in the diet, rather than absolute intakes. Furthermore, both MUN and MUN secretion (g/d) provided more accurate predictions of urinary N excretion (R-2 = 0.787 and 0.835, respectively) than measurements of the efficiency of N utilization for milk production (R-2 = 0.769). It is concluded that dietary CP content is the most important nutritional factor influencing MUN, and that measurements of MUN can be utilized as a diagnostic of protein feeding in the dairy cow and used to predict urinary N excretion.

The Evolution of Modularity in the Mammalian Skull II: Evolutionary Consequences

Changes in patterns and magnitudes of integration may influence the ability of a species to respond to selection. Consequently, modularity has often been linked to the concept of evolvability, but their relationship has rarely been tested empirically. One possible explanation is the lack of analytical tools to compare patterns and magnitudes of integration among diverse groups that explicitly relate these aspects to the quantitative genetics framework. We apply such framework here using the multivariate response to selection equation to simulate the evolutionary behavior of several mammalian orders in terms of their flexibility, evolvability and constraints in the skull. We interpreted these simulation results in light of the integration patterns and magnitudes of the same mammalian groups, described in a companion paper. We found that larger magnitudes of integration were associated with a blur of the modules in the skull and to larger portions of the total variation explained by size variation, which in turn can exert a strong evolutionary constraint, thus decreasing the evolutionary flexibility. Conversely, lower overall magnitudes of integration were associated with distinct modules in the skull, to smaller fraction of the total variation associated with size and, consequently, to weaker constraints and more evolutionary flexibility. Flexibility and constraints are, therefore, two sides of the same coin and we found them to be quite variable among mammals. Neither the overall magnitude of morphological integration, the modularity itself, nor its consequences in terms of constraints and flexibility, were associated with absolute size of the organisms, but were strongly associated with the proportion of the total variation in skull morphology captured by size. Therefore, the history of the mammalian skull is marked by a trade-off between modularity and evolvability. Our data provide evidence that, despite the stasis in integration patterns, the plasticity in the magnitude of integration in the skull had important consequences in terms of evolutionary flexibility of the mammalian lineages.

Identification of QTLs associated with citrus resistance to Phytophthora gummosis

Citrus gummosis, caused by Phylophthora spp., is an important citrus disease in Brazil. Almost all citrus rootstock varieties are susceptible to it to some degree, whereas resistance is present in Poncirus trifoliata, a closely related species. The objective of this study was to detect QTLs linked to citrus Phylophthora gummosis resistance. Eighty individuals of the F, progeny, obtained by controlled crosses between Sunki mandarin Citrus sunki (susceptible) and Poncirus trifoliata cv. Rubidoux (resistant), were evaluated. Resistance to Phytophthora parasitica was evaluated by inoculating stems of young plants with a disc of fungal mycelia and measuring lesion lengths a month later. Two QTLs linked to gummosis resistance were detected in linkage groups I and 5 of the P. trifoliala map, and one QTL in linkage group 2 of the C sunki map. The phenotypic variation explained by individual QTLs was 14% for C sunki and ranged from 16 to 24% for P. trifoliala. The low character heritability (h(2) = 18.7%) and the detection of more than one QTL associated with citrus Phytophthora gummosis resistance showed that inheritance of the resistance is quantitative.

Mapeamento de QTL e análises de espectroscopia para teor de óleo visando aplicação em programas de melhoramento de soja

Pós-graduação em Agronomia (Genética e Melhoramento de Plantas) - FCAV

Marker-trait association and epistasis for brown rust resistance in sugarcane

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The role of landuse and natural determinants for grassland vegetation composition in the Swiss Alps

The Alps provide a high habitat diversity for plant species, structured by broad- and fine-scale abiotic site conditions. In man-made grasslands, vegetation composition is additionally affected by the type of landuse. We recorded vegetation composition in 216 parcels of grassland in 12 municipalities representing an area of 170 x 70 km in the south-eastern part of the Swiss Alps. Each parcel was characterized by a combination of altitudinal level (valley, intermediate, alp). traditional landuse (mown. grazed), current management (mown, grazed, abandoned). and Fertilization (unfertilized, fertilized). For each parcel we also assessed the abiotic factors aspect, slope, pH value, and geographic coordinates, and for each municipality annual precipitation and its cultural tradition. We analysed vegetation composition using (i) variation partitioning in RDA. (ii) cover of graminoids. non-legume forbs, and legumes, and (iii) dominance and frequency of species. Species composition was determined by, in decreasing order of variation explained. landuse, broad-scale abiotic factors, fine-scale abiotic factors. and cultural tradition. Current socio-economically motivated landuse changes, such as grazing of unfertilized former meadows or their abandonment, strongly affect vegetation composition. In our study, the frequency of characteristic meadow species was significantly smaller in grazed and even smaller in abandoned parcels than in still mown ones, suggesting less severe consequences of grazing for vegetation composition than of abandonment. Therefore. low-intensity grazing and mowing every few years should be considered valuable conservation alternatives to abandonment. Furthermore. because each landuse type was characterized by different species. a high variety of landuse types should be promoted to preserve plant species diversity in Alpine grasslands. (C) 2007 Gesellschaft fur Okologie. Published by Elsevier GmbH. All rights reserved.

Piggery pond sludge as a nitrogen source for crops - 2. Assay of wet and stockpiled piggery pond sludge by successive cereal crops or direct measurement of soil available N

The appropriate use of wastes is a significant issue for the pig industry due to increasing pressure from regulatory authorities to protect the environment from pollution. Nitrogen contained in piggery pond sludge ( PPS) is a potential source of supplementary nutrient for crop production. Nitrogen contribution following the application of PPS to soil was obtained from 2 field experiments on the Darling Downs in southern Queensland on contrasting soil types, a cracking clay ( Vertosol) and a hardsetting sandy loam (Sodosol), and related to potentially mineralisable N from laboratory incubations conducted under controlled conditions and NO3- accumulation in the field. Piggery pond sludge was applied as-collected ( wet PPS) and following stockpiling to dry ( stockpiled PPS). Soil NO3- levels increased with increased application rates of wet and stockpiled PPS. Supplementary N supply from PPS estimated by fertiliser equivalence was generally unsatisfactory due to poor precision with this method, and also due to a high level of NO3- in the clay soil before the first assay crop. Also low recoveries of N by subsequent sorghum ( Sorghum bicolor) and wheat ( Triticum aestivum) assay crops at the 2 sites due to low in-crop rainfall in 1999 resulted in low apparent N availability. Over all, 29% ( range 12 - 47%) of total N from the wet PPS and 19% ( range 0 - 50%) from the stockpiled PPS were estimated to be plant-available N during the assay period. The high concentration of NO3- for the wet PPS application on sandy soil after the first assay crop ( 1998 barley, Hordeum vulgare) suggests that leaching of NO3- could be of concern when high rates of wet PPS are applied before infrequent periods of high precipitation, due primarily to the mineral N contained in wet PPS. Low yields, grain protein concentrations, and crop N uptake of the sorghum crop following the barley crop grown on the clay soil demonstrated a low residual value of N applied in PPS. NO3- in the sandy soil before sowing accounted for 79% of the variation in plant N uptake and was a better index than anaerobically mineralisable N ( 19% of variation explained). In clay soil, better prediction of crop N uptake was obtained when both anaerobically mineralisable N (39% of variation explained) and soil pro. le NO3- were used in combination (R-2 = 0.49).

Assembly processes of invertebrate communities in springs across different spatial scales

Il presente lavoro ha lo scopo di comprendere i processi sottesi ai pattern di coesistenza tra le specie di invertebrati sorgentizi, distinguendo tra dinamiche stocastiche e deterministiche. Le sorgenti sono ecosistemi complessi e alcune loro caratteristiche (ad esempio l’insularità, la stabilità termica, la struttura ecotonale “a mosaico”, la frequente presenza di specie rare ed endemiche, o l’elevata diversità in taxa) le rendono laboratori naturali utili allo studio dei processi ecologici, tra cui i processi di assembly. Al fine di studiare queste dinamiche è necessario un approccio multi-scala, per questo motivi sono state prese in considerazione tre scale spaziali. A scala locale è stato compiuto un campionamento stagionale su sette sorgenti (quattro temporanee e tre permanenti) del Monte Prinzera, un affioramento ofiolitico vicino alla città di Parma. In questa area sono stati valutati l’efficacia e l’impatto ambientale di diversi metodi di campionamento e sono stati analizzati i drivers ecologici che influenzano le comunità. A scala più ampia sono state campionate per due volte 15 sorgenti della regione Emilia Romagna, al fine di identificare il ruolo della dispersione e la possibile presenza di un effetto di niche-filtering. A scala continentale sono state raccolte informazioni di letteratura riguardanti sorgenti dell’area Paleartica occidentale, e sono stati studiati i pattern biogeografici e l’influenza dei fattori climatici sulle comunità. Sono stati presi in considerazione differenti taxa di invertebrati (macroinvertebrati, ostracodi, acari acquatici e copepodi), scegliendo tra quelli che si prestavano meglio allo studio dei diversi processi in base alle loro caratteristiche biologiche e all’approfondimento tassonomico raggiungibile. I campionamenti biologici in sorgente sono caratterizzati da diversi problemi metodologici e possono causare impatti sugli ambienti. In questo lavoro sono stati paragonati due diversi metodi: l’utilizzo del retino con un approccio multi-habitat proporzionale e l’uso combinato di trappole e lavaggio di campioni di vegetazione. Il retino fornisce dati più accurati e completi, ma anche significativi disturbi sulle componenti biotiche e abiotiche delle sorgenti. Questo metodo è quindi raccomandato solo se il campionamento ha come scopo un’approfondita analisi della biodiversità. D’altra parte l’uso delle trappole e il lavaggio della vegetazione sono metodi affidabili che presentano minori impatti sull’ecosistema, quindi sono adatti a studi ecologici finalizzati all’analisi della struttura delle comunità. Questo lavoro ha confermato che i processi niche-based sono determinanti nello strutturare le comunità di ambienti sorgentizi, e che i driver ambientali spiegano una rilevante percentuale della variabilità delle comunità. Infatti le comunità di invertebrati del Monte Prinzera sono influenzate da fattori legati al chimismo delle acque, alla composizione e all’eterogeneità dell’habitat, all’idroperiodo e alle fluttuazioni della portata. Le sorgenti permanenti mostrano variazioni stagionali per quanto riguarda le concentrazioni dei principali ioni, mentre la conduttività, il pH e la temperatura dell’acqua sono più stabili. È probabile che sia la stabilità termica di questi ambienti a spiegare l’assenza di variazioni stagionali nella struttura delle comunità di macroinvertebrati. L’azione di niche-filtering delle sorgenti è stata analizzata tramite lo studio della diversità funzionale delle comunità di ostracodi dell’Emilia-Romagna. Le sorgenti ospitano più del 50% del pool di specie regionale, e numerose specie sono state rinvenute esclusivamente in questi habitat. Questo è il primo studio che analizza la diversità funzionale degli ostracodi, è stato quindi necessario stilare una lista di tratti funzionali. Analizzando il pool di specie regionale, la diversità funzionale nelle sorgenti non è significativamente diversa da quella misurata in comunità assemblate in maniera casuale. Le sorgenti non limitano quindi la diversità funzionale tra specie coesistenti, ma si può concludere che, data la soddisfazione delle esigenze ecologiche delle diverse specie, i processi di assembly in sorgente potrebbero essere influenzati da fattori stocastici come la dispersione, la speciazione e le estinzioni locali. In aggiunta, tutte le comunità studiate presentano pattern spaziali riconoscibili, rivelando una limitazione della dispersione tra le sorgenti, almeno per alcuni taxa. Il caratteristico isolamento delle sorgenti potrebbe essere la causa di questa limitazione, influenzando maggiormente i taxa a dispersione passiva rispetto a quelli a dispersione attiva. In ogni caso nelle comunità emiliano-romagnole i fattori spaziali spiegano solo una ridotta percentuale della variabilità biologica totale, mentre tutte le comunità risultano influenzate maggiormente dalle variabili ambientali. Il controllo ambientale è quindi prevalente rispetto a quello attuato dai fattori spaziali. Questo risultato dimostra che, nonostante le dinamiche stocastiche siano importanti in tutte le comunità studiate, a questa scala spaziale i fattori deterministici ricoprono un ruolo prevalente. I processi stocastici diventano più influenti invece nei climi aridi, dove il disturbo collegato ai frequenti eventi di disseccamento delle sorgenti provoca una dinamica source-sink tra le diverse comunità. Si è infatti notato che la variabilità spiegata dai fattori ambientali diminuisce all’aumentare dell’aridità del clima. Disturbi frequenti potrebbero provocare estinzioni locali seguite da ricolonizzazioni di specie provenienti dai siti vicini, riducendo la corrispondenza tra gli organismi e le loro richieste ambientali e quindi diminuendo la quantità di variabilità spiegata dai fattori ambientali. Si può quindi concludere che processi deterministici e stocastici non si escludono mutualmente, ma contribuiscono contemporaneamente a strutturare le comunità di invertebrati sorgentizi. Infine, a scala continentale, le comunità di ostracodi sorgentizi mostrano chiari pattern biogeografici e sono organizzate lungo gradienti ambientali principalmente collegati altitudine, latitudine, temperatura dell’acqua e conducibilità. Anche la tipologia di sorgente (elocrena, reocrena o limnocrena) è influente sulla composizione delle comunità. La presenza di specie rare ed endemiche inoltre caratterizza specifiche regioni geografiche.

Periphyton light transmission relationships in Florida Bay and the Florida Keys, USA

Light transmission was measured through intact, submerged periphyton communities on artificial seagrass leaves. The periphyton communities were representative of the communities on Thalassia testudinum in subtropical seagrass meadows. The periphyton communities sampled were adhered carbonate sediment, coralline algae, and mixed algal assemblages. Crustose or film-forming periphyton assemblages were best prepared for light transmission measurements using artificial leaves fouled on both sides, while measurements through three-dimensional filamentous algae required the periphyton to be removed from one side. For one-sided samples, light transmission could be measured as the difference between fouled and reference artificial leaf samples. For two-sided samples, the percent periphyton light transmission to the leaf surface was calculated as the square root of the fraction of incident light. Linear, exponential, and hyperbolic equations were evaluated as descriptors of the periphyton dry weight versus light transmission relationship. Hyperbolic and exponential decay models were superior to linear models and exhibited the best fits for the observed relationships. Differences between the coefficients of determination (r2) of hyperbolic and exponential decay models were statistically insignificant. Constraining these models for 100% light transmission at zero periphyton load did not result in any statistically significant loss in the explanatory capability of the models. In most all cases, increasing model complexity using three-parameter models rather than two-parameter models did not significantly increase the amount of variation explained. Constrained two-parameter hyperbolic or exponential decay models were judged best for describing the periphyton dry weight versus light transmission relationship. On T. testudinum in Florida Bay and the Florida Keys, significant differences were not observed in the light transmission characteristics of the varying periphyton communities at different study sites. Using pooled data from the study sites, the hyperbolic decay coefficient for periphyton light transmission was estimated to be 4.36 mg dry wt. cm−2. For exponential models, the exponential decay coefficient was estimated to be 0.16 cm2 mg dry wt.−1.

Geographical patterns of the species richness of helminth parasites of moles (Talpa spp.) in Spain: Separating the effects of sampling effort from those of other conditioning factors

We analysed the viscera of 534 moles (Ta l p a spp.) from 30 of the 47 provinces of peninsular Spain, including 255 individuals of T. europaea from eight provinces, 154 individuals of T. occidentalis from 20 provinces, and 125 unidentified Ta l p a individuals from two provinces. We identified their helminth parasites and determined parasite species richness. We related parasite species richness with sampling effort using both a linear and a logarithmic function. We then performed stepwise linear regressions to predict mole parasite species richness from a small set of selected predictor variables that included sampling effort. We applied the resulting models to forecast T. euro p a e a, T. occidentalis, and Ta l p a spp. parasite species richness in all provinces with recorded host presence, assuming different levels of sampling eff o r t . F i n a l l y, we used partial regression analysis to partition the variation explained by each of the selected variables in the models. We found that mole parasite species richness is strongly conditioned by sampling effort, but that other factors such as cropland area and environmental disturbance have significant independent effects.

Relative importance of environment, human activity and spatial situation in determining the distribution of terrestrial mammal diversity in Argentina

We recorded the number of terrestrial mammal species in each Argentinian province, and the number of species belonging to particular groups (Marsupialia, Placentaria, and among the latter, Xenarthra, Carnivora, Ungulates and Rodentia). We performed multiple regressions of each group’s SR on environmental, human and spatial variables, to determine the amounts of variation explained by these factors. We then used a variance partitioning procedure to specify which proportion of the variation in SR is explained by each of the three factors exclusively and which proportions are attributable to interactions between factors.