Genome-wide association study of tick resistance in South African Nguni cattle

Genetic and genomic analysis of Belgian Blue's susceptibility for psoroptic mange

BACKGROUND

Psoroptic mange, caused by Psoroptes ovis mites, is affecting Belgian Blue cattle's welfare and production potential. The Belgian Blue cattle-known for its high degree of muscling, low feed conversion ratio and high beef quality-is highly susceptible for this disease.

RESULTS

In this study, we phenotyped 1975 Belgian Blue cattle from more than 100 different groups on commercial beef farms for their psoroptic mange susceptibility. Substantial individual differences were observed within these management groups, with lesion extent differences up to ± 15%. Animal models showed that estimated heritabilities were low for lesion extent and severe lesion extent (0.07 and 0.09, respectively) and 0.12 for the number of mites. A genome wide association study for mange susceptibility revealed signals on BTA6, BTA11, BTA15 and BTA24. In these regions, candidate genes GBA3, RAG2, and TRAF6 were identified.

CONCLUSIONS

Despite the challenges in phenotyping for psoroptic mange due to the timing of screening, the continuous evolution of lesions and different management conditions, we successfully conducted a study on the genetic susceptibility to psoroptic mange in Belgian Blue cattle. Our results clearly indicate that psoroptic mange is under polygenic control and the underlying candidate genes should be studied more thoroughly. This is the first study providing candidate genes for this complex disease. These results are already valuable for Belgian Blue breeding, however, further research is needed to unravel the architecture of this disease and to identify causal mutations.

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle

Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.

A Genome-Wide Association Study for Resistance to Tropical Theileriosis in Two Bovine Portuguese Autochthonous Breeds

The control of Tropical Theileriosis, a tick-borne disease with a strong impact on cattle breeding, can be facilitated using marker-assisted selection in breeding programs. Genome-wide association studies (GWAS) using high-density arrays are extremely important for the ongoing process of identifying genomic variants associated with resistance to Theileria annulata infection. In this work, single-nucleotide polymorphisms (SNPs) were analyzed in the Portuguese autochthonous cattle breeds Alentejana and Mertolenga. In total, 24 SNPs suggestive of significance (p ≤ 10-4) were identified for Alentejana cattle and 20 SNPs were identified for Mertolenga cattle. The genomic regions around these SNPs were further investigated for annotated genes and quantitative trait loci (QTLs) previously described by other authors. Regarding the Alentejana breed, the MAP3K1, CMTM7, SSFA2, and ATG13 genes are located near suggestive SNPs and appear as candidate genes for resistance to Tropical Theileriosis, considering its action in the immune response and resistance to other diseases. On the other hand, in the Mertolenga breed, the UOX gene is also a candidate gene due to its apparent link to the pathogenesis of the disease. These results may represent a first step toward the possibility of including genetic markers for resistance to Tropical Theileriosis in current breed selection programs.

Estimation of genetic parameters for the tick and hemoparasite burden in Angus cattle

The study was conducted with the objective of estimating genetic and phenotypic parameters for tick (CRM) and Babesia bigemina (IBBi), Babesia bovis (IBBo), and Anaplasma marginale (IAM) burden in Angus female breed in Brazil. The sample group was composed of Angus females raised in herds located in a region of endemic instability for cattle tick fever in the state of Rio Grande Sul (RS), Brazil. The variance components were estimated using Bayesian inference and Gibbs sampling algorithm, considering a multi-trait animal model. Heritability estimates showed values of low magnitude, ranging from 0.03 (IBBo) to 0.16 (CRM), while repeatability estimates ranged between 0.07 (IBBo) and 0.21 (CRM). Regarding the genetic correlation estimates, the values showed low (-0.01 for IBBo × IAM) to moderate (0.55 between IBBi × IAM) magnitudes. The results indicate that it is possible to use tick count and hemoparasite infection levels as selection criteria, with small genetic gains.

Selective versus strategic control against Rhipicephalus microplus in cattle: A comparative analysis of efficacy, animal health, productivity, cost, and resistance management

This study compared selective control versus strategic control against Rhipicephalus (Boophilus) microplus, under conditions of a real dairy farm. Strategic control requires that all animals be treated with acaricide at regular pre-determined intervals. Selective control only requires treatment of infested animals and only when they are at or above a pre-determined threshold. Tick counts on animals and in pasture were performed and the susceptibility of tick populations to the different treatment methods was evaluated at the beginning and end of the study using the Larval Packet Test, Larval Immersion Test, and Adult Immersion Test. Over the four years of the study strategic control was more advantageous than the selective control as the group experienced lower tick burden on animals and in pasture, absence of skin lesions and myiases, lower operating costs, treatments concentrated in fewer months of the year and lower resistance pressure.

Investigating volatile semiochemical production from Bos taurus and Bos indicus as a novel phenotype for breeding host resistance to ixodid ticks

Ticks and tick-borne diseases cause significant loss in livestock production with about 80% world's cattle at risk. The cost of chemical control is high and there is an ever-increasing tick resistance to chemical acaricides. Genetic selection as alternative long-term control strategy is constrained by laborious phenotyping using tick counts or scores. This study explored the use of host volatile semiochemicals that may be attractants or repellents to ticks as a phenotype for new tick resistance, with potential to be used as a proxy in selection programmes. Approximately 100 young cattle composed of Bos indicus and Bos taurus were artificially infested with 2,500 African blue tick, Rhipicephalus decoloratus larvae, with daily female tick (4.5 mm) counts taken from day 20 post-infestation. Volatile organic compounds were sampled from cattle before and after tick infestation by dynamic headspace collection, analysed by high-resolution gas chromatography (GC) and subjected to multivariate statistical analysis. Using 6-day repeated measure analysis, three pre-infestation GC peaks (BI938 - unknown, BI966 - 6-methyl-5-hepten-2-one and BI995 - hexyl acetate) and one post-infestation GC peak (AI933 - benzaldehyde / (E)-2-heptenal) were associated with tick resistance (P < 0.01 and P < 0.05 respectively). The high correlation coefficients (r = 0.66) between repeated records with all volatile compounds support the potential predictive value for volatile compounds in selective breeding programmes for tick resistance in cattle.

Tick Diversity and Distribution of Hard (Ixodidae) Cattle Ticks in South Africa

