Strategies for single nucleotide polymorphism (SNP) genotyping to enhance genotype imputation in Gyr (Bos indicus) dairy cattle: Comparison of commercially available SNP chips

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Synaptomenal complex analysis of four breeds of Bos taurus taurus x B. taurus indicus hybrids

The synaptonemal complex (SC) was analyzed in four F1 hybrids of Bos taurus taurus and B. taurus indicus including Gyr-Simmental (G-S), Nelore Simmental (N-S), Gyr-Holstein-Friesian (G-H) and Nelore-Piemontese (N-P). We analysed the frequency of various types of SC abnormalities and the frequency of cells with SC abnormalities. The results were compared with similar observations made on purebred animals. All the animals studied possessed 29 autosomal and one sex bivalent. The frequency of cells with abnormalities in the hybrids were 28.0% in the N-P, 29.1% in the G-S, 33.3% in the N-S and 40.0% in the G-H. The frequency of cells with abnormalities in the four hybrids was 31.5%; 57.9% of these abnormalities occurred in zygotene and 42.0% occurred in pachytene. The comparisons among the hybrids and among the hybrids and their parental breeds showed that the only significant difference was between Gyr and Gyr-Holstein-Friesian animals. Some aspects of the relationship between the frequency of cells with anomalies and the fertility of hybrids are discussed.

Genome-wide association study and ancestral origins of the slick-hair coat in tropically adapted cattle

The slick hair coat (SLICK) is a dominantly inherited trait typically associated with tropically adapted cattle that are from Criollo descent through Spanish colonization of cattle into the New World. The trait is of interest relative to climate change, due to its association with improved thermo-tolerance and subsequent increased productivity. Previous studies localized the SLICK locus to a 4 cM region on chromosome (BTA) 20 and identified signatures of selection in this region derived from Senepol cattle. The current study compares three slick-haired Criollo-derived breeds including Senepol, Carora, and Romosinuano and three additional slick-haired cross-bred lineages to non-slick ancestral breeds. Genome-wide association (GWA), haplotype analysis, signatures of selection, runs of homozygosity (ROH), and identity by state (IBS) calculations were used to identify a 0.8 Mb (37.7-38.5 Mb) consensus region for the SLICK locus on BTA20 in which contains SKP2 and SPEF2 as possible candidate genes. Three specific haplotype patterns are identified in slick individuals, all with zero frequency in non-slick individuals. Admixture analysis identified common genetic patterns between the three slick breeds at the SLICK locus. Principal component analysis (PCA) and admixture results show Senepol and Romosinuano sharing a higher degree of genetic similarity to one another with a much lesser degree of similarity to Carora. Variation in GWA, haplotype analysis, and IBS calculations with accompanying population structure information supports potentially two mutations, one common to Senepol and Romosinuano and another in Carora, effecting genes contained within our refined location for the SLICK locus.

Genomic analysis for managing small and endangered populations: a case study in Tyrol Grey cattle

