Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulae for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.

Meta-analysis of correlated traits via summary statistics from GWASs with an application in hypertension

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple-even distinct-traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10(-8)) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10(-7)) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.

Gene by environment QTL mapping through multiple trait analyses in blood pressure salt-sensitivity: identification of a novel QTL in rat chromosome 5

Background The genetic mechanisms underlying interindividual blood pressure variation reflect the complex interplay of both genetic and environmental variables. The current standard statistical methods for detecting genes involved in the regulation mechanisms of complex traits are based on univariate analysis. Few studies have focused on the search for and understanding of quantitative trait loci responsible for gene × environmental interactions or multiple trait analysis. Composite interval mapping has been extended to multiple traits and may be an interesting approach to such a problem. Methods We used multiple-trait analysis for quantitative trait locus mapping of loci having different effects on systolic blood pressure with NaCl exposure. Animals studied were 188 rats, the progenies of an F2 rat intercross between the hypertensive and normotensive strain, genotyped in 179 polymorphic markers across the rat genome. To accommodate the correlational structure from measurements taken in the same animals, we applied univariate and multivariate strategies for analyzing the data. Results We detected a new quantitative train locus on a region close to marker R589 in chromosome 5 of the rat genome, not previously identified through serial analysis of individual traits. In addition, we were able to justify analytically the parametric restrictions in terms of regression coefficients responsible for the gain in precision with the adopted analytical approach. Conclusion Future work should focus on fine mapping and the identification of the causative variant responsible for this quantitative trait locus signal. The multivariable strategy might be valuable in the study of genetic determinants of interindividual variation of antihypertensive drug effectiveness.

The multiple signals assessed by female satin bowerbirds: could they be used to narrow down females' choices of mates?

Female choice based on multiple male traits has been documented in many species but the functions of such multiple traits are still under debate. The satin bowerbird has a polygynous mating system in which males attract females to bowers for mating; females choose mates based on multiple aspects of males and their bowers. In this paper, we demonstrate that females use some cues to decide which males to examine closely and other cues to decide which males to mate with. Female visitation rates to bowers were significantly related to male size and the males' 'solitary' display rates, and, to a lesser extent, to the numbers of bower decorations. After controlling for female visitation rates, it was found that a male's mating success was significantly related to his size and the rate at which he 'painted' his bower with saliva and chewed up plant material.

A simple method to localise pleiotropic susceptibility loci using univariate linkage analyses of correlated traits

Univariate linkage analysis is used routinely to localise genes for human complex traits. Often, many traits are analysed but the significance of linkage for each trait is not corrected for multiple trait testing, which increases the experiment-wise type-I error rate. In addition, univariate analyses do not realise the full power provided by multivariate data sets. Multivariate linkage is the ideal solution but it is computationally intensive, so genome-wide analysis and evaluation of empirical significance are often prohibitive. We describe two simple methods that efficiently alleviate these caveats by combining P-values from multiple univariate linkage analyses. The first method estimates empirical pointwise and genome-wide significance between one trait and one marker when multiple traits have been tested. It is as robust as an appropriate Bonferroni adjustment, with the advantage that no assumptions are required about the number of independent tests performed. The second method estimates the significance of linkage between multiple traits and one marker and, therefore, it can be used to localise regions that harbour pleiotropic quantitative trait loci (QTL). We show that this method has greater power than individual univariate analyses to detect a pleiotropic QTL across different situations. In addition, when traits are moderately correlated and the QTL influences all traits, it can outperform formal multivariate VC analysis. This approach is computationally feasible for any number of traits and was not affected by the residual correlation between traits. We illustrate the utility of our approach with a genome scan of three asthma traits measured in families with a twin proband.

Developmental, metabolic and immunological costs of flea infestation in the common vole

Parasites use resources from their hosts, which can indirectly affect a number of host functions because of trade-offs in resource allocation. In order to get a comprehensive view of the costs imposed by blood sucking parasites to their hosts, it is important to monitor multiple components of the development and physiology of parasitized hosts over long time periods. The effect of infestation by fleas on body mass, body length growth, haematocrit, resistance to oxidative stress, resting metabolic rate and humoral immune response were experimentally evaluated. During a 3-month period, male common voles, Microtus arvalis, were either parasitized by rat fleas (Nosopsyllus fasciatus), which are naturally occurring generalist ectoparasites of voles, or reared without fleas. Then voles were challenged twice by injecting Keyhole Limpet Haemocyanin (KLH) to assess whether the presence of fleas affects the ability of voles to produce antibodies against a novel antigen. During the immune challenge we measured the evolution of body mass, haematocrit, resistance to oxidative stress and antibody production. Flea infestation negatively influenced the growth of voles. Moreover, parasitized voles had reduced haematocrit, higher resting metabolic rate and lower production of antibodies against the KLH. Resistance to oxidative stress was not influenced by the presence of fleas. During the immune challenge with KLH, body mass decreased in both groups, while the resistance to oxidative stress remained stable. In contrast, the haematocrit decreased only in parasitized voles. Our experiment shows that infestation by a haematophageous parasite negatively affects multiple traits like growth, energy consumption and immune response. Fleas may severely reduce the survival probability and reproductive success of their host in natural conditions.

Diversity of condensed tannin structures affects rumen in vitro methane production in sainfoin (Onobrychis viciifolia) accessions

