Genetic mapping of the belt pattern in Brown Swiss cattle to BTA3

Copy number variation analysis reveals variants associated with milk production traits in dairy goats

Copy number variation (CNV) is a major type of genomic structural variation. We investigated their impacts on goat dairy traits using the CaprineSNP50 array. From 120 samples of five dairy goat breeds, we totally identified 42 CNVs ranging from 56,044 bp to 4,337,625 bp. We found significant associations between two CNVs (CNV5 and CNV25) and two milk production traits (mean of milk fat yield and mean of milk protein yield) after false discovery rate (FDR) correction (P < 0.05). CNV5 overlaps the ADAMTS20 gene, which is involved in the differentiation of mammary cell and plays a crucial role in lactogenic activity of bovine mammary epithelial cells. CNV25 overlaps with PAPPA2, which has been found to be associated with bovine reproduction and milk production traits. Our results revealed that CNVs overlapped with ADAMTS20 and PAPPA2 could be involved in goat dairy traits and function as candidate markers for further genetic selection.

A non-coding regulatory variant in the 5'-region of the MITF gene is associated with white-spotted coat in Brown Swiss cattle

Recently, the Swiss breeding association reported an increasing number of white-spotted cattle in the Brown Swiss breed, which is normally solid brown coloured. A total of 60 Brown Swiss cattle with variably sized white abdominal spots, facial markings and depigmented claws were collected for this study. A genome-wide association study using 40k SNP genotypes of 20 cases and 1619 controls enabled us to identify an associated genome region on chromosome 22 containing the MITF gene, encoding the melanogenesis associated transcription factor. Variants at the MITF locus have been reported before to be associated with white or white-spotted phenotypes in other species such as horses, dogs and mice. Whole-genome sequencing of a single white-spotted cow and subsequent genotyping of 172 Brown Swiss cattle revealed two significantly associated completely linked single nucleotide variants (rs722765315 and rs719139527). Both variants are located in the 5'-regulatory region of the bovine MITF gene, and comparative sequence analysis showed that the variant rs722765315, located 139 kb upstream of the transcription start site of the bovine melanocyte-specific MITF transcript, is situated in a multi-species conserved sequence element which is supposed to be regulatory important. Therefore, we hypothesize that rs722765315 represents the most likely causative variant for the white-spotting phenotype observed in Brown Swiss cattle. Presence of the mutant allele in a heterozygous or homozygous state supports a dominant mode of inheritance with incomplete penetrance and results in a variable extent of coat colour depigmentation.

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.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0407-9) contains supplementary material, which is available to authorized users.

Genes related to the white face colour pattern in eight Russian cattle breeds

One of the major effects of domestication is change of animal coat colour to up to complete white colour of the whole body. It is possible that white colour of livestock animals had aesthetic significance for humans as well. The first step towards detection of genes and mutations controlling white colouring in animals is the genome-wide association studies. These studies, however, have not been done for the cattle breeds native to the Russian Federation. The aim of this study was therefore to identify genomic intervals and candidate genes that could be responsible for white face colouring in eight Russian cattle breeds. The data on genome-wide genotyping of 131,709 high-quality single nucleotide polymorphisms (SNPs) on 148 animas have been used in the program ­EMMAX. Association analysis has been performed using two related phenotypes: a) the white face with the rest of the body of any colour and b) white face with the rest of the body of different (non-white) colour. In the first case, the only statistically significant marker found was the SNP BovineHD0500019319 located on cattle chromosome (BTA) 5. The same SNP was the most significant within the cluster of three SNPs on BTA5: 68,803,879–69,365,854 associated also with the second phenotype. Five genes were found within this interval in the cattle genome, out of which the most likely functional candidate was SLC41A2, with the SNP BovineHD0500019319 found within its intronic sequence. SLC41A2 encodes a magnesium transporter protein. However, the function of this gene is not well established. Other members of this gene family are the key genes controlling differences in human skin and animal coat colour. Additional significant association signals with the second phenotype have been detected in BTA 1–4, 6–15, 18, 19, 24, 27, and 29. Overall, 37 genomic intervals have been detected associated with white face colouring in eight Russian native cattle breeds.

