Genomic regions associated with coat color in Gir cattle.

Indicine cattle breeds are adapted to the tropical climate, and their coat plays an important role in this process. Coat color influences thermoregulation and the adhesion of ectoparasites and may be associated with productive and reproductive traits. Furthermore, coat color is used for breed qualification, with breeders preferring certain colors. The Gir cattle is characterized by a wide variety of coat colors. Therefore, we performed genome-wide association studies to identify candidate genes for coat color in Gir cattle. Different phenotype scenarios were considered in the analyses and regions were identified on eight chromosomes. Some regions and many candidate genes are influencing coat color in the Gir cattle, which was found to be a polygenic trait. The candidate genes identified have been associated with white spotting patterns and base coat color in cattle and other species. In addition, a possible epistatic effect on coat color determination in the Gir cattle was suggested. This is the first published study that identified genomic regions and listed candidate genes associated with coat color in Gir cattle. The findings provided a better understanding of the genetic architecture of the trait in the breed and will allow to guide future fine-mapping studies for the development of genetic markers for selection.

Estimation of genetic parameters for reproductive indices in sheep.

This study aimed to estimate two reproductive efficiency indices in sheep based on the ratio between litter weight (at birth and weaning) and dam weight, as well as their genetic parameters. Phenotypic and pedigree data comprising the period from 1990 to 2018 were obtained from the Santa Inês sheep database of Embrapa Tabuleiros Costeiros. For estimation of the genetic parameters of the indices, a repeatability model was applied in single- and two-trait analyses by a Bayesian approach. The mean reproductive efficiency index was 0.069 ± 0.0163 and 0.43 ± 0.0955 at birth and weaning, respectively. These values indicate that, on average, ewes give birth to 69 g of lamb per kg body weight and wean 430 g of lamb per kg body weight. Described here for the first time, the heritability estimate obtained in single- and two-trait analyses was 0.24 for the index based on birth weights and ranged from 0.13 to 0.15 for the index based on weaning weights. The estimates indicate the possibility of genetic gain by selection and are similar to those reported for reproductive traits in sheep, representing an option for selection criterion. The genetic correlation between indices was positive and moderate (0.26). The repeatability estimates were high (0.49 for the birth weight index and 0.71 for the weaning weight index). These values indicate good prediction of future performance with few observations. The weaning weight index might be a good culling criterion of females.

Variability in the PRDM9 gene in Sindhi cattle.

BACKGROUND: Sindhi is a dual-purpose breed adapted to tropical environments. However, this breed has the smallest total population among indicine breeds in Brazil and the smallest effective number. In addition, the inbreeding coefficient is higher than 6.25% in ~ 60% of the population. Therefore, alternatives to increase genetic diversity are important. Within this context, the PRDM9 gene is particularly interesting since it is involved in meiotic recombination events, consequently enhancing genetic variability in the population by increasing the number of circulating haplotypes. Each allele of the gene induces recombination at a different hotspot. The larger the number of circulating alleles, the higher the recombination rate and the greater the genetic variability. METHODS: The aim of this study was to characterize alleles of the PRDM9 gene in Sindhi cattle. The region of the zinc finger domains of the gene was amplified by PCR, genotyped, and sequenced for allele identification in 50 Sindhi animals. RESULTS: Three alleles (A-cattle1, B-cattle14, and C-cattle19) and six genotypes (AA, BB, CC, AB, AC, and BC) were identified. CONCLUSION: The allele variation of the PRDM9 gene in the Sindhi breed enables to guide the mating of animals with different genotypes/alleles and to promote genetic variability by recombination if there is intralocus variability.

Guzerat indicine cattle and A2 milk production.

Cow milk might be associated to gastrointestinal disorders and abdominal pain in some people, which are, in part, due to the digestion of A1 beta-casein. Within this context, A2 milk has emerged as an alternative since it only contains A2 beta-casein that does not cause these complications. This milk is produced by cows with the A2A2 genotype. The aim of this study was to investigate the allele and genotype frequencies for the beta-casein gene in Guzerat cattle and to evaluate the feasibility of selection and production of A2 milk. The genotypes of 283 Guzerat cows from 10 herds were analyzed. The frequency of the A2A2 genotype was 0.80 and the frequency of the A2 allele was 0.90. These frequencies are slightly lower than those reported in previous studies involving populations of the same breed, but the number of animals, herein, genotyped was higher. Thus, the estimates are believed to be better and are also equally high. The Guzerat cattle has the potential for A2 milk production since most animals of the herd carry the favorable genotype for A2 milk production. Furthermore, the frequency of the A2A2 genotype can be rapidly increased by marker-assisted selection without compromising genetic variability.

Genomic association using principal components of morphometric traits in horses: identification of genes related to bone growth.

