During the last decade, microsatellites (short tandem repeats or STRs) have been successfully used for animal genetic identification, traceability and paternity, although in recent year single nucleotide polymorphisms (SNPs) have been increasingly used for this purpose. An efficient SNP identification system requires a marker set with enough power to identify individuals and their parents. Genetic diagnostics generally include the analysis of related animals. In this work, the degree of information provided by SNPs for a consanguineous herd of cattle was compared with that provided by STRs. Thirty-six closely related Angus cattle were genotyped for 18 STRs and 116 SNPs. Cumulative SNPs exclusion power values (Q) for paternity and sample matching probability (MP) yielded values greater than 0.9998 and 4.32E−42, respectively. Generally 2–3 SNPs per STR were needed to obtain an equivalent Q value. The MP showed that 24 SNPs were equivalent to the ISAG (International Society for Animal Genetics) minimal recommended set of 12 STRs (MP ∼ 10−11). These results provide valuable genetic data that support the consensus SNP panel for bovine genetic identification developed by the Parentage Recording Working Group of ICAR (International Committee for Animal Recording).

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome – or any mammalian transcriptome – appears significantly more feasible.

The analysis of runs of homozygosity (ROH), using high throughput genomic data, has become a valuable and frequently used methodology to characterize the genomic and inbreeding variation of livestock and wildlife animal populations. However, this methodology has been scarcely used in highly inbred domestic animals. Here, we analyzed and characterized the occurrence of ROH fragments in highly inbred (HI; average pedigree-based inbreeding coefficient FPED = 0.164; 0.103 to 0.306) and outbred Retinta bulls (LI; average FPED = 0.008; 0 to 0.025). We studied the length of the fragments, their abundance, and genome distribution using high-density microarray data. The number of ROH was significantly higher in the HI group, especially for long fragments (>8Mb). In the LI group, the number of ROH continuously decreased with fragment length. Genome-wide distribution of ROH was highly variable between samples. Some chromosomes presented a larger number of fragments (BTA1, BTA19, BTA29), others had longer fragments (BTA4, BTA12, BTA17), while other ones showed an increased ROH accumulation over specific loci (BTA2, BTA7, BTA23, BTA29). Similar differences were observed in the analysis of 12 individuals produced by a similar inbred event (FPED3 = 0.125). The correlation between the fraction of the genome covered by ROH (FROH) and FPED was high (0.79), suggesting that ROH-based estimations are indicative of inbreeding levels. On the other hand, the correlation between FPED and the microsatellite-based inbreeding coefficient (FMIC) was only moderate (r = 0.44), suggesting that STR-based inbreeding estimations should be avoided. Similarly, we found a very low correlation (r = -0.0132) between recombination rate and ROH abundance across the genome. Finally, we performed functional annotation analyses of genome regions with significantly enriched ROH abundance. Results revealed gene clusters related to pregnancy-associated proteins and immune reaction. The same analysis performed for regions enriched with recently formed ROH (> 8 Mb) showed gene clusters related to flagellum assembly. In both cases, the processes were related to male and female reproductive functions, which may partially explain the reduced fertility associated with inbred populations.

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F st ) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.

In vitro breeding constitutes a fast and intense selection strategy to genetically improve livestock populations.

No abstract

BackgroundIn bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Although the timing of puberty is likely to be a multigenic trait, previous studies indicate that there may also be single genes that exert major effects on the timing of puberty within the general population. Despite its economic importance, there are not many SNPs or genetic markers associated with the age of puberty in male cattle. In the present work, we selected three candidate genes, GNRHR, LHR and IGF1, and associated their polymorphisms with the age of puberty in Angus male cattle.ResultsAfter weaning, 276 Angus males were measured every month for weight (W), scrotal circumference (SC), sperm concentration (C) and percentage of motility (M). A total of 4 SNPs, two within GNRHR, one in LHR and one in IGF1 were genotyped using the pyrosequencing technique. IGF1-SnaBI SNP was significant associated (P < 0.01) with age at SC 28 cm, but it were not associated with age at M 10% and C 50 million. Genotype CC exhibited an average age at SC 28 cm of 7 and 11 days higher than CT (p = 0.037) and TT (p = 0.012), respectively. This SNP explained 1.5% of the genetic variance of age of puberty at SC28. LHR-I499L, GNRHR-SNP5 and GNRHR-SNP6 were not associated with any of the measurements. However, GNRHR haplotypes showed a suggestive association with age at SC 28 cm.ConclusionsThe findings presented here could support the hypothesis that IGF1 is a regulator of the arrival to puberty in male calves and is involved in the events that precede and initiate puberty in bull calves. Given that most studies in cattle, as well as in other mammals, were done in female, the present results are the first evidence of markers associated with age at puberty in male cattle.

FABP4 is a protein primarily expressed in adipocytes and macrophages that plays a key role in fatty acid trafficking and lipid hydrolysis. FABP4 gene polymorphisms have been associated with meat quality traits in cattle, mostly in Asian breeds under feedlot conditions. The objectives of this work were to characterize FABP4 genetic variation in several worldwide cattle breeds and evaluate possible genotype effects on fat content in a pasture-fed crossbred (Angus-Hereford-Limousin) population. We re-sequenced 43 unrelated animals from nine cattle breeds (Angus, Brahman, Creole, Hereford, Holstein, Limousin, Nelore, Shorthorn, and Wagyu) and obtained 22 single nucleotide polymorphisms (SNPs) over 3,164 bp, including four novel polymorphisms. Haplotypes and linkage disequilibrium analyses showed a high variability. Five SNPs were selected to perform validation and association studies in our crossbred population. Four SNPs showed well-balanced allele frequencies (minor frequency > 0.159), and three showed no significant deviations from Hardy-Weinberg proportions. SNPs showed significant effects on backfat thickness and fatty acid composition (P < 0.05). The protein structure of one of the missense SNPs was analyzed to elucidate its possible effect on fat content in our studied population. Our results revealed a possible blockage of the fatty acid binding site by the missense mutation.