The Complete Mitochondrial DNA Sequence of the Domestic Sheep (Ovis aries) and Comparison with the Other Major Ovine Haplotype

The complete mitochondrial DNA (mtDNA) molecule of the domestic sheep, Ovis aries, was sequenced, together with part of the mtDNA of a specimen representing the other major O. aries haplotype group. The length of the complete ovine mtDNA presented is 16,616 nucleotides (nt). This length is not absolute, however, due to heteroplasmy caused by the occurrence of different numbers of a 75-nt-long tandem repeat in the control region. The sequence data were included in analyses of intraspecific ovine molecular differences, molecular comparisons with bovine mtDNAs, and phylogenetic analyses based on complete mtDNAs. The comparisons with bovine mtDNAs were based on the central domains of the ovine control regions, representing both major ovine haplotype groups, and the corresponding domains of Bos taurus and B. indicus. The comparisons showed that the difference between the bovids was 1.4 times greater than the intraspecific ovine difference. These findings suggest that the strains of wild sheep from which domestic sheep originated were more closely related than were the B. primigenius subspecies which gave rise to B. indicus and B. taurus cattle. Datings based on complete mtDNAs suggest that the bovine and ovine lineages diverged about 30 million years before present. This dating is considerably earlier than that proposed previously. Received: 5 September 1997 / Accepted: 5 May 1998     

Occurrence of the Rumen Ciliate Oligoisotricha bubali in Domestic Cattle (Bos taurus)

Oligoisotricha bubali, previously observed twice in water buffalo, was detected in rumen contents of domestic cattle (Bos taurus) in two different areas of Tennessee. Concentrations ranged from <1 to 35% of the total protozoa in unweaned calves and up to 72% in older animals in feedlot. In contrast to the other genera of holotrichs, both total numbers and percent composition of O. bubali increased when animals were fed a corn silage-concentrate diet.     

Trace Element Distribution in Selected Edible Tissues of Zebu (Bos indicus) Cattle Slaughtered at Jimma,SW Ethiopia

The amount of trace elements present in edible bovine tissues is of importance for both animal health and human nutrition. This study presents data on trace element concentrations in semitendinosus and cardiac muscles, livers and kidneys of 60 zebu (Bos indicus) bulls, sampled at Jimma, Ethiopia. From 28 of these bulls, blood samples were also obtained. Deficient levels of copper were found in plasma, livers, kidneys and semitendinosus muscles. Suboptimal selenium concentrations were found in plasma and semitendinosus muscles. Semitendinosus muscles contained high iron concentrations. Trace elements were mainly stored in the liver, except for iron and selenium. Cardiac muscles generally contained higher concentrations of trace elements than semitendinous muscles except for zinc. A strong association was found between liver and kidney concentrations of copper, iron, cobalt and molybdenum. Liver storage was well correlated with storage in semitendinosus muscle for selenium and with cardiac muscle for cobalt and selenium. Plasma concentrations of copper, selenium, cobalt were well related with their respective liver concentrations and for cobalt and selenium, also with cardiac muscle concentrations. The data suggest multiple trace element deficiencies in zebu cattle in South-West Ethiopia, with lowered tissue concentrations as a consequence. Based on the comparison of our data with other literature, trace element concentrations in selected edible tissues of Bos indicus seem quite similar to those in Bos taurus. However, tissue threshold values for deficiency in Bos taurus cattle need to be refined and their applicability for Bos indicus cattle needs to be evaluated.     

Cattle demographic history modelled from autosomal sequence variation

The phylogeography of cattle genetic variants has been extensively described and has informed the history of domestication. However, there remains a dearth of demographic models inferred from such data. Here, we describe sequence diversity at 37 000 bp sampled from 17 genes in cattle from Africa, Europe and India. Clearly distinct population histories are suggested between Bos indicus and Bos taurus, with the former displaying higher diversity statistics. We compare the unfolded site frequency spectra in each to those simulated using a diffusion approximation method and build a best-fitting model of past demography. This implies an earlier, possibly glaciation-induced population bottleneck in B. taurus ancestry with a later, possibly domestication-associated demographic constriction in B. indicus. Strikingly, the modelled indicine history also requires a majority secondary admixture from the South Asian aurochs, indicating a complex, more diffuse domestication process. This perhaps involved multiple domestications and/or introgression from wild oxen to domestic herds; the latter is plausible from archaeological evidence of contemporaneous wild and domestic remains across different regions of South Asia.     

