Search

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is possible with models for inherited variability but may require a two‐step analysis. [ABSTRACT FROM AUTHOR]

The cerebellum has a highly conserved neural structure across species but varies widely in size. The wide variation in cerebellar size (both absolute and in proportion to the rest of the brain) among species and populations suggests that functional specialization is linked to its size. There is increasing recognition that the cerebellum contributes to cognitive processing and emotional control in addition to its role in motor coordination. However, to what extent cerebellum size reflects variation in these behavioral processes within species remains largely unknown. By using a unique intercross chicken population based on parental lines with high divergence in cerebellum size, we compared the behavior of individuals repeatedly exposed to the same fear test (emergence test) early in life and after sexual maturity (eight trials per age group) with proportional cerebellum size and cerebellum neural density. While proportional cerebellum size did not predict the initial fear response of the individuals (trial 1), it did increasingly predict adult individuals response as the trials progressed. Our results suggest that proportional cerebellum size does not necessarily predict an individual's fear response, but rather the habituation process to a fearful stimulus. Cerebellum neuronal density did not predict fear behavior in the individuals which suggests that these effects do not result from changes in neuronal density but due to other variables linked to proportional cerebellum size which might underlie fear habituation. [ABSTRACT FROM AUTHOR]

Introduction: Although breast and prostate cancers arise in different organs and are more frequent in the opposite sex, multiple studies have reported an association between their family history. Analysis of single nucleotide polymorphism data, based on distant relatives, has revealed a small positive genetic correlation between these cancers explained by common variants. The estimate of genetic correlation based on close relatives reveals the extent to which shared genetic risks are explained by both common and rare variants. This estimate is unknown for breast and prostate cancer. Method: We estimated the relative risks, heritability, and genetic correlation of breast cancer and prostate cancer based on the Minnesota Breast and Prostate Cancer Study, a family study of 141 families ascertained for breast cancer. Results: Heritability of breast cancer was 0.34 (95% credible interval: 0.23–0.49) and 0.65 (95% credible interval: 0.36–0.97) for prostate cancer, and the genetic correlation was 0.23. In terms of odds ratios, these values correspond to a 1.3 times higher odds of breast cancer among probands, given that the brother has prostate cancer. Conclusion: This study shows the inherent relation between prostate cancer and breast cancer; an incident of one in a family increases the risk of developing the other. The large difference between estimates of genetic correlation from distant and close relatives, if replicated, suggests that rare variants contribute to the shared genetic risk of breast and prostate cancer. However, the difference could stem from genotype-by-family effects shared between the two types of cancers. [ABSTRACT FROM AUTHOR]

We develop a Bayesian model (BayesRR-RC) that provides robust SNP-heritability estimation, an alternative to marker discovery, and accurate genomic prediction, taking 22 seconds per iteration to estimate 8.4 million SNP-effects and 78 SNP-heritability parameters in the UK Biobank. We find that only ≤10% of the genetic variation captured for height, body mass index, cardiovascular disease, and type 2 diabetes is attributable to proximal regulatory regions within 10kb upstream of genes, while 12-25% is attributed to coding regions, 32–44% to introns, and 22-28% to distal 10-500kb upstream regions. Up to 24% of all cis and coding regions of each chromosome are associated with each trait, with over 3,100 independent exonic and intronic regions and over 5,400 independent regulatory regions having ≥95% probability of contributing ≥0.001% to the genetic variance of these four traits. Our open-source software (GMRM) provides a scalable alternative to current approaches for biobank data. Improving inference in large-scale genetic data linked to electronic medical record data requires the development of novel computationally efficient regression methods. Here, the authors develop a Bayesian approach for association analyses to improve SNP-heritability estimation, discovery, fine-mapping and genomic prediction. [ABSTRACT FROM AUTHOR]

Statistical modeling with R: a dual frequentist and Bayesian approach for life scientists by Inchausti, Pablo. Oxford University Press. 2 November 2022. 480 p. £75. ISBN: 9780192859013By Lars RönnegårdReported by Author [Extracted from the article]

