Genetic basis of ear length in sheep breeds sampled across the region from the Middle East to the Alps.

Ear length in sheep (Ovis aries) shows a wide range of natural variation, from the absence of an outer ear structure (anotia), to small outer ears (microtia), to regular ear length. Up until now, the underlying genetics of this phenotype has been studied in four sheep breeds from China, Jordan and Italy. These studies revealed a broad range of genes significantly associated with ear length, potentially indicating genetic heterogeneity across breeds or geographic regions. In the current study, we performed genome-wide SNP genotyping and haplotype-based mapping, in a population of 340 individuals, to identify loci influencing ear length variation in additional sheep breeds from Slovenia, Croatia, Cyprus and Greece. Additionally, two previously described candidate variants were also genotyped in our mapping population. The mapping model without candidate variant genotypes revealed only one genome-wide significant signal, which was located next to HMX1 on OAR6. This region was previously described as being associated with ear length variation in the Altay and Awassi sheep breeds. The mapping model including the candidate duplication genotype near HMX1 as a fixed effect explained the phenotypic variance on OAR6 and revealed an additional genome-wide significant locus on OAR13 associated with ear length. Our results, combined with published evidence, suggest that a duplication in the evolutionarily conserved region near HMX1 is the major regulator of ear length in sheep breeds descended from a larger region from Central Asia, to the Middle East, Cyprus, Greece and to the Alps. This distribution suggests an ancient origin of the derived allele.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region.

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

Fine-mapping and identification of candidate causal genes for tail length in the Merinolandschaf breed.

Docking the tails of lambs in long-tailed sheep breeds is a common practice worldwide. But this practice is associated with pain. Breeding for a shorter tail could offer an alternative. Therefore, this study aimed to analyze the natural tail length variation in the Merinolandschaf and to identify causal alleles for the short tail phenotype segregating within long-tailed breeds. We used SNP-based association analysis and haplotype-based mapping in 362 genotyped (Illumina OvineSNP50) and phenotyped Merinolandschaf lambs. Genome-wide significant regions were capture sequenced in 48 lambs and comparatively analyzed in various long and short-tailed sheep breeds and wild sheep subspecies. Here we show a SNP located in the first exon of HOXB13 and a SINE element located in the promotor of HOXB13 as promising candidates. These results enable more precise breeding towards shorter tails, improve animal welfare by amplification of ancestral alleles and contribute to a better understanding of differential embryonic development.

Mitochondrial DNA Variation Contributes to the Aptitude for Dressage and Show Jumping Ability in the Holstein Horse Breed.

Maternal lineages are considered an important factor in breeding. Mitochondrial DNA (mtDNA) is maternally inherited and plays an important role in energy metabolism. It has already been associated with energy consumption and performances, e.g., stamina in humans and racehorses. For now, corresponding studies are lacking for sport performance of warmblood breeds. MtDNA sequences were available for 271 Holstein mares from 75 maternal lineages. As all mares within a lineage showed identical haplotypes regarding the non-synonymous variants, we expanded our data set by also including non-sequenced mares and assigning them to the lineage-specific haplotype. This sample consisting of 6334 to 16,447 mares was used to perform mitochondrial association analyses using breeding values (EBVs) estimated on behalf of the Fédération Équestre Nationale (FN) and on behalf of the Holstein Breeding Association (HOL). The association analyses revealed 20 mitochondrial SNPs (mtSNPs) significantly associated with FN-EBVs and partly overlapping 20 mtSNPs associated with HOL-EBVs. The results indicated that mtDNA contributes to performance differences between maternal lineages. Certain mitochondrial haplogroups were associated with special talents for dressage or show jumping. The findings encourage to set up innovative genetic evaluation models that also consider information on maternal lineages.

Exploring the Origin and Relatedness of Maternal Lineages Through Analysis of Mitochondrial DNA in the Holstein Horse.

