Heliodor Wierzbicki

A microsatellite study in the Łęgucki Młyn/Popielno hybrid zone reveals no genetic differentiation between two chromosome races of the common shrew (Sorex araneus)

This study investigated a chromosome hybrid zone between two chromosomal races of the common shrew (Sorex araneus). Gene flow and genetic structure of the hybrid zone, located in the northeast of Poland, were studied using seven polymorphic autosomal microsatellite loci (L9, L14, L33, L45, L67, L68, L97) and a Y-linked microsatellite locus (L8Y). Seventy-five animals (46 of the Łęgucki Młyn race and 29 of the Popielno race) from nine different localities were examined and the data were analyzed using hierarchical AMOVA and F-statistic. The studied microsatellite loci and races (divided into nine geographical populations) were characterized by observed heterozygosity (HO), expected heterozygosities within (HS), and between (HT) populations, inbreeding coefficient (FIS), fixation index (FST), and average allelic richness (A). We found that genetic structuring within and between the two chromosome races were weak and non-significant. This finding and unconstrained gene flow between the races indicates a high level of migration within the Łęgucki Młyn/Popielno hybrid zone, suggesting that evolutionarily important genetic structuring does not occur in interracial zones where races which are not genetically distinct come into contact.

Evolutionary history and phylogeographic relationships of shrews from Sorex araneus group

Shrews of the Sorex genus are an evolutionarily successful group that includes more than 77 species widely distributed in Eurasia and North America. The genus is one of the rare cases where karyotypic changes reflect well the evolutionary relationships among its species. The taxa showing the greatest variation in karyotype are usually classified into the Sorex araneus group. Its evolution was associated with chromosomal rearrangements, which could have promoted fast diversification of this group into many chromosomal races and species. These processes were additionally complicated by introgressions of mitochondrial DNA, which made the evolutionary history of this group quite complex and difficult to infer. To tackle the problem, we performed multi-method phylogenetic analyses based on mitochondrial cytochrome b that is considered a good molecular marker available for many representatives of Sorex. The results were compared with phylogenies based on chromosomal rearrangement data and put into temporal and spatial context using molecular dating and historical biogeography methods. We complemented the study with the estimation of diversification rates within the S. araneus group as well as comparing the results with paleontological records and climatic oscillations within the last 4 million years. Based on the gathered data, we proposed a hypothetical scenario for the evolution and geographic dispersion of species belonging to the S. araneus group. The shrews began to diversify about 2.7 million years ago in Eurasia and then migrated at least twice to North America. The evolution of shrews was driven by Pleistocene glacial and interglacial cycles, which increased their speciation rate and the emergence of new lineages. The migrations of populations were accompanied by introgressions of mitochondrial DNA into native shrews and occurred at least twice.

Microsatellite polymorphism and its association with body weight and selected morphometrics of farm red fox (Vulpes vulpes L.)

Polymorphism of 30 canine-derived microsatellites was studied in a group of 200 red foxes kept on 2 Polish farms. 22 out of 30 microsatellites were selected to study association between marker genotypes and body weight (BW), body length (BL), body circumference (BC), tail length (TL), ear height (EH), length of the right front limb (FRLL), length of the right rear limb (RRLL), length of the right front foot (FRFL) and length of the right rear foot (RRFL). A total of 112 alleles and 243 genotypes were found at 22 autosomal microsatellite loci. Three monomorphic loci deemed as uninformative were excluded from the study. The association between marker genotypes and the studied traits was analysed using general linear model (GLM) procedure and least squares means (LSM). Linkage disequilibrium (LD) was estimated to assess non-random association between microsatellite loci. Out of 19 microsatellites studied four markers showed no association with the studied traits, three markers had a significant effect on one trait, and another three markers had significant effect on two traits. Among ten microsatellites with significant effect on four economically important traits (BW, BL, BC, TL) four were associated with two characters: marker FH2613 with BW and BC, marker FH2097withBL and BC, marker ZUBECA6 with BW and BC, whereas marker REN75M10 was associated with BL and TL. The strongest LD (r2 ranged from 0.15 to 0.33) was estimated between nine loci with significant effect on economically important traits (BW, BL, BC, TL).

