Mixed model formulations for multi-environment trials

When studying genotype X environment interaction in multi-environment trials, plant breeders and geneticists often consider one of the effects, environments or genotypes, to be fixed and the other to be random. However, there are two main formulations for variance component estimation for the mixed model situation, referred to as the unconstrained-parameters (UP) and constrained-parameters (CP) formulations. These formulations give different estimates of genetic correlation and heritability as well as different tests of significance for the random effects factor. The definition of main effects and interactions and the consequences of such definitions should be clearly understood, and the selected formulation should be consistent for both fixed and random effects. A discussion of the practical outcomes of using the two formulations in the analysis of balanced data from multi-environment trials is presented. It is recommended that the CP formulation be used because of the meaning of its parameters and the corresponding variance components. When managed (fixed) environments are considered, users will have more confidence in prediction for them but will not be overconfident in prediction in the target (random) environments. Genetic gain (predicted response to selection in the target environments from the managed environments) is independent of formulation.

Phylogeography of the scarab beetle Antitrogus parvulus Britton (Coleoptera : Scarabaeidae) in south-eastern Queensland, Australia

This study surveys the population genetic structure of Childers canegrub, Antitrogus parvulus, to elucidate its population dynamics and gene flow. Antitrogus parvulus is a pest of sugarcane in the Bundaberg region and this knowledge can be used to optimise integrated pest management practices. Here, base-pair differences in the cytochrome oxidase II gene (COII) were used to characterise haplotypic diversity, infer levels of gene flow, and phylogenetic relationships of alleles and their phylogeographical structure. There were 28 unique haplotypes among the 70 sequenced individuals from the seven locations. All three variance components (among regions, among populations, within populations) are highly significant, with highest genetic diversity among regions and lowest among populations within regions. A positive correlation between migration rates and geographical distance and significant phylogeographical structure between four main geographical regions. The main implication of these findings for pest management is that if a grower can eliminate an existing infestation within a field, then reinvasion will be slow and further outbreaks within that field are unlikely to occur. The low dispersal ability of females also means that any resistance to insecticides that develops is likely to remain localised, but will rapidly become dominant within the affected population.

A genome scan for epidermal skin pattern in adolescent twins reveals suggestive linkage on 12p13.31

Deterioration in stratum corneum reticular patterning (skin pattern or skin wrinkling) has been associated with increased rates of solar keratoses and skin cancer. A previous analysis of data from the twin sample used in this investigation has shown that 86% of the variation in skin pattern is genetic at age 12 and 62% in an adult sample (mean age 47.5). Variation due to genetic influences is likely to be influenced by more than one locus. Here, we present results of a genome-wide linkage scan of skin pattern in adolescent twins and siblings from 428 nuclear twin families. Sib-pair linkage analysis was performed on skin pattern data collected from twins at age 12 (378 informative families) and 14 (316 families). Suggestive linkage was found at marker D12S397 (12p13.31, logarithm of the odds (lod) 1.94), when the effect of the trait locus was modelled to influence the skin pattern equally at both ages 12 and 14. In the same analysis, a peak was seen at 4q23 with a lod score of 1.55. A possible candidate for the peak at 12p13.31 is the protease inhibitor, alpha-2-macroglobulin.

Genome-wide scan of IQ finds significant linkage to a quantitative trait locus on 2q

A genome-wide linkage scan of 795 microsatellite markers (761 autosomal, 34 X chromosome) was performed on Multidimensional Aptitude Battery subtests and verbal, performance and full scale scores, the WAIS-R Digit Symbol subtest, and two word-recognition tests (Schonell Graded Word Reading Test, Cambridge Contextual Reading Test) highly predictive of IQ. The sample included 361 families comprising 2-5 siblings who ranged in age from 15.7 to 22.2 years; genotype, but not phenotype, data were available for 81% of parents. A variance components analysis which controlled for age and sex effects showed significant linkage for the Cambridge reading test and performance IQ to the same region on chromosome 2, with respective LOD scores of 4.15 and 3.68. Suggestive linkage (LOD score > 2.2) for various measures was further supported on chromosomes 6, 7, 11, 14, 21 and 22. Where location of linkage peaks converged for IQ subtests within the same scale, the overall scale score provided increased evidence for linkage to that region over any individual subtest. Association studies of candidate genes, particularly those involved in neural transmission and development, will be directed to genes located under the linkage peaks identified in this study.

Familial aggregation of human infection with Schistosoma japonicum in the Poyang Lake region, China

Despite the success of extensive control measures that have been implemented in China for over 50 years, the number of individuals infected with Schistosoma japonicum remains high in the existing endemic areas. A variance components analysis was undertaken to estimate the heritable and environmental components that contribute to S. japonicum infection in the Poyang Lake region of Jiangxi Province, PR China. The total target population was 3148 from four separate administrative villages. Two thousand seven hundred and five of these comprised 400 families ranging in size from 3 to 188. After adjustments were made for gender, water contact and past history of having had schistosomiasis, the heritable component was estimated to account for as much as 58% of the phenotype variation under the polygenic model. Household was not shown to be an important environmental factor. Incorporating village effects indicated that the results were valid for the total population. We conclude that genetic heritability in this region is high and plays an important role in determining risk of infection with S. japonicum. (c) 2005 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Global adaptation of spring bread and durum wheat lines near-isogenic for major reduced height genes

