Do non-bovine domestic animals produce A2 milk?: an in silico analysis.

A2 milk is an easily digestible product since it has only A2 beta-casein. In cattle, the A1 and A2 alleles are found in the population and the A2 milk is produced from A2A2 animals. Little is known about these alleles in other domestic dairy species. The present study aims to analyze sequence of genetic material available on public databases and quantify the animals genotyped. Eight domestic species were analyzed. There is strong evidence that domestic non-bovine species only carry A2 beta-casein. The data reported here for goats already confirm it due to the large number of animals genotyped as well as buffaloes. It means that they naturally produce A2 milk and no selection must be done. Thus, the fact that A2 milk is easier to digest can be used to add value to dairy product of these species. It helps to conquer new markets. It also improves people's health and breeder profitability.

Kappa and beta-casein alleles in dairy zebu cattle.

The kappa (CSN3) and beta-casein (CSN2) genes are intensively genotyped in dairy cattle for selection purposes. This information is also generated and disseminated for Zebu breeds adapted to tropical climates. The objective of this work was to gather information on the genotypes for the CSN3 and CSN2 genes in three breeds (Gyr, Guzerat and Sindhi), and to verify the genotypic frequencies in the populations. The genotype AA and allele A frequencies are high for the CSN3 gene, without changes in values over the years, possibly indicating a small gene participation in traits under selection. In addition, the A2A2 frequencies are high for the CSN2 gene (<∼0.80). It is recommended to verify the association and contribution of CSN3 genotypes in productive traits for these breeds. The potential of A2 milk production by these genetic groups is also confirmed.

Genome-wide association between single nucleotide polymorphisms with beef fatty acid profile in Nellore cattle using the single step procedure.

BACKGROUND: Saturated fatty acids can be detrimental to human health and have received considerable attention in recent years. Several studies using taurine breeds showed the existence of genetic variability and thus the possibility of genetic improvement of the fatty acid profile in beef. This study identified the regions of the genome associated with saturated, mono- and polyunsaturated fatty acids, and n-6 to n-3 ratios in the Longissimus thoracis of Nellore finished in feedlot, using the single-step method. RESULTS: The results showed that 115 windows explain more than 1 % of the additive genetic variance for the 22 studied fatty acids. Thirty-one genomic regions that explain more than 1 % of the additive genetic variance were observed for total saturated fatty acids, C12:0, C14:0, C16:0 and C18:0. Nineteen genomic regions, distributed in sixteen different chromosomes accounted for more than 1 % of the additive genetic variance for the monounsaturated fatty acids, such as the sum of monounsaturated fatty acids, C14:1 cis-9, C18:1 trans-11, C18:1 cis-9, and C18:1 trans-9. Forty genomic regions explained more than 1 % of the additive variance for the polyunsaturated fatty acids group, which are related to the total polyunsaturated fatty acids, C20:4 n-6, C18:2 cis-9 cis12 n-6, C18:3 n-3, C18:3 n-6, C22:6 n-3 and C20:3 n-6 cis-8 cis-11 cis-14. Twenty-one genomic regions accounted for more than 1 % of the genetic variance for the group of omega-3, omega-6 and the n-6:n-3 ratio. CONCLUSIONS: The identification of such regions and the respective candidate genes, such as ELOVL5, ESSRG, PCYT1A and genes of the ABC group (ABC5, ABC6 and ABC10), should contribute to form a genetic basis of the fatty acid profile of Nellore (Bos indicus) beef, contributing to better selection of the traits associated with improving human health.

Genome-wide association study of reproductive traits in Nellore heifers using Bayesian inference.

BACKGROUND: An important goal of Zebu breeding programs is to improve reproductive performance. A major problem faced with the genetic improvement of reproductive traits is that recording the time for an animal to reach sexual maturity is costly. Another issue is that accurate estimates of breeding values are obtained only a long time after the young bulls have gone through selection. An alternative to overcome these problems is to use traits that are indicators of the reproductive efficiency of the herd and are easier to measure, such as age at first calving. Another problem is that heifers that have conceived once may fail to conceive in the next breeding season, which increases production costs. Thus, increasing heifer's rebreeding rates should improve the economic efficiency of the herd. Response to selection for these traits tends to be slow, since they have a low heritability and phenotypic information is provided only later in the life of the animal. Genome-wide association studies (GWAS) are useful to investigate the genetic mechanisms that underlie these traits by identifying the genes and metabolic pathways involved. RESULTS: Data from 1853 females belonging to the Agricultural Jacarezinho LTDA were used. Genotyping was performed using the BovineHD BeadChip (777 962 single nucleotide polymorphisms (SNPs)) according to the protocol of Illumina - Infinium Assay II ® Multi-Sample HiScan with the unit SQ ™ System. After quality control, 305 348 SNPs were used for GWAS. Forty-two and 19 SNPs had a Bayes factor greater than 150 for heifer rebreeding and age at first calving, respectively. All significant SNPs for age at first calving were significant for heifer rebreeding. These 42 SNPs were next or within 35 genes that were distributed over 18 chromosomes and comprised 27 protein-encoding genes, six pseudogenes and two miscellaneous noncoding RNAs. CONCLUSIONS: The use of Bayes factor to determine the significance of SNPs allowed us to identify two sets of 42 and 19 significant SNPs for heifer rebreeding and age at first calving, respectively, which explain 11.35 % and 6.42 % of their phenotypic variance, respectively. These SNPs provide relevant information to help elucidate which genes affect these traits.

