Unlocking the bovine genome

Vitamin C Alleviates the Negative Effects of Heat Stress on Reproductive Processes by Regulating Amino Acid Metabolism in Granulosa Cells

Heat stress-induced biochemical alterations in ovarian follicles compromise the function of granulosa cells (GCs) and the developmental competence of oocytes. Summer heat stress can have a far-reaching negative impact on overall fertility and reproductive success. Together with the heat stress, the rise of assisted reproductive technologies (ART), potential confounding hazards of in vitro handling and the absence of systemic body support in ART makes it imperative to study the heat stress ameliorative effects of vitamin C under in vitro conditions. Using in vitro heat stress treatment of 43 °C for two hours in bovine GCs, we studied the effects of vitamin C on cell growth, oxidative stress, apoptosis and cell cycle progression together with a comprehensive metabolomics profiling. This study investigates the molecular milieu underlying the vitamin C (VC)-led alleviation of heat-related disruptions to metabolic processes in bovine GCs. The supplementation of VC ameliorated the detrimental effects of heat stress by reducing oxidative stress and apoptosis while restoring cell proliferation. Normal cell function restoration in treated GCs was demonstrated through the finding of significantly high levels of progesterone. We observed a shift in the metabolome from biosynthesis to catabolism, mostly dominated by the metabolism of amino acids (decreased tryptophan, methionine and tyrosine) and the active TCA cycle through increased Succinic acid. The Glutathione and tryptophan metabolism were important in ameliorating the inflammation and metabolism nexus under heat stress. Two significant enzymes were identified, namely tryptophan 2,3-dioxygenase (TDO2) and mitochondrial phenylalanyl-tRNA synthetase (FARS2). Furthermore, our findings provide insight into the significance of B-complex vitamins in the context of heat stress during VC supplementation. This study underscores the importance of VC supplementation in heat stress and designates multiple metabolic intervention faucets in the context of ameliorating heat stress and enhancing reproductive efficiency.

Progress and future perspectives of livestock genomics in India: a mini review

The field of genetics has evolved a lot after the emergence of molecular and advanced genomic technologies. The advent of Next Generation Sequencing, SNP genotyping platforms and simultaneous reduction in the cost of sequencing had opened the door to genomic research in farm animals. There are various applications of genomics in livestock, such as the use of genomic data: (i) to investigate genetic diversity and breed composition/population structure (ii) to identify genetic variants and QTLs related to economically important and ecological traits, genome-wide association studies (GWAS) and genomic signatures of selection; (iii) to enhance breeding programs by genomic selection. Compared to traditional methods, genomic selection is expected to improve selection response by increasing selection accuracy and reducing the generation interval due to early selection. Genomic selection (GS) in developed countries has led to rapid genetic gains, especially in dairy cattle, due to a well-established genetic evaluation system. Indian livestock system is still lagging behind developed nations in adopting these technologies. This review discusses the current status, challenges, and future perspectives of livestock genomics in India.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Investigating the impact of reference assembly choice on genomic analyses in a cattle breed

BACKGROUND

Reference-guided read alignment and variant genotyping are prone to reference allele bias, particularly for samples that are greatly divergent from the reference genome. A Hereford-based assembly is the widely accepted bovine reference genome. Haplotype-resolved genomes that exceed the current bovine reference genome in quality and continuity have been assembled for different breeds of cattle. Using whole genome sequencing data of 161 Brown Swiss cattle, we compared the accuracy of read mapping and sequence variant genotyping as well as downstream genomic analyses between the bovine reference genome (ARS-UCD1.2) and a highly continuous Angus-based assembly (UOA_Angus_1).

RESULTS

Read mapping accuracy did not differ notably between the ARS-UCD1.2 and UOA_Angus_1 assemblies. We discovered 22,744,517 and 22,559,675 high-quality variants from ARS-UCD1.2 and UOA_Angus_1, respectively. The concordance between sequence- and array-called genotypes was high and the number of variants deviating from Hardy-Weinberg proportions was low at segregating sites for both assemblies. More artefactual INDELs were genotyped from UOA_Angus_1 than ARS-UCD1.2 alignments. Using the composite likelihood ratio test, we detected 40 and 33 signatures of selection from ARS-UCD1.2 and UOA_Angus_1, respectively, but the overlap between both assemblies was low. Using the 161 sequenced Brown Swiss cattle as a reference panel, we imputed sequence variant genotypes into a mapping cohort of 30,499 cattle that had microarray-derived genotypes using a two-step imputation approach. The accuracy of imputation (Beagle R2) was very high (0.87) for both assemblies. Genome-wide association studies between imputed sequence variant genotypes and six dairy traits as well as stature produced almost identical results from both assemblies.

CONCLUSIONS