Ticks are amongst the important ectoparasites where livestock are concerned, as they adversely affect the animals through bloodsucking. In tropical and subtropical countries, they transmit pathogens such as babesiosis, theileriosis, ehrlichiosis, and anaplasmosis in cattle, causing a reduction in production rate and significant concomitant economic losses. Ticks affect 80% of the cattle population across the world, with an estimated economic loss of USD 20–30 billion per year. In South Africa, economic losses in the livestock industry caused by ticks and tick-borne diseases are estimated to exceed USD 33 million per year (ZAR 500 million). There are seven major genera of ixodid ticks in Southern Africa (i.e., Amblyomma, Dermacentor, Haemaphysalis, Hyalomma, Ixodes and Rhipicephalus). The environment in which a tick lives is made up of all the various biological and abiotic factors that are either necessary or unnecessary for its life. The areas where various ticks have been found have been documented in many publications. Using these data, maps of possible species’ habitats can be made. Historical records on tick distribution may be incorrect due to identification mistakes or a change in the tick’s name. All the sources used to generate the maps for this review were unpublished and came from a wide range of sources. To identify tick species and the accompanying microbial ecosystems, researchers are increasingly adopting tick identification methods including 16S and 18S rDNA gene sequencing. Indeed, little is known about the genetic alterations that give important traits, including the predilection for tick hosts, transmission, and acaricide resistance. Opportunities for exploring these changes in tick populations and subpopulations are provided by advancements in omics technologies. The literature on the variety of ixodid ticks, their direct and indirect effects, and control methods in South Africa is compiled in this review.

Genome-wide analysis identified candidate variants and genes associated with heat stress adaptation in Egyptian sheep breeds

Heat stress caused by climatic changes is one of the most significant stresses on livestock in hot and dry areas. It has particularly adverse effects on the ability of the breed to maintain homeothermy. Developing countries are advised to protect and prepare their animal resources in the face of potential threats such as climate change. The current study was conducted in Egypt’s three hot and dry agro-ecological zones. Three local sheep breeds (Saidi, Wahati, and Barki) were studied with a total of 206 ewes. The animals were exercised under natural heat stress. The heat tolerance index of the animals was calculated to identify animals with high and low heat tolerance based on their response to meteorological and physiological parameters. Genomic variation in these breeds was assessed using 64,756 single nucleotide polymorphic markers (SNPs). From the perspective of comparative adaptability to harsh conditions, our objective was to investigate the genomic structure that might control the adaptability of local sheep breeds to environmental stress under hot and dry conditions. In addition, indices of population structure and diversity of local breeds were examined. Measures of genetic diversity showed a significant influence of breed and location on populations. The standardized index of association (rbarD) ranged from 0.0012 (Dakhla) to 0.026 (Assuit), while for the breed, they ranged from 0.004 (Wahati) to 0.0103 (Saidi). The index of association analysis (I^a) ranged from 1.42 (Dakhla) to 35.88 (Assuit) by location and from 6.58 (Wahati) to 15.36 (Saidi) by breed. The most significant SNPs associated with heat tolerance were found in the MYO5A, PRKG1, GSTCD, and RTN1 genes (p ≤ 0.0001). MYO5A produces a protein widely distributed in the melanin-producing neural crest of the skin. Genetic association between genetic and phenotypic variations showed that OAR1_18300122.1, located in ST3GAL3, had the greatest positive effect on heat tolerance. Genome-wide association analysis identified SNPs associated with heat tolerance in the PLCB1, STEAP3, KSR2, UNC13C, PEBP4, and GPAT2 genes.

Identification and Characterization of Copy Number Variations Regions in West African Taurine Cattle

A total of 106 West African taurine cattle belonging to the Lagunaire breed of Benin (33), the N’Dama population of Burkina Faso (48), and N’Dama cattle sampled in Congo (25) were analyzed for Copy Number Variations (CNVs) using the BovineHDBeadChip of Illumina and two different CNV calling programs: PennCNV and QuantiSNP. Furthermore, 89 West African zebu samples (Bororo cattle of Mali and Zebu Peul sampled in Benin and Burkina Faso) were used as an outgroup to ensure that analyses reflect the taurine cattle genomic background. Analyses identified 307 taurine-specific CNV regions (CNVRs), covering about 56 Mb on all bovine autosomes. Gene annotation enrichment analysis identified a total of 840 candidate genes on 168 taurine-specific CNVRs. Three different statistically significant functional term annotation clusters (from ACt1 to ACt3) involved in the immune function were identified: ACt1 includes genes encoding lipocalins, proteins involved in the modulation of immune response and allergy; ACt2 includes genes encoding coding B-box-type zinc finger proteins and butyrophilins, involved in innate immune processes; and Act3 includes genes encoding lectin receptors, involved in the inflammatory responses to pathogens and B- and T-cell differentiation. The overlap between taurine-specific CNVRs and QTL regions associated with trypanotolerant response and tick-resistance was relatively low, suggesting that the mechanisms underlying such traits may not be determined by CNV alterations. However, four taurine-specific CNVRs overlapped with QTL regions associated with both traits on BTA23, therefore suggesting that CNV alterations in major histocompatibility complex (MHC) genes can partially explain the existence of genetic mechanisms shared between trypanotolerance and tick resistance in cattle. This research contributes to the understanding of the genomic features of West African taurine cattle.

Prevalence and tick loads in Nguni cattle reared in different environmental conditions across four provinces of South Africa

Background and Aim:

In tropical and subtropical countries, ixodid ticks are among livestock’s most economically important ectoparasites. Although Nguni cattle from South Africa have adapted to harsh environments, it is unknown whether they will be resistant to ticks, and the diseases carried by ticks under various climatic conditions. Therefore, this study aimed to compare tick load and estimate the prevalence of different tick species among Nguni cattle under different environmental conditions.

Materials and Methods:

Tick counts were conducted monthly under natural challenges over 2 years on 586 Nguni cattle located at ARC-Roodeplaat and Loskop farms (warmer climate), Mukhuthali Nguni Community and the University of Fort Hare farms (cooler climate). The generalized linear model procedure of the Statistical Analysis System was used to analyze the data. It fitted the location (farm), sex, year, month or season, and animal age as covariates.

Results:

The tick species (relative prevalence) observed were as follows: Amblyomma hebraeum (42%), Rhipicephalus evertsi (22%), Rhipicephalus (Boophilus) spp. (16%), Rhipicephalus appendiculatus (11%), Hyalomma marginatum (5%), and Rhipicephalus simus (4%). Tick infestation was significantly affected by location, season, year, month of the tick counting and age of the animal. Loskop farm had the highest tick count (m = 30.69) and showed the largest variation in tick count. Compared to the other seasons, higher tick counts were seen during the hot-dry (September–November) and hot-wet (December–February) seasons. A. hebraeum was the dominant tick species across all four farms, followed by R. evertsi. The perianal region (under the tail head), the perineum and the belly body locations were the most preferred tick attachment sites.

Conclusion:

These results provide useful information for developing appropriate control strategies for ticks and tick-borne diseases in these provinces of South Africa. Further work must investigate the feasibility of genetic improvement for tick resistance.

Sustained Effects of Muscle Calpain System Genotypes on Tenderness Phenotypes of South African Beef Bulls during Ageing up to 20 Days

Simple Summary