Maternal lineages are important for the breeding decision in the Holstein horse breed. To investigate the genetic diversity of the maternal lineages and the relationships between founder mares, the maternal inherited mitochondrial genome (except the repetitive part of the non-coding region) of 271 mares representing 75 lineages was sequenced. The sequencing predominantly revealed complete homology in the nucleotide sequences between mares from one lineage with exceptions in 13 lineages, where differences in one to three positions are probably caused by de novo mutations or alternate fixation of heteroplasmy. We found 78 distinct haplotypes that have not yet been described in other breeds. Six of these occurred in two or three different lineages indicating a common ancestry. Haplotypes can be divided into eight clusters with all mares from one lineage belonging to the same cluster. Within a cluster, the average number of pairwise differences ranged from zero to 16.49 suggesting close maternal relationships between these mares. The results showed that the current breeding population originated from at least eight ancestral founder mares.

A nonsense mutation of bone morphogenetic protein-15 (BMP15) causes both infertility and increased litter size in pigs.

BACKGROUND: Atypical external genitalia are often a sign of reproductive organ pathologies and infertility with both environmental or genetic causes, including karyotypic abnormalities. Genome-wide association studies (GWAS) provide a means for identifying chromosomal regions harboring deleterious DNA-variants causing such phenotypes. We performed a GWAS to unravel the causes of incidental cases of atypically small vulvae in German Landrace gilts. RESULTS: A case-control GWAS involving Illumina porcine SNP60 BeadChip-called genotypes of 17 gilts with atypically small vulvae and 1818 control animals (fertile German Landrace sows) identified a significantly associated region on the X-chromosome (P = 8.81 × 10- 43). Inspection of whole-genome sequencing data in the critical area allowed us to pinpoint a likely causal variant in the form of a nonsense mutation of bone morphogenetic protein-15 (BMP15; Sscrofa11.1_X:g.44618787C>T, BMP15:p.R212X). The mutant allele occurs at a frequency of 6.2% in the German Landrace breeding population. Homozygous gilts exhibit underdeveloped, most likely not functional ovaries and are not fertile. Male carriers do not seem to manifest defects. Heterozygous sows produce 0.41±0.02 (P=4.5 × 10-83) piglets more than wildtype animals. However, the mutant allele's positive effect on litter size accompanies a negative impact on lean meat growth. CONCLUSION: Our results provide an example for the power of GWAS in identifying the genetic causes of a fuzzy phenotype and add to the list of natural deleterious BMP15 mutations that affect fertility in a dosage-dependent manner, the first time in a poly-ovulatory species. We advise eradicating the mutant allele from the German Landrace breeding population since the adverse effects on the lean meat growth outweigh the larger litter size in heterozygous sows.

Genomic diversity and population structure of the indigenous Greek and Cypriot cattle populations.

BACKGROUND: The indigenous cattle populations from Greece and Cyprus have decreased to small numbers and are currently at risk of extinction due to socio-economic reasons, geographic isolation and crossbreeding with commercial breeds. This study represents the first comprehensive genome-wide analysis of 10 indigenous cattle populations from continental Greece and the Greek islands, and one from Cyprus, and compares them with 104 international breeds using more than 46,000 single nucleotide polymorphisms (SNPs). RESULTS: We estimated several parameters of genetic diversity (e.g. heterozygosity and allelic diversity) that indicated a severe loss of genetic diversity for the island populations compared to the mainland populations, which is mainly due to the declining size of their population in recent years and subsequent inbreeding. This high inbreeding status also resulted in higher genetic differentiation within the Greek and Cyprus cattle group compared to the remaining geographical breed groups. Supervised and unsupervised cluster analyses revealed that the phylogenetic patterns in the indigenous Greek breeds were consistent with their geographical origin and historical information regarding crosses with breeds of Anatolian or Balkan origin. Cyprus cattle showed a relatively high indicine ancestry. Greek island populations are placed close to the root of the tree as defined by Gir and the outgroup Yak, whereas the mainland breeds share a common historical origin with Busa. Unsupervised clustering and D-statistics analyses provided strong support for Bos indicus introgression in almost all the investigated local cattle breeds along the route from Anatolia up to the southern foothills of the Alps, as well as in most cattle breeds along the Apennine peninsula to the southern foothills of the Alps. CONCLUSIONS: All investigated Cyprus and Greek breeds present complex mosaic genomes as a result of historical and recent admixture events between neighbor and well-separated breeds. While the contribution of some mainland breeds to the genetic diversity pool seems important, some island and fragmented mainland breeds suffer from a severe decline of population size and loss of alleles due to genetic drift. Conservation programs that are a compromise between what is feasible and what is desirable should focus not only on the still highly diverse mainland breeds but also promote and explore the conservation possibilities for island breeds.