Analysis of genomic data is increasingly becoming part of the livestock industry. Therefore, the routine collection of genomic information would be an invaluable resource for effective management of breeding programs in small, endangered populations. The objective of the paper was to demonstrate how genomic data could be used to analyse (1) linkage disequlibrium (LD), LD decay and the effective population size (NeLD); (2) Inbreeding level and effective population size (NeROH) based on runs of homozygosity (ROH); (3) Prediction of genomic breeding values (GEBV) using small within-breed and genomic information from other breeds. The Tyrol Grey population was used as an example, with the goal to highlight the potential of genomic analyses for small breeds. In addition to our own results we discuss additional use of genomics to assess relatedness, admixture proportions, and inheritance of harmful variants. The example data set consisted of 218 Tyrol Grey bull genotypes, which were all available AI bulls in the population. After standard quality control restrictions 34,581 SNPs remained for the analysis. A separate quality control was applied to determine ROH levels based on Illumina GenCall and Illumina GenTrain scores, resulting into 211 bulls and 33,604 SNPs. LD was computed as the squared correlation coefficient between SNPs within a 10 mega base pair (Mb) region. ROHs were derived based on regions covering at least 4, 8, and 16 Mb, suggesting that animals had common ancestors approximately 12, 6, and 3 generations ago, respectively. The corresponding mean inbreeding coefficients (F ROH) were 4.0% for 4 Mb, 2.9% for 8 Mb and 1.6% for 16 Mb runs. With an average generation interval of 5.66 years, estimated NeROH was 125 (NeROH>16 Mb), 186 (NeROH>8 Mb) and 370 (NeROH>4 Mb) indicating strict avoidance of close inbreeding in the population. The LD was used as an alternative method to infer the population history and the Ne. The results show a continuous decrease in NeLD, to 780, 120, and 80 for 100, 10, and 5 generations ago, respectively. Genomic selection was developed for and is working well in large breeds. The same methodology was applied in Tyrol Grey cattle, using different reference populations. Contrary to the expectations, the accuracy of GEBVs with very small within breed reference populations were very high, between 0.13-0.91 and 0.12-0.63, when estimated breeding values and deregressed breeding values were used as pseudo-phenotypes, respectively. Subsequent analyses confirmed the high accuracies being a consequence of low reliabilities of pseudo-phenotypes in the validation set, thus being heavily influenced by parent averages. Multi-breed and across breed reference sets gave inconsistent and lower accuracies. Genomic information may have a crucial role in management of small breeds, even if its primary usage differs from that of large breeds. It allows to assess relatedness between individuals, trends in inbreeding and to take decisions accordingly. These decisions would be based on the real genome architecture, rather than conventional pedigree information, which can be missing or incomplete. We strongly suggest the routine genotyping of all individuals that belong to a small breed in order to facilitate the effective management of endangered livestock populations.

Imputation of microsatellite alleles from dense SNP genotypes for parentage verification across multiple Bos taurus and Bos indicus breeds

To assist cattle producers transition from microsatellite (MS) to single nucleotide polymorphism (SNP) genotyping for parental verification we previously devised an effective and inexpensive method to impute MS alleles from SNP haplotypes. While the reported method was verified with only a limited data set (N = 479) from Brown Swiss, Guernsey, Holstein, and Jersey cattle, some of the MS-SNP haplotype associations were concordant across these phylogenetically diverse breeds. This implied that some haplotypes predate modern breed formation and remain in strong linkage disequilibrium. To expand the utility of MS allele imputation across breeds, MS and SNP data from more than 8000 animals representing 39 breeds (Bos taurus and B. indicus) were used to predict 9410 SNP haplotypes, incorporating an average of 73 SNPs per haplotype, for which alleles from 12 MS markers could be accurately be imputed. Approximately 25% of the MS-SNP haplotypes were present in multiple breeds (N = 2 to 36 breeds). These shared haplotypes allowed for MS imputation in breeds that were not represented in the reference population with only a small increase in Mendelian inheritance inconsistancies. Our reported reference haplotypes can be used for any cattle breed and the reported methods can be applied to any species to aid the transition from MS to SNP genetic markers. While ~91% of the animals with imputed alleles for 12 MS markers had ≤1 Mendelian inheritance conflicts with their parents' reported MS genotypes, this figure was 96% for our reference animals, indicating potential errors in the reported MS genotypes. The workflow we suggest autocorrects for genotyping errors and rare haplotypes, by MS genotyping animals whose imputed MS alleles fail parentage verification, and then incorporating those animals into the reference dataset.

Genetic diversity of the Brazilian Creole cattle Pe-duro assessed by microsatellites and mitochondrial DNA

The objective of this study was to describe the genetic diversity and structure of the largest Pe-duro population by assessing variation at ten autosomal microsatellite (STR) loci and mitochondrial DNA (mtDNA) sequences. The mean expected heterozygosity was 0.755, the mean observed heterozygosity was 0.600 and significant inbreeding coefficient (Fis) and deviations from the Hardy-Weinberg equilibrium in most of analyzed loci demonstrate the impact of inbreeding and homozygosis on this population. A more in-depth genetic analysis could be achieved by expanding the STR list. The analysis of mtDNA provided evidence of ancestral African taurine haplotypes in Pe-duro and excluded maternal Zebuine introgression. In this report, the main Pe-duro population is genetically portrayed by sampling approximately 40% of it. As this herd represents the core of the Pe-duro conservation program, these findings are of outstanding value for the management and preservation of this Brazilian 'native' cattle breed.