Sainfoin is a non-bloating temperate forage legume with a moderate-to-high condensed tannin (CT) content. This study investigated whether the diversity of sainfoin accessions in terms of CT structures and contents could be related to rumen in vitro gas and methane (CH4) production and fermentation characteristics. The aim was to identify promising accessions for future investigations. Accessions differed (P < 0·0001) in terms of total gas and CH4 productions. Fermentation kinetics (i.e. parameters describing the shape of the gas production curve and half-time gas production) for CH4 production were influenced by accession (P ≤ 0·038), but not by PEG. Accession, PEG and time affected (P < 0·001) CH4 production, but accession and PEG interaction showed only a tendency (P = 0·08). Increase in CH4 due to PEG addition was not related to CT content. Further analysis of the relationships among multiple traits (nutritional composition, CT structure and CH4 production) using principal component analysis (PCA) based on optimally weighted variables revealed differences among accessions. The first two principal component axes, PC1 (57·6%) and PC2 (18·4%), explained 76·0% of the total variation among accessions. Loading of biplots derived from both PCAs made it possible to establish a relationship between the ratio of prodelphinidin:procyanidin (PD:PC) tannins and CH4 production in some accessions. The PD:PC ratio seems to be an important source of variation that is negatively related to CH4 production. These results suggested that sainfoin accessions collected from across the world exhibited substantial variation in terms of their effects on rumen in vitro CH4 production, revealing some promising accessions for future investigations.

Functional identity and diversity of animals predict ecosystem functioning better than species-based indices

Drastic biodiversity declines have raised concerns about the deterioration of ecosystem functions and have motivated much recent research on the relationship between species diversity and ecosystem functioning. A functional trait framework has been proposed to improve the mechanistic understanding of this relationship, but this has rarely been tested for organisms other than plants. We analysed eight datasets, including five animal groups, to examine how well a trait-based approach, compared with a more traditional taxonomic approach, predicts seven ecosystem functions below- and above-ground. Trait-based indices consistently provided greater explanatory power than species richness or abundance. The frequency distributions of single or multiple traits in the community were the best predictors of ecosystem functioning. This implies that the ecosystem functions we investigated were underpinned by the combination of trait identities (i.e. single-trait indices) and trait complementarity (i.e. multi-trait indices) in the communities. Our study provides new insights into the general mechanisms that link biodiversity to ecosystem functioning in natural animal communities and suggests that the observed responses were due to the identity and dominance patterns of the trait composition rather than the number or abundance of species per se.

Genetic associations of visual scores with subsequent rebreeding and days to first calving in Nellore cattle

Records of Nellore animals born from 1990 to 2006 were used to estimate genetic correlations of visual scores at yearling (conformation, C; finishing precocity, P; and muscling, M) with primiparous subsequent rebreeding (SR) and days to first calving (DC), because the magnitude of these associations is still unknown. Genetic parameters were estimated by multiple-traits Bayesian analysis, using a nonlinear (threshold) animal models for visual scores and SR and a linear animal models for weaning weight (WW) and DC. WW was included in the analysis to account for the effects of sequential selection. The posterior means of heritabilities estimated for C, P, M, SR and DC were 0.24 +/- 0.01, 0.31 +/- 0.01, 0.30 +/- 0.01, 0.18 +/- 0.02 and 0.06 +/- 0.02, respectively. The posterior means of genetic correlations estimated between SR and visual scores were low and positive, with values of 0.09 +/- 0.02 (C), 0.19 +/- 0.03 (P) and 0.18 +/- 0.05 (M). on the other hand, negative genetic correlations were found between DC and C (-0.11 +/- 0.09), P (-0.19 +/- 0.09) and M (-0.16 +/- 0.09). The primiparous rebreeding trait has genetic variability in Nellore cattle. The genetic correlations between visual scores, and SR and DC were low and favourable. The genetic changes in C, P and M were 0.02, 0.03 and 0.03/year, respectively. For SR and DC, genetic trends were 0.01/year and -0.01 days/year, respectively, indicating that the increase in genetic merit for reproductive traits was small over time. Direct selection for visual scores together with female reproductive traits is recommended to increase the fertility of beef cows.

Test-day model for milk yield of dairy buffaloes in colombia

The test-day model is the preferred method for genetic evaluations in dairy cattle. For this study, 28372 test-day records of 1220 lactations from 1997 to 2009 were used. The (co)variance components for milk in test-day were estimated using a Uni and multiple-traits repeated animal model with the Restricted Maximum Likelihood method (REML). The Contemporary Group (herd, year, and season of parity) and the age of parity (linear and quadratic) fixed effects, and the additive genetic, permanent environmental, and residual random effects were included in the model. The heritabilities ranged between 0.06 and 0.45 during lactation. The genetic correlations were greater than 0.93. In conclusion, the test-day model is appropriate for the genetic evaluation of dairy buffaloes in Colombia.

Estimativa de componentes de (co)variância de características de crescimento na raça Brahman utilizando inferência bayesiana

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estimativas de componentes de (CO) variância de características de crescimentos na raça Nelore, utilizando inferência bayesiana

Pós-graduação em Genética e Melhoramento Animal - FCAV

Abordagem multivariada na seleção de progênies de soja superiores e portadoras do gene RR

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

Estimativas de parâmetros genéticos em linhagens avançadas de soja

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

Prospecting major genes in dairy buffaloes

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