Detection of quantitative trait loci in Bos indicus and Bos taurus cattle using genome-wide association studies

Background

The apparent effect of a single nucleotide polymorphism (SNP) on phenotype depends on the linkage disequilibrium (LD) between the SNP and a quantitative trait locus (QTL). However, the phase of LD between a SNP and a QTL may differ between Bos indicus and Bos taurus because they diverged at least one hundred thousand years ago. Here, we test the hypothesis that the apparent effect of a SNP on a quantitative trait depends on whether the SNP allele is inherited from a Bos taurus or Bos indicus ancestor.

Methods

Phenotype data on one or more traits and SNP genotype data for 10 181 cattle from Bos taurus, Bos indicus and composite breeds were used. All animals had genotypes for 729 068 SNPs (real or imputed). Chromosome segments were classified as originating from B. indicus or B. taurus on the basis of the haplotype of SNP alleles they contained. Consequently, SNP alleles were classified according to their sub-species origin. Three models were used for the association study: (1) conventional GWAS (genome-wide association study), fitting a single SNP effect regardless of subspecies origin, (2) interaction GWAS, fitting an interaction between SNP and subspecies-origin, and (3) best variable GWAS, fitting the most significant combination of SNP and sub-species origin.

Results

Fitting an interaction between SNP and subspecies origin resulted in more significant SNPs (i.e. more power) than a conventional GWAS. Thus, the effect of a SNP depends on the subspecies that the allele originates from. Also, most QTL segregated in only one subspecies, suggesting that many mutations that affect the traits studied occurred after divergence of the subspecies or the mutation became fixed or was lost in one of the subspecies.

Conclusions

The results imply that GWAS and genomic selection could gain power by distinguishing SNP alleles based on their subspecies origin, and that only few QTL segregate in both B. indicus and B. taurus cattle. Thus, the QTL that segregate in current populations likely resulted from mutations that occurred in one of the subspecies and can have both positive and negative effects on the traits. There was no evidence that selection has increased the frequency of alleles that increase body weight.     

A novel polymorphism of GDF5 gene and its association with body measurement traits in <Emphasis Type="Italic">Bos taurus</Emphasis> and <Emphasis Type="Italic">Bos indicus</Emphasis> breeds

Body measurement traits, influenced by genes and environmental factors, play numerous important roles in the value assessment of productivity and economy. Growth differentiate factor 5 (GDF5), involved in the development and maintenance of bone and cartilage, is an important candidate gene for body measurement traits selection through marker-assisted selection (MAS). In this study, based on the PCR-RFLP technology, we discovered and evaluated the potential association of the single nucleotide polymorphism (SNP) (T586C in exon 1) of the bovine GDF5 gene with body measurement traits in 985 Bos taurus breed, 42 Bos indicus breed and 76 Bos indicus × Bos taurus individuals. As the SNP marker, there were the significant effects on the Body length (BL) in the Bos taurus (BT) and Bos indicus × Bos taurus (BMY) populations (P < 0.05). In BT population, animals with the genotype TT had lower mean values for BL and Hip width (HW) than these with the TC and CC genotype (P < 0.01). In BMY population, animals with the genotype TC had lower mean values for BL than these with the genotype CC (P < 0.05). These results suggest that the SNP of the GDF5 gene could be a very useful genetic marker for body measurement traits in the bovine reproduction and breeding.     

Expression of KCNQ1OT1, CDKN1C,H19, and PLAGL1 and the methylation patterns at the KvDMR1 and H19/IGF2 imprinting control regions is conserved between human and bovine

Background