A nonsense mutation in the optic atrophy 3 gene (OPA3) causes dilated cardiomyopathy in Red Holstein cattle

Cardiomyopathies are severe degenerative disorders of the myocardium that lead to heart failure. During the last three decades bovine dilated cardiomyopathy (BDCMP) was observed worldwide in cattle of Holstein-Friesian origin. In the Swiss cattle population BDCMP affects Fleckvieh and Red Holstein breeds. The heart of affected animals is enlarged due to dilation of both ventricles. Clinical signs are caused by systolic dysfunction and affected individuals die as a result of severe heart insufficiency. BDCMP follows an autosomal recessive pattern of inheritance and the disease-causing locus was mapped to bovine chromosome 18 (BTA18). In the present study we describe the successful identification of the causative mutation in the OPA3 gene located on BTA18 that was previously reported to cause 3-methylglutaconic aciduria type III in Iraqi-Jewish patients. We demonstrated conclusive genetic and functional evidence that the nonsense mutation c.343C>T in the bovine OPA3 gene causes the late-onset dilated cardiomyopathy in Red Holstein cattle.

Dual origins of dairy cattle farming - evidence from a comprehensive survey of European Y-chromosomal variation

BACKGROUND: Diversity patterns of livestock species are informative to the history of agriculture and indicate uniqueness of breeds as relevant for conservation. So far, most studies on cattle have focused on mitochondrial and autosomal DNA variation. Previous studies of Y-chromosomal variation, with limited breed panels, identified two Bos taurus (taurine) haplogroups (Y1 and Y2; both composed of several haplotypes) and one Bos indicus (indicine/zebu) haplogroup (Y3), as well as a strong phylogeographic structuring of paternal lineages. METHODOLOGY AND PRINCIPAL FINDINGS: Haplogroup data were collected for 2087 animals from 138 breeds. For 111 breeds, these were resolved further by genotyping microsatellites INRA189 (10 alleles) and BM861 (2 alleles). European cattle carry exclusively taurine haplotypes, with the zebu Y-chromosomes having appreciable frequencies in Southwest Asian populations. Y1 is predominant in northern and north-western Europe, but is also observed in several Iberian breeds, as well as in Southwest Asia. A single Y1 haplotype is predominant in north-central Europe and a single Y2 haplotype in central Europe. In contrast, we found both Y1 and Y2 haplotypes in Britain, the Nordic region and Russia, with the highest Y-chromosomal diversity seen in the Iberian Peninsula. CONCLUSIONS: We propose that the homogeneous Y1 and Y2 regions reflect founder effects associated with the development and expansion of two groups of dairy cattle, the pied or red breeds from the North Sea and Baltic coasts and the spotted, yellow or brown breeds from Switzerland, respectively. The present Y1-Y2 contrast in central Europe coincides with historic, linguistic, religious and cultural boundaries.

PRNP promoter polymorphisms are associated with BSE susceptibility in Swiss and German cattle

BACKGROUND: Non-synonymous polymorphisms within the prion protein gene (PRNP) influence the susceptibility and incubation time for transmissible spongiform encephalopathies (TSE) in some species such as sheep and humans. In cattle, none of the known polymorphisms within the PRNP coding region has a major influence on susceptibility to bovine spongiform encephalopathy (BSE). Recently, however, we demonstrated an association between susceptibility to BSE and a 23 bp insertion/deletion (indel) polymorphism and a 12 bp indel polymorphism within the putative PRNP promoter region using 43 German BSE cases and 48 German control cattle. The objective of this study was to extend this work by including a larger number of BSE cases and control cattle of German and Swiss origin. RESULTS: Allele, genotype and haplotype frequencies of the two indel polymorphisms were determined in 449 BSE cattle and 431 unaffected cattle from Switzerland and Germany including all 43 German BSE and 16 German control animals from the original study. When breeds with similar allele and genotype distributions were compared, the 23 bp indel polymorphism again showed a significant association with susceptibility to BSE. However, some additional breed-specific allele and genotype distributions were identified, mainly related to the Brown breeds. CONCLUSION: Our study corroborated earlier findings that polymorphisms in the PRNP promoter region have an influence on susceptibility to BSE. However, breed-specific differences exist that need to be accounted for when analyzing such data.