The effect of major dwarfing genes, Rht-B1 and Rht-D1, in bread (Triticum aestivum L.) and durum (Triticum turgidum L. var. durum) wheats varies with environment. Six reduced-height near-isogenic spring wheat lines, included in the International Adaptation Trial (IAT), were grown in 81 trials around the world. Of the 56 IAT trials yielding > 3 Mg ha(-1), the mean yield of semidwarfs was significantly greater than tails in 54% of trials; in the 27 trials yielding < 3 Mg ha-1, semidwarfs were superior in only 24%. Sixteen pairs of semidwarf-tall near-isolines were grown in six managed drought environment trials (DETs) in northwestern Mexico. In these trials, semidwarfs outyielded talls in all but the most droughted environment (2.5 Mg ha(-1)). The effect of the height alleles varied with genetic background and environment. For both yield and height, variance components for allele and environment by allele interaction were larger than those for genetic background and genetic background by environment. Pattern analysis showed that tall and semidwarf lines had similar adaptation to stressed environments (< 2.8 Mg ha(-1), low rainfall), while semidwarfs yielded more in less stressed environments (> 4.3 Mg ha(-1), high rainfall). The best adapted near-isogenic pair had a Kauz background, where the tall was only 16% taller than the dwarf. In the Kauz-derived pair, the semidwarf outyielded the tall in only 13% of trials with no differences in low yielding trials. This supports the idea that '' short talls '' may be useful in marginal environments (yield < 3 Mg ha(-1)).

Multi-level zero-inflated Poisson regression modelling of correlated count data with excess zeros

Count data with excess zeros relative to a Poisson distribution are common in many biomedical applications. A popular approach to the analysis of such data is to use a zero-inflated Poisson (ZIP) regression model. Often, because of the hierarchical Study design or the data collection procedure, zero-inflation and lack of independence may occur simultaneously, which tender the standard ZIP model inadequate. To account for the preponderance of zero counts and the inherent correlation of observations, a class of multi-level ZIP regression model with random effects is presented. Model fitting is facilitated using an expectation-maximization algorithm, whereas variance components are estimated via residual maximum likelihood estimating equations. A score test for zero-inflation is also presented. The multi-level ZIP model is then generalized to cope with a more complex correlation structure. Application to the analysis of correlated count data from a longitudinal infant feeding study illustrates the usefulness of the approach.

Linkage analyses of event-related potential slow wave phenotypes recorded in a working memory task

Working memory is an essential component of wide-ranging cognitive functions. It is a complex genetic trait probably influenced by numerous genes that individually have only a small influence. These genes may have an amplified influence on phenotypes closer to the gene action. In this study, event-related potential (ERP) phenotypes recorded during a working-memory task were collected from 656 adolescents from 299 families for whom genotypes were available. Univariate linkage analyses using the MERLIN variance-components method were conducted on slow wave phenotypes recorded at multiple sites while participants were required to remember the location of a target. Suggestive linkage (LOD > 2.2) was found on chromosomes 4, 5, 6, 10, 17, and 20. After correcting for multiple testing, suggestive linkage remained on chromosome 10. Empirical thresholds were computed for the most promising phenotypes. Those on chromosome 10 remained suggestive. A number of genes reported to regulate neural differentiation and function (i.e. NRP1, ANK3, and CHAT) were found under these linkage peaks and may influence the levels of neural activity occurring in individuals participating in a spatial working-memory task.

Comparison of Multivariate and Univariate Models for Genetic Evaluation of Milk Yield based on Test Day Data

2000 Mathematics Subject Classification: 62H12, 62P99

HERITABILITY FOR RESISTANCE TO ROSETTE DISEASE IN EXOTIC VALENCIA GROUNDNUTS

Groundnut rosette disease (GRD) is the most destructive virus disease of Valencia groundnuts ( Arachis hypogaea L.) in sub-Saharan Africa. Cultural, biological and chemical control measures have received limited success due to small scale farmers’ inability to use them. Use of host plant resistance provides the most effective and economically viable management option for the resource poor farmers. This study was conducted to determine heritability for resistance to GRD in Valencia groundnuts. Six crosses; Valencia C (P1) × ICGV-SM 90704 (P2), Valencia C (P1) × ICGV-SM 96801(P2), Valencia C (P1) × ICGV-SM 99566 (P2), NuMex-M3 (P1) × ICGV-SM 90704 (P2), NuMex-M3 × ICGV-SM 96801 (P2), and NuMex-M3 (P1) × ICGV-SM 99566 (P2), were made to generate F1, F2, BC1P1 and BC1P2 populations. Data on GRD severity were collected on a 1-9 score scale. Genetic Advance as a percentage of the mean (GAM) and heritability were estimated using variance components. Phenotypic Coefficient of Variation (PCV) and Genotypic Coefficient of Variation (GCV) estimates were high (20.04-70.1%) in the six crosses, except for Valencia C × ICGV-SM 96801(18.1%) and NuMex-M3 × ICGV-SM 96801(17.1%), which exhibited moderate GCV values. Broad and narrow sense heritability estimates for GRD disease score ranged from 64.1 to 73.7% and 31 to 41.9%, respectively, in all the crosses. GAM was high in all the crosses (21-50.7%), except for Valencia C x ICGV-SM 96801 (14.67), M3 x ICGV-SM 99566 (18%) and NuMex-M3 x ICGV-SM 96801 (13.5%) crosses that exhibited moderate GAM. The study revealed the presence of variability of GRD resistance, implying that genetic improvement of these exotic materials is possible.