Are scurs in heterozygous polled (Pp) cattle a complex quantitative trait?

BACKGROUND: Breeding genetically hornless, i.e. polled, cattle provides an animal welfare-friendly and non-invasive alternative to the dehorning of calves. However, the molecular regulation of the development of horns in cattle is still poorly understood. Studying genetic characters such as polledness and scurs, can provide valuable insights into this process. Scurs are hornlike formations that occur occasionally in a wide variety of sizes and forms as an unexpected phenotype when breeding polled cattle. METHODS: We present a unique dataset of 885 Holstein-Friesian cattle with polled parentage. The horn phenotype was carefully examined, and the phenotypic heterogeneity of the trait is described. Using a direct gene test for polledness, the polled genotype of the animals was determined. Subsequently, the existence of a putative scurs locus was investigated using high-density genotype data of a selected subset of 232 animals and two mapping approaches: mixed linear model-based association analyses and combined linkage disequilibrium and linkage analysis. RESULTS: The results of an exploratory data analysis indicated that the expression of scurs depends on age at phenotyping, sex and polled genotype. Scurs were more prevalent in males than in females. Moreover, homozygous polled animals did not express any pronounced scurs and we found that the Friesian polled allele suppresses the development of scurs more efficiently than the Celtic polled allele. Combined linkage and linkage disequilibrium mapping revealed four genome-wide significant loci that affect the development of scurs, one on BTA5 and three on BTA12. Moreover, suggestive associations were detected on BTA16, 18 and 23. The mixed linear model-based association analysis supports the results of the combined linkage and linkage disequilibrium analysis. None of the mapping approaches provided convincing evidence for a monogenic inheritance of scurs. CONCLUSIONS: Our results contradict the initial and still broadly accepted model for the inheritance of horns and scurs. We hypothesise an oligogenetic model to explain the development of scurs and polledness.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals.

Stature is affected by many polymorphisms of small effect in humans 1 . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

Short communication: Genomic prediction using imputed whole-genome sequence variants in Brown Swiss Cattle.

The accuracy of genomic prediction determines response to selection. It has been hypothesized that accuracy of genomic breeding values can be increased by a higher density of variants. We used imputed whole-genome sequence data and various single nucleotide polymorphism (SNP) selection criteria to estimate genomic breeding values in Brown Swiss cattle. The extreme scenarios were 50K SNP chip data and whole-genome sequence data with intermediate scenarios using linkage disequilibrium-pruned whole-genome sequence variants, only variants predicted to be missense, or the top 50K variants from genome-wide association studies. We estimated genomic breeding values for 3 traits (somatic cell score, nonreturn rate in heifers, and stature) and found differences in accuracy levels between traits. However, among different SNP sets, accuracy was very similar. In our analyses, sequence data led to a marginal increase in accuracy for 1 trait and was lower than 50K for the other traits. We concluded that the inclusion of imputed whole-genome sequence data does not lead to increased accuracy of genomic prediction with the methods.

Genome-wide association studies of fertility and calving traits in Brown Swiss cattle using imputed whole-genome sequences.

