Genome-wide association study for body weight in cattle populations from Siberia.

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome-wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed- and species-specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.

Some pitfalls in application of functional data analysis approach to association studies.

One of the most effective methods for gene-based mapping employs functional data analysis, which smoothes data using standard basis functions. The full functional linear model includes a functional representation of genotypes and their effects, while the beta-smooth only model smoothes the genotype effects only. Benefits and limitations of the beta-smooth only model should be studied before using it in practice. Here we analytically compare the full and beta-smooth only models under various scenarios. We show that when the full model employs two sets of basis functions equal in type and number, genotypes smoothing is eliminated from the model and it becomes analytically equivalent to the beta-smooth only model. If the basis functions differ only in type, genotypes smoothing is also eliminated from the full model, but the type of basis functions used for smoothing genotype effects becomes redefined. This leads to misinterpretation of the results and may reduce statistical power. When basis functions differ in number, no analytical comparison of the full and beta-smooth only models is possible. However, we show that the numbers of basis functions set unequal can become equal during the analysis, and the full model becomes disadvantageous.

[Functional linear models for region-based association analysis].

Regional association analysis is one of the most powerful tools for gene mapping because instead analysis of individual variants it simultaneously considers all variants in the region. Recent development of the models for regional association analysis involves functional data analysis approach. In the framework of this approach, genotypes of variants within region as well as their effects are described by continuous functions. Such approach allows us to use information about both linkage and linkage disequilibrium and reduce the influence of noise and/or observation errors. Here we define a functional linear mixed model to test association on independent and structured samples. We demonstrate how to test fixed and random effects of a set of genetic variants in the region on quantitative trait. Estimation of statistical properties of new methods shows that type I errors are in accordance with declared values and power is high especially for models with fixed effects of genotypes. We suppose that new functional regression linear models facilitate identification of rare genetic variants controlling complex human and animal traits. New methods are implemented in computer software FREGAT which is available for free download at http://mga.bionet.nsc.ru/soft/FREGAT/.

Linkage analysis of adult height in a large pedigree from a Dutch genetically isolated population.

Despite extensive research of genetic determinants of human adult height, the genes identified up until now allow to predict only a small proportion of the trait's variance. To identify new genes we analyzed 2,486 genotyped and phenotyped individuals in a large pedigree including 23,612 members in 18 generations. The pedigree was derived from a young genetically isolated Dutch population, where genetic heterogeneity is expected to be low and linkage disequilibrium has been shown to be increased. Complex segregation analysis confirmed high heritability of adult height, and suggested mixed model of height inheritance in this population. The estimates of the model parameters obtained from complex segregation analysis were used in parametric linkage analysis, which highlighted three genome-wide significant and additionally at least four suggestive loci involved in height. Significant peaks were located at the chromosomal regions 1p32 (LOD score = 3.35), 2p16 (LOD score = 3.29) and 16q24 (LOD score = 3.94). For the latter region, a strong association signal (FDR q < 0.05) was obtained for 19 SNPs, 17 of them were located in the CDH13 (cadherin 13) gene of which one (rs1035569) explained 1.5% of the total height variance.

Identification of all pachytene bivalents in the common shrew using DAPI-staining of synaptonemal complex spreads.

A major problem in studies of synaptonemal complexes (SC) is the difficulty in distinguishing individual chromosomes. This problem can be solved combining SC immunostaining with FISH of chromosome-specific sequences. However, this procedure is expensive, time-consuming and applicable only to a very limited number of species. In this paper we show how a combination of SC immunostaining and DAPI staining can allow identification of all chromosome arms in surface-spreads of the SC of the common shrew (Sorex araneus L.). Enhancement of brightness and contrast of the images with photo editing software allowed us to reveal clear DAPI-positive and negative bands with relative sizes and positions similar to DAPI landmarks on mitotic metaphase chromosomes. Using FISH with DNA probes prepared from chromosome arms m and n we demonstrated correct recognition of the chromosomes mp and hn on the basis of their DAPI pattern. We show that the approach we describe here may be applied to other species and can provide an important tool for identification of individual bivalents in pachytene surface-spreads.