The ARS-UCD1.2 and UOA_Angus_1 assemblies are suitable for reference-guided genome analyses in Brown Swiss cattle. Although differences in read mapping and genotyping accuracy between both assemblies are negligible, the choice of the reference genome has a large impact on detecting signatures of selection that already reached fixation using the composite likelihood ratio test. We developed a workflow that can be adapted and reused to compare the impact of reference genomes on genome analyses in various breeds, populations and species.

From Beef to Bees: High-Throughput Kinome Analysis to Understand Host Responses of Livestock Species to Infectious Diseases and Industry-Associated Stress

Within human health research, the remarkable utility of kinase inhibitors as therapeutics has motivated efforts to understand biology at the level of global cellular kinase activity (the kinome). In contrast, the diminished potential for using kinase inhibitors in food animals has dampened efforts to translate this research approach to livestock species. This, in our opinion, was a lost opportunity for livestock researchers given the unique potential of kinome analysis to offer insight into complex biology. To remedy this situation, our lab developed user-friendly, cost-effective approaches for kinome analysis that can be readily incorporated into most research programs but with a specific priority to enable the technology to livestock researchers. These contributions include the development of custom software programs for the creation of species-specific kinome arrays as well as comprehensive deconvolution and analysis of kinome array data. Presented in this review are examples of the application of kinome analysis to highlight the utility of the technology to further our understanding of two key complex biological events of priority to the livestock industry: host immune responses to infectious diseases and animal stress responses. These advances and examples of application aim to provide both mechanisms and motivation for researchers, particularly livestock researchers, to incorporate kinome analysis into their research programs.

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.

Comparison of long noncoding RNA between muscles and adipose tissues in Hanwoo beef cattle

Long noncoding RNAs (lncRNAs) regulate the expression of mRNA and can affect various biological processes and phenotypes. Currently, studies of lncRNAs in cattle are under way, but their exact function for several tissues has not yet been established. Hanwoo cattle (Bos taurus coreanae) have inhabited the Korean peninsula for about 6000 years and are one of the representative domesticated animals in Korea. As a result of intensive breeding, the meat of Hanwoo cattle is high in marbling content and is preferred by Koreans and other East Asian people. In this study, the expression of lncRNAs was identified in 36 samples from skeletal muscle and three adipose tissues (intramuscular, subcutaneous, and omental) of nine Hanwoo individuals. We identified 76 tissue-specific lncRNAs for each of the four tissues using the differences in expression levels. Through QTL information, we could identify 12 lncRNAs associated with shear force and six lncRNAs associated with body weight, which are two important traits in the Hanwoo population breeding strategy. By the physical position comparison of lncRNA and Bovine transcripts information, we could identify 11 lncRNAs that were in bovine transcripts, and four of the 11 genes related to transcripts of lncRNAs were biologically associated with muscle function. We believe this Hanwoo lncRNAs study will help reveal the lncRNA role in the physiological mechanisms of these four tissues.

Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled

Brahman cattle have a Bos indicus and Bos taurus mosaic genome, as a result of the process used to create the breed (repeat backcrossing of Bos taurus females to Bos indicus bulls). With the aim of identifying Bos taurus segments in the Brahman genome at sequence level resolution, we sequenced the genomes of 46 influential Brahman bulls. Using 36 million variants identified in the sequences, we searched for regions close to fixation for Bos indicus or Bos taurus segments that were longer than expected by chance (from simulation of the breed formation history of Brahman cattle). Regions close to fixation for Bos indicus content were enriched for protein synthesis genes, while regions of higher Bos taurus content included genes of the G-protein coupled receptor family (including genes implicated in puberty, such as THRS). The region with the most extreme Bos taurus enrichment was on chromosome 14 surrounding PLAG1. The introgressed Bos taurus allele at PLAG1 increases stature and the high frequency of the allele likely reflects strong selection for the trait. Finally, we provide evidence that the polled mutation in Brahmans, a desirable trait under very strong recent selection, is of Celtic origin and is introgressed from Bos taurus.

Transformation of animal genomics by next-generation sequencing technologies: a decade of challenges and their impact on genetic architecture

For more than a quarter of a century, sequencing technologies from Sanger's method to next-generation high-throughput techniques have provided fascinating opportunities in the life sciences. The continuing upward trajectory of sequencing technologies will improve livestock research and expedite the development of various new genomic and technological studies with farm animals. The use of high-throughput technologies in livestock research has increased interest in metagenomics, epigenetics, genome-wide association studies, and identification of single nucleotide polymorphisms and copy number variations. Such studies are beginning to provide revolutionary insights into biological and evolutionary processes. Farm animals, such as cattle, swine, and horses, have played a dual role as economically and agriculturally important animals as well as biomedical research models. The first part of this study explores the current state of sequencing methods, many of which are already used in animal genomic studies, and the second part summarizes the state of cattle, swine, horse, and chicken genome sequencing and illustrates its achievements during the last few years. Finally, we describe several high-throughput sequencing approaches for the improved detection of known, unknown, and emerging infectious agents, leading to better diagnosis of infectious diseases. The insights from viral metagenomics and the advancement of next-generation sequencing will strongly support specific and efficient vaccine development and provide strategies for controlling infectious disease transmission among animal populations and/or between animals and humans. However, prospective sequencing technologies will require further research and in-field testing before reaching the marketplace.