A structural variant in the 5'-flanking region of the TWIST2 gene affects melanocyte development in belted cattle

Belted cattle have a circular belt of unpigmented hair and skin around their midsection. The belt is inherited as a monogenic autosomal dominant trait. We mapped the causative variant to a 37 kb segment on bovine chromosome 3. Whole genome sequence data of 2 belted and 130 control cattle yielded only one private genetic variant in the critical interval in the two belted animals. The belt-associated variant was a copy number variant (CNV) involving the quadruplication of a 6 kb non-coding sequence located approximately 16 kb upstream of the TWIST2 gene. Increased copy numbers at this CNV were strongly associated with the belt phenotype in a cohort of 333 cases and 1322 controls. We hypothesized that the CNV causes aberrant expression of TWIST2 during neural crest development, which might negatively affect melanoblasts. Functional studies showed that ectopic expression of bovine TWIST2 in neural crest in transgenic zebrafish led to a decrease in melanocyte numbers. Our results thus implicate an unsuspected involvement of TWIST2 in regulating pigmentation and reveal a non-coding CNV underlying a captivating Mendelian character.

Serial translocation by means of circular intermediates underlies colour sidedness in cattle

Colour sidedness is a dominantly inherited phenotype of cattle characterized by the polarization of pigmented sectors on the flanks, snout and ear tips. It is also referred to as 'lineback' or 'witrik' (which means white back), as colour-sided animals typically display a white band along their spine. Colour sidedness is documented at least since the Middle Ages and is presently segregating in several cattle breeds around the globe, including in Belgian blue and brown Swiss. Here we report that colour sidedness is determined by a first allele on chromosome 29 (Cs(29)), which results from the translocation of a 492-kilobase chromosome 6 segment encompassing KIT to chromosome 29, and a second allele on chromosome 6 (Cs(6)), derived from the first by repatriation of fused 575-kilobase chromosome 6 and 29 sequences to the KIT locus. We provide evidence that both translocation events involved circular intermediates. This is the first example, to our knowledge, of a phenotype determined by homologous yet non-syntenic alleles that result from a novel copy-number-variant-generating mechanism.

Genetic heterogeneity at the bovine KIT gene in cattle breeds carrying different putative alleles at the spotting locus

According to classical genetic studies, piebaldism in cattle is largely influenced by the allelic series at the spotting locus (S), which includes the S(H) (Hereford pattern), S(+) (non-spotted) and s (spotted) alleles. The S locus was mapped on bovine chromosome 6 in the region containing the KIT gene. We investigated the KIT gene, analysing its variability and haplotype distribution in cattle of three breeds (Angus, Hereford and Holstein) with different putative alleles (S(+), S(H) and s respectively) at the S locus. Resequencing of a whole of 0.485 Mb revealed 111 polymorphisms. The global nucleotide diversity was 0.087%. Tajima's D-values were negative for all breeds, indicating putative directional selection. Of the 28 inferred haplotypes, only five were observed in the Hereford breed, in which one was the most frequent. Coalescent simulation showed that it is highly unlikely (P < 10E-6) to obtain this low number of haplotypes conditionally on the observed number of segregating SNPs. Therefore, the neutral model could be rejected for the Hereford breed, suggesting that a selection sweep occurred at the KIT locus. Twelve haplotypes were inferred in Holstein and Angus. For these two breeds, the neutral model could not be rejected. High heterogeneity of the KIT gene was confirmed from a phylogenetic analysis. Our results suggest a role of the KIT gene in determining the S(H) allele(s) in the Hereford, but no evidence of selective sweep was obtained in Holstein, suggesting that complex mechanisms (or other genes) might be the cause of the spotted phenotype in this breed.