Background and Objectives: Levodopa concentration in patients with Parkinson's disease is frequently modelled with ordinary differential equations (ODEs). Here, we investigate a pharmacokinetic model of plasma levodopa concentration in patients with Parkinson's disease by introducing stochasticity to separate the intra-individual variability into measurement and system noise, and to account for auto-correlated errors. We also investigate whether the induced stochasticity provides a better fit than the ODE approach. Methods: In this study, a system noise variable is added to the pharmacokinetic model for duodenal levodopa/carbidopa gel (LCIG) infusion described by three ODEs through a standard Wiener process, leading to a stochastic differential equations (SDE) model. The R package population stochastic modelling (PSM) was used for model fitting with data from previous studies for modelling plasma levodopa concentration and parameter estimation. First, the diffusion scale parameter (σw), measurement noise variance, and bioavailability are estimated with the SDE model. Second, σw is fixed to certain values from 0 to 1 and bioavailability is estimated. Cross-validation was performed to compare the average root mean square errors (RMSE) of predicted plasma levodopa concentration. Results: Both the ODE and the SDE models estimated bioavailability to be approximately 75%. The SDE model converged at different values of σw that were significantly different from zero. The average RMSE for the ODE model was 0.313, and the lowest average RMSE for the SDE model was 0.297 when σw was fixed to 0.9, and these two values are significantly different. Conclusions: The SDE model provided a better fit for LCIG plasma levodopa concentration by approximately 5.5% in terms of mean percentage change of RMSE. [ABSTRACT FROM AUTHOR]

The ability of a population to adapt to changes in their living conditions, whether in nature or captivity, often depends on polymorphisms in multiple genes across the genome. In-depth studies of such polygenic adaptations are difficult in natural populations, but can be approached using the resources provided by artificial selection experiments. Here, we dissect the genetic mechanisms involved in long-term selection responses of the Virginia chicken lines, populations that after 40 generations of divergent selection for 56-day body weight display a 9-fold difference in the selected trait. In the F15 generation of an intercross between the divergent lines, 20 loci explained >60% of the additive genetic variance for the selected trait. We focused particularly on fine-mapping seven major QTL that replicated in this population and found that only two fine-mapped to single, bi-allelic loci; the other five contained linked loci, multiple alleles or were epistatic. This detailed dissection of the polygenic adaptations in the Virginia lines provides a deeper understanding of the range of different genome-wide mechanisms that have been involved in these long-term selection responses. The results illustrate that the genetic architecture of a highly polygenic trait can involve a broad range of genetic mechanisms, and that this can be the case even in a small population bred from founders with limited genetic diversity. [ABSTRACT FROM AUTHOR]

Individuals involved in a social interaction exhibit different behavioral traits that, in combination, form the individual's behavioral responses. Selectively bred strains of silver foxes ( Vulpes vulpes) demonstrate markedly different behaviors in their response to humans. To identify the genetic basis of these behavioral differences we constructed a large F population including 537 individuals by cross-breeding tame and aggressive fox strains. 98 fox behavioral traits were recorded during social interaction with a human experimenter in a standard four-step test. Patterns of fox behaviors during the test were evaluated using principal component (PC) analysis. Genetic mapping identified eight unique significant and suggestive QTL. Mapping results for the PC phenotypes from different test steps showed little overlap suggesting that different QTL are involved in regulation of behaviors exhibited in different behavioral contexts. Many individual behavioral traits mapped to the same genomic regions as PC phenotypes. This provides additional information about specific behaviors regulated by these loci. Further, three pairs of epistatic loci were also identified for PC phenotypes suggesting more complex genetic architecture of the behavioral differences between the two strains than what has previously been observed. [ABSTRACT FROM AUTHOR]