BACKGROUND: The detection of quantitative trait loci has accelerated with recent developments in genomics. The introduction of genomic selection in combination with sequencing efforts has made a large amount of genotypic data available. Functional traits such as fertility and calving traits have been included in routine genomic estimation of breeding values making large quantities of phenotypic data available for these traits. This data was used to investigate the genetics underlying fertility and calving traits and to identify potentially causative genomic regions and variants. We performed genome-wide association studies for 13 functional traits related to female fertility as well as for direct and maternal calving ease based on imputed whole-genome sequences. Deregressed breeding values from ~1000-5000 bulls per trait were used to test for associations with approximately 10 million imputed sequence SNPs. RESULTS: We identified a QTL on BTA17 associated with non-return rate at 56 days and with interval from first to last insemination. We found two significantly associated non-synonymous SNPs within this QTL region. Two more QTL for fertility traits were identified on BTA25 and 29. A single QTL was identified for maternal calving traits on BTA13 whereas three QTL on BTA19, 21 and 25 were identified for direct calving traits. The QTL on BTA19 co-localizes with the reported BH2 haplotype. The QTL on BTA25 is concordant for fertility and calving traits and co-localizes with a QTL previously reported to influence stature and related traits in Brown Swiss dairy cattle. CONCLUSION: The detection of QTL and their causative variants remains challenging. Combining comprehensive phenotypic data with imputed whole genome sequences seems promising. We present a QTL on BTA17 for female fertility in dairy cattle with two significantly associated non-synonymous SNPs, along with five additional QTL for fertility traits and calving traits. For all of these we fine mapped the regions and suggest candidate genes and candidate variants.

Whole-genome analysis of introgressive hybridization and characterization of the bovine legacy of Mongolian yaks.

The yak is remarkable for its adaptation to high altitude and occupies a central place in the economies of the mountainous regions of Asia. At lower elevations, it is common to hybridize yaks with cattle to combine the yak's hardiness with the productivity of cattle. Hybrid males are sterile, however, preventing the establishment of stable hybrid populations, but not a limited introgression after backcrossing several generations of female hybrids to male yaks. Here we inferred bovine haplotypes in the genomes of 76 Mongolian yaks using high-density SNP genotyping and whole-genome sequencing. These yaks inherited ∼1.3% of their genome from bovine ancestors after nearly continuous admixture over at least the last 1,500 years. The introgressed regions are enriched in genes involved in nervous system development and function, and particularly in glutamate metabolism and neurotransmission. We also identified a novel mutation associated with a polled (hornless) phenotype originating from Mongolian Turano cattle. Our results suggest that introgressive hybridization contributed to the improvement of yak management and breeding.

The 1.78-kb insertion in the 3'-untranslated region of RXFP2 does not segregate with horn status in sheep breeds with variable horn status.

BACKGROUND: The mode of inheritance of horn status in sheep is far more complex than a superficial analysis might suggest. Observations, which were mostly based on crossbreeding experiments, indicated that the allele that results in horns is dominant in males and recessive in females, and some authors even speculated about the involvement of more than two alleles. However, all recent genome-wide association analyses point towards a very strong effect of a single autosomal locus on ovine chromosome 10, which was narrowed down to a putatively causal insertion polymorphism in the 3'-untranslated region of the relaxin/insulin-like family peptide receptor 2 gene (RXFP2). The main objective of this study was to test this insertion polymorphism as the causal mutation in diverse sheep breeds, including breeds with a variable and/or sex-dependent horn status. RESULTS: After re-sequencing a region of about 246 kb that covered the RFXP2 gene and its flanking regions for 24 sheep from six completely horned and six completely polled breeds, we identified the same insertion polymorphism that was previously published as segregating with horn status in these breeds. Multiplex PCR genotyping of 489 sheep from 34 breeds and some crosses between sheep breeds showed a nearly perfect segregation of the insertion polymorphism with horn status in sheep breeds of Central and Western European origin. In these breeds and their crossings, heterozygous males were horned and heterozygous females were polled. However, this segregation pattern was not, or at least not completely, reproducible in breeds with sex-dependent and/or variable horn status, especially in sheep that originated from even more southern European regions and from Africa. In such breeds, we observed almost all possible combinations of genotype, sex and horn status phenotype. CONCLUSIONS: The 1.78-kb insertion polymorphism in the 3'-untranslated region of RXFP2 and SNPs in the 3'-UTR, exon 14 and intron 11 of this gene that we analyzed in this study cannot be considered as the only cause of polledness in sheep and are not useful as a universal marker to define the genetic horn status in sheep.

