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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Methods</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>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.</p

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

Bovine polledness

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

A heuristic two-dimensional presentation of microsatellite-based data applied to dogs and wolves

Methods based on genetic distance matrices usually lose information during the process of tree-building by converting a multi-dimensional matrix into a phylogenetic tree. We applied a heuristic method of two-dimensional presentation to achieve a better resolution of the relationship between breeds and individuals investigated. Four hundred and nine individuals from nine German dog breed populations and one free-living wolf population were analysed with a marker set of 23 microsatellites. The result of the two-dimensional presentation was partly comparable with and complemented a model-based analysis that uses genotype patterns. The assignment test and the neighbour-joining tree based on allele sharing estimate allocated 99% and 97% of the individuals according to their breed, respectively. The application of the two-dimensional presentation to distances on the basis of the proportion of shared alleles resulted in comparable and further complementary insight into inferred population structure by multilocus genotype data. We expect that the inference of population structure in domesticated species with complex breeding histories can be strongly supported by the two-dimensional presentation based on the described heuristic method

First Insight Into Genetic Diversity of Alpine ibex (Capra ibex) in Slovenia

In Europe, the Alpine ibex (Capra ibex) was on the brink of extinction in the 19th century. Therefore, different conservation measures were implemented, and several reintroductions were made in the Alpine arc, starting from the only surviving population in Gran Paradiso, Italy. An extreme historical bottleneck and additional reintroductions have strongly shaped the genetic make-up of recent populations, resulting in significant genetic drift and profound inbreeding across the species range. To support science-based conservation actions, molecular methods have been increasingly used. However, such analyses did not include populations in Slovenia. We analysed neutral loci (partial fragment of mitochondrial cytochrome b, mtDNA) and the adaptive major histocompatibility complex (MHC DRB exon 2) of the Alpine ibex from both Slovenian populations (Julian and Kamnik-Savinja Alps) to understand how past reintroductions and recent management have affected the genetic diversity of the species. Results showed that both populations are genetically severely depleted, carrying only one mtDNA haplotype and one functional allele for MHC DRB exon 2, Caib-DRB*01. This calls for further conservation actions, including the reintroduction of individuals with different genetic background. However, the Alpine ibex is currently considered a non-native species in Slovenia, which makes conservation actions extremely difficult and threatens the long-term survival of the species. Therefore, scientists and population managers are urging policy/decision makers to change the status of the species to the native one and consequently to allow reintroductions. These appeals are supported by previous archaeological data on the existence of bones assigned to Alpine ibex in the Julian Alps, and evidence of severe genetic depletion in current ibex populations confirmed in this study. PRVI VPOGLED V GENETSKO RAZNOLIKOST ALPSKEGA KOZOROGA (Capra ibex) V SLOVENIJI Izvleček: V Evropi je bil alpski kozorog (Capra ibex) v 19. stoletju na robu izumrtja. Izvajali so se različni ukrepi za njegovo ohranjanje. V alpskem loku je bilo izvedenih več ponovnih naselitev, najprej z edino ohranjeno populacijo v kraju Gran Paradiso v Italiji. Izredna okrnjenost vrste v preteklosti in dodatne ponovne naselitve so močno vplivale na genetsko sestavo populacije, kar je povzročilo znaten genetski zdrs in parjenje v sorodstvu na celotnem območju vrste. V podporo znanstveno utemeljenim ukrepom ohranjanja se vse pogosteje uporabljajo molekularne metode, vendar takšne analize niso vključevale populacij v Sloveniji. Da bi razumeli, kako je ponovno naseljevanje in nedavno upravljanje vplivalo na genetsko raznolikost vrste, smo analizirali nevtralne lokuse (delni fragment mitohondrijskega citokroma b, mtDNA) in adaptivni poglavitni histokompatibilnostni kompleks (MHC DRB ekson 2) alpskega kozoroga iz obeh slovenskih populacij (Julijske in Kamniško-Savinjske Alpe). Rezultati so pokazali, da sta obe populaciji genetsko zelo osiromašeni, saj nosita le en haplotip mtDNA in en funkcionalni alel za MHC DRB ekson 2, Caib-DRB*01. Zato so potrebni nadaljnji ukrepi za ohranjanje, vključno s ponovno naselitvijo živali z drugačnim genetskim ozadjem. Vendar alpski kozorog v Sloveniji trenutno velja za tujerodno vrsto, kar zelo otežuje ukrepe za njegovo ohranitev in ogroža dolgoročno preživetje vrste. Znanstveniki in upravljavci populacij zato pozivajo politike/odločevalce, naj spremenijo status vrste v avtohtono in posledično omogočijo ponovno naselitev. Ti pozivi so podprti s predhodnimi arheološkimi podatki o obstoju kosti alpskega kozoroga v Julijskih Alpah in z dokazi o izraziti genetski osiromašenosti sedanjih populacij kozoroga, potrjenimi v tej študiji. Ključne besede: Capra ibex; mitohondrijska DNA; MHC DRB exon2; ponovna naselitev; upravljanj

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

Recovery of Native Genetic Background in Admixed Populations Using Haplotypes, Phenotypes, and Pedigree Information - Using Cika Cattle as a Case Breed