The measurement of morphometric traits in horses is important for determining breed qualification and is one of the main selection criteria for the species. The development of an index (HPC) that consists of principal components weighted by additive genetic values allows to explore the most relevant relationships using a reduced number of variables that explain the greatest amount of variation in the data. Genome-wide association studies (GWAS) using HPC are a relatively new approach that permits to identify regions related to a set of traits. The aim of this study was to perform GWAS using HPC for 15 linear measurements as the explanatory variable in order to identify associated genomic regions and to elucidate the biological mechanisms linked to this index in Campolina horses. For GWAS, weighted single-step GBLUP was applied to HPC. The eight genomic windows that explained the highest proportion of additive genetic variance were identified. The sum of the additive variance explained by the eight windows was 95.89%. Genes involved in bone and cartilage development were identified (SPRY2, COL9A2, MIR30C, HEYL, BMP8B, LTBP1, FAM98A, and CRIM1). They represent potential positional candidates for the HPC of the linear measurements evaluated. The HPC is an efficient alternative to reduce the 15 usually measured traits in Campolina horses. Moreover, candidate genes inserted in region that explained high additive variance of the HPC were identified and might be fine-mapped for searching putative mutation/markers.

Economic considerations of breeding for polledness versus disbudding in beef cattle.

The domestication of animals has rendered horns less necessary for survival. Moreover, the use of polled and disbudded animals is interesting in order to avoid injuries of animals and handlers, among other advantages. We therefore conducted a comparative economic analysis of different traditional disbudding techniques versus selective breeding for polledness in Nelore cattle, the main beef breed of tropical systems in Brazil. The cost to obtain animals without horns was estimated in three different scenarios: disbudding with hot iron, disbudding with caustic paste, and phenotypic selection for polled animals. Price quotations of the materials were obtained in different states of the country and averaged. An initial frequency of horned animals of 92.16% was obtained based on the records of the Brazilian Association of Zebu Breeders. Selective breeding was found to be the best cost-effective scenario. This result differs from intensive production systems of dairy cattle in which traditional disbudding continues to be the best cost-effective scenario. The main explanation is the lack of difference in the price of Nelore semen from polled and horned bulls. Phenotypic selection for polled animals is the best cost-effective method, and it is in accordance with welfare practices. Care should be taken regarding the intensive use of few polled breeding animals in order to avoid inbreeding depression in other traits.

AFFINITY OF BRAZILIAN WILD MAMMAL IMMUNOGLOBULINS TO BACTERIAL PROTEINS A AND G.

Staphylococcal A and streptococcal G proteins are widely used in immunoassays when specific immunological reagents are unavailable, such as for wild animals. The affinity of bacterial proteins A and G to the immunoglobulins of seven Brazilian mammals were tested, including black-tufted marmoset (Callithrix penicillata, n = 5), golden-bellied capuchin (Sapajus xanthosternos, n = 13), woolly mouse opossum (Micoureus demerarae, n = 6), long-nosed armadillo (Dasypus novemcinctus, n = 5), collared anteater (Tamandua tetradactyla, n = 5), ocelot (Leopardus pardalis, n = 6), and vampire bat (Desmodus rotundus, n = 5). Blood samples were collected from animals that were rescued in peri-urban rainforest fragments. Sera pools of each species were tested by ELISA to determine the intensity of each bacterial protein affinity to the immunoglobulins. When comparing the affinity to both proteins, immunoglobulins from D. rotundus, S. xanthosternos, and T. tetradactyla presented a higher affinity to protein G, whereas a higher affinity to protein A was found for immunoglobulins of C. penicillata and L. pardalis. The only species that presented a very low affinity to both bacterial proteins was M. demerarae. This study can be used as a reference for further studies on the development of sensitive and specific immunodiagnostic assays to be used for the monitoring of the health of these wild mammals.

Genetic origin of donkeys in Brazil.

INTRODUCTION: The genetic groups of native donkeys in Brazil are characterized by adaptation to the local environment. However, the donkey population in the country is declining, mainly because of agricultural mechanization and transportation that has led to the abandonment and the consequent indiscriminate slaughter of these animals. There are three local genetic groups of distinct geographic and temporal formation. However, analyses of their origin, phylogenetic relationship, and population structure are scarce. Within this context, molecular markers such as the mitochondrial control region (D-loop) are useful for these analyses. MATERIAL AND METHODS: This study aimed to evaluate the variation and origin of maternal lineages of groups of naturalized donkeys in Brazil (Brazilian, Nordestino, and Pêga). RESULTS AND DISCUSSION: We detected five mitochondrial haplotypes with 19 polymorphic sites, two of them exclusively found in the Nordestino donkey; this group is in fact more distant from the others. Phylogenetic analysis indicates maternal contributions of two clades (Nubian and Somali) to the formation of the genetic groups of donkeys, a fact that explains the high diversity, structure, and distances of the groups, reported here for the first time. CONCLUSION : This analysis contributes production and conservation of native donkey breeds. It also gives clues about the formation of the Iberian breeds from which Brazilian donkeys originated.