Confirmation of a non-synonymous SNP in PNPLA8 as a candidate causal mutation for Weaver syndrome in Brown Swiss cattle.

BACKGROUND: Bovine progressive degenerative myeloencephalopathy (Weaver syndrome) is a neurodegenerative disorder in Brown Swiss cattle that is characterized by progressive hind leg weakness and ataxia, while sensorium and spinal reflexes remain unaffected. Although the causal mutation has not been identified yet, an indirect genetic test based on six microsatellite markers and consequent exclusion of Weaver carriers from breeding have led to the complete absence of new cases for over two decades. Evaluation of disease status by imputation of 41 diagnostic single nucleotide polymorphisms (SNPs) and a common haplotype published in 2013 identified several suspected carriers in the current breeding population, which suggests a higher frequency of the Weaver allele than anticipated. In order to prevent the reemergence of the disease, this study aimed at mapping the gene that underlies Weaver syndrome and thus at providing the basis for direct genetic testing and monitoring of today's Braunvieh/Brown Swiss herds. RESULTS: Combined linkage/linkage disequilibrium mapping on Bos taurus chromosome (BTA) 4 based on Illumina Bovine SNP50 genotypes of 43 Weaver-affected, 31 Weaver carrier and 86 Weaver-free animals resulted in a maximum likelihood ratio test statistic value at position 49,812,384 bp. The confidence interval (0.853 Mb) determined by the 2-LOD drop-off method was contained within a 1.72-Mb segment of extended homozygosity. Exploitation of whole-genome sequence data from two official Weaver carriers and 1145 other bulls that were sequenced in Run4 of the 1000 bull genomes project showed that only a non-synonymous SNP (rs800397662) within the PNPLA8 gene at position 49,878,773 bp was concordant with the Weaver carrier status. Targeted SNP genotyping confirmed this SNP as a candidate causal mutation for Weaver syndrome. Genotyping for the candidate causal mutation in a random sample of 2334 current Braunvieh animals suggested a frequency of the Weaver allele of 0.26 %. CONCLUSIONS: Through combined use of exhaustive sequencing data and SNP genotyping results, we were able to provide evidence that supports the non-synonymous mutation at position 49,878,773 bp as the most likely causal mutation for Weaver syndrome. Further studies are needed to uncover the exact mechanisms that underlie this syndrome.

The 80-kb DNA duplication on BTA1 is the only remaining candidate mutation for the polled phenotype of Friesian origin.

BACKGROUND: The absence of horns, called polled phenotype, is the favored trait in modern cattle husbandry. To date, polled cattle are obtained primarily by dehorning calves. Dehorning is a practice that raises animal welfare issues, which can be addressed by selecting for genetically hornless cattle. In the past 20 years, there have been many studies worldwide to identify unique genetic markers in complete association with the polled trait in cattle and recently, two different alleles at the POLLED locus, both resulting in the absence of horns, were reported: (1) the Celtic allele, which is responsible for the polled phenotype in most breeds and for which a single candidate mutation was detected and (2) the Friesian allele, which is responsible for the polled phenotype predominantly in the Holstein-Friesian breed and in a few other breeds, but for which five candidate mutations were identified in a 260-kb haplotype. Further studies based on genome-wide sequencing and high-density SNP (single nucleotide polymorphism) genotyping confirmed the existence of the Celtic and Friesian variants and narrowed down the causal Friesian haplotype to an interval of 145 kb. RESULTS: Almost 6000 animals were genetically tested for the polled trait and we detected a recombinant animal which enabled us to reduce the Friesian POLLED haplotype to a single causal mutation, namely a 80-kb duplication. Moreover, our results clearly disagree with the recently reported perfect co-segregation of the POLLED mutation and a SNP at position 1 390 292 bp on bovine chromosome 1 in the Holstein-Friesian population. CONCLUSION: We conclude that the 80-kb duplication, as the only remaining variant within the shortened Friesian haplotype, represents the most likely causal mutation for the polled phenotype of Friesian origin.

Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae.