DNA-based diagnosis of rare diseases in veterinary medicine: a 4.4 kb deletion of ITGB4 is associated with epidermolysis bullosa in Charolais cattle.

BACKGROUND Rare diseases in livestock animals are traditionally poorly diagnosed. Other than clinical description and pathological examination, the underlying causes have, for the most part, remained unknown. A single case of congenital skin fragility in cattle was observed, necropsy, histological and ultrastructural examinations were carried out and whole genome sequencing was utilized to identify the causative mutation. RESULTS A single purebred female Charolais calf with severe skin lesions was delivered full-term and died spontaneously after birth. The clinical and pathological findings exactly matched the gross description given by previous reports on epitheliogenesis imperfecta and epidermolysis bullosa (EB) in cattle. Histological and ultrastructural changes were consistent with EB junctionalis (EBJ). Genetic analysis revealed a previously unpublished ITGB4 loss-of-function mutation; the affected calf was homozygous for a 4.4 kb deletion involving exons 17 to 22, and the dam carried a single copy of the deletion indicating recessive inheritance. The homozygous mutant genotype did not occur in healthy controls of various breeds but some heterozygous carriers were found among Charolais cattle belonging to the affected herd. The mutant allele was absent in a representative sample of unrelated sires of the German Charolais population. CONCLUSION This is the first time in which a recessively inherited ITGB4 associated EBJ has been reported in cattle. The identification of heterozygous carriers is of importance in avoiding the transmission of this defect in future. Current DNA sequencing methods offer a powerful tool for understanding the genetic background of rare diseases in domestic animals having a reference genome sequence available.

Whole Genome Analysis of Pinkeye Scores in Angus Cattle

Infectious Bovine Keratoconjunctivitis (IBK), known as pinkeye, is a common infectious disease affecting the eyes of cattle. It is characterized by excessive tearing, inflammation of the conjunctiva, and ulceration of the cornea. Although pinkeye is non-fatal, it has a marked economic impact on the cattle industry, due to the decreased performance of infected individuals. Genetic effects on the susceptibility of IBK have been studied and Hereford, Jersey, and Holstein breeds were found to be more susceptible to IBK than Bos Indicus breeds. The objectives of our study were: 1) to estimate genetic parameters of IBK scored in different categories by using genomic threshold model, and 2) to detect markers in linkage disequilibrium with quantitative tract loci (QTL) associated with IBK.

Double muscling in cattle due to mutations in the myostatin gene

Myostatin (GDF-8) is a member of the transforming growth factor β superfamily of secreted growth and differentiation factors that is essential for proper regulation of skeletal muscle mass in mice. Here we report the myostatin sequences of nine other vertebrate species and the identification of mutations in the coding sequence of bovine myostatin in two breeds of double-muscled cattle, Belgian Blue and Piedmontese, which are known to have an increase in muscle mass relative to conventional cattle. The Belgian Blue myostatin sequence contains an 11-nucleotide deletion in the third exon which causes a frameshift that eliminates virtually all of the mature, active region of the molecule. The Piedmontese myostatin sequence contains a missense mutation in exon 3, resulting in a substitution of tyrosine for an invariant cysteine in the mature region of the protein. The similarity in phenotypes of double-muscled cattle and myostatin null mice suggests that myostatin performs the same biological function in these two species and is a potentially useful target for genetic manipulation in other farm animals.

Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1

Background Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation. Results Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue. The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen. Conclusions The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.

Cow-keeping in India : a simple and practical book on their care and treatment, their various breeds, and the means of rendering them profitable /

Includes index.

Observations on live stock: : containing hints for choosing and improving the best breeds of the most useful kinds of domestic animals. /

"Appendix."--p. 93-103.