When searching for genetic markers for the selection of more tender beef, it is important to maintain minimal environmental variation from pre-slaughter, right through to the ageing process, to ensure the accuracy of the obtained phenotypes. This is because beef quality traits have a large environmental component that can greatly alter the characteristics of the meat, which would not reflect a true genetic effect. We propose that variable ageing times are especially important in determining whether markers are associated with tenderization or not. Our analyses included candidate genes for the protein degrading enzyme system for calpains, because they contribute the most to tenderization. We were able to validate these markers in South African beef cattle, where they could be useful for selection. The timing of the collection of tenderness data was critical, as only a few (6/134) genetic markers sustained their association with tenderization over ageing to 20 days. A larger tenderization effect earlier in ageing, as shown here for the capn1_187 and capn1_4751 markers, would decrease the length of ageing. This would not only increase profits, but also decrease the energy needed during the storage and refrigeration of aged beef, decreasing the carbon footprint of beef production.

Abstract

The most important factor that determines beef tenderness is its proteolytic activity, and the balance between calpain-1 protease activity and calpastatin inhibition is especially important, while contributions can also arise from calpain-2 and, possibly, calpain-3. The meat ageing process itself affects these processes. To determine whether genotypes in the calpain–calpastatin system can enhance tenderness through a 20-day ageing period, South African purebred beef bulls (n = 166) were genotyped using the Illumina BovineHD SNP BeadChip through a gene-based association analysis targeting the cast, capn3, capn2 and capn1 genes. The Warner–Bratzler shear force (WBSF) and myofibril fragment length (MFL) of Longissimus thoracis et lumborum (LTL) steaks were evaluated between d 3 and d 20 of ageing, with protease enzyme activity in the first 20 h post-mortem. Although several of the 134 SNPs are associated with tenderness, only seven SNP in the cast, capn2 and capn1 genes sustained genetic associations, additive to the ageing-associated increases in tenderness for at least three of the four ageing periods. While most genomic associations were relatively stable over time, some genotypes within the SNP responded differently to ageing, resulting in altered genomic effects over time. The level of ageing at which genomic associations are performed is an important factor that determines whether SNPs affect tenderness phenotypes.

Myostatin mutation causing double muscling could affect increased psoroptic mange sensitivity in dual purpose Belgian Blue cattle

Belgian Blue cattle are known for their high degree of muscling and good carcass qualities. This high degree of muscling is mainly caused by a mutation in the myostatin gene (MSTN). Although the MSTN mutation is considered as fixed in the Belgian Blue breed, segregation is occurring in a sub-population bred for dual purpose. In the latter population, we observed an association between the mutation in MSTN and susceptibility to psoroptic mange, a skin disease caused by Psoroptes ovis mites that heavily plagues Belgian Blue cattle. In total, 291 animals were sampled and screened for their susceptibility for mange lesions and their MSTN genotype. Via linear mixed modelling, we observed that homozygous mutant animals had a significant increase in the size of mange lesions (+2.51% lesion extent) compared to homozygous wild type. These findings were confirmed with zero-inflated modelling, an animal model and odds analysis. Risk ratios for developing severe mange lesions were 5.9 times as high for homozygous mutant animals. All analyses confirmed an association between the MSTN genotype and psoroptic mange lesion size.

Suitability of GWAS as a Tool to Discover SNPs Associated with Tick Resistance in Cattle: A Review

Understanding the biological mechanisms underlying tick resistance in cattle holds the potential to facilitate genetic improvement through selective breeding. Genome wide association studies (GWAS) are popular in research on unraveling genetic determinants underlying complex traits such as tick resistance. To date, various studies have been published on single nucleotide polymorphisms (SNPs) associated with tick resistance in cattle. The discovery of SNPs related to tick resistance has led to the mapping of associated candidate genes. Despite the success of these studies, information on genetic determinants associated with tick resistance in cattle is still limited. This warrants the need for more studies to be conducted. In Africa, the cost of genotyping is still relatively expensive; thus, conducting GWAS is a challenge, as the minimum number of animals recommended cannot be genotyped. These population size and genotype cost challenges may be overcome through the establishment of collaborations. Thus, the current review discusses GWAS as a tool to uncover SNPs associated with tick resistance, by focusing on the study design, association analysis, factors influencing the success of GWAS, and the progress on cattle tick resistance studies.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

Genetic Analyses and Genome-Wide Association Studies on Pathogen Resistance of Bos taurus and Bos indicus Cattle Breeds in Cameroon

Autochthonous taurine and later introduced zebu cattle from Cameroon differ considerably in their resistance to endemic pathogens with little to no reports of the underlying genetic make-up. Breed history and habitat variations are reported to contribute significantly to this diversity worldwide, presumably in Cameroon as well, where locations diverge in climate, pasture, and prevalence of infectious agents. In order to investigate the genetic background, the genotypes of 685 individuals of different Cameroonian breeds were analysed by using the BovineSNP50v3 BeadChip. The variance components including heritability were estimated and genome-wide association studies (GWAS) were performed. Phenotypes were obtained by parasitological screening and categorised in Tick-borne pathogens (TBP), gastrointestinal nematodes (GIN), and onchocercosis (ONC). Estimated heritabilities were low for GIN and TBP (0.079 (se = 0.084) and 0.109 (se = 0.103) respectively) and moderate for ONC (0.216 (se = 0.094)). Further than revealing the quantitative nature of the traits, GWAS identified putative trait-associated genomic regions on five chromosomes, including the chromosomes 11 and 18 for GIN, 20 and 24 for TBP, and 12 for ONC. The results imply that breeding for resistant animals in the cattle population from Northern Cameroon might be possible for the studied pathogens; however, further research in this field using larger datasets will be required to improve the resistance towards pathogen infections, propose candidate genes or to infer biological pathways, as well as the genetic structures of African multi-breed populations.

Multiple Country and Breed Genomic Prediction of Tick Resistance in Beef Cattle