Modelos linear e não linear em análises genéticas para sobrevivência de crias de ovinos da raça Santa Inês.

Resumo: Registros de sobrevivência do nascimento ao desmame de 3846 crias de ovinos da raça Santa Inês foram analisados por modelos de reprodutor linear e não linear (modelo de limiar), para estimar componentes de variância e herdabilidade. Os modelos usados para sobrevivência, analisada como característica da cria, incluíram os efeitos fixos de sexo, da combinação tipo de nascimento-criação da cria e da idade da ovelha ao parto, efeito da covariável peso da cria ao nascer e efeitos aleatórios de reprodutor, da classe rebanho-ano-estação e do resíduo. Componentes de variância para o modelo linear foram estimados pelo método da máxima verossimilhança restrita (REML) e para o modelo não linear por uma aproximação da máxima verossimilhança marginal (MML), pelo programa CMMAT2. O coeficiente de herdabilidade (h2) estimado pelo modelo de limiar foi de 0,29, e pelo modelo linear, 0,14. A correlação de ordem de Spearman entre as capacidades de transmissão dos reprodutores, com base nos dois modelos foi de 0,96. As estimativas de h2 obtidas indicam a possibilidade de se obter, por seleção, ganho genético para sobrevivência. [Linear and nonlinear models in genetic analyses of lamb survival in the Santa Inês hair sheep breed]. Abstract: Records of 3,846 lambs survival from birth to weaning of Santa Inês hair sheep breed, were analyzed by linear and non linear sire models (threshold model) to estimate variance components and heritability (h2). The models that were used to analyze survival, considered in this study as a lamb trait, included the fixed effects of sex of the lamb, combination of type of birth-rearing of lamb, and age of ewe, birth weight of lamb as covariate, and random effects of sire, herd-year-season and residual. Variance components were obtained using restricted maximum likelihood (REML), in linear model and marginal maximum likelihood in threshold model through CMMAT2 program. Estimate of heritability (h2) obtained by threshold model was 0.29 and by linear model was 0.14. Rank correlation of Spearman, between sire solutions based on the two models was 0.96. The obtained estimates in this study indicate that it is possible to acquire genetic gain to survival by selection.

Estimates of genetic parameters using RELM/BLUP for intrapopulational genetic breeding of Platonia insignis Mart.

The aim of this study was to estimate genetic parameters to support the selection of bacuri progenies for a first cycle of recurrent selection, using the REML/BLUP (restricted maximum likelihood/best linear unbiased prediction) procedure to estimate the variance components and genotypic values. Were evaluated twelve variables in a total of 210 fruits from 39 different seed trees, from a field trial with an experimental design of incomplete blocks with clonal replies among subplots. The three variables related with the fruit development (weight, diameter, length) showed strong correlation, and where fruit length showed higher heritability and potential to be used for indirect selection. Among the 39 progenies evaluated in this study, five present potential to compose the next cycle of recurrent selection, due they hold good selection differential either to agrotechnological variables as to development of bacuri fruit.

Effect of inclusion or non-inclusion of short lactations and cow and/or dam genetic group on genetic evaluation of Girolando dairy cattle.

The objective of this study was to evaluate the effects of inclusion or non-inclusion of short lactations and cow (CGG) and/or dam (DGG) genetic group on the genetic evaluation of 305-day milk yield (MY305), age at first calving (AFC), and first calving interval (FCI) of Girolando cows. Covariance components were estimated by the restricted maximum likelihood method in an animal model of single trait analyses. The heritability estimates for MY305, AFC, and FCI ranged from 0.23 to 0.29, 0.40 to 0.44, and 0.13 to 0.14, respectively, when short lactations were not included, and from 0.23 to 0.28, 0.39 to 0.43, and 0.13 to 0.14, respectively, when short lactations were included. The inclusion of short lactations caused little variation in the variance components and heritability estimates of traits, but their non-inclusion resulted in the re-ranking of animals. Models with CGG or DGG fixed effects had higher heritability estimates for all traits compared with models that consider these two effects simultaneously. We recommend using the model with fixed effects of CGG and inclusion of short lactations for the genetic evaluation of Girolando cattle.

Components of variance method and partitioning method of genetic analysis applied to weight per fruit of tomato hybrid and parental populations /

Includes bibliographical references (p. 56-57).