The depredation of semi-domesticated reindeer by large carnivores reflects an important human–wildlife conflict in Fennoscandia. Recent studies have revealed that brown bears (Ursus arctos) may kill substantial numbers of reindeer calves (Rangifer tarandus tarandus) in forest areas in Sweden. Several authors have suggested that predation risk is an important driver of habitat selection in wild Rangifer populations where predation is a limiting factor, but little is known about these mechanisms in semi-domesticated populations. We examined the habitat selection of female reindeer in relation to spatial and temporal variations in brown bear predation risk on the reindeer calving grounds and evaluated the simultaneous responses of brown bears and reindeer to landscape characteristics. We used GPS data from 110 reindeer years (97 individuals) and 29 brown bear years (19 individuals), from two reindeer herding districts in the forest area of northern Sweden. Our results did not indicate that reindeer alter their behavior in response to spatiotemporal variation in brown bear predation risk, on the scale of the calving range. Instead, we suggest that spatiotemporal behavioral adjustments by brown bears were the main driver of prey-predator interactions in our study system. Contrasting responses by brown bears and reindeer to clear-cuts and young forest indicate that forestry can influence species interactions and possibly yield negative consequences for the reindeer herd. Even if clear-cuts may be beneficial in terms of calf survival, logging activity will eventually cause greater abundance of young regenerating forest, reducing available reindeer habitats and increasing habitat preferred by brown bears. Domestication may have made semi-domesticated reindeer in Fennoscandia less adapted to cope with predators. Areal restrictions, limiting the opportunity for dispersion and escape, possibly make the calves more susceptible to predation. Also, a generally higher population density in semi-domesticated herds compared to wild populations can make dispersion a less efficient strategy and the reindeer calves easier prey. Overall, the lack of ability of the reindeer females to reduce brown bear encounter risk on the scale of the calving range is probably an important reason for the high brown bear predation rates on reindeer calves documented in our study areas. [ABSTRACT FROM AUTHOR]

We analyze a real data set pertaining to reindeer fecal pellet-group counts obtained from a survey conducted in a forest area in northern Sweden. In the data set, over 70% of counts are zeros, and there is high spatial correlation. We use conditionally autoregressive random effects for modeling of spatial correlation in a Poisson generalized linear mixed model ( GLMM), quasi-Poisson hierarchical generalized linear model ( HGLM), zero-inflated Poisson ( ZIP), and hurdle models. The quasi-Poisson HGLM allows for both under- and overdispersion with excessive zeros, while the ZIP and hurdle models allow only for overdispersion. In analyzing the real data set, we see that the quasi-Poisson HGLMs can perform better than the other commonly used models, for example, ordinary Poisson HGLMs, spatial ZIP, and spatial hurdle models, and that the underdispersed Poisson HGLMs with spatial correlation fit the reindeer data best. We develop R codes for fitting these models using a unified algorithm for the HGLMs. Spatial count response with an extremely high proportion of zeros, and underdispersion can be successfully modeled using the quasi-Poisson HGLM with spatial random effects. [ABSTRACT FROM AUTHOR]

Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long-term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single-nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful., Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long-term study of collared flycatchers in Sweden., The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response., The method is available in the R package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long-term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations. [ABSTRACT FROM AUTHOR]

This study covers a period when society changed from a pre-industrial agricultural society to a postindustrial service-producing society. Parallel with this social transformation, major population changes took place. In this study, we analyse to what extent local population change is affected by neighbouring populations. To do this, we focused on the last 190 years of local population change that redistributed population in Sweden. We used literature to identify several different processes in the population redistribution. The different processes implied different spatial dependencies between local population change and the surrounding populations. The analysis is based on an unchanged historical parish division, and we used an index of local spatial correlation to describe different types of spatial dependencies that influenced the redistribution of the population. To control inherent time dependencies, we introduced a non-separable spatial-temporal correlation model into the analysis of population redistribution. Hereby, several different spatial dependencies could be simultaneously observed over time. The main conclusions are that while local population changes have been highly dependent on neighbouring populations in the 19th century, this spatial dependence became insignificant already when two parishes are separated by 5km in the late 20th century. It is argued that the only process that significantly redistributed the population at the end of the 20th century is the immigration to Sweden. [ABSTRACT FROM AUTHOR]