Beckwith-Wiedemann syndrome (BWS) is a loss-of-imprinting pediatric overgrowth syndrome. The primary features of BWS include macrosomia, macroglossia, and abdominal wall defects. Secondary features that are frequently observed in BWS patients are hypoglycemia, nevus flammeus, polyhydramnios, visceromegaly, hemihyperplasia, cardiac malformations, and difficulty breathing. BWS is speculated to occur primarily as the result of the misregulation of imprinted genes associated with two clusters on chromosome 11p15.5, namely the KvDMR1 and H19/IGF2. A similar overgrowth phenotype is observed in bovine and ovine as a result of embryo culture. In ruminants this syndrome is known as large offspring syndrome (LOS). The phenotypes associated with LOS are increased birth weight, visceromegaly, skeletal defects, hypoglycemia, polyhydramnios, and breathing difficulties. Even though phenotypic similarities exist between the two syndromes, whether the two syndromes are epigenetically similar is unknown. In this study we use control Bos taurus indicus X Bos taurus taurus F1 hybrid bovine concepti to characterize baseline imprinted gene expression and DNA methylation status of imprinted domains known to be misregulated in BWS. This work is intended to be the first step in a series of experiments aimed at determining if LOS will serve as an appropriate animal model to study BWS.

Results

The use of F1 B. t. indicus x B. t. taurus tissues provided us with a tool to unequivocally determine imprinted status of the regions of interest in our study. We found that imprinting is conserved between the bovine and human in imprinted genes known to be associated with BWS. KCNQ1OT1 and PLAGL1 were paternally-expressed while CDKN1C and H19 were maternally-expressed in B. t. indicus x B. t. taurus F1 concepti. We also show that in bovids, differential methylation exists at the KvDMR1 and H19/IGF2 ICRs.

Conclusions

Based on these findings we conclude that the imprinted gene expression of KCNQ1OT1, CDKN1C, H19, and PLAGL1 and the methylation patterns at the KvDMR1 and H19/IGF2 ICRs are conserved between human and bovine. Future work will determine if LOS is associated with misregulation at these imprinted loci, similarly to what has been observed for BWS.     

A whole-genome assembly of the domestic cow, Bos taurus

Background

The genome of the domestic cow, Bos taurus, was sequenced using a mixture of hierarchical and whole-genome shotgun sequencing methods.

Results

We have assembled the 35 million sequence reads and applied a variety of assembly improvement techniques, creating an assembly of 2.86 billion base pairs that has multiple improvements over previous assemblies: it is more complete, covering more of the genome; thousands of gaps have been closed; many erroneous inversions, deletions, and translocations have been corrected; and thousands of single-nucleotide errors have been corrected. Our evaluation using independent metrics demonstrates that the resulting assembly is substantially more accurate and complete than alternative versions.

Conclusions

By using independent mapping data and conserved synteny between the cow and human genomes, we were able to construct an assembly with excellent large-scale contiguity in which a large majority (approximately 91%) of the genome has been placed onto the 30 B. taurus chromosomes. We constructed a new cow-human synteny map that expands upon previous maps. We also identified for the first time a portion of the B. taurus Y chromosome.     

Estimate of the population of preantral follicles in the ovaries of Bos taurus indicus and Bos taurus taurus cattle

The number of oocytes recovered from Bos taurus indicus females subjected to ovum pick-up averaged two to four times greater compared to Bos taurus taurus females. The objective of the present study was to test the hypothesis that this difference in oocyte yield was due to more preantral follicles in the ovaries of Bos indicus females. Ovaries (n = 64) from Nelore (Bos indicus) fetuses (n = 10), heifers (n = 12), and cows (n = 10), and Aberdeen Angus (Bos taurus) fetuses (n = 10), heifers (n = 12), and cows (n = 10) were cut longitudinally into halves, fixed, and processed for histological evaluation. The number of preantral follicles was estimated by counting them in each histological section, using the oocyte nucleus as a marker and employing a correction factor. The average number of preantral follicles in the ovaries of Bos indicus vs Bos taurus was (mean ± SD) 143,929 ± 64,028 vs 285,155 ± 325,195 for fetuses, 76,851 ± 78,605 vs 109,673 ± 86,078 for heifers, and 39,438 ± 31,017 vs 89,577 ± 86,315 for cows (P > 0.05). The number of preantral follicles varied greatly among individual animals within the same category, as well as between breeds. In conclusion, we inferred that the higher oocyte yield from Bos indicus females was not due to a greater ovarian reserve of preantral follicles. Therefore, mechanisms controlling follicle development after the preantral stage likely accounted for differences between Bos indicus and Bos taurus females in number of oocytes retrieved at ovum pick-up.     

Sequence comparison of mitochondrial tRNA genes and origin of light strand replication in Bos taurus and Nellore (Bos indicus) breeds