Vitamin K2 improves proliferation and migration of bovine skeletal muscle cells in vitro

Skeletal muscle function is highly dependent on the ability to regenerate, however, during ageing or disease, the proliferative capacity is reduced, leading to loss of muscle function. We have previously demonstrated the presence of vitamin K2 in bovine skeletal muscles, but whether vitamin K has a role in muscle regulation and function is unknown. In this study, we used primary bovine skeletal muscle cells, cultured in monolayers in vitro, to assess a potential effect of vitamin K2 (MK-4) during myogenesis of muscle cells. Cell viability experiments demonstrate that the amount of ATP produced by the cells was unchanged when MK-4 was added, indicating viable cells. Cytotoxicity analysis show that MK-4 reduced the lactate dehydrogenase (LDH) released into the media, suggesting that MK-4 was beneficial to the muscle cells. Cell migration, proliferation and differentiation was characterised after MK-4 incubation using wound scratch analysis, immunocytochemistry and real-time PCR analysis. Adding MK-4 to the cells led to an increased muscle proliferation, increased gene expression of the myogenic transcription factor myod as well as increased cell migration. In addition, we observed a reduction in the fusion index and relative gene expression of muscle differentiation markers, with fewer complex myotubes formed in MK-4 stimulated cells compared to control cells, indicating that the MK-4 plays a significant role during the early phases of muscle proliferation. Likewise, we see the same pattern for the relative gene expression of collagen 1A, showing increased gene expression in proliferating cells, and reduced expression in differentiating cells. Our results also suggest that MK-4 incubation affect low density lipoprotein receptor-related protein 1 (LRP1) and the low-density lipoprotein receptor (LDLR) with a peak in gene expression after 45 min of MK-4 incubation. Altogether, our experiments show that MK-4 has a positive effect on muscle cell migration and proliferation, which are two important steps during early myogenesis.

Applying the One Health Concept to Mycobacterial Research - Overcoming Parochialism

Mycobacterial infections remain a public health problem. Historically important, globally ubiquitous and with a wide host range, we are still struggling to control mycobacterial infections in humans and animals. While previous reviews have focused on individual mycobacterial infections in either humans or animals, a comprehensive review of the zoonotic aspect of mycobacteria in the context of the One Health initiative is lacking. With the purpose of providing a concise and comprehensive resource, we have collated literature to address the zoonotic potential of different mycobacterial species and elaborate on the necessity for an inter-sectorial approach to attain a new vision to combat mycobacterial infections.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

In vitro secretion of zymogens by bovine pancreatic acini and ultra-structural analysis of exocytosis

The aim of this study is to establish a bovine pancreatic acinar cell culture model with longer viability and functionality. The cells could be maintained in a functional state for upto 20 days with normal morphology. Cells were positive for amylase as observed by immunofluorescence staining. Acinar cells are spherical and range about 2–3 µm in diameter. The porosome formed by exocytosis and heterogenous enzyme granules of size ranging 100–300 nm were seen on the surface of cells by electron microscopy. The activity of the enzymes was high on day 15 and the activity profile of the enzymes is in the order: protease>lipase>amylase and the enzymes were identified by SDS-PAGE. Long-term culture of bovine pancreatic acini could be useful in studying the pathogenesis of pancreatitis. Since the bovine genome shares about 80% identity with the human genome, the cells derived from bovine pancreas can be engineered and used as a potential xenotransplant to treat conditions like pancreatitis as the tissue source is abundantly available.

•

Acinar cells are functional unit of pancreas.

•

In vitro release of enzyme granules by exocytosis are well demonstrated.

•

Amylase shares 86% protein identity with that of human.

•

Acinar culture showed longer viability upto 20 days with normal cellular function.

•

Useful culture model in understanding the acinar cell biology and physiology.

Genome-wide association study for claw disorders and trimming status in dairy cattle