We present a new version (⩾ 2.0) of the hglm package for fitting hierarchical generalized linear models (HGLMs) with spatially correlated random effects. CAR() and SAR() families for conditional and simultaneous autoregressive random effects were implemented. Eigen decomposition of the matrix describing the spatial structure (e.g., the neighborhood matrix) was used to transform the CAR/SAR random effects into an independent, but heteroscedastic, Gaussian random effect. A linear predictor is fitted for the random effect variance to estimate the parameters in the CAR and SAR models. This gives a computationally efficient algorithm for moderately sized problems. [ABSTRACT FROM AUTHOR]

Context: Over the last decade, we have seen a massive increase in the construction of wind farms in northern Fennoscandia. Wind farms comprising hundreds of wind turbines are being built, with little knowledge of the possible cumulative adverse effects on the habitat use and migration of semi-domesticated free-ranging reindeer. Objectives: We assessed how reindeer responded to wind farm construction in an already fragmented landscape, with specific reference to the effects on use of movement corridors and reindeer habitat selection. Methods: We used GPS-data from reindeer during calving and post-calving in the Malå reindeer herding community in Sweden. We analysed data from the pre-development years compared to the construction years of two relatively small wind farms. Results: During construction of the wind farms, use of original migration routes and movement corridors within 2 km of development declined by 76 %. This decline in use corresponded to an increase in activity of the reindeer measured by increased step lengths within 0-5 km. The step length was highest nearest the development and declining with distance, as animals moved towards migration corridors and turned around or were observed in holding patterns while not crossing. During construction, reindeer avoided the wind farms at both regional and landscape scale of selection. Conclusions: The combined construction activities associated with even a few wind turbines combined with power lines and roads in or close to central movement corridors caused a reduction in the use of such corridors and grazing habitat and increased the fragmentation of the reindeer calving ranges. [ABSTRACT FROM AUTHOR]

Reindeer herding in Sweden is a form of pastoralism practised by the indigenous Sámi population. The economy is mainly based on meat production. Herd size is generally regulated by harvest in order not to overuse grazing ranges and keep a productive herd. Nonetheless, herd growth and room for harvest is currently small in many areas. Negative herd growth and low harvest rate were observed in one of two herds in a reindeer herding community in Central Sweden. The herds (A and B) used the same ranges from April until the autumn gathering in October–December, but were separated on different ranges over winter. Analyses of capture-recapture for 723 adult female reindeer over five years (2007–2012) revealed high annual losses (7.1% and 18.4%, for herd A and B respectively). A continuing decline in the total reindeer number in herd B demonstrated an inability to maintain the herd size in spite of a very small harvest. An estimated breakpoint for when herd size cannot be kept stable confirmed that the observed female mortality rate in herd B represented a state of herd collapse. Lower calving success in herd B compared to A indicated differences in winter foraging conditions. However, we found only minor differences in animal body condition between the herds in autumn. We found no evidence that a lower autumn body mass generally increased the risk for a female of dying from one autumn to the next. We conclude that the prime driver of the on-going collapse of herd B is not high animal density or poor body condition. Accidents or disease seem unlikely as major causes of mortality. Predation, primarily by lynx and wolverine, appears to be the most plausible reason for the high female mortality and state of collapse in the studied reindeer herding community. [ABSTRACT FROM AUTHOR]