Polymorphisms in TOX and NCOA2 genes and their associations with reproductive traits in cattle.

Reproductive traits are an important component of the economic selection index for beef cattle in the tropics. Phenotypic expression of these traits occurs late because they are measured when the animals reach reproductive age. Association studies using high-density markers have been conducted to identify genes that influence certain traits. The identification of causal mutations in these genes permits the inclusion of these single nucleotide polymorphisms (SNPs) in customised DNA chips to increase efficiency and validity. Therefore, the aim of the present study was to detect causal mutations in the TOX and NCOA2 genes, previously identified by genome-wide association studies of zebu cattle. DNA was extracted from 385 Nellore females and polymorphisms were investigated by polymerase chain reaction sequencing. Five polymorphisms were detected in the NCOA2 gene and four in the TOX gene that were associated with reproductive traits. Analysis of variance showed that SNP 1718 in the NCOA2 gene was significant for early pregnancy probability (P=0.02) and age at first calving (P=0.03), and SNP 2038 in the same gene was significant for days to calving (P=0.03). Studies investigating polymorphisms in other regions of the gene and in other genes should be conducted to identify causal mutations.

Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes.

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal.

Milk fatty acid characterization and genetic parameter estimates for milk conjugated linoleic acid in buffaloes.

The objectives of this study were to analyse buffalo milk fat composition, to verify the activity of Delta(9)-desaturase enzyme in the mammary gland, as well as to estimate additive genetic variances for milk, fat and protein yield, and milk cis-9,trans-11 conjugated linoleic acid percentage (cis-9,trans-11 CLA%). A total of 3929 lactation milk yields (MY) records from 2130 buffaloes and 1598 lactation fat (FY) and protein (PY) yield records from 914 buffaloes were analysed. For cis-9,trans-11 CLA%percentage, a total of 661 milk samples from 225 buffaloes, daughters of 8 sires, belonging to 4 herds and calving in 2003 and 2004, were used. The genetic parameters and variance components were estimated by Restricted Maximum Likelihood applying an animal model. The fixed effects considered in the model were: contemporary group (herd, year, calving season) and age at calving (linear and quadratic effects) and lactation length (linear and quadratic effects) as covariables. Additive genetic and permanent environment effects were considered as random. The MY, FY, PY and CLA% means were 1482±355 kg, 90.1±24.6 kg, 56.9±15.2 kg and 0.69±0.16%, respectively. Heritability estimates for MY, FY, PY and CLA% were 0.28±0.05, 0.26±0.11, 0.25±0.11 and 0.35±0.14, respectively. There is enough additive genetic variation for buffalo milk, protein and fat yield to improve these traits through selection. The cis-9,trans-11 CLA% can be enhanced by selection in buffaloes and will contribute to improving human health. The activity and efficiency of Delta(9)-desaturase in the mammary was measured and confirmed.

Genetic parameters for stayability in Murrah buffaloes.

In order to contribute to the breeding programmes of Asian water buffalo, the aim of this study was to analyse the influence of genetic effects in the stayability of Murrah dairy buffaloes. The stayability trait (ST) was defined as the female's ability to stay in the herd for one (ST1), two (ST2), three (ST3), four (ST4), five (ST5) or six years (ST6) after the first calving. The same trait was also considered as continuous and was designated stayability in days up to one (STD1), two (STD2), three (STD3), four (STD4), five (STD5) or six years (STD6) after the first calving. Data from 1016 females reared in nine herds located in the State of São Paulo, Brazil, were analysed. Statistical models included the additive genetic effect of the animal and the fixed effects of the buffalo breeding herd, birth year and birth season. Additive effects for ST were estimated by approximate restricted maximum likelihood using a threshold model, while for STD, the additive effects were estimated by restricted maximum likelihood. Heritability estimates were lower for ST, except for ST1, (0.11+/-0.07, 0.17+/-0.06, 0.23+/-0.06, 0.16+/-0.08, 0.14+/-0.09 and 0.16+/-0.10 for ST1, ST2, ST3, ST4, ST5 and ST6, respectively) when compared with STD (0.05+/-0.06, 0.18+/-0.08, 0.40+/-0.10, 0.49+/-0.11, 0.41+/-0.11 and 0.30+/-0.13, for STD1, STD2, STD3, STD4, STD5 and STD6, respectively). Considering the values of heritability and owing to the serial nature of STD to a specific age, selection for STD3 should have a favourable influence on STD to other ages.

