High-throughput genotyping methods have increased the analytical power to study complex traits but high cost has remained a barrier for large scale use in animal improvement. We have adapted ...genotyping-by-sequencing (GBS) used in plants for genotyping 47 animals representing 7 taurine and indicine breeds of cattle from the US and Africa. Genomic DNA was digested with different enzymes, ligated to adapters containing one of 48 unique bar codes and sequenced by the Illumina HiSeq 2000. PstI was the best enzyme producing 1.4 million unique reads per animal and initially identifying a total of 63,697 SNPs. After removal of SNPs with call rates of less than 70%, 51,414 SNPs were detected throughout all autosomes with an average distance of 48.1 kb, and 1,143 SNPs on the X chromosome at an average distance of 130.3 kb, as well as 191 on unmapped contigs. If we consider only the SNPs with call rates of 90% and over, we identified 39,751 on autosomes, 850 on the X chromosome and 124 on unmapped contigs. Of these SNPs, 28,843 were not tightly linked to other SNPs. Average marker density per autosome was highly correlated with chromosome size (coefficient of correlation = -0.798, r(2) = 0.637) with higher density in smaller chromosomes. Average SNP call rate was 86.5% for all loci, with 53.0% of the loci having call rates >90% and the average minor allele frequency being 0.212. Average observed heterozygosity ranged from 0.046-0.294 among individuals, and from 0.064-0.197 among breeds, with Brangus showing the highest diversity as expected. GBS technique is novel, flexible, sufficiently high-throughput, and capable of providing acceptable marker density for genomic selection or genome-wide association studies at roughly one third of the cost of currently available genotyping technologies.
High-throughput sequencing technologies have increased the ability to detect sequence variations for complex trait improvement. A high throughput genome wide genotyping-by-sequencing (GBS) method was ...used to generate 515,787 single nucleotide polymorphisms (SNPs), from which 76,355 SNPs with call rates >85% and minor allele frequency ≥1.5% were used in genome wide association study (GWAS) of 44 milk traits in 1,246 Canadian Holstein cows. GWAS was accomplished with a mixed linear model procedure implementing the additive and dominant models. A strong signal within the centromeric region of bovine chromosome 14 was associated with test day fat percentage. Several SNPs were associated with eicosapentaenoic acid, docosapentaenoic acid, arachidonic acid, CLA:9c11t and gamma linolenic acid. Most of the significant SNPs for 44 traits studied are novel and located in intergenic regions or introns of genes. Novel potential candidate genes for milk traits or mammary gland functions include ERCC6, TONSL, NPAS2, ACER3, ITGB4, GGT6, ACOX3, MECR, ADAM12, ACHE, LRRC14, FUK, NPRL3, EVL, SLCO3A1, PSMA4, FTO, ADCK5, PP1R16A and TEP1. Our study further demonstrates the utility of the GBS approach for identifying population-specific SNPs for use in improvement of complex dairy traits.
We examined the effects of two direct-fed microbials (DFM) containing multiple microbial species and their fermentation products on energy status, nutrient digestibility, and ruminal fermentation, ...bacterial community, and metabolome of beef steers. Nine ruminally cannulated Holstein steers (mean ± SD body weight: 243 ± 12.4 kg) were assigned to three treatments arranged in a triplicated 3 × 3 Latin square design with three 21-d periods. Dietary treatments were 1) control (CON; basal diet), 2) Commence (PROB; basal diet plus 19 g/d of Commence), and 3) RX3 (SYNB; basal diet plus 28 g/d of RX3). Commence and RX3 are both multispecies DFM products. From day 16 to 20 of each period, feed and fecal samples were collected daily to determine the apparent total tract digestibilities of nutrients using indigestible neutral detergent fiber method. On day 21 of each period, blood samples were collected for analysis of plasma glucose and nonesterified fatty acid. Ruminal contents were collected at approximately 1, 3, 6, 9, 12, and 18 h after feeding on day 21 for analysis of volatile fatty acids (VFA), lactate, ammonia-N concentrations, bacterial community, and metabolome profile. Total tract digestibilities of nutrients did not differ (P > 0.05) among treatments. Compared with CON, steers fed either supplemental PROB or SYNB had greater (P = 0.04) plasma glucose concentrations. Compared with CON, total ruminal VFA, propionate, isovalerate, and valerate concentrations increased (P ≤ 0.05) or tended to increase (P ≤ 0.10) with either supplemental PROB or SYNB, but were not different (P > 0.05) between PROB and SYNB. Compared with CON, PROB reduced (P ≤ 0.05) the relative abundance of Prevotella 1 and Prevotellaceae UCG-001 but increased (P ≤ 0.05) the relative abundance of Rikenellaceae RC9, Succinivibrionaceae UCG-001, Succiniclasticum, and Ruminococcaceae UCG-002. Supplemental SYNB decreased (P ≤ 0.05) the relative abundance of Prevotella 1 and Prevotellaceae UCG-001 but increased (P ≤ 0.05) the relative abundance of Prevotella 7, Succinivibrio, Succiniclasticum, and Ruminococcaceae UCG-014. Compared with CON, metabolome analysis revealed that some amino acids were increased (P ≤ 0.05) in steers fed PROB. This study demonstrated that, compared with CON, supplementation of either PROB or SYNB altered the ruminal bacterial community and metabolome differently; however, their effects on the ruminal VFA profile and energy status of the steers were not different from each other.
Genomic selection methodologies and genome-wide association studies use powerful statistical procedures that correlate large amounts of high-density SNP genotypes and phenotypic data. Actual 305-d ...milk (MY), fat (FY), and protein (PY) yield data on 695 cows and 76,355 genotyping-by-sequencing-generated SNP marker genotypes from Canadian Holstein dairy cows were used to characterize linkage disequilibrium (LD) structure of Canadian Holstein cows. Also, the comparison of pedigree-based BLUP, genomic BLUP (GBLUP), and Bayesian (BayesB) statistical methods in the genomic selection methodologies and the comparison of Bayesian ridge regression and BayesB statistical methods in the genome-wide association studies were carried out for MY, FY, and PY. Results from LD analysis revealed that as marker distance decreases, LD increases through chromosomes. However, unexpected high peaks in LD were observed between marker pairs with larger marker distances on all chromosomes. The GBLUP and BayesB models resulted in similar heritability estimates through 10-fold cross-validation for MY and PY; however, the GBLUP model resulted in higher heritability estimates than BayesB model for FY. The predictive ability of GBLUP model was significantly lower than that of BayesB for MY, FY, and PY. Association analyses indicated that 28 high-effect markers and markers on Bos taurus autosome 14 located within 6 genes (DOP1B, TONSL, CPSF1, ADCK5, PARP10, and GRINA) associated significantly with FY.
We examined the effects of two direct-fed microbials (DFM) containing multiple microbial species and their fermentation products on energy status, nutrient digestibility, and ruminal fermentation, ...bacterial community, and metabolome of beef steers. Nine ruminally cannulated Holstein steers (mean ± SD body weight: 243 ± 12.4 kg) were assigned to three treatments arranged in a triplicated 3 x 3 Latin square design with three 21-d periods. Dietary treatments were 1) control (CON; basal diet), 2) Commence (PROB; basal diet plus 19 g/d of Commence), and 3) RX3 (SYNB; basal diet plus 28 g/d of RX3). Commence and RX3 are both multispecies DFM products. From day 16 to 20 of each period, feed and fecal samples were collected daily to determine the apparent total tract digestibilities of nutrients using indigestible neutral detergent fiber method. On day 21 of each period, blood samples were collected for analysis of plasma glucose and nonesterified fatty acid. Ruminal contents were collected at approximately 1, 3, 6, 9, 12, and 18 h after feeding on day 21 for analysis of volatile fatty acids (VFA), lactate, ammonia-N concentrations, bacterial community, and metabolome profile. Total tract digestibilities of nutrients did not differ (P > 0.05) among treatments. Compared with CON, steers fed either supplemental PROB or SYNB had greater (P = 0.04) plasma glucose concentrations. Compared with CON, total ruminal VFA, propionate, isovalerate, and valerate concentrations increased (P < 0.05) or tended to increase (P < 0.10) with either supplemental PROB or SYNB, but were not different (P > 0.05) between PROB and SYNB. Compared with CON, PROB reduced (P < 0.05) the relative abundance of Prevotella 1 and Prevotellaceae UCG-001 but increased (P < 0.05) the relative abundance of Rikenellaceae RC9, Succinivibrionaceae UCG-001, Succiniclasticum, and Ruminococcaceae UCG-002. Supplemental SYNB decreased (P < 0.05) the relative abundance of Prevotella 1 and Prevotellaceae UCG-001 but increased (P < 0.05) the relative abundance of Prevotella 7, Succinivibrio, Succiniclasticum, and Ruminococcaceae UCG-014. Compared with CON, metabolome analysis revealed that some amino acids were increased (P < 0.05) in steers fed PROB. This study demonstrated that, compared with CON, supplementation of either PROB or SYNB altered the ruminal bacterial community and metabolome differently; however, their effects on the ruminal VFA profile and energy status of the steers were not different from each other.
The present study was carried out to estimate the genetic parameters for direct and maternal influences on Mecheri sheep (
Ovis aries
) growth traits using Bayesian multi-trait animal model. The ...genetic parameters were calculated using data from 2825 Mecheri lambs born between 2010 and 2020 that were kept in semi-arid tropical climate. Mecheri sheep body weight (mean ± SE) at various stages, viz. BW, WW, and BW12, were 2.6 ± 0.01, 11.1 ± 0.05, and 20.7 ± 0.13 kg, respectively. The Mecheri sheep gained 71.5 percent of their body weight at the age of 6 months. With the exception of birth weight, the weights of the animals varied considerably (
P
< 0.01) by the year of birth. The fixed effect of sex significantly (
P
< 0.01) influenced all the growth traits examined. The direct estimates of heritability (± SD) for BW, WW, BW6, and BW12 was 0.21 ± 0.041, 0.21 ± 0.041, 0.12 ± 0.052, and 0.13 ± 0.053, respectively, and the maternal heritability for BW, WW, BW6, and BW12 was 0.18 ± 0.021, 0.08 ± 0.023, 0.11 ± 0.022, and 0.13 ± 0.033, respectively. Significant variance was indicated by moderately larger direct heritability estimates for BWT and WWT, indicating that there will be more opportunities for selection response during the genetic improvement programme. For the majority of the variables examined, direct heritability values were higher than maternal heritability values. The additive genetic correlation between WW and BW6, BW9, and BW12 was 0.70 ± 0.145, 0.57 ± 0.171, and 0.50 ± 0.194, respectively. The maternal genetic correlations ranged from 0.06 ± 0.152 (BW–BW12) to 0.86 ± 0.046 (BW6–BW9), and the residual correlation varied from 0.18 ± 0.034 (BW–WW) to 0.85 ± 0.013 (BW9–BW12). The BW had a stronger genetic relationship with WW and a weaker relationship with other growth traits. The WW and BW6 showed a favourable genetic relationship, even if the tendency was decreasing in the latter stages of maturation. The BW6 and BW9 demonstrated the strongest positive genetic relationship (0.90 ± 0.052) of all the variables analysed. After examining the heritabilities and genetic correlation between WW and BW6, it was established that selection based on at WW was beneficial, since it manifests early in life and would result in moderate genetic progress through selection.
We investigated whole blood and hepatic mRNA expressions of immune genes and rumen microbiome of crossbred beef steers with divergent residual feed intake phenotype to identify relevant biological ...processes underpinning feed efficiency in beef cattle. Low-RFI beef steers (n = 20; RFI = - 1.83 kg/d) and high-RFI beef steers (n = 20; RFI = + 2.12 kg/d) were identified from a group of 108 growing crossbred beef steers (average BW = 282 ± 30.4 kg) fed a high-forage total mixed ration after a 70-d performance testing period. At the end of the 70-d testing period, liver biopsies and blood samples were collected for total RNA extraction and cDNA synthesis. Rumen fluid samples were also collected for analysis of the rumen microbial community. The mRNA expression of 84 genes related to innate and adaptive immunity was analyzed using pathway-focused PCR-based arrays. Differentially expressed genes were determined using P-value ≤ 0.05 and fold change (FC) ≥ 1.5 (in whole blood) or ≥ 2.0 (in the liver). Gene ontology analysis of the differentially expressed genes revealed that pathways related to pattern recognition receptor activity, positive regulation of phagocytosis, positive regulation of vitamin metabolic process, vascular endothelial growth factor production, positive regulation of epithelial tube formation and T-helper cell differentiation were significantly enriched (FDR < 0.05) in low-RFI steers. In the rumen, the relative abundance of PeH15, Arthrobacter, Moryella, Weissella, and Muribaculaceae was enriched in low-RFI steers, while Methanobrevibacter, Bacteroidales_BS11_gut_group, Bacteroides and Clostridium_sensu_stricto_1 were reduced. In conclusion, our study found that low-RFI beef steers exhibit increased mRNA expression of genes related to immune cell functions in whole blood and liver tissues, specifically those involved in pathogen recognition and phagocytosis regulation. Additionally, these low-RFI steers showed differences in the relative abundance of some microbial taxa which may partially account for their improved feed efficiency compared to high-RFI steers.
The African continent is home to diverse populations of livestock breeds adapted to harsh environmental conditions with more than 70% under traditional systems of management. Animal productivity is ...less than optimal in most cases and is faced with numerous challenges including limited access to adequate nutrition and disease management, poor institutional capacities and lack of adequate government policies and funding to develop the livestock sector. Africa is home to about 1.3 billion people and with increasing demand for animal proteins by an ever growing human population, the current state of livestock productivity creates a significant yield gap for animal products. Although a greater section of the population, especially those living in rural areas depend largely on livestock for their livelihoods; the potential of the sector remains underutilized and therefore unable to contribute significantly to economic development and social wellbeing of the people. With current advances in livestock management practices, breeding technologies and health management, and with inclusion of all stakeholders, African livestock populations can be sustainably developed to close the animal protein gap that exists in the continent. In particular, advances in gene technologies, and application of genomic breeding in many Western countries has resulted in tremendous gains in traits like milk production with the potential that, implementation of genomic selection and other improved practices (nutrition, healthcare, etc.) can lead to rapid improvement in traits of economic importance in African livestock populations. The African livestock populations in the context of this review are limited to cattle, goat, pig, poultry, and sheep, which are mainly exploited for meat, milk, and eggs. This review examines the current state of livestock productivity in Africa, the main challenges faced by the sector, the role of various stakeholders and discusses in-depth strategies that can enable the application of genomic technologies for rapid improvement of livestock traits of economic importance.
The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The ...interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats' populations. Fu and Li tests were significantly positive but Tajima's D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model dN/dS ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω0 (dN/dS) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.