Performing a genome-wide association study (GWAS) might add to a better understanding of the development of claw disorders and the need for trimming. Therefore, the aim of the current study was to perform a GWAS on claw disorders and trimming status and to validate the results for claw disorders based on an independent data set. Data consisted of 20,474 cows with phenotypes for claw disorders and 50,238 cows with phenotypes for trimming status. Recorded claw disorders used in the current study were double sole (DS), interdigital hyperplasia (IH), sole hemorrhage (SH), sole ulcer (SU), white line separation (WLS), a combination of infectious claw disorders consisting of (inter-)digital dermatitis and heel erosion, and a combination of laminitis-related claw disorders (DS, SH, SU, and WLS). Of the cows with phenotypes for claw disorders, 1,771 cows were genotyped and these cow data were used for the GWAS on claw disorders. A SNP was considered significant when the false discovery rate≤0.05 and suggestive when the false discovery rate≤0.20. An independent data set of 185 genotyped bulls having at least 5 daughters with phenotypes (6,824 daughters in total) for claw disorders was used to validate significant and suggestive SNP detected based on the cow data. To analyze the trait "trimming status" (i.e., the need for claw trimming), a data set with 327 genotyped bulls having at least 5 daughters with phenotypes (18,525 daughters in total) was used. Based on the cow data, in total 10 significant and 45 suggestive SNP were detected for claw disorders. The 10 significant SNP were associated with SU, and mainly located on BTA8. The suggestive SNP were associated with DS, IH, SU, and laminitis-related claw disorders. Three of the suggestive SNP were validated in the data set of 185 bulls, and were located on BTA13, BTA14, and BTA17. For infectious claw disorders, SH, and WLS, no significant or suggestive SNP associations were detected. For trimming status, 1 significant and 1 suggestive SNP were detected, both located close to each other on BTA15. Some significant and suggestive SNP were located close to SNP detected in studies on feet and leg conformation traits. Genes with major effects could not be detected and SNP associations were spread across the genome, indicating that many SNP, each explaining a small proportion of the genetic variance, influence claw disorders. Therefore, to reduce the incidence of claw disorders by breeding, genomic selection is a promising approach.

The bovine lung in biomedical research: visually guided bronchoscopy, intrabronchial inoculation and in vivo sampling techniques

There is an ongoing search for alternative animal models in research of respiratory medicine. Depending on the goal of the research, large animals as models of pulmonary disease often resemble the situation of the human lung much better than mice do. Working with large animals also offers the opportunity to sample the same animal repeatedly over a certain course of time, which allows long-term studies without sacrificing the animals.

The aim was to establish in vivo sampling methods for the use in a bovine model of a respiratory Chlamydia psittaci infection. Sampling should be performed at various time points in each animal during the study, and the samples should be suitable to study the host response, as well as the pathogen under experimental conditions.

Bronchoscopy is a valuable diagnostic tool in human and veterinary medicine. It is a safe and minimally invasive procedure. This article describes the intrabronchial inoculation of calves as well as sampling methods for the lower respiratory tract. Videoendoscopic, intrabronchial inoculation leads to very consistent clinical and pathological findings in all inoculated animals and is, therefore, well-suited for use in models of infectious lung disease. The sampling methods described are bronchoalveolar lavage, bronchial brushing and transbronchial lung biopsy. All of these are valuable diagnostic tools in human medicine and could be adapted for experimental purposes to calves aged 6-8 weeks. The samples obtained were suitable for both pathogen detection and characterization of the severity of lung inflammation in the host.

Genomic architecture of histone 3 lysine 27 trimethylation during late ovine skeletal muscle development

The ruminant developmental transition from late foetus to lamb is associated with marked changes in skeletal muscle structure and function that reflect programming for new physiological demands following birth. To determine whether epigenetic changes are involved in this transition, we investigated the genomic architecture of the chromatin modification, histone 3 lysine 27 trimethylation (H3K27me3), which typically regulates early life developmental processes; however, its role in later life processes is unclear. Chromatin immunoprecipitation coupled with next‐generation sequencing was used to map H3K27me3 nucleosomes in ovine longissimus lumborum skeletal muscle at 100 days of gestation and 12 weeks post‐partum. In both states, H3K27me3 modification was associated with genes, transcription start sites and CpG islands and with transcriptional silencing. The H3K27me3 peaks consisted of two major categories, promoter specific and regional, with the latter the dominant feature. Genes encoding homeobox transcription factors regulating early life development and genes involved in neural functions, particularly gated ion channels, were strongly modified by H3K27me3. Gene promoters differentially modified by H3K27me3 in the foetus and lamb were enriched for gated ion channels, which may reflect changes in neuromuscular function. However, most modified genes showed no changes, indicating that H3K27me3 does not have a large role in late muscle maturation. Notably, promyogenic transcription factors were strongly modified with H3K27me3 but showed no differences between the late gestation foetus and lamb, likely reflecting their lack of involvement in the myofibre fusion process occurring in this transition. H3K27me3 is a major architectural feature of the epigenetic landscape of ruminant skeletal muscle, and it comments on gene transcription and gene function in the context of late skeletal muscle development.

The combination of glycosaminoglycans and fibrous proteins improves cell proliferation and early differentiation of bovine primary skeletal muscle cells