Statistical aspects of genetic association testing in small samples, based on selective DNA pooling data in the arctic fox

We analysed data from a selective DNA pooling experiment with 130 individuals of the arctic fox (Alopex lagopus), which originated from 2 different types regarding body size. The association between alleles of 6 selected unlinked molecular markers and body size was tested by using univariate and multinomial logistic regression models, applying odds ratio and test statistics from the power divergence family. Due to the small sample size and the resulting sparseness of the data table, in hypothesis testing we could not rely on the asymptotic distributions of the tests. Instead, we tried to account for data sparseness by (i) modifying confidence intervals of odds ratio; (ii) using a normal approximation of the asymptotic distribution of the power divergence tests with different approaches for calculating moments of the statistics; and (iii) assessingP values empirically, based on bootstrap samples. As a result, a significant association was observed for 3 markers. Furthermore, we used simulations to assess the validity of the normal approximation of the asymptotic distribution of the test statistics under the conditions of small and sparse samples.

Linear models for breeding values prediction in haplotype-assisted selection - an analysis of QTL-MAS Workshop 2011 Data

BackgroundThe aim of this study was to estimate haplotype effects and then to predict breeding values using linear models. The haplotype based analysis enables avoidance of loosing information due to linkage disequilibrium between single markers. There are also less explanatory variables in the linear model which makes the estimation more reliable.MethodsDifferent methods and criteria for marker and haplotype selection were considered. First, markers with MAF lower than 5% where excluded from the data set. Then, SNPs in complete linkage disequilibrium where selected. Next step was to construct haplotypes and to estimate their frequencies basing on selected SNPs. The haplotypes with a frequency lower than 1% were not considered in further analysis. Chosen haplotypes were used as the explanatory variables in the linear models for breeding values prediction. Linear models with fixed and random haplotype effects as well as animal model were tested.ResultsThe number of markers was limited to 1206, 1189, 1249, 1288 and 1167 for chromosome 1, 2, 3, 4 and 5, respectively due to MAF criterion. In total 409 subsets of SNPs with r2=1 were found. 1476 haplotypes with different lengths were inferred. The frequencies of 817 haplotypes were higher than 1% - 184 for the first chromosome, 172 for the second, 131 for the third, 146 for the forth and 184 haplotypes for the fifth chromosome. The haplotype effects estimated using random models were comparable and more precise in prediction for individuals with unknown phenotypes. A few haplotypes with large effects were found when their effects were defined as fixed in the linear model . The correlations of the predicted breeding values with true breeding values were not that high. This could be brought about by selection criteria imposed on the genotype data which led to substantial reduction of number of markers.ConclusionsAlthough not many markers were considered in the study, the results obtained show that the implemented approach can be considered as quite promising. The haplotype approach let to avoid high dimensional models as compared with single SNPs models.

Comparative analysis of morphometrics of wild and farm foxes (Vulpes vulpes L.) – preliminary results

Nowadays, an interesting aspect of animal farming is the possibility to analyze genetic, morphometric and behavioural changes in farm animals, which have taken place as a result of selective breeding. Then, the comparative analysis with wild animals can be done. Fur animals (fox, mink, raccoon dog) have relatively short history of domestication, so an interesting point is if this process affected body measurements of those species. The aim of the present study was to analyze and compare several body measurements of wild and farm foxes. The analyzed traits were: body weight, body length and circumference, tail length, ears height, length of front and hind limbs, length of front and hind feet. Farm animals were characterized by higher values of all traits (except tail length and ear height), than wild ones and these differences were statistically significant. Smaller differences between wild and captive populations were found for traits which were not crucial for species survival (e.g. tail length and ears height). Comparison between sexes showed that they differed in many size measurements, but measurements of females from the farm population and males from the wild population were comparable.