Weak genotype x environment interaction suggests that measuring scrotal circumference at 12 and 18 mo of age is helpful to select precocious Brahman cattle.

The aim of this study was to evaluate the genotype x environment interaction (GxE) for scrotal circumference (SC) measured at different ages using pedigree-based (A-1) and pedigree and genomic-based (H-1) relationship matrices. Data from 1,515 Brahman bulls, from the Cooperative Research Centre for Beef Genetic Technologies (Beef CRC) experimental dataset, were used in this study. SC was adjusted to age and body weight measured at 6 mo (SC6), 12 mo (SC12), 18 mo (SC18), and 24 mo of age (SC24). Body weight (BW) measured at 6 mo (BW6), 12 mo (BW12), 18 mo (BW18), and 24 mo of age (BW24) was used as criteria to describe the environment for SC in each age. All the animals measured were genotyped using medium-density SNP chips ("50k" or "70k" SNP) and their genotype were imputed using a reference panel with 729,068 SNP. The environment gradient (EG) was obtained by standardizing the solutions of the contemporary groups obtained by Animal Model with BW as the dependent variable. Then, the reaction norms (RN) were determined through a Random Regression Model. The breeding values (EBV) were estimated using either A-1 or H-1. The rank correlation was obtained using Spearman's correlation among the EBV estimated for the traits in analysis. For SC6 and SC24, higher estimates of heritability (h²) were obtained using A-1, when compared with those observed with H-1. In those ages, the improvement of the environment decreases the h² coefficient. On the other hand, the h² for SC12 and SC18 increased as the environment became more favorable, regardless of the matrix used. The RN for SC6 and SC24 estimated using A-1 and H-1 showed a decrease of variance from the worst to the best environment, an indication of existence of GxE. On the other hand, for SC12 and SC18, there were no significant differences between the EBV estimated in the lower and in the higher environments, regardless of the kinship matrix used, suggesting absence of GxE on those ages. Spearman's correlation among EBV estimated using A-1 and H-1 in different EG was practically equal to unity for all traits evaluated. In our study, there was weak evidence of GxE effect on SC in ages suitable for selection for sexual precocity. So, the absence of GxE at 12 and 18 mo means that these ages are advantageous for measuring SC to selection for sexual precocity. The advantage is that no changes in classification were observed when the sires were evaluated in different environments.

Beef production systems rely on efficient cow-calf operations, that is, farms where the cow herd has a high level of fertility and pregnancies are common. Bull fertility also plays an important role in terms of pregnancy rates. To increase herd fertility, cattle breeders and genetic selection programs use some indicator traits that are related to fertility. A common indicator trait used is scrotal circumference (SC), which like most reproduction traits are influenced by the animal's genetics and its environment. For some traits, when the environment has a large effect and it interacts with the genetics of the animals, selection might be less successful. Therefore, it is important to investigate genotype by environment interactions and their effect on reproduction traits used for selection. In this study, SC was measured at four different ages in Brahman cattle. We found weak evidence of genotype by environment effect on SC measured at 12 and 18 mo. In short, SC measured at these ages can be a good indicator of sexual precocity. No changes in sire rankings were observed when SC was measured at those ages, meaning that selecting the best sire is more straightforward than if the environment was playing a bigger role.

Network Analyses Predict Small RNAs That Might Modulate Gene Expression in the Testis and Epididymis of Bos indicus Bulls.

Spermatogenesis relies on complex molecular mechanisms, essential for the genesis and differentiation of the male gamete. Germ cell differentiation starts at the testicular parenchyma and finishes in the epididymis, which has three main regions: head, body, and tail. RNA-sequencing data of the testicular parenchyma (TP), head epididymis (HE), and tail epididymis (TE) from four bulls (three biopsies per bull: 12 samples) were subjected to differential expression analyses, functional enrichment analyses, and co-expression analyses. The aim was to investigate the co-expression and infer possible regulatory roles for transcripts involved in the spermatogenesis of Bos indicus bulls. Across the three pairwise comparisons, 3,826 differentially expressed (DE) transcripts were identified, of which 384 are small RNAs. Functional enrichment analysis pointed to gene ontology (GO) terms related to ion channel activity, detoxification of copper, neuroactive receptors, and spermatogenesis. Using the regulatory impact factor (RIF) algorithm, we detected 70 DE small RNAs likely to regulate the DE transcripts considering all pairwise comparisons among tissues. The pattern of small RNA co-expression suggested that these elements are involved in spermatogenesis regulation. The 3,826 DE transcripts (mRNAs and small RNAs) were further subjected to co-expression analyses using the partial correlation and information theory (PCIT) algorithm for network prediction. Significant correlations underpinned the co-expression network, which had 2,216 transcripts connected by 158,807 predicted interactions. The larger network cluster was enriched for male gamete generation and had 15 miRNAs with significant RIF. The miRNA bta-mir-2886 showed the highest number of connections (601) and was predicted to down-regulate ELOVL3, FEZF2, and HOXA13 (negative co-expression correlations and confirmed with TargetScan). In short, we suggest that bta-mir-2886 and other small RNAs might modulate gene expression in the testis and epididymis, in Bos indicus cattle.

Polymorphisms in major histocompatibility complex genes and its associations with milk quality in Murrah buffaloes.

