BackgroundThe existence of moderate to high levels of linkage disequilibrium (LD) between genetic markers and quantitative trait loci (QTL) affecting traits of interest is fundamental for the success of genome-wide association (GWAS) and genomic selection (GS) studies. Knowledge about the extent and the pattern of LD in livestock populations is essential to determine the density of single nucleotide polymorphisms (SNP) required for accurate GWAS and GS. Moreover, observed LD is related to historical effective population sizes (Ne), and can provide insights into the genetic diversity history of populations. Estimates of the consistency of linkage phase across breeds (RH,B) can be used to determine if there is sufficient relationship to use pooled reference populations in multi-breed GS programs. The objective of this study was to estimate LD levels, persistence of phase and effective population size in Hereford and Braford cattle populations sampled in Brazil.ResultsMean LD estimates, measured using the squared correlation of alleles at two loci (r2), obtained between adjacent SNP across all chromosomes were 0.21 ± 0.27 for Herefords (391 samples with 41,241 SNP) and 0.16 ± 0.22 for Brafords (2044 samples and 41,207 SNP). Estimated r2 was > 0.2 and 0.3, respectively, for 34 and 25 % of adjacent markers in Herefords, and 26 and 17 % in Brafords. Estimated Ne for Brafords and Herefords at the current generation was 220 and 153 individuals, respectively. The two breeds demonstrated moderate to strong persistence of phase at all distances (RH,B = 0.53 to 0.97). The largest phase correlations were found in the 0 to 50 Kb bins (RH,B = 0.92 to 0.97). Estimated LD decreased rapidly with increasing distance between SNP, however, useful linkage for GWAS and GS (r2 > 0.2) was found spanning to ~50 Kb.ConclusionsPanels containing about 50,000 and 150,000 SNP markers are necessary to detect minimal levels of LD between adjacent markers that would be useful for GWAS and GS studies to Hereford and Braford breeds, respectively. Markers are expected to be linked to the same QTL alleles in distances < 50 Kb in both populations due to observed high persistence of phase levels.

BackgroundGenetic 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.ResultsIntergroup 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.ConclusionThis 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.

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.

This study aimed to verify the influence of bovine genetic resistance on biological traits of the Rhipicephalus (Boophilus) microplus tick. Genetic resistance or susceptibility was determined according to breeding values for tick counts, predicted using a dataset of 9007 Hereford and Braford (Hereford×Zebu) bovines naturally infested and raised under extensive production systems in southern Brazil. From a total of 974 Braford heifers born in 2008, 20 were classified as genetically tick-resistant and 20 classified as genetically tick-susceptible, and used to obtain the ticks samples used in this study. The 40 heifers were exposed to four subsequent artificial infestations with approximately 20,000 larvae at 14-day intervals. From the 19th to 23rd day of each infestation tick counts were performed on the left body side of the heifers. Engorged ticks were manually collected on the day of highest observed burden after each infestation. Tick counts on susceptible heifers were 5.5, 10.5, 11.1 and 6.9 times larger than on resistant heifers, respectively, after the first, second, third and fourth artificial infestations. In the third infestation, ticks from resistant heifers showed lower egg production index (P<0.0001) than ticks from susceptible heifers. In the fourth infestation, ticks from susceptible group showed higher egg mass weight (P<0.05) and nutrient index (P<0.0001) than ticks from resistant heifers. Tick initial weights showed a positive association with egg production index in susceptible heifers (P<0.05) and a negative association in the resistant group (P<0.05), suggesting a host defense mechanism that reduces the conversion efficiency of ingested blood to eggs in engorged ticks from resistant cattle. This shows that bovine genetic tick resistance, in addition to affecting the number of ticks carried by the animals, also affected the egg mass weight, egg production and nutrient indexes of ticks. The results of the present study imply that the selection of resistant animals could be used as a strategic tool for tick control in production systems, reducing infestation levels on cattle and environment.

No abstract