Background Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. Results We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. Conclusions Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Navel injuries caused by friction against the pasture can promote infection, reproductive problems and costly treatments in beef cattle raised in extensive systems. A haplotype-based genome-wide association study (GWAS) was performed for visual scores of navel length at yearling in Nellore cattle (Bos indicus) using data from 2,016 animals and 503,088 single nucleotide polymorphism (SNP) markers. The strongest signal (p = 1.01 × 10-9) was found on chromosome 5 spanning positions 47.9–48.2 Mbp. This region contains introns 3 and 4 and exons 4 and 5 of the high mobility group AT-hook 2 gene (HMGA2). Further inspection of the region with whole genome sequence data of 21 Nellore bulls revealed correlations between counts of the significant haplotype and copy number gains of a ∼6.2 kbp segment of intron 3 of HMGA2. Analysis of genome sequences from five African B. indicus and four European Bos taurus breeds revealed that the copy number variant (CNV) is indicine-specific. This intronic CNV was then validated through quantitative polymerase chain reaction (qPCR) using Angus animals as copy neutral controls. Importantly, the CNV was not detectable by means of conventional SNP-based GWAS or SNP probe intensity analyses. Given that HMGA2 affects the expression of the insulin-like growth factor 2 gene (IGF2) together with the pleomorphic adenoma gene 1 (PLAG1), and that the latter has been repeatedly shown to be associated with quantitative traits of economic importance in cattle, these findings highlight the emerging role of variants impacting the insulin-like growth factor pathway to cattle breeding.

No abstract

SummaryIn the past two decades, average litter size (ALS) in Entlebucher Mountain dogs decreased by approximately 0.8 puppies. We conducted a GWAS for ALS using the single‐step methodology to take advantage of 1632 pedigree records, 892 phenotypes and 372 genotypes (173 662 markers) for which only 12% of the dogs had both phenotypes and genotypes available. Our analysis revealed associations towards the growth differentiation factor 9 gene (GDF9), which is known to regulate oocyte maturation. The trait heritability was estimated at 43.1%, from which approximately 15% was accountable by the GDF9 locus alone. Therefore, markers flanking GDF9 explained approximately 6.5% of the variance in ALS. Analysis of WGSs revealed two missense substitutions in GDF9, one of which (g.11:21147009G>A) affected a highly conserved nucleotide in vertebrates. The derived allele A was validated in 111 dogs and shown to be associated with decreased ALS (−0.75 ± 0.22 puppies per litter). The variant was further predicted to cause a proline to serine substitution. The affected residue was immediately followed by a six‐residue deletion that is fixed in the canine species but absent in non‐canids. We further confirmed that the deletion is prevalent in the Canidae family by sequencing three species of wild canids. Since canids uniquely ovulate oocytes at the prophase stage of the first meiotic division, requiring maturation in the oviduct, we conjecture that the amino acid substitution and the six‐residue deletion of GDF9 may serve as a model for insights into the dynamics of oocyte maturation in canids.

Brazil is the largest exporter of beef in the world, and most of that beef derives from Nellore cattle. Although considered a zebu breed (Bos indicus), the history of Nellore cattle in Brazil is marked by the importation of bulls from India, the use of a Creole taurine (Bos taurus) maternal lineage to quickly expand the herds and backcrossing to Nellore bulls to recover zebu ancestry. As a consequence, the current Brazilian Nellore population carries an average taurine ancestry of approximately 1%. Although that percentage seems small, some taurine variants deviate substantially from that average, with the better-known cases being the PLAG1-Q haplotype involved with body size variation and the Guarani (P G ) polled variant producing hornless animals. Here, we report taurine haplotypes in 9074 Nellore animals genotyped for 539 657 imputed SNP markers. Apart from PLAG1-Q and P G , our analysis further revealed common taurine haplotypes (>3%) spanning genes related to immunity, growth, reproduction and hair and skin phenotypes. Using data from 22 economically important traits, we showed that many of the major QTL previously reported in the breed are at least partially driven by taurine haplotypes. As B. taurus and B. indicus haplotypes are highly divergent, presenting widely different sets of functional variants, our results provide promising targets for future scrutiny in Nellore cattle.