Ticks cause substantial production losses for beef and dairy cattle. Cattle resistance to ticks is one of the most important factors affecting tick control, but largely neglected due to the challenge of phenotyping. In this study, we evaluate the pooling of tick resistance phenotyped reference populations from multi-country beef cattle breeds to assess the possibility of improving host resistance through multi-trait genomic selection. Data consisted of tick counts or scores assessing the number of female ticks at least 4.5 mm length and derived from seven populations, with breed, country, number of records and genotyped/phenotyped animals being respectively: Angus (AN), Brazil, 2,263, 921/1,156, Hereford (HH), Brazil, 6,615, 1,910/2,802, Brangus (BN), Brazil, 2,441, 851/851, Braford (BO), Brazil, 9,523, 3,062/4,095, Tropical Composite (TC), Australia, 229, 229/229, Brahman (BR), Australia, 675, 675/675, and Nguni (NG), South Africa, 490, 490/490. All populations were genotyped using medium density Illumina SNP BeadChips and imputed to a common high-density panel of 332,468 markers. The mean linkage disequilibrium (LD) between adjacent SNPs varied from 0.24 to 0.37 across populations and so was sufficient to allow genomic breeding values (GEBV) prediction. Correlations of LD phase between breeds were higher between composites and their founder breeds (0.81 to 0.95) and lower between NG and the other breeds (0.27 and 0.35). There was wide range of estimated heritability (0.05 and 0.42) and genetic correlation (-0.01 and 0.87) for tick resistance across the studied populations, with the largest genetic correlation observed between BN and BO. Predictive ability was improved under the old-young validation for three of the seven populations using a multi-trait approach compared to a single trait within-population prediction, while whole and partial data GEBV correlations increased in all cases, with relative improvements ranging from 3% for BO to 64% for TC. Moreover, the multi-trait analysis was useful to correct typical over-dispersion of the GEBV. Results from this study indicate that a joint genomic evaluation of AN, HH, BN, BO and BR can be readily implemented to improve tick resistance of these populations using selection on GEBV. For NG and TC additional phenotyping will be required to obtain accurate GEBV.

Genomic Loci Affecting Milk Production in German Black Pied Cattle (DSN)

German Black Pied cattle (DSN) is an endangered population of about 2,550 dual-purpose cattle in Germany. Having a milk yield of about 2,500 kg less than the predominant dairy breed Holstein, the preservation of DSN is supported by the German government and the EU. The identification of the genomic loci affecting milk production in DSN can provide a basis for selection decisions for genetic improvement of DSN in order to increase market chances through the improvement of milk yield. A genome-wide association analysis of 30 milk traits was conducted in different lactation periods and numbers. Association using multiple linear regression models in R was performed on 1,490 DSN cattle genotyped with BovineSNP50 SNP-chip. 41 significant and 20 suggestive SNPs affecting milk production traits in DSN were identified, as well as 15 additional SNPs for protein content which are less reliable due to high inflation. The most significant effects on milk yield in DSN were detected on chromosomes 1, 6, and 20. The region on chromosome 6 was located nearby the casein gene cluster and the corresponding haplotype overlapped the CSN3 gene (casein kappa). Associations for fat and protein yield and content were also detected. High correlation between traits of the same lactation period or number led to some SNPs being significant for multiple investigated traits. Half of all identified SNPs have been reported in other studies, previously. 15 SNPs were associated with the same traits in other breeds. The other associated SNPs have been reported previously for traits such as exterior, health, meat and carcass, production, and reproduction traits. No association could be detected between DGAT1 and other known milk genes with milk production traits despite the close relationship between DSN and Holstein. The results of this study confirmed that many SNPs identified in other breeds as associated with milk traits also affect milk traits in dual-purpose DSN cattle and can be used for further genetic analysis to identify genes and causal variants that affect milk production in DSN cattle.

Genetic characterisation of non-descript cattle populations in communal areas of South Africa

Context Indigenous cattle breeds represent an important genetic resource for livelihood of communal-area inhabitants. Indigenous breeds have the ability to withstand harsh climatic conditions, can adapt genetically to poor-quality forages and are resistant to parasites and diseases. These unique traits possessed by indigenous breeds are under threat because of unrestrained crossing with exotic commercial breeds, and this can lead to total loss of a breed. Aims The study was conducted to assess the genetic diversity and population structure of South African non-descript communal beef cattle populations by using 25 microsatellite markers. Methods Unrelated and non-descript animals (n = 150) were sampled from communal areas from five (5) provinces of South Africa, namely, Eastern Cape, KwaZulu–Natal, Limpopo, Mpumalanga and the North West, with 30 samples per breed taken. Six (6) known cattle breeds (n = 180) were used as a reference population. This included Angus, Afrikaner, Bonsmara, Brahman, Drakensberger and the Nguni, with 30 samples per breed. Key results High level of genetic diversity was found across the five non-descript populations, with an average heterozygosity of 75%. The Limpopo population was found to be the most diverse population, with the highest average number of alleles (8.5) and heterozygosity (ranging between observed heterozygosity of 70% and expected heterozygosity of 79%). STRUCTURE software assigned populations (2 ≤ K ≤ 20), with the most probable cluster being at K = 7. The Eastern Cape, KwaZulu–Natal and Limpopo populations had genetic material similar to those possessed by the Nguni and Bonsmara reference populations. Conclusions Results from the study showed that most genetic differentiation occurred within populations rather than among populations, and this might be due to the fact that there is no selection for or against any specific production trait expressed in the populations. Implications The obtained information will serve as a baseline for the development and implementation of sound breeding programs that will assist in controlling the gene flow, so as to lower the possible genetic dilution of the currently available genetic material.

A multi-breed GWAS for morphometric traits in four Beninese indigenous cattle breeds reveals loci associated with conformation, carcass and adaptive traits

BACKGROUND

Specific adaptive features including disease resistance and growth abilities in harsh environments are attributed to indigenous cattle breeds of Benin, but these breeds are endangered due to crossbreeding. So far, there is a lack of systematic trait recording, being the basis for breed characterizations, and for structured breeding program designs aiming on conservation. Bridging this gap, own phenotyping for morphological traits considered measurements for height at withers (HAW), sacrum height (SH), heart girth (HG), hip width (HW), body length (BL) and ear length (EL), including 449 cattle from the four indigenous Benin breeds Lagune, Somba, Borgou and Pabli. In order to utilize recent genomic tools for breed characterizations and genetic evaluations, phenotypes for novel traits were merged with high-density SNP marker data. Multi-breed genetic parameter estimations and genome-wide association studies (GWAS) for the six morphometric traits were carried out. Continuatively, we aimed on inferring genomic regions and functional loci potentially associated with conformation, carcass and adaptive traits.

RESULTS

SNP-based heritability estimates for the morphometric traits ranged between 0.46 ± 0.14 (HG) and 0.74 ± 0.13 (HW). Phenotypic and genetic correlations ranged from 0.25 ± 0.05 (HW-BL) to 0.89 ± 0.01 (HAW-SH), and from 0.14 ± 0.10 (HW-BL) to 0.85 ± 0.02 (HAW-SH), respectively. Three genome-wide and 25 chromosome-wide significant SNP positioned on different chromosomes were detected, located in very close chromosomal distance (±25 kb) to 15 genes (or located within the genes). The genes PIK3R6 and PIK3R1 showed direct functional associations with height and body size. We inferred the potential candidate genes VEPH1, CNTNAP5, GYPC for conformation, growth and carcass traits including body weight and body fat deposition. According to their functional annotations, detected potential candidate genes were associated with stress or immune response (genes PTAFR, PBRM1, ADAMTS12) and with feed efficiency (genes MEGF11 SLC16A4, CCDC117).

CONCLUSIONS

