Project description:The Toll-like receptor (TLR) and peptidoglycan recognition protein 1 (PGLYRP1) genes play key roles in the innate immune systems of mammals. While the TLRs recognize a variety of invading pathogens and induce innate immune responses, PGLYRP1 is directly microbicidal. We used custom allele-specific assays to genotype and validate 220 diallelic variants, including 54 nonsynonymous SNPs in 11 bovine innate immune genes (TLR1-TLR10, PGLYRP1) for 37 cattle breeds. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and we were unable to differentiate between the specialized B. t. taurus beef and dairy breeds, despite an average polymorphism density of one locus per 219 bp. Ninety-nine tagSNPs and one tag insertion-deletion polymorphism were sufficient to predict 100% of the variation at all 11 innate immune loci in both subspecies and their hybrids, whereas 58 tagSNPs captured 100% of the variation at 172 loci in B. t. taurus. PolyPhen and SIFT analyses of nonsynonymous SNPs encoding amino acid replacements indicated that the majority of these substitutions were benign, but up to 31% were expected to potentially impact protein function. Several diversity-based tests provided support for strong purifying selection acting on TLR10 in B. t. taurus cattle. These results will broadly impact efforts related to bovine translational genomics.
Project description:BackgroundWe present here the assembly of the bovine genome. The assembly method combines the BAC plus WGS local assembly used for the rat and sea urchin with the whole genome shotgun (WGS) only assembly used for many other animal genomes including the rhesus macaque.ResultsThe assembly process consisted of multiple phases: First, BACs were assembled with BAC generated sequence, then subsequently in combination with the individual overlapping WGS reads. Different assembly parameters were tested to separately optimize the performance for each BAC assembly of the BAC and WGS reads. In parallel, a second assembly was produced using only the WGS sequences and a global whole genome assembly method. The two assemblies were combined to create a more complete genome representation that retained the high quality BAC-based local assembly information, but with gaps between BACs filled in with the WGS-only assembly. Finally, the entire assembly was placed on chromosomes using the available map information.Over 90% of the assembly is now placed on chromosomes. The estimated genome size is 2.87 Gb which represents a high degree of completeness, with 95% of the available EST sequences found in assembled contigs. The quality of the assembly was evaluated by comparison to 73 finished BACs, where the draft assembly covers between 92.5 and 100% (average 98.5%) of the finished BACs. The assembly contigs and scaffolds align linearly to the finished BACs, suggesting that misassemblies are rare. Genotyping and genetic mapping of 17,482 SNPs revealed that more than 99.2% were correctly positioned within the Btau_4.0 assembly, confirming the accuracy of the assembly.ConclusionThe biological analysis of this bovine genome assembly is being published, and the sequence data is available to support future bovine research.
Project description:With regulatory roles in development, cell proliferation and disease, micro-RNA (miRNA) biology is of great importance and a potential key to novel RNA-based therapeutic regimens. Biochemically based sequencing approaches have provided robust means of uncovering miRNA binding landscapes on transcriptomes of various species. However, a current limitation to the therapeutic potential of miRNA biology in cattle is the lack of validated miRNAs targets. Here, we use cross-linking immunoprecipitation (CLIP) of the Argonaute (AGO) proteins and unambiguous miRNA-target identification through RNA chimeras to define a regulatory map of miRNA interactions in the cow (Bos taurus). The resulting interactome is the deepest reported to date for any species, demonstrating that comprehensive maps can be empirically obtained. We observe that bovine miRNA targeting principles are consistent with those observed in other mammals. Motif and structural analyses define expanded pairing rules with most interactions combining seed-based pairing with distinct, miRNA-specific patterns of auxiliary pairing. Further, miRNA-target chimeras had predictive value in evaluating true regulatory sites of the miR-17 family. Finally, we define miRNA-specific targeting for >5000 mRNAs and determine gene ontologies (GO) for these targets. This confirmed repression of genes important for embryonic development and cell cycle progress by the let-7 family, and repression of those involved in cell cycle arrest by the miR-17 family, but it also suggested a number of unappreciated miRNA functions. Our results provide a significant resource for transcriptomic understanding of bovine miRNA regulation, and demonstrate the power of experimental methods for establishing comprehensive interaction maps.
Project description:Proteins are the major constituents of muscle and are key molecules regulating the metabolic changes during conversion of muscle to meat. Brazil is one of the largest exporters of beef and most Brazilian cattle are composed by zebu (Nellore) genotype. Bos indicus beef is generally leaner and tougher than Bos taurus such as Angus. The aim of this study was to compare the muscle proteomic and phosphoproteomic profile of Angus and Nellore. Seven animals of each breed previously subjected the same growth management were confined for 84 days. Proteins were extracted from Longissimus lumborum samples collected immediately after slaughter and separated by two-dimensional electrophoresis. Pro-Q Diamond stain was used in phosphoproteomics. Proteins identification was performed using matrix assisted laser desorption/ionization time-of-flight mass spectrometry. Tropomyosin alpha-1 chain, troponin-T, myosin light chain-1 fragment, cytoplasmic malate dehydrogenase, alpha-enolase and 78 kDa glucose-regulated protein were more abundant in Nellore, while myosin light chain 3, prohibitin, mitochondrial stress-70 protein and heat shock 70 kDa protein 6 were more abundant in Angus (P<0.05). Nellore had higher phosphorylation of myosin regulatory light chain-2, alpha actin-1, triosephosphate isomerase and 14-3-3 protein epsilon. However, Angus had greater phosphorylation of phosphoglucomutase-1 and troponin-T (P<0.05). Therefore, proteins involved in contraction and muscle organization, myofilaments expressed in fast or slow-twitch fibers and heat shock proteins localized in mitochondria or sarcoplasmic reticulum and involved in cell flux of calcium and apoptosis might be associated with differences in beef quality between Angus and Nellore. Furthermore, prohibitin appears to be a potential biomarker of intramuscular fat in cattle. Additionally, differences in phosphorylation of myofilaments and glycolytic enzymes could be involved with differences in muscle contraction force, susceptibility to calpain, apoptosis and postmortem glycolysis, which might also be related to differences in beef quality among Angus and Nellore.
Project description:Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.
Project description:BackgroundReference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus-B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus-B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references.ResultsThe alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity.ConclusionsWe report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.