Primary muscle cell model systems from farm animals are widely used to acquire knowledge about muscle development, muscle pathologies, overweight issues and tissue regeneration. The morphological properties of a bovine primary muscle cell model system, in addition to cell proliferation and differentiation features, were characterized using immunocytochemistry, western blotting and real-time PCR. We observed a reorganization of the Golgi complex in differentiated cells. The Golgi complex transformed to a highly fragmented network of small stacks of cisternae positioned throughout the myotubes as well as around the nucleus. Different extracellular matrix (ECM) components were used as surface coatings in order to improve cell culture conditions. Our experiments demonstrated improved proliferation and early differentiation for cells grown on surface coatings containing a mixture of both glycosaminoglycans (GAGs) and fibrous proteins. We suggest that GAGs and fibrous proteins mixed together into a composite biomaterial can mimic a natural ECM, and this could improve myogenesis for in vitro cell cultures.

Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity

Background

Hanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo.

Results

The Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. ‘Regulation of actin filament length’ was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus.

Conclusions

We identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.

Mis-assembled "segmental duplications" in two versions of the Bos taurus genome

We analyzed the whole genome sequence coverage in two versions of the Bos taurus genome and identified all regions longer than five kilobases (Kbp) that are duplicated within chromosomes with >99% sequence fidelity in both copies. We call these regions High Fidelity Duplications (HFDs). The two assemblies were Btau 4.2, produced by the Human Genome Sequencing Center at Baylor College of Medicine, and UMD Bos taurus 3.1 (UMD 3.1), produced by our group at the University of Maryland. We found that Btau 4.2 has a far greater number of HFDs, 3111 versus only 69 in UMD 3.1. Read coverage analysis shows that 39 million base pairs (Mbp) of sequence in HFDs in Btau 4.2 appear to be a result of a mis-assembly and therefore cannot be qualified as segmental duplications. UMD 3.1 has only 0.41 Mbp of sequence in HFDs that are due to a mis-assembly.

Identification and functional analysis of three isoforms of bovine BST-2

Human BST-2 (hBST-2) has been identified as a cellular antiviral factor that blocks the release of various enveloped viruses. Orthologues of BST-2 have been identified in several species, including human, monkeys, pig, mouse, cat and sheep. All have been reported to possess antiviral activity. Duplication of the BST-2 gene has been observed in sheep and the paralogues are referred to as ovine BST-2A and BST2-B, although only a single gene corresponding to BST-2 has been identified in most species. In this study, we identified three isoforms of bovine BST-2, named bBST-2A1, bBST-2A2 and bBST-2B, in bovine cells treated with type I interferon, but not in untreated cells. Both bBST-2A1 and bBST-2A2 are posttranslationally modified by N-linked glycosylation and a GPI-anchor as well as hBST-2, while bBST-2B has neither of these modifications. Exogenous expression of bBST-2A1 or bBST-2A2 markedly reduced the production of bovine leukemia virus and vesicular stomatitis virus from cells, while the antiviral activity of bBST-2B was much weaker than those of bBST-2A1 and bBST-2A2. Our data suggest that bBST-2A1 and bBST-2A2 function as part of IFN-induced innate immunity against virus infection. On the other hand, bBST-2B may have a different physiological function from bBST-2A1 and bBST-2A2.

A bovine model of respiratory Chlamydia psittaci infection: challenge dose titration

This study aimed to establish and evaluate a bovine respiratory model of experimentally induced acute C. psittaci infection. Calves are natural hosts and pathogenesis may resemble the situation in humans. Intrabronchial inoculation of C. psittaci strain DC15 was performed in calves aged 2–3 months via bronchoscope at four different challenge doses from 10⁶ to 10⁹ inclusion-forming units (ifu) per animal. Control groups received either UV-inactivated C. psittaci or cell culture medium. While 10⁶ ifu/calf resulted in a mild respiratory infection only, the doses of 10⁷ and 10⁸ induced fever, tachypnea, dry cough, and tachycardia that became apparent 2–3 days post inoculation (dpi) and lasted for about one week. In calves exposed to 10⁹ ifu C. psittaci, the respiratory disease was accompanied by severe systemic illness (apathy, tremor, markedly reduced appetite). At the time point of most pronounced clinical signs (3 dpi) the extent of lung lesions was below 10% of pulmonary tissue in calves inoculated with 10⁶ and 10⁷ ifu, about 15% in calves inoculated with 10⁸ and more than 30% in calves inoculated with 10⁹ ifu C. psittaci. Beside clinical signs and pathologic lesions, the bacterial load of lung tissue and markers of pulmonary inflammation (i.e., cell counts, concentration of proteins and eicosanoids in broncho-alveolar lavage fluid) were positively associated with ifu of viable C. psittaci. While any effect of endotoxin has been ruled out, all effects could be attributed to infection by the replicating bacteria. In conclusion, the calf represents a suitable model of respiratory chlamydial infection. Dose titration revealed that both clinically latent and clinically manifest infection can be reproduced experimentally by either 10⁶ or 10⁸ ifu/calf of C. psittaci DC15 while doses above 10⁸ ifu C. psittaci cannot be recommended for further studies for ethical reasons. This defined model of different clinical expressions of chlamydial infection allows studying host-pathogen interactions.