Accurate measurements contributed to large SNP heritabilities for some morphological traits, even for a small mixed-breed sample size. Multi-breed GWAS detected different loci associated with conformation or carcass traits. The identified potential candidate genes for immune response or feed efficiency indicators reflect the evolutionary development and adaptability features of the breeds.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data

BACKGROUND

The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics.

RESULTS

A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds.

CONCLUSION

The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Genomic Study of Babesia bovis Infection Level and Its Association With Tick Count in Hereford and Braford Cattle

Bovine babesiosis is a tick-borne disease caused by intraerythrocytic protozoa and leads to substantial economic losses for the livestock industry throughout the world. Babesia bovis is considered the most pathogenic species, which causes bovine babesiosis in Brazil. Genomic data could be used to evaluate the viability of improving resistance against B. bovis infection level (IB) through genomic selection, and, for that, knowledge of genetic parameters is needed. Furthermore, genome-wide association studies (GWAS) could be conducted to provide a better understanding of the genetic basis of the host response to B. bovis infection. No previous work in quantitative genetics of B. bovis infection was found. Thus, the objective of this study was to estimate the genetic correlation between IB and tick count (TC), evaluate predictive ability and applicability of genomic selection, and perform GWAS in Hereford and Braford cattle. The single-step genomic best linear unbiased prediction method was used, which allows the estimation of both breeding values and marker effects. Standard phenotyping was conducted for both traits. IB quantifications from the blood of 1,858 animals were carried using quantitative PCR assays. For TC, one to three subsequent tick counts were performed by manually counting adult female ticks on one side of each animal's body that was naturally exposed to ticks. Animals were genotyped using the Illumina BovineSNP50 panel. The posterior mean of IB heritability, estimated by the Bayesian animal model in a bivariate analysis, was low (0.10), and the estimations of genetic correlation between IB and TC were also low (0.15). The cross-validation genomic prediction accuracy for IB ranged from 0.18 to 0.35 and from 0.29 to 0.32 using k-means and random clustering, respectively, suggesting that genomic predictions could be used as a tool to improve genetics for IB, especially if a larger training population is developed. The top 10 single nucleotide polymorphisms from the GWAS explained 5.04% of total genetic variance for IB, which were located on chromosomes 1, 2, 5, 6, 12, 17, 18, 16, 24, and 26. Some candidate genes participate in immunity system pathways indicating that those genes are involved in resistance to B. bovis in cattle. Although the genetic correlation between IB and TC was weak, some candidate genes for IB were also reported in tick infestation studies, and they were also involved in biological resistance processes. This study contributes to improving genetic knowledge regarding infection by B. bovis in cattle.

Identification of Whole-Genome Significant Single Nucleotide Polymorphisms in Candidate Genes Associated With Serum Biochemical Traits in Chinese Holstein Cattle

A genome-wide association study (GWAS) was conducted on 23 serum biochemical traits in Chinese Holstein cattle. The experimental population consisted of 399 cattle, each genotyped by a commercial bovine 50K SNP chip, which had 49,663 SNPs. After data cleaning, 41,092 SNPs from 361 Holstein cattle were retained for GWAS. The phenotypes were measured values of serum measurements of these animals that were taken at 11 days after parturition. Two statistical models, a fixed-effect linear regression model (FLM) and a mixed-effect linear model (MLM), were used to estimate the association effects of SNPs. Genome-wide significant and suggestive thresholds were set up to be 1.22E-06 and 2.43E-06, respectively. In the Chinese Holstein population, FLM identified 81 genome-wide significant (0.05/41,092 = 1.22E-06) SNPs associated with 11 serum traits. Among these SNPs, five SNPs (BovineHD0100005950, ARS-BFGL-NGS-115158, BovineHD1500021175, BovineHD0800028900, and BTB-00442438) were also identified by the MLM to have genome-wide suggestive effects on CHE, DBIL, and LDL. Both statistical models pinpointed two SNPs that had significant effects on the Holstein population. The SNP BovineHD0800028900 (located near the gene LOC101903458 on chromosome 8) was identified to be significantly associated with serum high- and low-density lipoprotein (HDL and LDL), whereas BovineHD1500021175 (located in 73.4Mb on chromosome 15) was an SNP significantly associated with total bilirubin and direct bilirubin (TBIL and DBIL). Further analyses are needed to identify the causal mutations affecting serum traits and to investigate the correlation of effects for loci associated with fatty liver disease in dairy cattle.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin

Background

Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance.

Results

Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively.

Conclusion

This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Single-plant GWAS coupled with bulk segregant analysis allows rapid identification and corroboration of plant-height candidate SNPs

BACKGROUND

Genome wide association studies (GWAS) are a powerful tool for identifying quantitative trait loci (QTL) and causal single nucleotide polymorphisms (SNPs)/genes associated with various important traits in crop species. Typically, GWAS in crops are performed using a panel of inbred lines, where multiple replicates of the same inbred are measured and the average phenotype is taken as the response variable. Here we describe and evaluate single plant GWAS (sp-GWAS) for performing a GWAS on individual plants, which does not require an association panel of inbreds. Instead sp-GWAS relies on the phenotypes and genotypes from individual plants sampled from a randomly mating population. Importantly, we demonstrate how sp-GWAS can be efficiently combined with a bulk segregant analysis (BSA) experiment to rapidly corroborate evidence for significant SNPs.

RESULTS

In this study we used the Shoepeg maize landrace, collected as an open pollinating variety from a farm in Southern Missouri in the 1960's, to evaluate whether sp-GWAS coupled with BSA can efficiently and powerfully used to detect significant association of SNPs for plant height (PH). Plant were grown in 8 locations across two years and in total 768 individuals were genotyped and phenotyped for sp-GWAS. A total of 306 k polymorphic markers in 768 individuals evaluated via association analysis detected 25 significant SNPs (P ≤ 0.00001) for PH. The results from our single-plant GWAS were further validated by bulk segregant analysis (BSA) for PH. BSA sequencing was performed on the same population by selecting tall and short plants as separate bulks. This approach identified 37 genomic regions for plant height. Of the 25 significant SNPs from GWAS, the three most significant SNPs co-localize with regions identified by BSA.

CONCLUSION

Overall, this study demonstrates that sp-GWAS coupled with BSA can be a useful tool for detecting significant SNPs and identifying candidate genes. This result is particularly useful for species/populations where association panels are not readily available.

Genetic diversity and relationships among three Southern African Nguni cattle populations

