Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle.

BACKGROUND: Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (FROH), genomic relationship matrix approach (FGRM) and based on the observed versus expected number of homozygous genotypes (FHOM), and from pedigree-based coefficient (FPED). RESULTS: ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH1-2 Mb) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of FROH, FGRM, FHOM, and FPED approaches. FPED estimates ranged from 0.00 to 0.327 and FROH from 0.001 to 0.201. Low to moderate correlations were observed between FPED-FROH and FGRM-FROH, with values ranging from -0.11 to 0.51. Low to high correlations were observed between FROH-FHOM and moderate between FPED-FHOM and FGRM-FHOM. Correlations between FROH from different lengths and FPED gradually increased with ROH length. CONCLUSIONS: Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low FPED-FROH correlations for small segments indicate that FPED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits.

BACKGROUND: The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (FROH), genomic relationship matrix (FGRM), and pedigree-based coefficient (FPED). RESULTS: The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of FPED and FROH, and their correlations (- 0.05 to 0.26) were low. Estimates of correlation between FGRM-FPED was zero, while the correlation (- 0.01 to - 0.07) between FGRM-FROH decreased as a function of ROH length, except for FROH > 8Mb (- 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands. CONCLUSIONS: Low FPED-FROH correlation estimates indicate that FPED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and FROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Genomic selection in dairy cattle simulated populations.

Genomic selection is arguably the most promising tool for improving genetic gain in domestic animals to emerge in the last few decades, but is an expensive process. The aim of this study was to evaluate the economic impact related to the implementation of genomic selection in a simulated dairy cattle population. The software QMSim was used to simulate genomic and phenotypic data. The simulated genome contained 30 chromosomes with 100 cm each, 1666 SNPs markers equally spread and 266 QTLs randomly designated for each chromosome. The numbers of markers and QTLs were designated according to information available from Animal QTL (http://www.animalgenome.org/QTLdb) and Bovine QTL (http://bovineqtl.tamu.edu/). The allelic frequency changes were assigned in a gamma distribution with alpha parameters equal to 0·4. Recurrent mutation rates of 1·0e-4 were assumed to apply to markers and QTLs. A historic population of 1000 individuals was generated and the total number of animals was reduced gradually along 850 generations until we obtained a number of 200 animals in the last generation, characterizing a bottleneck effect. Progenies were created along generations from random mating of the male and female gametes, assuming the same proportion of both genders. Than the population was extended for another 150 generations until we obtained 17 000 animals, with only 320 male individuals in the last generation. After this period a 25 year of selection was simulated taking into account a trait limited by sex with heritability of 0·30 (i.e. milk yield), one progeny/cow/year and variance equal to 1·0. Annually, 320 bulls were mated with 16 000 dams, assuming a replacement rate of 60 and 40% for males and females, respectively. Selection and discard criteria were based in four strategies to obtain the EBVs assuming as breeding objective to maximize milk yield. The progeny replaced the discarded animals creating an overlapping generation structure. The selection strategies were: RS is selection based on random values; PS is selection based on phenotypic values; Blup is selection based on EBVs estimated by BLUP; and GEBV is selection based on genomic estimated breeding values in one step, using high (GBlup) and low (GBlupi) density panels. Results indicated that the genetic evaluation using the aid of genomic information could provide better genetic gain rates in dairy cattle breeding programs as well as reduce the average inbreeding coefficient in the population. The economic viability indicators showed that only Blup and GBlup/GBlupi strategies, the ones that used milk control and genetic evaluation were economic viable, considering a discount rate of 6·32% per year.

Genetic analysis on accumulated productivity and calving intervals in Nelore cattle.

The aim of this study was to estimate genetic parameters for accumulated productivity (ACP), first calving interval (CI1), second calving interval (CI2), and mean calving interval (MCI) in Nelore beef cattle. The ACP trait is a reproduction index and comprises the total number of calves born per dam, weight of weaned calves, and age of the dam at calving. Genetic parameters were estimated by the average information restricted maximum likelihood method in two-trait analyses. The average heritability estimate for ACP was 0.17 (0.03). For CI1, CI2, and MCI, the heritability estimates were 0.02, 0.02, and 0.06, respectively. Genetic correlations between ACP with CI1, CI2, and MCI were -0.16 ± 0.47, -0.29 ± 0.53, and -0.40 ± 0.27, respectively. Despite of the low heritability estimates obtained in our study, reproduction traits should be further studied and their inclusion in the selection criteria must be evaluated in order to improve the performance of females because these traits are of great economic importance in beef cattle. Accumulated productivity could contribute to decrease the mean calving interval in Nelore cattle.

Effect of IGF1, GH, and PIT1 markers on the genetic parameters of growth and reproduction traits in Canchim cattle.

The availability of dense genomic information has increased genome-wide association studies for the bovine species; however research to assess the effect of single genes on production traits is still important to elucidate the genes functions. On this study the association of IGF1, GH, and PIT1 markers with growth and reproductive traits (birth weight, weaning weight, weight at 12 and 18 months of age, preweaning average daily weight gain, age and weight at first calving, and scrotal circumference at 12 and 18 months of age) were assessed by means of the variance component approach. The phenotypes were adjusted and then analyzed under two animal models, one which considered the polygenic and genotype (IGF1, GH or PIT1 markers) effects (Model 1), and the other which considers only the polygenic effect (Model 2). When the likelihood ratio test and the Bonferroni correction was applied at 5 % significance level, the genetic markers for the IGF1, GH, and PIT1 genes did not influence significantly the traits (p > 0.002). However, evidence of association of IGF1 with birth weight (p = 0.06) and GH with weight at first calving (p = 0.03) and with weight at 12 months of age (p = 0.08) was observed. In conclusion we could not confirm the associations between IGF1, GH, and PIT1 and growth traits that were previously reported in Canchim cattle, and no association was observed between these genes and reproductive traits. Future studies involving functional markers of IGF1, GH and PIT1 genes may help to clarify the role of these genes in growth and reproductive processes.

Genetic parameter estimation for pre- and post-weaning traits in Brahman cattle in Brazil.

Beef cattle producers in Brazil use body weight traits as breeding program selection criteria due to their great economic importance. The objectives of this study were to evaluate different animal models, estimate genetic parameters, and define the most fitting model for Brahman cattle body weight standardized at 120 (BW120), 210 (BW210), 365 (BW365), 450 (BW450), and 550 (BW550) days of age. To estimate genetic parameters, single-, two-, and multi-trait analyses were performed using the animal model. The likelihood ratio test was verified between all models. For BW120 and BW210, additive direct genetic, maternal genetic, maternal permanent environment, and residual effects were considered, while for BW365 and BW450, additive direct genetic, maternal genetic, and residual effects were considered. Finally, for BW550, additive direct genetic and residual effects were considered. Estimates of direct heritability for BW120 were similar in all analyses; however, for the other traits, multi-trait analysis resulted in higher estimates. The maternal heritability and proportion of maternal permanent environmental variance to total variance were minimal in multi-trait analyses. Genetic, environmental, and phenotypic correlations were of high magnitude between all traits. Multi-trait analyses would aid in the parameter estimation for body weight at older ages because they are usually affected by a lower number of animals with phenotypic information due to culling and mortality.

Modern livestock farming under tropical conditions using sensors in grazing systems.

The aim of this study was to evaluate a commercial sensor-a three-axis accelerometer-to predict animal behavior with a variety of conditions in tropical grazing systems. The sensor was positioned on the underjaw of young bulls to detect the animals' movements. A total of 22 animals were monitored in a grazing system, during both seasons (wet and dry), with different quality and quantity forage allowance. The machine learning (ML) methods used were random forest (RF), convolutional neural net and linear discriminant analysis; the metrics used to determine the best method were accuracy, Kappa coefficient, and a confusion matrix. After predicting animal behavior using the best ML method, a forecast for animal performance was developed using a mechanistic model: multiple linear regression to correlate intermediate average daily gain (iADG) observed versus iADG predicted. The best ML method yielded accuracy of 0.821 and Kappa coefficient of 0.704, was RF. From the forecast for animal performance, the Pearson correlation was 0.795 and the mean square error was 0.062. Hence, the commercial Ovi-bovi sensor, which is a three-axis accelerometer, can act as a powerful tool for predicting animal behavior in beef cattle production developed under a variety tropical grazing condition.

Cross-population selection signatures in Canchim composite beef cattle.

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Genome-Wide Association Study Provides Insights into Important Genes for Reproductive Traits in Nelore Cattle.

The identification of genomic regions associated with reproductive traits as well as their biological processes allows a better understanding of the phenotypic variability of these traits. This information could be applied to animal breeding programs to accelerate genetic gain. The aim of this study was to evaluate the association between single nucleotide polymorphisms (SNP) with a scrotal circumference at 365 days of age (SC365) and at 450 days of age (SC450), gestation length (GL) as a calf trait, age at first calving (AFC), accumulated productivity (ACP), heifer early calving until 30 months (HC30), and stayability (STAY) traits, in order to identify candidate genes and biological pathways associated with reproductive traits in Nelore cattle. The data set consisted of pedigree, phenotypes, and genotypes of Nelore cattle from the "Associação Nacional de Criadores e Pesquisadores" (ANCP). The association analyses were performed using the Weighted Single-Step Genome-Wide Association method; the regions, consisting of 10 consecutive SNP, which explained more than 0.5% of additive genetic variance, were considered as a significant association. A total of 3, 6, 7, 5, 10, 25, and 12 windows were associated with SC355, SC450, GL, AFC, ACP, HC30, and STAY, respectively. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, among them CAMK1D, TASP1, ACOXL, RAB11FIP5, and SFXN5. Moreover, the genes were enriched in functional terms, like negative regulation of fat cell differentiation, fatty acid alpha-oxidation, and sphingolipids signaling pathway. The identification of the genes associated with the traits, as well as genes enriched in the terms and pathway mentioned above, should contribute to future biological validation studies and may be used as candidate genes in Nelore breeding programs.

Genetic evaluation for reproductive and productive traits in Brahman cattle.

Our objectives were to estimate genetic parameters for male and female reproductive traits and their genetic correlations with body weight and carcass traits, evaluate the genetic trends over the years, and verify the effect of inbreeding on the phenotypes of Brahman cattle. The traits evaluated were body weights at 120, 210, 365, and 450 days of age (W120, W210, W365, and W450); scrotal circumference at 365 and 450 days of age (SC365 and SC450), age at first calving (AFC), gestation length (GL), stayability (STAY), ribeye area (REA), backfat thickness (BFT), and rump fat thickness (RFT). Direct heritability estimates ranged from 0.10 ± 0.03 (AFC) to 0.43 ± 0.06 (GL). Maternal heritability estimates for body weights, scrotal circumferences, and GL ranged from 0.07 ± 0.02 to 0.15 ± 0.03. The proportion of the maternal permanent environment for W120 and W210 was equal to 0.11 ± 0.02. Genetic correlations varied between -0.60 ± 0.25 (STAY and BFT) to 0.97 ± 0.01 (W365 and W450). Except for AFC, all genetic trends were significant (p < 0.05) and presented favorable annual genetic gains. Unfavorable effects due to the increase of inbreeding coefficients were observed for body weights and AFC, suggesting greater attention be paid to the applied mating systems to control inbreeding. Reproductive traits, such as AFC and STAY, could be assisted indirectly by scrotal circumference selection. The emphasis applied to body weight selection, especially at W210, may assist REA. The BFT and RFT traits presented genetic variability and have responded to selection, although not included in the Brahman selection index.

Exploring the genetic architecture of feed efficiency traits in chickens.

Chicken feed efficiency (FE) traits are the most important economic traits in broiler production. Several studies evaluating genetic factors affecting food consumption in chickens are available. However, most of these studies identified genomic regions containing putative quantitative trait loci for each trait separately. It is still a challenge to find common gene networks related to these traits. Therefore, here, a genome-wide association study (GWAS) was conducted to explore candidate genomic regions responsible for Feed Intake (FI), Body Weight Gain (BWG) and Feed Conversion Ratio (FCR) traits and their gene networks. A total of 1430 broilers from an experimental population was genotyped with the high density Affymetrix 600K SNP array. A total of 119 associated SNPs located in 20 chromosomes were identified, where some of them were common in more than one FE trait. In addition, novel genomic regions were prospected considering the SNPs dominance effects and sex interaction, identifying putative candidate genes only when these effects were fit in the model. Relevant candidate genes such as ATRNL1, PIK3C2A, PTPRN2, SORCS3 and gga-mir-1759 were highlighted in this study helping to elucidate the genomic architecture of feed efficiency traits. These results provide new insights on the mechanisms underlying the consumption and utilization of food in chickens.

Immune response pattern of the popliteal lymph nodes of dogs with visceral leishmaniasis.

The present study aimed to estimate the cell response and parasite load in the popliteal lymph nodes of dogs with visceral leishmaniasis (VL), comparing these findings with the clinical staging of the disease. From the necropsy, 33 dogs were classified as symptomatic (S), asymptomatic (A), or oligosymptomatic (O). Cytology and histopathology were used to determine any presence of microscopic lesions and immunohistochemistry, for parasite load. Dog hyperimmune serum was used as the primary antibody. The inflammatory infiltrate in lymph nodes consisted of macrophages and plasmocytes. The granulomas invaded the trabecular and sinusoid regions and sometimes compressed the lymphocytes of the cortical region (atrophy) and medullary cord cells. Parasite load intensity was unrelated to the density of the macrophages infiltrating the lymph node. Significant differences in parasite load (P < 0.05) were observed between the three groups of infected dogs. Follicular hyperplasia of the cortical region occurred among A and O, while follicular atrophy predominated among S. The parasite load was the greatest among S, followed by O. It can be concluded that, regardless of clinical condition, the most evident cell response consisted of macrophages and plasmocytes. Lymphoid atrophy was observed among animals with intense granulomatous reaction and high parasite load, such as among the symptomatic dogs (P < 0.05). Likewise, the oligosymptomatic dogs also presented high density of parasites in the lymph nodes. Thus, we can confirm that dogs with clinical manifestations of VL have an immune system that is less effective for controlling infection by Leishmania chagasi, thereby favoring parasite multiplication.

Evaluation of imputation accuracy using the combination of two high-density panels in Nelore beef cattle.

This study compared imputation from lower-density commercial and customized panels to high-density panels and a combined panel (Illumina and Affymetrix) in Nelore beef cattle. Additionally, linkage disequilibrium and haplotype block conformation were estimated in individual high-density panels and compared with corresponding values in the combined panel after imputation. Overall, 814 animals were genotyped using BovineHD BeadChip (IllumHD), and 93 of these animals were also genotyped using the Axion Genome-Wide BOS 1 Array Plate (AffyHD). In general, customization considering linkage disequilibrium and minor allele frequency had the highest accuracies. The IllumHD panel had higher values of linkage disequilibrium for short distances between SNPs than AffyHD and the combined panel. The combined panel had an increased number of small haplotype blocks. The use of a combined panel is recommended due to its increased density and number of haplotype blocks, which in turn increase the probability of a marker being close to a quantitative trait locus of interest. Considering common SNPs between IllumHD and AffyHD for the customization of a low-density panel increases the imputation accuracy for IllumHD, AffyHD and the combined panel.

Genomic structure of a crossbred Landrace pig population.

Single nucleotide polymorphism (SNP) markers are used to study population structure and conservation genetics, which permits assessing similarities regarding the linkage disequilibrium and information about the relationship among individuals. To investigate the population genomic structure of 300 females and 25 males from a commercial maternal pig line we analyzed linkage disequilibrium extent, inbreeding coefficients using genomic and conventional pedigree data, and population stratification. The average linkage disequilibrium (r2) was 0.291 ± 0.312 for all adjacent SNPs, distancing less than 100 Kb (kilobase) between markers. The average inbreeding coefficients obtained from runs of homozygosity (ROH) and pedigree analyses were 0.119 and 0.0001, respectively. Low correlation was observed between the inbreeding coefficients possibly as a result of genetic recombination effect accounted for the ROH estimates or caused by pedigree identification errors. A large number of long ROHs might indicate recent inbreeding events in the studied population. A total of 36 homozygous segments were found in more than 30% of the population and these ROH harbor genes associated with reproductive traits. The population stratification analysis indicated that this population was possibly originated from two distinct populations, which is a result from crossings between the eastern and western breeds used in the formation of the line. Our findings provide support to understand the genetic structure of swine populations and may assist breeding companies to avoid a high level of inbreeding coefficients to maintain genetic diversity, showing the effectiveness of using genome-wide SNP information for quantifying inbreeding when the pedigree was incomplete or incorrect.

Genome-wide association study provides insights into genes related with horn development in Nelore beef cattle.

The causal mutation for polledness in Nelore (Bos taurus indicus) breed seems to have appeared first in Brazil in 1957. The expression of the polled trait is known to be ruled by a few groups of alleles in taurine breeds; however, the genetic basis of this trait in indicine cattle is still unclear. The aim of this study was to identify genomic regions associated with the hornless trait in a commercial Nelore population. A total of 107,294 animals had phenotypes recorded and 2,238 were genotyped/imputed for 777k SNP. The weighted single-step approach for genome-wide association study (WssGWAS) was used to estimate the SNP effects and variances accounted for by 1 Mb sliding SNP windows. A centromeric region of chromosome 1 with 3.11 Mb size (BTA1: 878,631-3,987,104 bp) was found to be associated with hornless in the studied population. A total of 28 protein-coding genes are mapped in this region, including the taurine Polled locus and the IFNAR1, IFNAR2, IFNGR2, KRTAP11-1, MIS18A, OLIG1, OLIG2, and SOD1 genes, which expression can be related to the horn formation as described in literature. The functional enrichment analysis by DAVID tool revealed cytokine-cytokine receptor interaction, JAK-STAT signaling, natural killer cell mediated cytotoxicity, and osteoclast differentiation pathways as significant (P < 0.05). In addition, a runs of homozygosity (ROH) analysis identified a ROH island in polled animals with 2.47 Mb inside the region identified by WssGWAS. Polledness in Nelore cattle is associated with one region in the genome with 3.1 Mb size in chromosome 1. Several genes are harbored in this region, and they may act together in the determination of the polled/horned phenotype. Fine mapping the locus responsible for polled trait in Nelore breed and the identification of the molecular mechanisms regulating the horn growth deserve further investigation.

Evaluation of transformation growth factor beta1, interleukin-10, and interferon-gamma in male symptomatic and asymptomatic dogs naturally infected by Leishmania (Leishmania) chagasi.

The aims of this study were to evaluate the immunomodulatory role of TGF-beta1, IL-10, and INF-gamma in spleen and liver extracts and supernatant cultures of white spleen cells from male symptomatic and asymptomatic dogs, naturally infected by Leishmania (Leishmania) chagasi. Thirty dogs from Araçatuba, São Paulo, Brazil, an endemic leishmaniosis area, were selected by positive ELISA serological reaction for Leishmania sp. and divided into two groups: asymptomatic (n=15) and symptomatic (n=15) consisting of animals with at least three characteristic signs (fever, dermatitis, lymphoadenopathy, onychogryphosis, weight loss, cachexia, locomotion problems, conjunctivitis, epistaxis, hepatosplenomegaly, edema, and apathy). After euthanasia, spleen and liver fragments were collected for ex vivo quantification of TGF-beta1, IL-10, and INF-gamma. Naturally active in vitro produced TGF-beta1 was also evaluated in spleen cell culture supernatant. Spleen and liver extract of asymptomatic dogs had higher mean TGF-beta1 levels than symptomatic dogs. High concentrations of IL-10 were found in spleen, and mainly in liver extract of both groups. Higher INF-gamma concentrations were found in spleen extracts of symptomatic dogs, and in liver extracts of asymptomatic dogs. Extract of this cytokine was lower in spleen extract. Although INF-gamma is being produced in canine infection, mean levels of TGF-beta1 and IL-10 from spleen and liver extracts were quantitatively much higher; suggesting that immune response in both asymptomatic and symptomatic dogs was predominantly type Th2.

An investigation of Leishmania spp. in Didelphis spp. from urban and peri-urban areas in Bauru (São Paulo, Brazil).

Blood and bone marrow samples were taken from 112 Didelphis spp., collected between March 2005 and February 2006, from urban and peri-urban areas of Bauru, São Paulo State, Brazil, to evaluate the hypothesis that these animals might constitute a reservoir of Leishmania spp. Anti-Leishmania ssp. antibodies were screened in the serum samples using an enzyme-linked immuno-sorbent assay (ELISA) and the polymerase chain reaction (PCR). PCR was performed on fragments of DNA samples from Leishmania spp. using primers 13A and 13B, and showed a positive outcome in 91.6% of the 112 samples tested. Of the 107 samples analyzed by ELISA, 71% were positive. Evidence of epidemiological risk factors such as a circulating parasite and freely moving vectors suggests that Didelphis spp. may participate in the transmission cycle of Leishmania spp. in Bauru.

Estimation of genetic parameters for longevity considering the cow's age at last calving.

The aim of this study was to estimate heritability and predict breeding values for longevity among cows in herds of Nellore breed, considering the trait cow's age at last calving (ALC), by means of survival analysis methodology. The records of 11,791 animals from 22 farms were used. The variable ALC has been used by a criterion that made it possible to include cows not only at their first calving but also at their ninth calving. The criterion used was the difference between the date of each cow's last calving and the date of the last calving on each farm. If this difference was greater than 36 months, the cow was considered to have failed and uncensored. If not, this cow was censored, thus indicating that future calving remained possible for this cow. The survival model used for the analyses was the proportional hazards model, and the base risk was given by a Weibull distribution. The heritability estimate obtained was equal to 0.25. It was found that the ALC variable had the capacity to respond to selection for the purpose of increasing the longevity of the cows in the herds.

Study on the introgression of beef breeds in Canchim cattle using single nucleotide polymorphism markers.

The aim of this study was to evaluate the level of introgression of breeds in the Canchim (CA: 62.5% Charolais-37.5% Zebu) and MA genetic group (MA: 65.6% Charolais-34.4% Zebu) cattle using genomic information on Charolais (CH), Nelore (NE), and Indubrasil (IB) breeds. The number of animals used was 395 (CA and MA), 763 (NE), 338 (CH), and 37 (IB). The Bovine50SNP BeadChip from Illumina panel was used to estimate the levels of introgression of breeds considering the Maximum likelihood, Bayesian, and Single Regression method. After genotype quality control, 32,308 SNPs were considered in the analysis. Furthermore, three thresholds to prune out SNPs in linkage disequilibrium higher than 0.10, 0.05, and 0.01 were considered, resulting in 15,286, 7,652, and 1,582 SNPs, respectively. For k = 2, the proportion of taurine and indicine varied from the expected proportion based on pedigree for all methods studied. For k = 3, the Regression method was able to differentiate the animals in three main clusters assigned to each purebred breed, showing more reasonable according to its biological viewpoint. Analyzing the data considering k = 2 seems to be more appropriate for Canchim-MA animals due to its biological interpretation. The usage of 32,308 SNPs in the analyses resulted in similar findings between the estimated and expected breed proportions. Using the Regression approach, a contribution of Indubrasil was observed in Canchim-MA when k = 3 was considered. Genetic parameter estimation could account for this breed composition information as a source of variation in order to improve the accuracy of genetic models. Our findings may help assemble appropriate reference populations for genomic prediction for Canchim-MA in order to improve prediction accuracy. Using the information on the level of introgression in each individual could also be useful in breeding or crossing design to improve individual heterosis in crossbred cattle.

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds.

Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.