Despite massive research efforts, the molecular etiology of bovine polledness and the developmental pathways involved in horn ontogenesis are still poorly understood. In a recent article, we provided evidence for the existence of at least two different alleles at the Polled locus and identified candidate mutations for each of them. None of these mutations was located in known coding or regulatory regions, thus adding to the complexity of understanding the molecular basis of polledness. We confirm previous results here and exhaustively identify the causative mutation for the Celtic allele (PC) and four candidate mutations for the Friesian allele (PF). We describe a previously unreported eyelash-and-eyelid phenotype associated with regular polledness, and present unique histological and gene expression data on bovine horn bud differentiation in fetuses affected by three different horn defect syndromes, as well as in wild-type controls. We propose the ectopic expression of a lincRNA in PC/p horn buds as a probable cause of horn bud agenesis. In addition, we provide evidence for an involvement of OLIG2, FOXL2 and RXFP2 in horn bud differentiation, and draw a first link between bovine, ovine and caprine Polled loci. Our results represent a first and important step in understanding the genetic pathways and key process involved in horn bud differentiation in Bovidae.

Bovine polledness--an autosomal dominant trait with allelic heterogeneity.

The persistent horns are an important trait of speciation for the family Bovidae with complex morphogenesis taking place briefly after birth. The polledness is highly favourable in modern cattle breeding systems but serious animal welfare issues urge for a solution in the production of hornless cattle other than dehorning. Although the dominant inhibition of horn morphogenesis was discovered more than 70 years ago, and the causative mutation was mapped almost 20 years ago, its molecular nature remained unknown. Here, we report allelic heterogeneity of the POLLED locus. First, we mapped the POLLED locus to a ∼381-kb interval in a multi-breed case-control design. Targeted re-sequencing of an enlarged candidate interval (547 kb) in 16 sires with known POLLED genotype did not detect a common allele associated with polled status. In eight sires of Alpine and Scottish origin (four polled versus four horned), we identified a single candidate mutation, a complex 202 bp insertion-deletion event that showed perfect association to the polled phenotype in various European cattle breeds, except Holstein-Friesian. The analysis of the same candidate interval in eight Holsteins identified five candidate variants which segregate as a 260 kb haplotype also perfectly associated with the POLLED gene without recombination or interference with the 202 bp insertion-deletion. We further identified bulls which are progeny tested as homozygous polled but bearing both, 202 bp insertion-deletion and Friesian haplotype. The distribution of genotypes of the two putative POLLED alleles in large semi-random sample (1,261 animals) supports the hypothesis of two independent mutations.

Mapping of a milk production quantitative trait locus to a 1.056 Mb region on bovine chromosome 5 in the Fleckvieh dual purpose cattle breed.