Genomic association for sexual precocity in beef heifers using pre-selection of genes and haplotype reconstruction.

Reproductive traits are of the utmost importance for any livestock farming, but are difficult to measure and to interpret since they are influenced by various factors. The objective of this study was to detect associations between known polymorphisms in candidate genes related to sexual precocity in Nellore heifers, which could be used in breeding programs. Records of 1,689 precocious and non-precocious heifers from farms participating in the Conexão Delta G breeding program were analyzed. A subset of single nucleotide polymorphisms (SNP) located in the region of the candidate genes at a distance of up to 5 kb from the boundaries of each gene, were selected from the panel of 777,000 SNPs of the High-Density Bovine SNP BeadChip. Linear mixed models were used for statistical analysis of early heifer pregnancy, relating the trait with isolated SNPs or with haplotype groups. The model included the contemporary group (year and month of birth) as fixed effect and parent of the animal (sire effect) as random effect. The fastPHASE® and GenomeStudio® were used for reconstruction of the haplotypes and for analysis of linkage disequilibrium based on r2 statistics. A total of 125 candidate genes and 2,024 SNPs forming haplotypes were analyzed. Statistical analysis after Bonferroni correction showed that nine haplotypes exerted a significant effect (p<0.05) on sexual precocity. Four of these haplotypes were located in the Pregnancy-associated plasma protein-A2 gene (PAPP-A2), two in the Estrogen-related receptor gamma gene (ESRRG), and one each in the Pregnancy-associated plasma protein-A gene (PAPP-A), Kell blood group complex subunit-related family (XKR4) and mannose-binding lectin genes (MBL-1) genes. Although the present results indicate that the PAPP-A2, PAPP-A, XKR4, MBL-1 and ESRRG genes influence sexual precocity in Nellore heifers, further studies are needed to evaluate their possible use in breeding programs.

Genome-Wide Association Study of Meat Quality Traits in Nellore Cattle.

The objective of this study was to identify genomic regions that are associated with meat quality traits in the Nellore breed. Nellore steers were finished in feedlots and slaughtered at a commercial slaughterhouse. This analysis included 1,822 phenotypic records of tenderness and 1,873 marbling records. After quality control, 1,630 animals genotyped for tenderness, 1,633 animals genotyped for marbling, and 369,722 SNPs remained. The results are reported as the proportion of variance explained by windows of 150 adjacent SNPs. Only windows with largest effects were considered. The genomic regions were located on chromosomes 5, 15, 16 and 25 for marbling and on chromosomes 5, 7, 10, 14 and 21 for tenderness. These windows explained 3,89% and 3,80% of the additive genetic variance for marbling and tenderness, respectively. The genes associated with the traits are related to growth, muscle development and lipid metabolism. The study of these genes in Nellore cattle is the first step in the identification of causal mutations that will contribute to the genetic evaluation of the breed.

Low frequency of Y anomaly detected in Australian Brahman cow-herds.

Indicine cattle have lower reproductive performance in comparison to taurine. A chromosomal anomaly characterized by the presence Y markers in females was reported and associated with infertility in cattle. The aim of this study was to investigate the occurrence of the anomaly in Brahman cows. Brahman cows (n = 929) were genotyped for a Y chromosome specific region using real time-PCR. Only six out of 929 cows had the anomaly (0.6%). The anomaly frequency was much lower in Brahman cows than in the crossbred population, in which it was first detected. It also seems that the anomaly doesn't affect pregnancy in the population. Due to the low frequency, association analyses couldn't be executed. Further, SNP signal of the pseudoautosomal boundary region of the Y chromosome was investigated using HD SNP chip. Pooled DNA of "non-pregnant" and "pregnant" cows were compared and no difference in SNP allele frequency was observed. Results suggest that the anomaly had a very low frequency in this Australian Brahman population and had no effect on reproduction. Further studies comparing pregnant cows and cows that failed to conceive should be executed after better assembly and annotation of the Y chromosome in cattle.

Association of ADIPOQ, OLR1 and PPARGC1A gene polymorphisms with growth and carcass traits in Nelore cattle.

In beef cattle farming, growth and carcass traits are important for genetic breeding programs. Molecular markers can be used to assist selection and increase genetic gain. The ADIPOQ, OLR1 and PPARGC1A genes are involved in lipid synthesis and fat accumulation in adipose tissue. The objective of this study was to identify polymorphisms in these genes and to assess the association with growth and carcass traits in Nelore cattle. A total of 639 animals were genotyped by PCR-RFLP for rs208549452, rs109019599 and rs109163366 in ADIPOQ, OLR1 and PPARGC1A gene, respectively. We analyzed the association of SNPs identified with birth weight, weaning weight, female yearling weight, female hip height, male yearling weight, male hip height, loin eye area, rump fat thickness, and backfat thickness. The OLR1 marker was associated with rump fat thickness and weaning weight (P < 0.05) and the PPARGC1 marker was associated with female yearling weight.