Transcriptome and proteome research in veterinary science: what is possible and what questions can be asked?

In recent years several technologies for the complete analysis of the transcriptome and proteome have reached a technological level which allows their routine application as scientific tools. The principle of these methods is the identification and quantification of up to ten thousands of RNA and proteins species in a tissue, in contrast to the sequential analysis of conventional methods such as PCR and Western blotting. Due to their technical progress transcriptome and proteome analyses are becoming increasingly relevant in all fields of biological research. They are mainly used for the explorative identification of disease associated complex gene expression patterns and thereby set the stage for hypothesis-driven studies. This review gives an overview on the methods currently available for transcriptome analysis, that is, microarrays, Ref-Seq, quantitative PCR arrays and discusses their potentials and limitations. Second, the most powerful current approaches to proteome analysis are introduced, that is, 2D-gel electrophoresis, shotgun proteomics, MudPIT and the diverse technological concepts are reviewed. Finally, experimental strategies for biomarker discovery, experimental settings for the identification of prognostic gene sets and explorative versus hypothesis driven approaches for the elucidation of diseases associated genes and molecular pathways are described and their potential for studies in veterinary research is highlighted.

Genetic architecture of gene expression in ovine skeletal muscle

Background

In livestock populations the genetic contribution to muscling is intensively monitored in the progeny of industry sires and used as a tool in selective breeding programs. The genes and pathways conferring this genetic merit are largely undefined. Genetic variation within a population has potential, amongst other mechanisms, to alter gene expression via cis- or trans-acting mechanisms in a manner that impacts the functional activities of specific pathways that contribute to muscling traits. By integrating sire-based genetic merit information for a muscling trait with progeny-based gene expression data we directly tested the hypothesis that there is genetic structure in the gene expression program in ovine skeletal muscle.

Results

The genetic performance of six sires for a well defined muscling trait, longissimus lumborum muscle depth, was measured using extensive progeny testing and expressed as an Estimated Breeding Value by comparison with contemporary sires. Microarray gene expression data were obtained for longissimus lumborum samples taken from forty progeny of the six sires (4-8 progeny/sire). Initial unsupervised hierarchical clustering analysis revealed strong genetic architecture to the gene expression data, which also discriminated the sire-based Estimated Breeding Value for the trait. An integrated systems biology approach was then used to identify the major functional pathways contributing to the genetics of enhanced muscling by using both Estimated Breeding Value weighted gene co-expression network analysis and a differential gene co-expression network analysis. The modules of genes revealed by these analyses were enriched for a number of functional terms summarised as muscle sarcomere organisation and development, protein catabolism (proteosome), RNA processing, mitochondrial function and transcriptional regulation.

Conclusions

This study has revealed strong genetic structure in the gene expression program within ovine longissimus lumborum muscle. The balance between muscle protein synthesis, at the levels of both transcription and translation control, and protein catabolism mediated by regulated proteolysis is likely to be the primary determinant of the genetic merit for the muscling trait in this sheep population. There is also evidence that high genetic merit for muscling is associated with a fibre type shift toward fast glycolytic fibres. This study provides insight into mechanisms, presumably subject to strong artificial selection, that underpin enhanced muscling in sheep populations.

Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping

Background

Integration of genomic variation with phenotypic information is an effective approach for uncovering genotype-phenotype associations. This requires an accurate identification of the different types of variation in individual genomes.

Results

We report the integration of the whole genome sequence of a single Holstein Friesian bull with data from single nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) array technologies to determine a comprehensive spectrum of genomic variation. The performance of resequencing SNP detection was assessed by combining SNPs that were identified to be either in identity by descent (IBD) or in copy number variation (CNV) with results from SNP array genotyping. Coding insertions and deletions (indels) were found to be enriched for size in multiples of 3 and were located near the N- and C-termini of proteins. For larger indels, a combination of split-read and read-pair approaches proved to be complementary in finding different signatures. CNVs were identified on the basis of the depth of sequenced reads, and by using SNP and CGH arrays.

Conclusions

Our results provide high resolution mapping of diverse classes of genomic variation in an individual bovine genome and demonstrate that structural variation surpasses sequence variation as the main component of genomic variability. Better accuracy of SNP detection was achieved with little loss of sensitivity when algorithms that implemented mapping quality were used. IBD regions were found to be instrumental for calculating resequencing SNP accuracy, while SNP detection within CNVs tended to be less reliable. CNV discovery was affected dramatically by platform resolution and coverage biases. The combined data for this study showed that at a moderate level of sequencing coverage, an ensemble of platforms and tools can be applied together to maximize the accurate detection of sequence and structural variants.