Animal breeding programs have used molecular genetic tools as an auxiliary method to identify and select animals with superior genetic merit for milk production and milk quality traits as well as disease resistance. Genes of the major histocompatibility complex (MHC) are important molecular markers for disease resistance that could be applied for genetic selection. The aim of this study was to identify single nucleotide polymorphisms (SNPs) and haplotypes in DRB2, DRB3, DMA, and DMB genes in Murrah breed and to analyze the association between molecular markers and milk, fat, protein and mozzarella production, fat and protein percentage, and somatic cell count. Two hundred DNA samples from Murrah buffaloes were used. The target regions of candidate genes were amplified by polymerase chain reaction (PCR) followed by sequencing and identification of polymorphisms. Allele and genotype frequencies, as well as linkage disequilibrium between SNPs, were calculated. Genotypes were used in association analyses with milk production and quality traits. Except for the DMA gene, identified as monomorphic, the other genes presented several polymorphisms. The DMB, DRB2, and DRB3 genes presented two, six, and seven SNPs, respectively. Fifty-seven haplotype blocks were constructed from 15 SNPs identified, which was used in association analyses. All the studied traits had at least one associated haplotype. In conclusion, it is suggested that the haplotypes found herein can be associated with important traits related to milk production and quality.

RAPID COMMUNICATION: Multi-breed validation study unraveled genomic regions associated with puberty traits segregating across tropically adapted breeds1.

An efficient strategy to improve QTL detection power is performing across-breed validation studies. Variants segregating across breeds are expected to be in high linkage disequilibrium (LD) with causal mutations affecting economically important traits. The aim of this study was to validate, in a Tropical Composite cattle (TC) population, QTL associations identified for sexual precocity traits in a Nellore and Brahman meta-analysis genome-wide association study. In total, 2,816 TC, 8,001 Nellore, and 2,210 Brahman animals were available for the analysis. For that, genomic regions significantly associated with puberty traits in the meta-analysis study were validated for the following sexual precocity traits in TC: age at first corpus luteum (AGECL), first postpartum anestrus interval (PPAI), and scrotal circumference at 18 months of age (SC). We considered validated QTL those underpinned by significant markers from the Nellore and Brahman meta-analysis (P ≤ 10-4) that were also significant for a TC trait, i.e., presenting a P-value of ≤10-3 for AGECL, PPAI, or SC. We also considered as validated QTL those regions where significant markers in the reference population were at ±250 kb from significant markers in the validation population. Using this criteria, 49 SNP were validated for AGECL, 4 for PPAI, and 14 for SC, from which 5 were in common with AGECL, totaling 62 validated SNP for these traits and 30 candidate genes surrounding them. Considering just candidate genes closest to the top SNP of each chromosome, for AGECL 8 candidate genes were identified: COL8A1, PENK, ENSBTAG00000047425, BPNT1, ADAMTS17, CCHCR1, SUFU, and ENSBTAG00000046374. For PPAI, 3 genes emerged as candidates (PCBP3, KCNK10, and MRPS5), and for SC 8 candidate genes were identified (SNORA70, TRAC, ASS1, BPNT1, LRRK1, PKHD1, PTPRM, and ENSBTAG00000045690). Several candidate regions presented here were previously associated with puberty traits in cattle. The majority of emerging candidate genes are related to biological processes involved in reproductive events, such as maintenance of gestation, and some are known to be expressed in reproductive tissues. Our results suggested that some QTL controlling early puberty seem to be segregating across cattle breeds adapted to tropical conditions.

Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle.

Multitrait meta-analyses are a strategy to produce more accurate genome-wide association studies, especially for complex phenotypes. We carried out a meta-analysis study for traits related to sexual precocity in tropical beef cattle (Nellore and Brahman) aiming to identify important genomic regions affecting these traits. The traits included in the analyses were age at first calving (AFC), early pregnancy (EP), age at first corpus luteum (AGECL), first postpartum anoestrus interval (PPAI), and scrotal circumference (SC). The traits AFC, EP, and SCN were measured in Nellore cattle, while AGECL, PPAI, and SCB were measured in Brahman cattle. Meta-analysis resulted in 108 significant single-nucleotide polymorphisms (SNPs), at an empirical threshold P-value of 1.39 × 10-5 (false discovery rate [FDR] < 0.05). Within 0.5 Mb of the significant SNP, candidate genes were annotated and analyzed for functional enrichment. Most of the closest genes to the SNP with higher significance in each chromosome have been associated with important roles in reproductive function. They are TSC22D2, KLF7, ARHGAP29, 7SK, MAP3K5, TLE3, WDR5, TAF3, TMEM68, PPP1R15B, NR2F2, GALR1, SUFU, and KCNU1. We did not observe any significant SNP in BTA5, BTA12, BTA17, BTA18, BTA19, BTA20, BTA22, BTA23, BTA25, and BTA28. Although the majority of significant SNPs are in BTA14, it was identified significant associations in multiple chromosomes (19 out of 29 autosomes), which is consistent with the postulation that reproductive traits are complex polygenic phenotypes. Five proposed association regions harbor the majority of the significant SNP (76%) and were distributed over four chromosomes (P < 1.39 × 10-5, FDR < 0.05): BTA2 (5.55%) from 95 to 96 Mb, BTA4 (5.55%) from 94.1 to 94.8 Mb, BTA14 (59.26%) from 24 to 25 Mb and 29 to 30 Mb, and BTA21 (5.55%) from 6.7 Mb to 11.4 Mb. These regions harbored key genes related to reproductive function. Moreover, these genes were enriched for functional groups associated with immune response, maternal-fetal tolerance, pregnancy maintenance, embryo development, fertility, and response to stress. Further studies including other breeds and precocity traits could confirm the importance of these regions and identify new candidate regions for sexual precocity in beef cattle.

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle.

BACKGROUND: Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits. RESULTS: Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits. CONCLUSIONS: The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

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.

Bovine gene polymorphisms related to fat deposition and meat tenderness.

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.

Bovine gene polymorphisms related to fat deposition and meat tenderness

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.