The rapid development of new technologies for large-scale analysis of genetic variation in the genomes of individuals and populations has presented statistical geneticists with a grand challenge to develop efficient methods for identifying the small proportion of all identified genetic polymorphisms that have effects on traits of interest. To address such a 'large p small n' problem, we have developed a heteroscedastic effects model (HEM) that has been shown to be powerful in high-throughput genetic analyses. Here, we describe how this whole-genome model can also be utilized in chemometric analysis. As a proof of concept, we use HEM to predict analyte concentrations in silage using Fourier transform infrared spectroscopy signals. The results show that HEM often outperforms the classic methods and in addition to this presents a substantial computational advantage in the analyses of such high-dimensional data. The results thus show the value of taking an interdisciplinary approach to chemometric analysis and indicate that large-scale genomic models can be a promising new approach for chemometric analysis that deserve to be evaluated more by experts in the field. The software used for our analyses is freely available as an R package at . Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Background: Canalization is defined as the stability of a genotype against minor variations in both environment and genetics. Genetic variation in degree of canalization causes heterogeneity of within-family variance. The aims of this study are twofold: (1) quantify genetic heterogeneity of (within-family) residual variance in Atlantic salmon and (2) test whether the observed heterogeneity of (within-family) residual variance can be explained by simple scaling effects. Results: Analysis of body weight in Atlantic salmon using a double hierarchical generalized linear model (DHGLM) revealed substantial heterogeneity of within-family variance. The 95% prediction interval for within-family variance ranged from ~0.4 to 1.2 kg2, implying that the within-family variance of the most extreme high families is expected to be approximately three times larger than the extreme low families. For cross-sectional data, DHGLM with an animal mean sub-model resulted in severe bias, while a corresponding sire-dam model was appropriate. Heterogeneity of variance was not sensitive to Box-Cox transformations of phenotypes, which implies that heterogeneity of variance exists beyond what would be expected from simple scaling effects. Conclusions: Substantial heterogeneity of within-family variance was found for body weight in Atlantic salmon. A tendency towards higher variance with higher means (scaling effects) was observed, but heterogeneity of within-family variance existed beyond what could be explained by simple scaling effects. For cross-sectional data, using the animal mean sub-model in the DHGLM resulted in biased estimates of variance components, which differed substantially both from a standard linear mean animal model and a sire-dam DHGLM model. Although genetic differences in canalization were observed, selection for increased canalization is difficult, because there is limited individual information for the variance sub-model, especially when based on cross-sectional data. Furthermore, potential macro-environmental changes (diet, climatic region, etc.) may make genetic heterogeneity of variance a less stable trait over time and space. [ABSTRACT FROM AUTHOR]

Background: Genetic variation for environmental sensitivity indicates that animals are genetically different in their response to environmental factors. Environmental factors are either identifiable (e.g. temperature) and called macro-environmental or unknown and called micro-environmental. The objectives of this study were to develop a statistical method to estimate genetic parameters for macro- and micro-environmental sensitivities simultaneously, to investigate bias and precision of resulting estimates of genetic parameters and to develop and evaluate use of Akaike's information criterion using h-likelihood to select the best fitting model. Methods: We assumed that genetic variation in macro- and micro-environmental sensitivities is expressed as genetic variance in the slope of a linear reaction norm and environmental variance, respectively. A reaction norm model to estimate genetic variance for macro-environmental sensitivity was combined with a structural model for residual variance to estimate genetic variance for micro-environmental sensitivity using a double hierarchical generalized linear model in ASReml. Akaike's information criterion was constructed as model selection criterion using approximated h-likelihood. Populations of sires with large half-sib offspring groups were simulated to investigate bias and precision of estimated genetic parameters. Results: Designs with 100 sires, each with at least 100 offspring, are required to have standard deviations of estimated variances lower than 50% of the true value. When the number of offspring increased, standard deviations of estimates across replicates decreased substantially, especially for genetic variances of macro- and micro-environmental sensitivities. Standard deviations of estimated genetic correlations across replicates were quite large (between 0.1 and 0.4), especially when sires had few offspring. Practically, no bias was observed for estimates of any of the parameters. Using Akaike's information criterion the true genetic model was selected as the best statistical model in at least 90% of 100 replicates when the number of offspring per sire was 100. Application of the model to lactation milk yield in dairy cattle showed that genetic variance for micro- and macro-environmental sensitivities existed. Conclusion: The algorithm and model selection criterion presented here can contribute to better understand genetic control of macro- and micro-environmental sensitivities. Designs or datasets should have at least 100 sires each with 100 offspring. [ABSTRACT FROM AUTHOR]