BACKGROUND: In a previous study in the Fleckvieh dual purpose cattle breed, we mapped a quantitative trait locus (QTL) affecting milk yield (MY1), milk protein yield (PY1) and milk fat yield (FY1) during first lactation to the distal part of bovine chromosome 5 (BTA5), but the confidence interval was too large for positional cloning of the causal gene. Our objective here was to refine the position of this QTL and to define the candidate region for high-throughput sequencing. METHODS: In addition to those previously studied, new Fleckvieh families were genotyped, in order to increase the number of recombination events. Twelve new microsatellites and 240 SNP markers covering the most likely QTL region on BTA5 were analysed. Based on haplotype analysis performed in this complex pedigree, families segregating for the low frequency allele of this QTL (minor allele) were selected. Single- and multiple-QTL analyses using combined linkage and linkage disequilibrium methods were performed. RESULTS: Single nucleotide polymorphism haplotype analyses on representative family sires and their ancestors revealed that the haplotype carrying the minor QTL allele is rare and most probably originates from a unique ancestor in the mapping population. Analyses of different subsets of families, created according to the results of haplotype analysis and availability of SNP and microsatellite data, refined the previously detected QTL affecting MY1 and PY1 to a region ranging from 117.962 Mb to 119.018 Mb (1.056 Mb) on BTA5. However, the possibility of a second QTL affecting only PY1 at 122.115 Mb was not ruled out. CONCLUSION: This study demonstrates that targeting families segregating for a less frequent QTL allele is a useful method. It improves the mapping resolution of the QTL, which is due to the division of the mapping population based on the results of the haplotype analysis and to the increased frequency of the minor allele in the families. Consequently, we succeeded in refining the region containing the previously detected QTL to 1 Mb on BTA5. This candidate region contains 27 genes with unknown or partially known function(s) and is small enough for high-throughput sequencing, which will allow future detailed analyses of candidate genes.

Genetic diversity of European cattle breeds highlights the conservation value of traditional unselected breeds with high effective population size.

In times of rapid global and unforeseeable environmental changes, there is an urgent need for a sustainable cattle breeding policy, based on a global view. Most of the indigenous breeds are specialized in a particular habitat or production system but are rapidly disappearing. Thus, they represent an important resource to meet present and future breeding objectives. Based on 105 microsatellites, we obtained thorough information on genetic diversity and population structure of 16 cattle breeds that cover a geographical area from the domestication centre near Anatolia, through the Balkan and alpine regions, to the North-West of Europe. Breeds under strict artificial selection and indigenous breeds under traditional breeding schemes were included. The overall results showed that the genetic diversity is widespread in Busa breeds in the Anatolian and Balkan areas, when compared with the alpine and north-western European breeds. Our results reflect long-term evolutionary and short-term breeding events very well. The regular pattern of allele frequency distribution in the entire cattle population studied clearly suggests conservation of rare alleles by conservation of preferably unselected traditional breeds with large effective population sizes. From a global and long-term conservation genetics point of view, the native and highly variable breeds closer to the domestication centre could serve as valuable sources of genes for future needs, not only for cattle but also for other farm animals.

Fine-mapping and candidate gene analysis of bovine spinal muscular atrophy.

Bovine spinal muscular atrophy (SMA), an autosomal recessive neurodegenerative disease, has been mapped at moderate resolution to the distal part of Chromosome 24. In this article we confirm this location and fine-map the SMA locus to an interval of approximately 0.8 cM at the very distal end of BTA24. Despite remarkable similarity to human SMA, the causative gene SMN can be excluded in bovine SMA. However, the interval where the disease now has been mapped contains BCL2, like SMN an antiapoptotic factor, and shown to bind to SMN. Moreover, knockout mice lacking the BCL2 gene show rapid motor neuron degeneration with early postnatal onset, as observed in bovine SMA. A comparative cattle/human map of the distal end of BTA24, based on the emerging bovine genome sequencing data, shows conserved synteny to HSA18 with hints of a segmental duplication and pericentric inversion just after the last available bovine marker DIK4971. This synteny lets us conclude that SMA is in immediate vicinity of the telomere. Candidate gene analysis of BCL2, however, excludes most of this gene, except its promoter region, and draws attention to the neighboring gene VPS4B, part of the endosomal protein-sorting machinery ESCRT-III which is involved in several neurodegenerative diseases.