Tracing cattle breeds with principal components analysis ancestry informative SNPs

The recent release of the Bovine HapMap dataset represents the most detailed survey of bovine genetic diversity to date, providing an important resource for the design and development of livestock production. We studied this dataset, comprising more than 30,000 Single Nucleotide Polymorphisms (SNPs) for 19 breeds (13 taurine, three zebu, and three hybrid breeds), seeking to identify small panels of genetic markers that can be used to trace the breed of unknown cattle samples. Taking advantage of the power of Principal Components Analysis and algorithms that we have recently described for the selection of Ancestry Informative Markers from genomewide datasets, we present a decision-tree which can be used to accurately infer the origin of individual cattle. In doing so, we present a thorough examination of population genetic structure in modern bovine breeds. Performing extensive cross-validation experiments, we demonstrate that 250-500 carefully selected SNPs suffice in order to achieve close to 100% prediction accuracy of individual ancestry, when this particular set of 19 breeds is considered. Our methods, coupled with the dense genotypic data that is becoming increasingly available, have the potential to become a valuable tool and have considerable impact in worldwide livestock production. They can be used to inform the design of studies of the genetic basis of economically important traits in cattle, as well as breeding programs and efforts to conserve biodiversity. Furthermore, the SNPs that we have identified can provide a reliable solution for the traceability of breed-specific branded products.

Bovine proteins containing poly-glutamine repeats are often polymorphic and enriched for components of transcriptional regulatory complexes

Background

About forty human diseases are caused by repeat instability mutations. A distinct subset of these diseases is the result of extreme expansions of polymorphic trinucleotide repeats; typically CAG repeats encoding poly-glutamine (poly-Q) tracts in proteins. Polymorphic repeat length variation is also apparent in human poly-Q encoding genes from normal individuals. As these coding sequence repeats are subject to selection in mammals, it has been suggested that normal variations in some of these typically highly conserved genes are implicated in morphological differences between species and phenotypic variations within species. At present, poly-Q encoding genes in non-human mammalian species are poorly documented, as are their functions and propensities for polymorphic variation.

Results

The current investigation identified 178 bovine poly-Q encoding genes (Q ≥ 5) and within this group, 26 genes with orthologs in both human and mouse that did not contain poly-Q repeats. The bovine poly-Q encoding genes typically had ubiquitous expression patterns although there was bias towards expression in epithelia, brain and testes. They were also characterised by unusually large sizes. Analysis of gene ontology terms revealed that the encoded proteins were strongly enriched for functions associated with transcriptional regulation and many contributed to physical interaction networks in the nucleus where they presumably act cooperatively in transcriptional regulatory complexes. In addition, the coding sequence CAG repeats in some bovine genes impacted mRNA splicing thereby generating unusual transcriptional diversity, which in at least one instance was tissue-specific. The poly-Q encoding genes were prioritised using multiple criteria for their likelihood of being polymorphic and then the highest ranking group was experimentally tested for polymorphic variation within a cattle diversity panel. Extensive and meiotically stable variation was identified.

Conclusions

Transcriptional diversity can potentially be generated in poly-Q encoding genes by the impact of CAG repeat tracts on mRNA alternative splicing. This effect, combined with the physical interactions of the encoded proteins in large transcriptional regulatory complexes suggests that polymorphic variations of proteins in these complexes have strong potential to affect phenotype.

Extremophiles and their application to veterinary medicine

Extremophiles can be defined as organisms that can survive in extreme environments that cannot support mammalian life. They include microorganisms that can tolerate temperature extremes, extremes of pH, salinity, hydrostatic pressure and ionizing radiation, as well as low oxygen tension, desiccation and the presence of heavy metals. Psychrophilic organisms also include fish in polar waters and animals that withstand freezing. Rare examples of thermophilic pathogens exist, and the main category of extremophilic animal pathogens comprises psychrophilic and psychrotrophic microorganisms that cause fish diseases, e.g. Flavobacterium psychrophilum, Moritella viscosa, Aliivibrio wodanis and Aliivibrio salmonicida. The most widely known application of an extremophile product in veterinary medicine is DNA polymerase from thermophiles, which is a mainstay of PCR-based diagnostics for an extensive range of animal pathogens. DNA polymerases and other extremophile enzymes are also used in many molecular biology applications and animal genomics. Other extremophile products may find application in veterinary medicine in the future. These include enzymes in biosensors, compatible solutes in skin care products, drug excipients, treatments for respiratory disease, radioprotectants, peptide antibiotics, archaeal lipids for drug delivery and anti-cancer therapeutics.

Copy number variation in the bovine genome

Background

