Genetic diversity and genomic signatures of selection among cattle breeds from Siberia, eastern and northern Europe.

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P-value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.

Genome-wide analysis reveals the ancient and recent admixture history of East African Shorthorn Zebu from Western Kenya.

The Kenyan East African zebu cattle are valuable and widely used genetic resources. Previous studies using microsatellite loci revealed the complex history of these populations with the presence of taurine and zebu genetic backgrounds. Here, we estimate at genome-wide level the genetic composition and population structure of the East African Shorthorn Zebu (EASZ) of western Kenya. A total of 548 EASZ from 20 sub-locations were genotyped using the Illumina BovineSNP50 v. 1 beadchip. STRUCTURE analysis reveals admixture with Asian zebu, African and European taurine cattle. The EASZ were separated into three categories: substantial (⩾12.5%), moderate (1.56%

Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium.

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)-pairs reflects distant breeding history while LD from distal SNP-pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD-based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD-based actual effective population size was estimated in the range of 87-149, whereas pedigree-based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.

Influence of edentulousness on gonial angle, ramus and condylar height.

It has been shown that the shape of the mandible correlates with occlusal condition and the function of the masticatory muscles. Edentulous subjects have a wider gonial angle than dentate subjects, and a gender difference has also been shown. However, some studies have reported differing results. Less is known about the effect of dental status and gender on the ramus and condylar height. The aim of this study was to evaluate the association of tooth loss on the shape of mandible (i.e., gonial angle, ramus height and condylar height) in subjects aged 60 years and older. A total of 1036 subjects (667 dentate, 389 edentulous; 554 women and 482 men) were included in the study. Interviews and clinical and panoramic radiographic examinations were carried out. The gonial angle of the mandible and the mandibular and condylar height were measured using panoramic radiographs. In edentulous subjects, the gonial angle was significantly larger, while the ramus and condylar height was significantly smaller on both sides compared with dentate subjects. Women had a significantly larger gonial angle and smaller ramus and condylar height on both sides compared with men. In conclusion, the morphology of the mandible changes as a consequence of tooth loss, which can be expressed as a widening of the gonial angle and shortening of the ramus and condylar height. These findings highlight the importance of rehabilitation of the masticatory system to maintain good functioning of the masticatory muscles for as long as possible.

Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds.

The extent of linkage disequilibrium (LD) and effective population size in Finnish Landrace and Finnish Yorkshire pig populations were studied using a whole genome SNP panel (Illumina PorcineSNP60 BeadChip) and pedigree data. Genotypic data included 86 Finnish Landrace and 32 Finnish Yorkshire boars. Pedigree data included 608,138 Finnish Landrace 554,237 and Finnish Yorkshire pigs, and on average 15 ancestral generations were known for the reference animals, born in 2005 to 2009. The breeding animals of the 2 populations have been kept separate in the breeding programs. Based on the pedigree data, the current effective population size for Finnish Landrace is 91 and for Finnish Yorkshire 61. Linkage disequilibrium measures (D' and r(2)) were estimated for over 1.5 million pairs of SNP. Average r(2) for SNP 30 kb apart was 0.47 and 0.49 and for SNP 5 Mb apart 0.09 and 0.12 for Finnish Landrace and Finnish Yorkshire, respectively. Average LD (r(2)) between adjacent SNP in the Illumina PorcineSNP60 BeadChip was 0.43 (57% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Landrace and 0.46 (60% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Yorkshire, and average r(2) > 0.2 extended to 1.0 and 1.5 Mb for Finnish Landrace and Finnish Yorkshire, respectively. Effective population size estimates based on the decay of r(2) with distance were similar to those based on the pedigree data: 80 and 55 for Finnish Landrace and Finnish Yorkshire, respectively. Thus, the results indicate that the effective population size of Finnish Yorkshire is smaller than of Finnish Landrace and has a clear effect on the extent of LD. The current effective population size of both breeds is above the recommended minimum of 50 but may get smaller than that in the near future, if no action is taken to balance the inbreeding rate and selection response. Because a moderate level of LD extends over a long distance, selection based on whole genome SNP markers (genomic selection) is expected to be efficient for both breeds.

Mitochondrial sequence reveals high levels of gene flow between breeds of domestic sheep from Asia and Europe.

Sequence variation present within the mitochondrial genome was used to investigate genetic diversity within sheep breeds from Asia and Europe. Comparison of 2027 bp of sequence from 121 animals revealed 44 phylogenetically informative nucleotide positions and a single insertion/deletion. A total of 57 haplotypes were observed which formed two distinct clades. Type A haplotypes were found in breeds from Asia (India, Indonesia, Mongolia, and Tibet), while type B haplotypes were observed at the highest frequency in breeds sourced from Europe (nine breeds from Austria, Aland, Finland, Spain, and northwestern Russia). The distribution of haplotypes indicates sheep appear to have the weakest population structure and the highest rate of intercontinental dispersal of any domestic animal reported to date. Only 2.7% of the sequence variation observed was partitioned between continents, which is lower than both goat (approximately 10%) and cattle (approximately 50%). Diagnostic restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) tests which distinguish type A and B haplotypes were used to test an additional 223 animals from 17 breeds of European and Asian origin. A mixture of the two lineages was found in every breed except Suffolk and the Indian Garole, indicating introgression has played a major part during breed development and subsequent selection.