The phenotypic effect of a gene is normally described by the mean-difference between alternative genotypes. A gene may, however, also influence the phenotype by causing a difference in variance between genotypes. Here, we reanalyze a publicly available Arabidopsis thaliana dataset [1] and show that genetic variance heterogeneity appears to be as common as normal additive effects on a genomewide scale. The study also develops theory to estimate the contributions of variance differences between genotypes to the phenotypic variance, and this is used to show that individual loci can explain more than 20% of the phenotypic variance. Two well-studied systems, cellular control of molybdenum level by the ion-transporter MOT1 and flowering-time regulation by the FRI-FLC expression network, and a novel association for Leaf serration are used to illustrate the contribution of major individual loci, expression pathways, and gene-by-environment interactions to the genetic variance heterogeneity. [ABSTRACT FROM AUTHOR]

A number of recent works have introduced statistical methods for detecting genetic loci that affect phenotypic variability, which we refer to as variability-controlling quantitative trait loci (vQTL). These are genetic variants whose allelic state predicts how much phenotype values will vary about their expected means. Such loci are of great potential interest in both human and non-human genetic studies, one reason being that a detected vQTL could represent a previously undetected interaction with other genes or environmental factors. The simultaneouspublication of these new methods in different journals has in many cases precluded opportunity for comparison. We survey some of these methods, the respective trade-offs they imply, and the connections between them. The methods fall into three main groups: classical non-parametric, fully parametric, and semi-parametric two-stage approximations. Choosing between alternatives involves balancing the need for robustness, flexibility, and speed. For each method, we identify important assumptions and limitations, including those of practical importance, such as their scope for including covariates and random effects. We show in simulations that both parametric methods and their semi-parametric approximations can give elevated false positive rates when they ignore mean-variancerelationships intrinsic to the data generation process. We conclude that choice of method depends on the trait distribution, the need to include non-genetic covariates, and the population size and structure, coupled with a critical evaluation of how these fit with the assumptions of the statistical model. [ABSTRACT FROM AUTHOR]

Dealing with genotype uncertainty is an ongoing issue in genetic analyses of complex traits. Here we consider genotype uncertainty in quantitative trait loci (QTL) analyses for large crosses in variance component models, where the genetic information is included in identity-by-descent (IBD) matrices. An IBD matrix is one realization from a distribution of potential IBD matrices given available marker information. In QTL analyses, its expectation is normally used resulting in potentially reduced accuracy and loss of power. Previously, IBD distributions have been included in models for small human full-sib families. We develop an Expectation–Maximization (EM) algorithm for estimating a full model based on Monte Carlo imputation for applications in large animal pedigrees. Our simulations show that the bias of variance component estimates using traditional expected IBD matrix can be adjusted by accounting for the distribution and that the calculations are computationally feasible for large pedigrees. [ABSTRACT FROM PUBLISHER]

Background: Genome-wide dense markers have been used to detect genes and estimate relative genetic values. Among many methods, Bayesian techniques have been widely used and shown to be powerful in genome-wide breeding value estimation and association studies. However, computation is known to be intensive under the Bayesian framework, and specifying a prior distribution for each parameter is always required for Bayesian computation. We propose the use of hierarchical likelihood to solve such problems. Results: Using double hierarchical generalized linear models, we analyzed the simulated dataset provided by the QTLMAS 2010 workshop. Marker-specific variances estimated by double hierarchical generalized linear models identified the QTL with large effects for both the quantitative and binary traits. The QTL positions were detected with very high accuracy. For young individuals without phenotypic records, the true and estimated breeding values had Pearson correlation of 0.60 for the quantitative trait and 0.72 for the binary trait, where the quantitative trait had a more complicated genetic architecture involving imprinting and epistatic QTL. Conclusions: Hierarchical likelihood enables estimation of marker-specific variances under the likelihoodist framework. Double hierarchical generalized linear models are powerful in localizing major QTL and computationally fast. [ABSTRACT FROM AUTHOR]