The aim of this study was to obtain unbiased estimates of the diversity parameters, the population history, and the degree of admixture in Cika cattle which represents the local admixed breeds at risk of extinction undergoing challenging conservation programs. Genetic analyses were performed on the genome-wide Single Nucleotide Polymorphism (SNP) Illumina Bovine SNP50 array data of 76 Cika animals and 531 animals from 14 reference populations. To obtain unbiased estimates we used short haplotypes spanning four markers instead of single SNPs to avoid an ascertainment bias of the BovineSNP50 array. Genome-wide haplotypes combined with partial pedigree and type trait classification show the potential to improve identification of purebred animals with a low degree of admixture. Phylogenetic analyses demonstrated unique genetic identity of Cika animals. Genetic distance matrix presented by rooted Neighbour-Net suggested long and broad phylogenetic connection between Cika and Pinzgauer. Unsupervised clustering performed by the admixture analysis and two-dimensional presentation of the genetic distances between individuals also suggest Cika is a distinct breed despite being similar in appearance to Pinzgauer. Animals identified as the most purebred could be used as a nucleus for a recovery of the native genetic background in the current admixed population. The results show that local well-adapted strains, which have never been intensively managed and differentiated into specific breeds, exhibit large haplotype diversity. They suggest a conservation and recovery approach that does not rely exclusively on the search for the original native genetic background but rather on the identification and removal of common introgressed haplotypes would be more powerful. Successful implementation of such an approach should be based on combining phenotype, pedigree, and genome-wide haplotype data of the breed of interest and a spectrum of reference breeds which potentially have had direct or indirect historical contribution to the genetic makeup of the breed of interest

Tissue-specific and minor inter-individual variation in imprinting of <i>IGF2R</i> is a common feature of <i>Bos taurus</i> concepti and not correlated with fetal weight

The insulin-like growth factor 2 receptor (IGF2R) is essential for prenatal growth regulation and shows gene dosage effects on fetal weight that can be affected by in-vitro embryo culture. Imprinted maternal expression of murine Igf2r is well documented for all fetal tissues excluding brain, but polymorphic imprinting and biallelic expression were reported for IGF2R in human. These differences have been attributed to evolutionary changes correlated with specific reproductive strategies. However, data from species suitable for testing this hypothesis are lacking. The domestic cow (Bos taurus) carries a single conceptus with a similar gestation length as human. We identified 12 heterozygous concepti informative for imprinting studies among 68 Bos taurus fetuses at Day 80 of gestation (28% term) and found predominantly maternal IGF2R expression in all fetal tissues but brain, which escapes imprinting. Inter-individual variation in allelic expression bias, i.e. expression of the repressed paternal allele relative to the maternal allele, ranged from 4.6−8.9% in heart, 4.3−10.2% in kidney, 6.1−11.2% in liver, 4.6−15.8% in lung and 3.2−12.2% in skeletal muscle. Allelic bias for mesodermal tissues (heart, skeletal muscle) differed significantly (P&lt;0.05) from endodermal tissues (liver, lung). The placenta showed partial imprinting with allelic bias of 22.9−34.7% and differed significantly (P&lt;0.001) from all other tissues. Four informative fetuses were generated by in-vitro fertilization (IVF) with embryo culture and two individuals displayed fetal overgrowth. However, there was no evidence for changes in imprinting or DNA methylation after IVF, or correlations between allelic bias and fetal weight. In conclusion, imprinting of Bos taurus IGF2R is similar to mouse except in placenta, which could indicate an effect of reproductive strategy. Common minor inter-individual variation in allelic bias and absence of imprinting abnormalities in IVF fetuses suggest changes in IGF2R expression in overgrown fetuses could be modulated through other mechanisms than changes in imprinting

Remapping of the belted phenotype in cattle on BTA3 identifies a multiplication event as the candidate causal mutation

Background: It has been known for almost a century that the belted phenotype in cattle follows a pattern of dominant inheritance. In 2009, the approximate position of the belt locus in Brown Swiss cattle was mapped to a 922-kb interval on bovine chromosome 3 and, subsequently, assigned to a 336-kb haplotype block based on an animal set that included, Brown Swiss, Dutch Belted (Lakenvelder) and Belted Galloway individuals. A possible candidate gene in this region i.e. HES6 was investigated but the causal mutation remains unknown. Thus, to elucidate the causal mutation of this prominent coat color phenotype, we decided to remap the belted phenotype in an independent animal set of several European bovine breeds, i.e. Gurtenvieh (belted Brown Swiss), Dutch Belted and Belted Galloway and to systematically scan the candidate region. We also checked the presence of the detected causal mutation in the genome of belted individuals from a Siberian cattle breed. Results: A combined linkage disequilibrium and linkage analysis based on 110 belted and non-belted animals identified a candidate interval of 2.5 Mb. Manual inspection of the haplotypes in this region identified four candidate haplotypes that consisted of five to eight consecutive SNPs. One of these haplotypes overlapped with the initial 922-kb interval, whereas two were positioned proximal and one was positioned distal to this region. Next-generation sequencing of one heterozygous and two homozygous belted animals identified only one private belted candidate allele, i.e. a multiplication event that is located between 118,608,000 and 118,614,000 bp. Targeted locus amplification and quantitative real-time PCR confirmed an increase in copy number of this region in the genomes of both European (Belted Galloway, Dutch Belted and Gurtenvieh) and Siberian (Yakutian cattle) breeds. Finally, using nanopore sequencing, the exact breakpoints were determined at 118,608,362 and 118,614,132 bp. The closest gene to the candidate causal mutation (16 kb distal) is TWIST2. Conclusions: Based on our findings and those of a previously published study that identified the same multiplication event, a quadruplication on bovine chromosome 3 between positions 118,608,362 and 118,614,132 bp is the most likely candidate causal mutation for the belted phenotype in cattle

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 Bu\va. 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