Copy number variations (CNVs), which represent a significant source of genetic diversity in mammals, have been shown to be associated with phenotypes of clinical relevance and to be causative of disease. Notwithstanding, little is known about the extent to which CNV contributes to genetic variation in cattle.

Results

We designed and used a set of NimbleGen CGH arrays that tile across the assayable portion of the cattle genome with approximately 6.3 million probes, at a median probe spacing of 301 bp. This study reports the highest resolution map of copy number variation in the cattle genome, with 304 CNV regions (CNVRs) being identified among the genomes of 20 bovine samples from 4 dairy and beef breeds. The CNVRs identified covered 0.68% (22 Mb) of the genome, and ranged in size from 1.7 to 2,031 kb (median size 16.7 kb). About 20% of the CNVs co-localized with segmental duplications, while 30% encompass genes, of which the majority is involved in environmental response. About 10% of the human orthologous of these genes are associated with human disease susceptibility and, hence, may have important phenotypic consequences.

Conclusions

Together, this analysis provides a useful resource for assessment of the impact of CNVs regarding variation in bovine health and production traits.

Epigenetic regulation of milk production in dairy cows

It is well established that milk production of the dairy cow is a function of mammary epithelial cell (MEC) number and activity and that these factors can be influenced by diverse environmental influences and management practises (nutrition, milk frequency, photoperiod, udder health, hormonal and local effectors). Thus, understanding how the mammary gland is able to respond to these environmental cues provides a huge potential to enhance milk production of the dairy cow. In recent years our understanding of molecular events within the MEC underlying bovine lactation has been advanced through mammary microarray studies and will be further advanced through the recent availability of the bovine genome sequence. In addition, the potential of epigenetic regulation (non-sequence inheritable chemical changes in chromatin, such as DNA methylation and histone modifications, which affect gene expression) to manipulate mammary function is emerging. We propose that a substantial proportion of unexplained phenotypic variation in the dairy cow is due to epigenetic regulation. Heritability of epigenetic marks also highlights the potential to modify lactation performance of offspring. Understanding the response of the MEC (cell signaling pathways and epigenetic mechanisms) to external stimuli will be an important prerequisite to devising new technologies for maximising their activity and, hence, milk production in the dairy cow.

Annotation of novel transcripts putatively relevant for bovine fat metabolism

Two bovine transcripts encoded by the interleukin-1 receptor-associated kinase 1 (IRAK1) gene and the locus LOC618944 predicted as similar to human chromosome 6 open reading frame 52 (C6orf52) gene had indicated divergent expression in bovine skeletal muscle containing different amount of intramuscular fat in a pilot screening experiment. However, for both loci any role in the regulation of energy or fat metabolism is not yet described. In this study, we validated and refined gene structure, screened for mRNA splice variants and analyzed the tissue-specific gene expression patterns of both loci as a prerequisite to elucidate their potential physiological function. Based on comparative sequence analysis, a new full-length gene model for the bovine IRAK1 gene was developed and confirmed experimentally. Expression of IRAK1 mRNA was found in a variety of tissues, and a splice variant was identified in skeletal muscle caused by an in-frame deleted segment of 210 bp affecting regions of intrinsic disorder in the respective protein. For the locus LOC618944, our data contributed to a revised gene model and its assignment to BTA23 (bovine chromosome 23) on the current bovine genome assembly supported by comparative similarity analysis between the bovine and human genomes and experimental data. Furthermore, we identified several splice variants in mammary gland, fat and skeletal muscle tissue and detected a highly similar processed pseudogene on BTA26. All transcript variants of LOC618944 detected in the analyzed tissues represent noncoding RNAs. For both loci, our results suggest yet undetected physiological functions in tissues relevant for fat or energy metabolism in cattle.

Analysis of structure and gene expression of bovine CCDC3 gene indicates a function in fat metabolism

Our study reports the molecular analysis of the bovine gene encoding the coiled-coil domain-containing protein 3 (CCDC3). Based on comparative sequence analysis and in silico sequence merging of predicted gene models, a new full-length gene model for the bovine CCDC3 gene was predicted and confirmed experimentally. The CCDC3 gene was assigned to bovine chromosome 13. It consists of three exons comprising 2599bp encoding for a respective protein of 274 amino acids. The strong CCDC3 sequence homology on amino acid level between species suggests a conserved universal function of this protein. In mice, the CCDC3 gene had been found to be highly expressed in adipocytes and regulated by hormonal-nutritional alternations and in obesity. The tissue expression pattern of bovine CCDC3 mRNA indicates a ubiquitous physiological function of the gene. Significant differences in CCDC3 mRNA expression in skeletal muscle between individuals characterized by divergent intramuscular fat deposition support the potential function of the gene in fat or energy metabolism, which possibly could also be inferred for other mammalian species. This first report of structural analysis and molecular characterization of the CCDC3 gene in cattle will contribute to a better understanding of the yet unknown physiological role of the respective protein in mammals.