Selection effect for growth traits on energy requirements in beef Nelore steers.

Growth data of 77,372 Nelore steers were used to estimate the selection effect on energy requirements considering two beef production systems: cow-calf and slaughter cycles. All the animals had measures from 120 days to 7 years old. The parameters necessary to evaluate the selection effect on energy requirements were obtained by random regression analysis using Legendre polynomials. The models included additive direct and maternal effects, and animal and maternal permanent environmental effects as random terms. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of animal age (cubic regression) were considered random covariables. The coefficients from the model M3353_5 were used to calculate the genetic gains necessary to predict the increase in phenotypes. The selection was simulated for body weight (BW) and weight gain (WG) at different ages and energy requirements were calculated using NRC equations. The cost of feed was calculated for a cow-calf and slaughter cycle of production considering a system of Brachiaria decumbens pasture without supplementation. In slaughter system, the selection for weight of 365 days of age is the best option. In cow-calf systems, the selection W120 is the best choice.

The effect of mitochondrial DNA polymorphisms on cattle reproduction.

The aim of this study was to identify SNPs located in mitochondrial DNA that are associated with reproductive traits in beef cows. A total of 1999 Nelore females genotyped with the high-density Illumina BovineHD BeadChip (Illumina Inc., San Diego, CA, USA) were used to study the association of mitochondrial DNA variants with reproductive traits using a single-step procedure. In a preliminary analysis, the present results indicate a small participation of the mitogenome in the expression of reproductive traits in beef cattle. However, possible difficulties related to the biological characteristics of mitochondrial DNA and its inheritance, genotyping, and annotation of the phenotypes studied may also explain the results.

Comparison of methods for predicting genomic breeding values for growth traits in Nellore cattle.

The objective of this study was to evaluate the accuracy of genomic predictions of growth traits in Nellore cattle. Data from 5064 animals belonging to farms that participate in the Conexão DeltaGen and PAINT breeding programs were used. Genotyping was performed with the Illumina BovineHD BeadChip (777,962 SNPs). After quality control of the genomic data, 412,993 SNPs were used. Deregressed EBVs (DEBVs) were calculated using the estimated breeding values (EBVs) and accuracies of birth weight (BW), weight gain from birth to weaning (GBW), postweaning weight gain (PWG), yearling height (YH), and cow weight (CW) provided by GenSys. Three models were used to estimate marker effects: genomic best linear unbiased prediction (GBLUP), BayesCπ, and improved Bayesian least absolute shrinkage and selection operator (IBLASSO). The prediction ability of genomic estimated breeding value (GEBVs) was estimated by the average Pearson correlation between DEBVs and GEBVs, predicted with the different methodologies in the validation populations. The regression coefficients of DEBVs on GEBVs in the validation population were calculated and used as indicators of prediction bias of GEBV. In general, the Bayesian methods provided slightly more accurate predictions of genomic breeding values than GBLUP. The BayesCπ and IBLASSO were similar for all traits (BW, GBW, PWG, and YH), except for CW. Thus, there does not seem to be a more suitable method for the estimation of SNP effects and genomic breeding values. Bayesian regression models are of interest for future applications of genomic selection in this population, but further improvements are needed to reduce deflation of their predictions.

Water buffaloes (Bubalus bubalis) only have A2A2 genotype for beta-casein.

Beta-casein is a milk protein that has two variants: A1 and A2. Some individuals have difficulties digesting beta-casein A1, which can cause gastrointestinal disorders. A2 milk has emerged as an alternative. This milk only contains beta-casein A2 and is obtained from females carrying the A2A2 genotype of the gene. In cattle, allele and genotype frequencies vary according to breed and marker-assisted selection is performed to obtain A2A2 animals and the consequent production of A2 milk that is easier to digest. This study aimed to evaluate the alleles of beta-casein in buffaloes. A total of 657 buffaloes of four different breeds were genotyped and all animals carried the A2A2 genotype, i.e., allele A1 does not exist in the buffalo species. Thus, all milk products of buffaloes are naturally A2. This result adds value to products derived from buffalo milk.

Investigating the accuracy of imputing autosomal variants in Nellore cattle using the ARS-UCD1.2 assembly of the bovine genome.