The Nguni cattle breed has distinct populations that are adapted to the different ecological zones of Southern Africa. This study was carried out to assess genetic diversity and establish the relationships among South African (SA), Mozambican (Landim), and Swazi Nguni cattle populations, using 25 microsatellite markers. Genotypic data were generated from deoxyribonucleic acid (DNA) samples of 90 unrelated individuals of the three cattle populations, collected from government conservations and stud herds. DNA profiles of five local beef breeds were used as the reference populations. Most of the 25 microsatellite markers were highly polymorphic across the studied populations, with an overall polymorphic information content (PIC) mean of 0.676. Genetic diversity within populations was high with expected heterozygosity varying from 0.705 ± 0.024 (Landim) to 0.748 ± 0.021 (SA Nguni) and mean number of alleles being highest in the SA Nguni (7.52 ± 0.42). Average observed heterozygosity (0.597 ± 0.046) compared to the expected heterozygosity (0.719 ± 0.022) was lowest for the Swazi Nguni, which also had a high number of Hardy-Weinberg Equilibrium (HWE) deviated loci (13), confirming the relatively high level of inbreeding (0.158 ± 0.058) in that population. Analysis of molecular variance revealed only 9.61% of the total variation between the populations and 90.39% within populations. A short genetic distance (0.299) was observed between Landim and Swazi Nguni, with the SA Nguni (> 0.500) being the most genetically distant population. The distant relationship between SA Nguni and the other two Nguni cattle populations was further confirmed by a principal coordinates analysis. The three Nguni populations clustered independently from each other, despite some evidence of admixture. Therefore, it can be concluded that SA Nguni, Landim, and Swazi Nguni populations in Southern Africa exhibit high levels of genetic diversity and are genetically distant; with the two latter populations being less genetically apart. These results present useful information for the development of strategies for regional management of animal genetic resources, through conservation and sustainable utilisation.

Genomic clues of the evolutionary history of Bos indicus cattle

Together with their sister subspecies Bos taurus, zebu cattle (Bos indicus) have contributed to important socioeconomic changes that have shaped modern civilizations. Zebu cattle were domesticated in the Indus Valley 8000 years before present (YBP). From the domestication site, they expanded to Africa, East Asia, southwestern Asia and Europe between 4000 and 1300 YBP, intercrossing with B. taurus to form clinal variations of zebu ancestry across the landmass of Afro-Eurasia. In the past 150 years, zebu cattle reached the Americas and Oceania, where they have contributed to the prosperity of emerging economies. The zebu genome is characterized by two mitochondrial haplogroups (I1 and I2), one Y chromosome haplogroup (Y3) and three major autosomal ancestral groups (Indian-Pakistani, African and Chinese). Phenotypically, zebu animals are recognized by their hump, large ears and excess skin. They are rustic, resilient to parasites and capable of bearing the hot and humid climates of the tropics. Many resources are available to study the zebu genome, including commercial arrays of SNP, reference assemblies and publicly available genotypes and whole-genome sequences. Nevertheless, many of these resources were initially developed to support research and subsidize industrial applications in B. taurus, and therefore they can produce bias in data analysis. The combination of genomics with precision agriculture holds great promise for the identification of genetic variants affecting economically important traits such as tick resistance and heat tolerance, which were naturally selected for millennia and played a major role in the evolution of B. indicus cattle.

Genome Analysis Reveals Genetic Admixture and Signature of Selection for Productivity and Environmental Traits in Iraqi Cattle

The Near East cattle are adapted to different agro-ecological zones including desert areas, mountains habitats, and humid regions along the Tigris and Euphrates rivers system. The region was one of the earliest and most significant areas of cattle husbandry. Currently, four main breeds of Iraqi cattle are recognized. Among these, the Jenoubi is found in the southern more humid part of Iraq, while the Rustaqi is found in the middle and drier region of the country. Despite their importance, Iraqi cattle have up to now been poorly characterized at the genome level. Here, we report at a genome-wide level the diversity and signature of positive selection in these two breeds. Thirty-five unrelated Jenoubi cattle, sampled in the Maysan and Basra regions, and 60 Rustaqi cattle, from around Baghdad and Babylon, were genotyped using the Illumina Bovine HD BeadChip (700K). Genetic population structure and diversity level were studied using principal component analysis (PCA), expected heterozygosity (He), observed heterozygosity (Ho), and admixture. Signatures of selection were studied using extended haplotype homozygosity (EHH) (iHS and Rsb) and inter-population Wright's Fst. The results of PCA and admixture analysis, including European taurine, Asian indicine, African indicine, and taurine indicate that the two breeds are crossbreed zebu × taurine, with more zebu background in Jenoubi cattle compared with Rustaqi. The Rustaqi has the greatest mean heterozygosity (He = 0.37) among all breeds. iHS and Rsb signatures of selection analyses identify 68 candidate genes under positive selection in the two Iraqi breeds, while Fst analysis identifies 220 candidate genes including genes related to the innate and acquired immunity responses, different environmental selection pressures (e.g., tick resistance and heat stress), and genes of commercial interest (e.g., marbling score).

Identifying Loci Under Positive Selection in Yellow Korean Cattle (Hanwoo)

Jeju Black cattle is one of the aboriginal Korean cattle breeds that has been isolated in Jeju island for a long time, while Yellow Hanwoo cattle has been extensively selected for beef production traits for the last several decades. Aiming to investigate broader patterns of selection, we genotyped 352 Yellow Hanwoo and 169 Jeju Black cattle using a customized 150 K bovine chip. Our composite selection signals’ analysis to identify selection signatures (cross-population extended haplotype homozygosity [XP-EHH], ΔSAF, and F_ST) identified recent and strong signature of selection near many loci with mutations affecting the traits under strong selection as outlier in Yellow Hanwoo, including SCP2 (P = 8.41 × 10⁻¹⁰) that may be involved in the meat quality. We found nine candidate regions with significant clusters of selection signals, and further bioinformatics analyses of the genes located within these regions revealed mainly genes involved in G-protein coupled receptor signaling pathway (GO:0007186) or olfactory transduction (bta04740), which may be due to adaptation to natural environments in Jeju island. Based on the stronger correlation of Ne₁₀/Ne₁₀₀ ratio between Yellow Hanwoo (0.61) and Jeju Black (0.66) cattle, our results suggest that the difference of chromosomal regions of selection signature between the 2 cattle breeds was due to a consequence of selection processes to adapt to environmental differences between Jeju island and the main inland, Korean peninsula.

Breeding objectives of Brangus cattle in Brazil

Our goal was to define a breeding objective for Brangus beef cattle in Brazil. Bioeconomic models were produced and used to estimate economic values (EVs). The scenarios simulated were typical full-cycle beef production systems that are used in tropical and subtropical regions. The breeding objective contained pregnancy rate (PR), warm carcass weight (WCW), mature cow weight (MCW), number of nematode eggs per gram of faeces (EPG) and tick count (TICK). Two models were used in series to estimate the EV. A deterministic model was used to simulate effects of PR, WCW and MCW on profitability with a constant parasite load. Subsequently, stochastic models were used to estimate economic values for TICK and EPG as consequences of their environmental effects on weight gains, mortality and health costs. The EV of PR, WCW, MCW, EPG and TICK, was US$1.59, US$2.11, -US$0.24, -US$5.35 and -US$20.88, respectively. Results indicate positive emphasis should be placed on PR (12.49%) and WCW (65.07%) with negative emphasis on MCW (13.92%), EPG (2.77%) and TICK (5.75%). In comparison with the indexes usually used, these results suggest a reformulation in the selection indexes of the beef production system in tropical and subtropical regions in order to obtain greater profitability.