Background: Linkage mapping is used to identify genomic regions affecting the expression of complex traits. However, when experimental crosses such as F2 populations or backcrosses are used to map regions containing a Quantitative Trait Locus (QTL), the size of the regions identified remains quite large, i.e. 10 or more Mb. Thus, other experimental strategies are needed to refine the QTL locations. Advanced Intercross Lines (AIL) are produced by repeated intercrossing of F2 animals and successive generations, which decrease linkage disequilibrium in a controlled manner. Although this approach is seen as promising, both to replicate QTL analyses and fine-map QTL, only a few AIL datasets, all originating from inbred founders, have been reported in the literature. Methods: We have produced a nine-generation AIL pedigree (n = 1529) from two outbred chicken lines divergently selected for body weight at eight weeks of age. All animals were weighed at eight weeks of age and genotyped for SNP located in nine genomic regions where significant or suggestive QTL had previously been detected in the F2 population. In parallel, we have developed a novel strategy to analyse the data that uses both genotype and pedigree information of all AIL individuals to replicate the detection of and fine-map QTL affecting juvenile body weight. Results: Five of the nine QTL detected with the original F2 population were confirmed and fine-mapped with the AIL, while for the remaining four, only suggestive evidence of their existence was obtained. All original QTL were confirmed as a single locus, except for one, which split into two linked QTL. Conclusions: Our results indicate that many of the QTL, which are genome-wide significant or suggestive in the analyses of large intercross populations, are true effects that can be replicated and fine-mapped using AIL. Key factors for success are the use of large populations and powerful statistical tools. Moreover, we believe that the statistical methods we have developed to efficiently study outbred AIL populations will increase the number of organisms for which in-depth complex traits can be analyzed. [ABSTRACT FROM AUTHOR]

We present the hglm package for fitting hierarchical generalized linear models. It can be used for linear mixed models and generalized linear mixed models with random effects for a variety of links and a variety of distributions for both the outcomes and the random effects. Fixed effects can also be fitted in the dispersion part of the model. [ABSTRACT FROM AUTHOR]

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM. [ABSTRACT FROM AUTHOR]

Variance component models are commonly used to detect quantitative trait loci (QTL) in general pedigrees. The variance-covariance structure of the random QTL effect is given by the identity by descent (IBD) between genotypes. Epistatic effects have previously been modeled, both for unlinked and linked loci, as a random effect with a variance-covariance structure given by the Hadamard product between the IBD matrices of the direct QTL effects. In the original papers, the model was given but not derived. Here, we identify the underlying assumptions of this previously proposed model. It assumes that either an unlinked QTL or a fully informative marker (i.e., all marker alleles are unique in the base generation) is located between the loci. We discuss the need of developing a general algorithm to estimate the variance-covariance structure of the random epistatic effect for linked loci. [ABSTRACT FROM AUTHOR]

Robust and efficient optimization methods for variance component estimation using Restricted Maximum Likelihood (REML) models for genetic mapping of quantitative traits are considered. We show that the standard Newton-AI scheme may fail when the optimum is located at one of the constraint boundaries, and we introduce different approaches to remedy this by taking the constraints into account. We approximate the Hessian of the objective function using the average information matrix and also by using an inverse BFGS formula. The robustness and efficiency is evaluated for problems derived from two experimental data from the same animal populations. [ABSTRACT FROM AUTHOR]

The article offers information on quality assurance in education and research using a peer‐review system at the university. It mentions views on several research topics including humanities, social science, health, engineering and natural sciences; minor political scandal at the national level in Sweden; and importance of traditional peer reviewed journals.

Understanding the genetic basis of traits involved in adaptation is a major challenge in evolutionary biology but remains poorly understood. Here, we use genome-wide association mapping using a custom 50 k single nucleotide polymorphism (SNP) array in a natural population of collared flycatchers to examine the genetic basis of clutch size, an important life-history trait in many animal species. We found evidence for an association on chromosome 18 where one SNP significant at the genome-wide level explained 3.9% of the phenotypic variance. We also detected two suggestive quantitative trait loci (QTLs) on chromosomes 9 and 26. Fitness differences among genotypes were generally weak and not significant, although there was some indication of a sex-by-genotype interaction for lifetime reproductive success at the suggestive QTL on chromosome 26. This implies that sexual antagonism may play a role in maintaining genetic variation at this QTL. Our findings provide candidate regions for a classic avian life-history trait that will be useful for future studies examining the molecular and cellular function of, as well as evolutionary mechanisms operating at, these loci. [ABSTRACT FROM AUTHOR]