Although the complete bovine mitochondrial DNA molecule has been previously sequencedand sequence comparisons of the mitochondrial displacement loop have been performed, detailed sequence information is limited on coding regions of mitochondrial DNA within and among breeds of Bos taurus and Bos indicus. This study analysed polymorphism of the mitochondrial DNA transfer RNA genes for trypto-phan, alanine, asparagine, cysteine, tyrosine and the origin of light strand replication among Ayrshire, Canadian, Belgium Blue, Brown Swiss, Hereford, Jersey, Limousine, Piedmon-taise, Red Angus, Simmental (Bos taurus) and a Nellore (Bos indicus). Nucleotide sequence analysis of a 420-bp fragment of mitochondrial DNA comprising the five transfer RNA genes showed 100% homology among single individuals of the Bos taurus breeds. The Nellore breed showed guanine to adenine substitutions in the DHU arm of asparagine tRNA and in the origin of light-strand replication. This equates to a 0.5% sequence difference between the Nellore andBos taurus breeds and may reflect an independent evolutionary origin of the species.     

Tissue-Specific and Minor Inter-Individual Variation in Imprinting of IGF2R Is a Common Feature of Bos taurus 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<0.05) from endodermal tissues (liver, lung). The placenta showed partial imprinting with allelic bias of 22.9−34.7% and differed significantly (P<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.     

A catalogue of validated single nucleotide polymorphisms in bovine orthologs of mammalian imprinted genes and associations with beef production traits

Genetic (or ‘genomic’) imprinting, a feature of approximately 100 mammalian genes, results in monoallelic expression from one of the two parentally inherited chromosomes. To date, most studies have been directed on imprinted genes in murine or human models; however, there is burgeoning interest in the effects of imprinted genes in domestic livestock species. In particular, attention has focused on imprinted genes that influence foetal growth and development and that are associated with several economically important production traits in cattle, sheep and pigs. We have re-sequenced regions in 20 candidate bovine imprinted genes in order to validate single nucleotide polymorphisms (SNPs) that may influence important production traits in cattle. Putative SNPs detected via re-sequencing were subsequently re-formatted for high-throughput SNP genotyping in 185 cattle samples comprising 138 performance-tested European Bos taurus (all Limousin bulls), 29 African B. taurus and 18 Indian B. indicus samples. Analysis of the resulting genotypic data identified 117 validated SNPs. Preliminary genotype–phenotype association analyses using 83 SNPs that were polymorphic in the Limousin samples with minor allele frequencies ⩾0.05 revealed significant associations between two candidate bovine imprinted genes and a range of important beef production traits: average daily gain, average feed intake, live weight, feed conversion ratio, residual feed intake and residual gain. These genes were the Ras protein-specific guanine nucleotide releasing factor gene (RASGRF1) and the zinc finger, imprinted 2 gene (ZIM2). Despite the relatively small sample size used in these analyses, the observed associations with production traits are supported by the purported biological function of the RASGRF1 and ZIM2 gene products. These results support the hypothesis that imprinted genes contribute significantly to important complex production traits in cattle. Furthermore, these SNPs may be usefully incorporated into future marker-assisted and genomic selection breeding schemes.     

Mass spectrometric measurements of the apolipoproteins of bovine (Bos taurus) HDL

As is the case in most mammals, high density lipoproteins (HDL) also comprise the major group of lipid carriers that circulate in bovine (Bos taurus) blood. As a continuation of our proteogenomic studies of mammalian apolipoproteins, we have obtained molecular masses for several of the apolipoproteins associated with bovine HDL. The major apolipoprotein on the HDL surface is apoA-I, but other apolipoproteins were also detected. Using electrospray-ionization mass spectrometry (ESI-MS), we report on values for apolipoproteins, A-I, proA-I and A-II, as well as post-translationally modified apoA-I. Analyses of tryptic fragments did reveal the presence of apoA-IV and apoC-III. However, in contrast to our previous studies of other mammalian HDL, we did not detect apoC-I. Interestingly, examination of the current assembly for the bovine genome does not show any evidence for an apoC-I gene.     

The production of mixed-species Bos taurus-Bos indicus twin calves

Bos taurus-Bos indicus twin calves were produced by transferring Day 6–8 Bos indicus or Bos taurus embryos to previously inseminated Bos taurus or Bos indicus cows. Embryos were transferred to 112 recipients of which 64 (57.1%) were diagnosed as pregnant by rectal palpation at 7–8 weeks gestation. Twin foetuses were diagnosed in 38 (65.5%) of 58 cows in which the location and number of foetuses was recorded. Substantial foetal losses occurred after 7–8 weeks and only 44 cows calved giving 21 sets of twins and 23 singletons. Six of the single calves were from transferred embryos. The duration of gestation of the twin pregnancies (282 days) was similar to that of single Bos taurus pregnancies (282 days for Jerseys and 286 days for Friesians) and 14 days shorter than single Bos indicus pregnancies (296 days). Consequently Bos indicus members of mixed-species twins were premature at birth and required intensive care to ensure their survival.     

Somatic cell mapping of conglutinin (CGN1) to cattle syntenic group U29 and fluorescence in situ localization to Chromosome 28

A 260-bp genomic PstI fragment, which encodes a portion of the carbohydrate recognition domain, was used along with hybrid somatic cells to map the conglutinin gene (CGN1) to domestic cow (Bos taurus) syntenic group U29. In turn, a cosmid containing the entire bovine CGN1 was used with fluorescence in situ hybridization to sublocalize this gene to cattle chromosome (Chr) (BTA) 28 band 18. Since BTA 28 and several of the other small acrocentric autosomes of cattle are difficult to discriminate, we have also chromosomally sublocalized CGN1 to the p arm of the lone biarmed autosome of the gaur (Bos gaurus). The use of the gaur 2/28 Robertsonian as a marker chromosome and our assignment of CGN1 to BTA 28 should help resolve some of the nomenclatural questions involving this cattle chromosome.     

Serosurvey of Smooth Brucella,Leptospira spp. and Toxoplasma gondii in Free-Ranging Jaguars (Panthera onca) and Domestic Animals from Brazil

This study investigated the exposure of jaguar populations and domestic animals to smooth Brucella, Leptospira spp. and Toxoplasma gondii in the Cerrado, Pantanal and Amazon biomes of Brazil. Between February 2000 and January 2010, serum samples from 31 jaguars (Panthera onca), 1,245 cattle (Bos taurus), 168 domestic dogs (Canis lupus familiaris) and 29 domestic cats (Felis catus) were collected and analysed by rose bengal test for smooth Brucella, microscopic agglutination test for Leptospira spp. and modified agglutination test for T. gondii. Cattle populations from all sites (9.88%) were exposed to smooth Brucella, but only one jaguar from Cerrado was exposed to this agent. Jaguars captured in the Cerrado (60.0%) and in the Pantanal (45.5%) were seropositive for different serovars of Leptospira spp., cattle (72.18%) and domestic dogs (13.1%) from the three sites and one domestic cat from Pantanal were also seropositive for the agent. The most prevalent serotype of Leptospira spp. identified in jaguars from the Cerrado (Grippotyphosa) and the Pantanal (Pomona) biomes were distinct from those found in the domestic animals sampled. Jaguars (100%), domestic dogs (38.28%) and domestic cats (82.76%) from the three areas were exposed to T. gondii. Our results show that brucellosis and leptospirosis could have been transmitted to jaguars by domestic animals; and jaguars probably play an important role in the maintenance of T. gondii in nature.     

Introgression,admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.     

Sequence and evolution of related bovine and caprine satellite DNAs: Identification of a short DNA sequence potentially involved in satellite DNA amplification

The satellite II DNAs of the domestic ox Bos taurus and sheep Ovis aries have been sequenced, and that of the domestic goat Capra hircus partially sequenced. All three are related, and consist of repeat units of about 700 base-pairs. There is no evidence of internal repetition within these repeat units. When matched for maximum homology, the goat and sheep sequences show 83% homology, whereas the ox and sheep sequences share only 70% homology. Factors contributing to the uncertainty of the exact homology between these sequences are discussed, but the results are nevertheless consistent with their progenitor sequence being present in the common ancestor of cattle and sheep. Goat satellite II DNA is shown to contain another, unrelated, tandemly repeated sequence, which is composed of 22 base-pair repeat units. Both this sequence and a region of ox satellite II share good homology with the 11 base-pair progenitor sequence of ox 1.706 g/cm³ satellite DNA. It is suggested that this shared sequence could play a role in bovine satellite DNA amplification.