Livestock Genomics for Developing Countries - African Examples in Practice

African livestock breeds are numerous and diverse, and typically well adapted to the harsh environment conditions under which they perform. They have been used over centuries to provide livelihoods as well as food and nutritional security. However, African livestock systems are dynamic, with many small- and medium-scale systems transforming, to varying degrees, to become more profitable. In these systems the women and men livestock keepers are often seeking new livestock breeds or genotypes - typically those that increase household income through having enhanced productivity in comparison to traditional breeds while maintaining adaptedness. In recent years genomic approaches have started to be utilized in the identification and development of such breeds, and in this article we describe a number of examples to this end from sub-Saharan Africa. These comprise case studies on: (a) dairy cattle in Kenya and Senegal, as well as sheep in Ethiopia, where genomic approaches aided the identification of the most appropriate breed-type for the local productions systems; (b) a cross-breeding program for dairy cattle in East Africa incorporating genomic selection as well as other applications of genomics; (c) ongoing work toward creating a new cattle breed for East Africa that is both productive and resistant to trypanosomiasis; and (d) the use of African cattle as resource populations to identify genomic variants of economic or ecological significance, including a specific case where the discovery data was from a community based breeding program for small ruminants in Ethiopia. Lessons learnt from the various case studies are highlighted, and the concluding section of the paper gives recommendations for African livestock systems to increasingly capitalize on genomic technologies.

Genome-wide associations and functional gene analyses for endoparasite resistance in an endangered population of native German Black Pied cattle

Background

Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections.

Results

A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (p_Bonf = 4.47 × 10^− 7), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10^− 6 on BTA 1 to 2.18 × 10^− 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections.

Conclusions

We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5659-4) contains supplementary material, which is available to authorized users.

Towards a new phenotype for tick resistance in beef and dairy cattle: a review

About 80% of the world’s cattle are affected by ticks and tick-borne diseases, both of which cause significant production losses. Cattle host resistance to ticks is the most important factor affecting the economics of tick control, but it is largely neglected in tick-control programs due to technical difficulties and costs associated with identifying individual-animal variation in resistance. The present paper reviews the scientific literature to identify factors affecting resistance of cattle to ticks and the biological mechanisms of host tick resistance, to develop alternative phenotype(s) for tick resistance. If new cost-effective phenotype(s) can be developed and validated, then tick resistance of cattle could be genetically improved using genomic selection, and incorporated into breeding objectives to simultaneously improve cattle productive attributes and tick resistance. The phenotype(s) could also be used to improve tick control by using cattle management. On the basis of the present review, it is recommended that three possible phenotypes (haemolytic analysis; measures of skin hypersensitivity reactions; simplified artificial tick infestations) be further developed to determine their practical feasibility for consistently, cost-effectively and reliably measuring cattle tick resistance in thousands of individual animals in commercial and smallholder farmer herds in tropical and subtropical areas globally. During evaluation of these potential new phenotypes, additional measurements should be included to determine the possibility of developing a volatile-based resistance phenotype, to simultaneously improve cattle resistance to both ticks and biting flies. Because the current measurements of volatile chemistry do not satisfy the requirements of a simple, cost-effective phenotype for use in commercial cattle herds, consideration should also be given to inclusion of potentially simpler measures to enable indirect genetic selection for volatile-based resistance to ticks.

Combining Landscape Genomics and Ecological Modelling to Investigate Local Adaptation of Indigenous Ugandan Cattle to East Coast Fever

East Coast fever (ECF) is a fatal sickness affecting cattle populations of eastern, central, and southern Africa. The disease is transmitted by the tick Rhipicephalus appendiculatus, and caused by the protozoan Theileria parva parva, which invades host lymphocytes and promotes their clonal expansion. Importantly, indigenous cattle show tolerance to infection in ECF-endemically stable areas. Here, the putative genetic bases underlying ECF-tolerance were investigated using molecular data and epidemiological information from 823 indigenous cattle from Uganda. Vector distribution and host infection risk were estimated over the study area and subsequently tested as triggers of local adaptation by means of landscape genomics analysis. We identified 41 and seven candidate adaptive loci for tick resistance and infection tolerance, respectively. Among the genes associated with the candidate adaptive loci are PRKG1 and SLA2. PRKG1 was already described as associated with tick resistance in indigenous South African cattle, due to its role into inflammatory response. SLA2 is part of the regulatory pathways involved into lymphocytes' proliferation. Additionally, local ancestry analysis suggested the zebuine origin of the genomic region candidate for tick resistance.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F₂ population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F₂ population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

Genetic Improvement in South African Livestock: Can Genomics Bridge the Gap Between the Developed and Developing Sectors?

South Africa (SA) holds a unique position on the African continent with a rich diversity in terms of available livestock resources, vegetation, climatic regions and cultures. The livestock sector has been characterized by a dual system of a highly developed commercial sector using modern technology vs. a developing sector including emerging and smallholder farmers. Emerging farmers typically aim to join the commercial sector, but lag behind with regard to the use of modern genetic technologies, while smallholder farmers use traditional practices aimed at subsistence. Several factors influence potential application of genomics by the livestock industries, which include available research funding, socio-economic constraints and extension services. State funded Beef and Dairy genomic programs have been established with the aim of building reference populations for genomic selection with most of the potential beneficiaries in the well-developed commercial sector. The structure of the beef, dairy and small stock industries is fragmented and the outcomes of selection strategies are not perceived as an advantage by the processing industry or the consumer. The indigenous and local composites represent approximately 40% of the total beef and sheep populations and present valuable genetic resources. Genomic research has mostly provided insight on genetic biodiversity of these resources, with limited attention to novel phenotypes associated with adaptation or disease tolerance. Genetic improvement of livestock through genomic technology needs to address the role of adapted breeds in challenging environments, increasing reproductive and growth efficiency. National animal recording schemes contributed significantly to progress in the developed sector with regard to genetic evaluations and estimated breeding values (EBV) as a selection tool over the past three decades. The challenge remains on moving the focus to novel traits for increasing efficiency and addressing welfare and environmental issues. Genetic research programs are required that will be directed to bridge the gap between the elite breeders and the developing livestock sector. The aim of this review was to provide a perspective on the dichotomy in the South African livestock sector arguing that a realistic approach to the use of genomics in beef, dairy and small stock is required to ensure sustainable long term genetic progress.

Cattle Tick Rhipicephalus microplus-Host Interface: A Review of Resistant and Susceptible Host Responses