Native breeds demonstrate high contributions to the molecular variation in northern European sheep.

Population contribution to genetic diversity can be estimated using neutral variation. However, population expansion or hybridization of diverged ancestries may weaken correlation between neutral and non-neutral variation. Microsatellite variation was studied at 25 loci in 20 native and 12 modern or imported northern European sheep breeds. Breed contributions to total gene diversity, allelic richness and mean allele-sharing distance between individuals were measured. Indications of changes in population size and admixtures of divergent ancestries were investigated and the extent of inbreeding was estimated. The northern European sheep demonstrated signs of reduction in effective population size. Many old, small populations made a substantial positive contribution to total molecular variation, but populations with several divergent major ancestries did not contribute substantially to molecular variation, with the exception of the Norwegian Rygja sheep. However, several diverged major ancestries may cause it to contribute less to non-neutral variation than expected from the microsatellite data. Breed uniqueness and within-breed variability generally had opposite effects on breed contributions to molecular diversity. The degree of inbreeding did not reflect the breed contribution to total gene diversity or allelic richness, but inbred populations increased the mean allele-sharing distance between individuals. Our study indicates breed conservation to be especially important in maintaining allelic variation in northern European sheep and supports the evolutionary importance of peripheral populations.

Unfolding of population structure in Baltic sheep breeds using microsatellite analysis.

Studies of domestic animals are performed on breeds, but a breed does not necessarily equate to a genetically defined population. The division of sheep from three native and four modern Baltic sheep breeds was studied using 21 microsatellite loci and applying a Bayesian clustering method. A traditional breed-wise approach was compared to that relying on the pattern of molecular diversity. In this study, a breed was found to be inconsistent with a distinct genetic population for three reasons: (i) a lack of differentiation between modern Baltic breeds, since the majority of the studied sheep formed a single population; (ii) the presence of individuals of foreign ancestry within the breed; and (iii) an undefined local Saaremaa sheep was referred to as a breed, but was shown to consist of separate populations. In the breed-wise approach, only the clearly distinct Ruhnu sheep demonstrated low within-breed variation, although the newly identified Saaremaa populations also have low variability. Providing adequate management recommendations for the Saaremaa sheep is not possible without further studies, but the potential harmful effects of inbreeding in the Ruhnu sheep could be reduced through the use of two genetically related Saaremaa populations. In other breeds, excessive crossing appears to be a larger concern than inbreeding. Assigning individuals into populations based on the pattern of genetic diversity offers potentially unbiased means of elucidating the genetic population structure of species. Combining these genetic populations with phenotypic and aetiological data will enable formulation of the most informed recommendations for gene resource management.

Comparison of microsatellite and blood protein diversity in sheep: inconsistencies in fragmented breeds.

Finnsheep, Romanov, Oxford Down and three local breeds from Finland or northwestern Russia were assessed at 15 microsatellite and 7 protein loci. A novel albumin allele was identified. Diversity patterns were mostly concordant between marker types, but discrepancies appeared for the local Viena and Vepsia sheep, both demonstrating frequent linkage disequilibria for both marker types and excess of homozygotes for microsatellites, and in the case of Vepsia also for proteins as signs of breed fragmentation. On the basis of microsatellite data, the neighbour-joining tree and two-dimensional map constructed from DA distances suggested that difference in longitude of breed origin would relate to breed relationship, whereas on the basis of protein data latitude would have this quality. These different impressions resulted because genetic distances involving Vepsia sheep were relatively low for protein variation compared with microsatellites. Microsatellite variation correlated positively with protein variation, but for the local Viena sheep protein variation was comparatively low. Populations had significant differences in allelic richness, but not in genetic diversity. Analysis implied that at least 30 polymorphic loci were needed to detect a difference in diversity between populations using a paired t-test, if the true mean diversity difference was 0.2. In the total sample, proteins demonstrated larger theta-values, but this was reversed for Finnsheep, for which model-based clustering of microsatellite genotypes revealed a structure associated with coat colour. Imported and rare sheep exhibited lowered allelic variability and increased frequency of pairwise disequilibria between unlinked markers. Our results emphasize that more loci are required for studying fragmented breeds.