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.

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.

Simultaneous prediction of acidity parameters (pH and titratable acidity) in Kefir using near infrared reflectance spectroscopy.

Acidity in terms of pH and titratable acids influences the texture and flavour of fermented dairy products, such as Kefir. However, the methods for determining pH and titratable acidity (TA) are time consuming. Near infrared (NIR) spectroscopy is a non-destructive method, which simultaneously predicts multiple traits from a single scan and can be used to predict pH and TA. The best pH NIR calibration model was obtained with no spectral pre-treatment applied, whereas smoothing was found to be the best pre-treatment to develop the TA calibration model. Using cross-validation, the prediction results were found acceptable for both pH and TA. With external validation, similar results were found for pH and TA, and both models were found to be acceptable for screening purposes.

A genome-wide linkage scan provides evidence for both new and previously reported loci influencing common migraine

Latent class analysis was performed on migraine symptom data collected in a Dutch population sample (N = 12,210, 59% female) in order to obtain empirical groupings of individuals suffering from symptoms of migraine headache. Based on these heritable groupings (h(2) = 0.49, 95% CI: 0.41-0.57) individuals were classified as affected (migrainous headache) or unaffected. Genome-wide linkage analysis was performed using genotype data from 105 families with at least 2 affected siblings. In addition to this primary phenotype, linkage analyses were performed for the individual migraine symptoms. Significance levels, corrected for the analysis of multiple traits, were determined empirically via a novel simulation approach. Suggestive linkage for migrainous headache was found on chromosomes 1 (LOD = 1.63; pointwise P = 0.0031), 13 (LOD = 1.63; P = 0.0031), and 20 (LOD = 1.85; P = 0.0018). Interestingly, the chromosome 1 peak was located close to the ATP1A2 gene, associated with familial hemiplegic migraine type 2 (FHM2). Individual symptom analysis produced a LOD score of 1.97 (P = 0.0013) on chromosome 5 (photo/phonophobia), a LOD score of 2.13 (P = 0.0009) on chromosome 10 (moderate/severe pain intensity) and a near significant LOD score of 3.31 (P = 0.00005) on chromosome 13 (pulsating headache). These peaks were all located near regions previously reported in migraine linkage studies. Our results provide important replication and support for the presence of migraine susceptibility genes within these regions, and further support the utility of an LCA-based phenotyping approach and analysis of individual symptoms in migraine genetic research. Additionally, our novel "2-step" analysis and simulation approach provides a powerful means to investigate linkage to individual trait components.

Cannibalism can drive the evolution of behavioural phase polyphenism in locusts

During outbreaks, locust swarms can contain millions of insects travelling thousands of kilometers while devastating vegetation and crops. Such large-scale spatial organization is preceded locally by a dramatic density-dependent phenotypic transition in multiple traits. Behaviourally, low-density solitarious individuals avoid contact with one another; above a critical local density, they undergo a rapid behavioural transition to the gregarious phase whereby they exhibit mutual attraction. Although proximate causes of this phase polyphenism have been widely studied, the ultimate driving factors remain unclear. Using an individual-based evolutionary model, we reveal that cannibalism, a striking feature of locust ecology, could lead to the evolution of density-dependent behavioural phase-change in juvenile locusts. We show that this behavioural strategy minimizes risk associated with cannibalistic interactions and may account for the empirically observed persistence of locust groups during outbreaks. Our results provide a parsimonious explanation for the evolution of behavioural plasticity in locusts.

Manuel Bandeira nas cartas e nas crônicas: percurso literário e arte poética

Este trabalho é um exercício reflexivo sobre o percurso literário e a arte poética de Manuel Carneiro de Souza Bandeira Filho nos decênios de 1920 e 1930. Na leitura da correspondência trocada por ele e Mário Raul de Morais Andrade, e de seus textos escritos para a imprensa, buscamos analisar as múltiplas feições assumidas pelo poeta pernambucano e procuramos desentranhar traços gerais de sua concepção de vida e de poesia. Deste modo, na presente dissertação, o gênero epistolar e a crônica são considerados como notáveis objetos de estudo para a literatura brasileira. A correspondência com Mário de Andrade adquire relevo na produção intelectual de Manuel Bandeira, pois ela possibilita um maior entendimento do poeta a partir de sua escrita de si, assim como permite observar a preocupação do escritor com a memória da cultura brasileira. Da mesma forma, a crônica torna-se importante porque traz o testemunho do cronista sobre o tempo circundante. Ao percorremos uma parte significativa da prosa de Manuel Bandeira, estabelecemos o cotejo com o seu texto memorialístico Itinerário de Pasárgada, rastreando, nessas fontes, os escritos que denunciam o seu posicionamento em relação às inovações da arte moderna

Padrões de assimetria associados à ecologia de cinco espécies de lagartos tropicais Cnemidophorus (Squamata, Teiidae) em áreas de restinga ao longo da costa leste brasileira

O estudo da instabilidade do desenvolvimento (ID) em organismos pode ser útil para identificar os primeiros sinais de estresse ambiental. Alterações na ID também influenciam a aptidão dos organismos. A variação entre a simetria dos lados esquerdo e direito, manifestada como assimetria flutuante (AF), é um indicador da ID. Ao quantificarmos a literatura científica referente aos estudos de assimetria em répteis identificamos as seguintes tendências: perturbações ambientais aumentam a AF nas populações, fêmeas tendem a escolher machos com características sexuais simétricas, a instabilidade térmica durante o desenvolvimento do ovo pode produzir indivíduos com maior assimetria, a ocorrência de lesões é mais frequente em indivíduos assimétricos para um determinado lado e há uma correlação negativa entre a assimetria e o desempenho locomotor. Considerando os nossos critério de busca, não encontramos estudos sobre assimetria em répteis na América do Sul. Neste estudo, utilizamos diferentes caracteres merísticos e métricos bilaterais para identificar a AF em diferentes populações de cinco espécies de lagartos do gênero Cnemidophorus (C. abaetensis, C. lacertoides, C. littoralis, C. nativo, C. ocellifer) em 15 áreas de restinga ao longo da costa leste do Brasil. A AF foi maior em populações de Cnemidophorus com maior densidade, menor em populações de restingas com maior degradação ambiental e a sua relação com a concentração de mercúrio variou dependendo da espécie. Provavelmente os agentes estressores ambientais avaliados atuam com intensidades diferentes. O nível de degradação ambiental e a concentração de mercúrio causariam a deterioração precoce dos indivíduos com maiores valores de AF. Isso resulta em uma relação inversa entre esses estresses e a AF das populações. Indivíduos com vestígios de autotomia caudal tiveram maiores valores de AF, o que pode ser um indicativo de que indivíduos com menor aptidão efetuam mais autotomia caudal do que aqueles com maior aptidão. Lagartos com evidências de autotomia também tiveram maior tamanho corpóreo, provavelmente devido a estes estarem expostos a mais tempo aos riscos de predação. Não encontramos diferenciação entre as espécies na frequência de autotomia caudal sendo esta talvez melhor explicada por fatores ambientais como a intensidade de predação. Os índices de AF para caracteres múltiplos se mostraram mais eficientes do que a AF de caracteres individuais nas comparações entre indivíduos vivendo em populações com diferentes níveis de estresse e entre indivíduos com e sem autotomia caudal

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.

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)