BACKGROUND: Imputation accuracy among other things depends on the size of the reference panel, the marker's minor allele frequency (MAF), and the correct placement of single nucleotide polymorphism (SNP) on the reference genome assembly. Using high-density genotypes of 3938 Nellore cattle from Brazil, we investigated the accuracy of imputation from 50 K to 777 K SNP density using Minimac3, when map positions were determined according to the bovine genome assemblies UMD3.1 and ARS-UCD1.2. We assessed the effect of reference and target panel sizes on the pre-phasing based imputation quality using ten-fold cross-validation. Further, we compared the reliability of the model-based imputation quality score (Rsq) from Minimac3 to the empirical imputation accuracy. RESULTS: The overall accuracy of imputation measured as the squared correlation between true and imputed allele dosages (R2dose) was almost identical using either the UMD3.1 or ARS-UCD1.2 genome assembly. When the size of the reference panel increased from 250 to 2000, R2dose increased from 0.845 to 0.917, and the number of polymorphic markers in the imputed data set increased from 586,701 to 618,660. Advantages in both accuracy and marker density were also observed when larger target panels were imputed, likely resulting from more accurate haplotype inference. Imputation accuracy increased from 0.903 to 0.913, and the marker density in the imputed data increased from 593,239 to 595,570 when haplotypes were inferred in 500 and 2900 target animals. The model-based imputation quality scores from Minimac3 (Rsq) were systematically higher than empirically estimated accuracies. However, both metrics were positively correlated and the correlation increased with the size of the reference panel and MAF of imputed variants. CONCLUSIONS: Accurate imputation of BovineHD BeadChip markers is possible in Nellore cattle using the new bovine reference genome assembly ARS-UCD1.2. The use of large reference and target panels improves the accuracy of the imputed genotypes and provides genotypes for more markers segregating at low frequency for downstream genomic analyses. The model-based imputation quality score from Minimac3 (Rsq) can be used to detect poorly imputed variants but its reliability depends on the size of the reference panel and MAF of the imputed variants.

Mapping genomic regions for reproductive traits in beef cattle: Inclusion of the X chromosome.

Although the second largest chromosome of the genome, the X chromosome is usually excluded from genome-wide association studies (GWAS). Considering the presence and importance of genes on this chromosome that are involved in reproduction, the aim of this study was to evaluate the effect of its inclusion in GWAS on reproductive traits (scrotal circumference [SC], early pregnancy [P16] and age at first calving [AFC]) in a Nelore herd. Genotype data from 3,263 animals with the above-mentioned phenotypes were used. The results showed an increase in the variances explained by the autosomal markers for all traits when the X chromosome was not included. For SC, there was an increase of more than 10% for the windows on chromosomes 2 and 6. For P16, the effect was increased by almost 20% for windows on chromosome 5. The same pattern was found for AFC, with an increase of more than 10% for the most important windows. The results indicate that the noninclusion of the X chromosome can overestimate the effects of autosomes on SC, P16 and AFC not only because of the additive effect of the X chromosome itself but also because of its epistatic effect on autosomal genes.

Inconsistencies in horse coat color registration: A case study.

Grullo is a dun dilution on a black coat that is common in the Campolina horse: an autochthonous Brazilian breed. The aims of this case study were to evaluate inconsistencies in grullo coat color registration and to explain their possible causes. A total of 3,270 grullo Campolina horses were evaluated. To confirm the genetic possibility of having grullo animals, the coat color genotypes of parents were inferred by phenotype and compared with those of progeny. A total of 242 horses that were registered as grullos could not have this coat based on their parents' information. Possible explanations for incorrect registration are errors of paternity and in coat color identification. We suggest maintaining obligatory paternity testing and enhancing training in coat color identification.

Twinning rate in buffaloes: A case report.

The information of twin pregnancy frequencies in buffaloes is scarce in the literature. Evidence even indicates that twin pregnancies are not carried to term, but it does not seem to be the case. In this case report, the event of twin pregnancy in buffaloes is reported as well as its frequency. Twin birth is the presence of two or more embryos in a single gestation. Multiple births are uncommon in buffaloes since the species is uniparous. We studied a population of Murrah buffaloes in which slightly more than 13,000 births were observed. In this study, the frequency of occurrence of twin births in Murrah buffaloes was 0.14%. This low incidence impaired the estimation of genetic parameters. Twins can have advantages and disadvantages depending on the production system, which are discussed in this case report.

Consumption of rabbit meat in Brazil: Potential and limitations.

Rabbit production is still small in Brazil and organization of its production chain is therefore necessary. The aim of this study was to investigate consumer perceptions of rabbit meat in the country. An online survey was conducted with 2614 participants. The responses were analyzed using the nonparametric Mann-Whitney U test. The survey revealed that the greatest acceptance of rabbit meat occurred in the population of adult men with a higher income and educational level who live in capital cities. Regional differences exist in Brazil. Consumption is more consolidated in the south region. The main limiting factor for rabbit meat consumption is the unavailability of the product. The limited availability of rabbit meat results in low consumption. Producers ensuring a regular supply of rabbit meat and verticalization by the industry are essential. These proposed actions were found to be even more important than demystify the image of the rabbit as a pet.