Ticks are able to transmit tick-borne infectious agents to vertebrate hosts which cause major constraints to public and livestock health. The costs associated with mortality, relapse, treatments, and decreased production yields are economically significant. Ticks adapted to a hematophagous existence after the vertebrate hemostatic system evolved into a multi-layered defense system against foreign invasion (pathogens and ectoparasites), blood loss, and immune responses. Subsequently, ticks evolved by developing an ability to suppress the vertebrate host immune system with a devastating impact particularly for exotic and crossbred cattle. Host genetics defines the immune responsiveness against ticks and tick-borne pathogens. To gain an insight into the naturally acquired resistant and susceptible cattle breed against ticks, studies have been conducted comparing the incidence of tick infestation on bovine hosts from divergent genetic backgrounds. It is well-documented that purebred and crossbred Bos taurus indicus cattle are more resistant to ticks and tick-borne pathogens compared to purebred European Bos taurus taurus cattle. Genetic studies identifying Quantitative Trait Loci markers using microsatellites and SNPs have been inconsistent with very low percentages relating phenotypic variation with tick infestation. Several skin gene expression and immunological studies have been undertaken using different breeds, different samples (peripheral blood, skin with tick feeding), infestation protocols and geographic environments. Susceptible breeds were commonly found to be associated with the increased expression of toll like receptors, MHC Class II, calcium binding proteins, and complement factors with an increased presence of neutrophils in the skin following tick feeding. Resistant breeds had higher levels of T cells present in the skin prior to tick infestation and thus seem to respond to ticks more efficiently. The skin of resistant breeds also contained higher numbers of eosinophils, mast cells and basophils with up-regulated proteases, cathepsins, keratins, collagens and extracellular matrix proteins in response to feeding ticks. Here we review immunological and molecular determinants that explore the cattle tick Rhipicephalus microplus-host resistance phenomenon as well as contemplating new insights and future directions to study tick resistance and susceptibility, in order to facilitate interventions for tick control.

Analyses of reaction norms reveal new chromosome regions associated with tick resistance in cattle

Despite single nucleotide polymorphism (SNP) availability and frequent cost reduction has allowed genome-wide association studies even in complex traits as tick resistance, the use of this information source in SNP by environment interaction context is unknown for many economically important traits in cattle. We aimed at identifying putative genomic regions explaining differences in tick resistance in Hereford and Braford cattle under SNP by environment point of view as well as to identify candidate genes derived from outliers/significant markers. The environment was defined as contemporary group means of tick counts, since they seemed to be the most appropriate entities to describe the environmental gradient in beef cattle. A total of 4363 animals having tick counts (n=10 673) originated from 197 sires and 3966 dams were used. Genotypes were acquired on 3591 of these cattle. From top 1% SNPs (410) having the greatest effects in each environment, 75 were consistently relevant in all environments, which indicated SNP by environment interaction. The outliers/significant SNPs were mapped on chromosomes 1, 2, 5, 6, 7, 9, 11, 13, 14, 15, 16, 18, 21, 23, 24, 26 and 28, and potential candidate genes were detected across environments. The presence of SNP by environment interaction for tick resistance indicates that genetic expression of resistance depends upon tick burden. Markers with major portion of genetic variance explained across environments appeared to be close to genes with different direct or indirect functions related to immune system, inflammatory process and mechanisms of tissue destruction/repair, such as energy metabolism and cell differentiation.

Genetic parameters for tick counts across months for different tick species and anatomical locations in South African Nguni cattle

The objective of the study was to characterise genetic parameters across months for different tick species and anatomical locations in South African Nguni cattle. Tick counts were conducted monthly, over a 2-year period, on 586 Nguni cattle under natural infestation, from four herds located in different provinces of South Africa. The counts were recorded for six species of ticks (Amblyomma hebraeum, Rhipicephalus evertsi evertsi, Rhipicephalus decoleratus and microplus (Boofilids), Rhipicephalus appendiculatus, Rhipicephalus simus and Hyalomma marginatum) attached on eight anatomical locations on the animals and were summed by species and anatomical location. Heritability estimates, phenotypic and genetic correlations were estimated on a monthly basis using mixed linear models, fitting univariate and bivariate sire models. Fixed effects considered were location, sex, year and age as a covariate. Tick counts were higher in the hot months, and A. hebraeum was the most dominant tick species. Heritability estimates for tick count varied by month and trait and ranged from 0 to 0.89. Genetic correlations were mostly positive, and low to high, with some negative correlations with high standard error. Phenotypic correlations were low to moderate. In general, high genetic correlations were observed between whole body count and the anatomical location counts, suggesting that it may not be necessary to conduct whole body counts. Counts from the belly and perineum appeared to be the most suitable surrogate traits for whole body count. These findings provide useful information for developing strategies for the practical implementation of genetic selection, as a supplement to the traditional tick control measures.

Tag SNP selection for prediction of tick resistance in Brazilian Braford and Hereford cattle breeds using Bayesian methods

BACKGROUND

Cattle resistance to ticks is known to be under genetic control with a complex biological mechanism within and among breeds. Our aim was to identify genomic segments and tag single nucleotide polymorphisms (SNPs) associated with tick-resistance in Hereford and Braford cattle. The predictive performance of a very low-density tag SNP panel was estimated and compared with results obtained with a 50 K SNP dataset.

RESULTS

BayesB (π = 0.99) was initially applied in a genome-wide association study (GWAS) for this complex trait by using deregressed estimated breeding values for tick counts and 41,045 SNP genotypes from 3455 animals raised in southern Brazil. To estimate the combined effect of a genomic region that is potentially associated with quantitative trait loci (QTL), 2519 non-overlapping 1-Mb windows that varied in SNP number were defined, with the top 48 windows including 914 SNPs and explaining more than 20% of the estimated genetic variance for tick resistance. Subsequently, the most informative SNPs were selected based on Bayesian parameters (model frequency and t-like statistics), linkage disequilibrium and minor allele frequency to propose a very low-density 58-SNP panel. Some of these tag SNPs mapped close to or within genes and pseudogenes that are functionally related to tick resistance. Prediction ability of this SNP panel was investigated by cross-validation using K-means and random clustering and a BayesA model to predict direct genomic values. Accuracies from these cross-validations were 0.27 ± 0.09 and 0.30 ± 0.09 for the K-means and random clustering groups, respectively, compared to respective values of 0.37 ± 0.08 and 0.43 ± 0.08 when using all 41,045 SNPs and BayesB with π = 0.99, or of 0.28 ± 0.07 and 0.40 ± 0.08 with π = 0.999.

CONCLUSIONS

Bayesian GWAS model parameters can be used to select tag SNPs for a very low-density panel, which will include SNPs that are potentially linked to functional genes. It can be useful for cost-effective genomic selection tools, when one or a few key complex traits are of interest.