Genomic analysis of genotype-environment interaction in age at first calving of Murrah buffaloes.

Age at first calving (AFC) is a measure of sexual maturity associated with the start of productive life of dairy animals. Additionally, a lower AFC reduces the generation interval and early culling of females. However, AFC has low heritability, making it a trait highly influenced by environmental factors. In this scenario, one way to improve the reproductive performance of buffalo cows is to select robust animals according to estimated breeding value (EBV) using models that include genotype-environment interaction (GEI) with the application of reaction norm models (RNMs). This can be achieved by understanding the genomic basis related to GEI of AFC. Thus, in this study, we aimed to predict EBV considering GEI via the RNM and identify candidate genes related to this component in dairy buffaloes through genome-wide association studies (GWAS). We used 1795 AFC records from three Murrah buffalo herds and formed environmental gradients (EGs) from contemporary group solutions obtained from genetic analysis of 270-day cumulative milk yield. Heritability estimates ranged from 0.15 to 0.39 along the EG. GWAS of the RNM slope parameter identified important genomic regions. The genomic window that explained the highest percentage of genetic variance of the slope (0.67%) was located on BBU1. After functional analysis, five candidate genes were detected, involved in two biological processes. The results suggested the existence of a GEI for AFC in Murrah buffaloes, with reclassification of animals when different environmental conditions were considered. The inclusion of genomic information increased the accuracy of breeding values for the intercept and slope of the reaction norm. GWAS analysis suggested that important genes associated with the AFC reaction norm slope were possibly also involved in biological processes related to lipid metabolism and immunity.

Effects of aeration and season of the year on fish waste composting and compost quality.

Employing forced aeration (FA) in composting static windrows (SW) from fish waste (FW) has the potential to enhance the development of process and, organic fertiliser quality. However, due to the impact of season, the FA may lead to excessive drying of SW and, difficulty in thermophilic temperature maintenance. The aim of this study was to assess the effects of passive aeration (PA) and FA on the composting of FW in SW during the summer and winter seasons. The temperatures of the windrows remained within the thermophilic range for most of the composting period, with peak temperatures observed shortly after starting and turning the windrows (at 50 and 70 days). The aeration benefited the initial TS degradations, resulting in 86.66 and 45.99% of the TS total reduced to FA and PA piles, at 50 days during the winter. The C organic reduction was 77.77 and 76.33% in summer and winter to FA piles, respectively, but this reduction was 59.24% and 67.82% for winter and summer, respectively, in PA windrows. At 50 days, the N reduction in FA piles was already at 70.32% and 71.87% for winter and summer. The volatile solids reductions were significantly higher (p < 0.01) in FA piles during the summer. Although the FA has been shown to enhance the organic constituents' degradation during the composting of FW, its adoption was not enough to improve the compost composition. Thus, by conducting piles on a small scale, with the perforated wall, as described in this study, the FA could be dispensed.

Identification of genomic regions related to age at first calving and first calving interval in water buffalo using single-step GBLUP.

In Brazil, water buffaloes have been used to produce milk for mozzarella cheese production. Consequently, the main selection criterion applied for the buffalo genetic improvement is the estimated mozzarella yield as a function of milk, fat and protein production. However, given the importance of reproductive traits in production systems, this study aimed to use techniques for identifying genomic regions that affect the age at first calving (AFC) and first calving interval (FCI) in buffalo cows and to select candidate genes for the identification of QTL and gene expression studies. The single-step GBLUP method was used for the identification of genomic regions. Windows of 1 Mb containing single-nucleotide polymorphisms were constructed and the 10 windows that explained the greatest proportion of genetic variance were considered candidate regions for each trait. Genes present into the selected windows were identified using the UOA_WB_1 assembly as the reference, and their ontology was defined with the Panther tool. Candidate regions for both traits were identified on BBU 3, 12, 21 and 22; for AFC, candidates were detected on BBU 6, 7, 8, 9 and 15 and for first calving interval on BBU 4, 14 and 19. This study identified regions with great contribution to the additive genetic variance of age at first calving and first calving interval in the population of buffalo cows studied. The ROCK2, PMVK, ADCY2, MAP2K6, BMP10 and GFPT1 genes are main candidates for reproductive traits in water dairy buffaloes, and these results may have future applications in animal breeding programs or in gene expression studies of the species.

Estimation of (co)variance components for birth weight in Nellore cattle using alternative data editing criteria.

The aim of this study was to estimate the (co)variance components and breeding values for birthweight (BW) in Nellore cattle by considering or not identical weights that exhibit a high frequency within the contemporary group (CG). A total of 175,258 BW records of Nellore cattle born between 2002 and 2018 were used. The CG was formed by farm, year of birth, sex and feeding regime at birth. CGs with more than 16% of identical BW values were eliminated, generating a data file called BWd. Another file was created without removing these animals (BWt). A mixed linear model was used for statistical analysis, which included fixed and random effects. In both data files analysed, single-trait analysis was performed by Bayesian inference. The mean direct and maternal heritability for BW and the correlation between direct and maternal effects were 0.27, 0.07 and -0.07 for BWt, respectively, and 0.30, 0.093 and -0.07 for BWd. This method should affect the estimation of genetic merits of animals for BW, providing greater safety in the choice of sires.

Analysis of the population structure of buffaloes in Brazil.

The aim of this study was to evaluate the population structure of Brazilian buffaloes using of pedigree information. The pedigree used in the analyses included records of 16,915 animals. The population parameters were obtained through the ENDOG software. The estimates of mean of inbreeding (F) and average relatedness (AR) coefficients were 3.22% and 5.99%, respectively. The average generation interval was 6.39 years. The effective number of founders and ancestors was 28 and 22, respectively. In this study, we concluded that the selection of individuals with lower AR is necessary to avoid an increase in matings between inbreeding individuals in this population, in order to obtain greater genetic gain by selection.

Linkage Disequilibrium-Based Inference of Genome Homology and Chromosomal Rearrangements Between Species.

The aim of this study was to analyze the genomic homology between cattle (Bos taurus) and buffaloes (Bubalus bubalis) and to propose a rearrangement of the buffalo genome through linkage disequilibrium analyses of buffalo SNP markers referenced in the cattle genome assembly and also compare it to the buffalo genome assembly. A panel of bovine SNPs (single nucleotide polymorphisms) was used for hierarchical, non-hierarchical and admixture cluster analyses. Thus, the linkage disequilibrium information between markers of a specific panel of buffalo was used to infer chromosomal rearrangement. Haplotype diversity and imputation accuracy of the submetacentric chromosomes were also analyzed. The genomic homology between the species enabled us to use the bovine genome assembly to recreate a buffalo genomic reference by rearranging the submetacentric chromosomes. The centromere of the submetacentric chromosomes exhibited high linkage disequilibrium and low haplotype diversity. It allowed hypothesizing about chromosome evolution. It indicated that buffalo submetacentric chromosomes are a centric fusion of ancestral acrocentric chromosomes. The chronology of fusions was also suggested. Moreover, a linear regression between buffalo and cattle rearranged assembly and the imputation accuracy indicated that the rearrangement of the chromosomes was adequate. When using the bovine reference genome assembly, the rearrangement of the buffalo submetacentric chromosomes could be done by SNP BTA (chromosome of Bos taurus) calculations: shorter BTA (shorter arm of buffalo chromosome) was given as [(shorter BTA length - SNP position in shorter BTA)] and larger BTA length as [shorter BTA length + (larger BTA length - SNP position in larger BTA)]. Finally, the proposed linkage disequilibrium-based method can be applied to elucidate other chromosomal rearrangement events in other species with the possibility of better understanding the evolutionary relationship between their genomes.

Genetic parameters for tonic immobility, body weight, and morphological traits of the red-winged tinamou (Rhynchotus rufescens).

This study was carried out to estimate genetic parameters for morphology, body weight, and tonic immobility traits in the red-winged tinamou (Rhynchotus rufescens). Information on 690 birds was used and genetic parameters were estimated using Bayesian methods under a multi-trait animal model. The following traits were considered in this study: tarsal length (TL), bill length (BL), wing length (WL), head width (HW), bill width (BW), mature weight (MW), weight at 90 days (W90), and tonic immobility (TI). The heritability showed estimates between 0.15 for wing length and 0.56 for bill length. Positive and negative genetic correlations were estimated, ranging from - 0.33 to 0.81. All the morphological, production, and behavioral traits studied will have moderate to high response to selection. The body weight at 90 days is a better alternative for use in breeding programs and its selection would not lead to an increase in the time of tonic immobility. Both the selection for weight gain and for reduction of tonic immobility time would lead to an increase in the size of the legs of the red-winged tinamou, which could be advantageous for thermal control of these birds in tropical systems.

Genome-wide association studies for growth traits in buffaloes using the single step genomic BLUP.

Growth traits are important for the profitability of buffalo breeding systems, since in general, these animals are raised both for meat and milk. In this study, the single-step genomic BLUP method was employed to prospect the genomic regions' associated with weight at standard ages of 100, 210, 365, and 550 days in a buffalo population, aiming to identify genes with stronger expression for those characteristics. We found 6, 1, 2, and 5 SNPs significantly associated (p value < 10-5) with weight at 100, 210, 365, and 550 days of age, respectively, where those SNPs respectively explained 0.164, 0.040, 0.044, and 0.213% of the additive variance of each trait. SNP AX-85099682 (BBU24) was significant for weight at 100, 210, and 365 days, indicating the existence of a possible QTL affecting the initial growth rate of buffaloes. All told, eight genes (CBLB, TRNAG-UCC, GADD45B, LOC112583811, MGAT4C, KCNMA1, SLC5A2, and TGFB1I1) were identified as candidates for the growth traits of buffaloes. However, molecular and gene expression studies are necessary to validate these genes for subsequent use in programs for genetic improvement of the species.

Genotype-environment interaction for age at first calving in buffaloes, using the reaction norm model.

The objective of this study was to evaluate the genotype-environment interaction effect on age at first calving in buffaloes. The records were analysed using two approaches: (a) standard animal model and (b) reaction norm model. For the reaction norm analysis, two environmental gradients were formed, using age of first calving or milk yield group contemporary average. The results showed differences in the heritability estimates when using the two approaches. The reaction norm model indicated high heritability in more favourable environments and low magnitude genetic correlations between extreme environments. Based on our findings, we verified the significance of the genotype-environment interaction effect on age at first calving in buffaloes.

Use of single-step genome-wide association studies for prospecting genomic regions related to milk production and milk quality of buffalo.

The aim of this research communication was to identify chromosome regions and genes that could be related to milk yield (MY), milk fat (%F) and protein percentage (%P) in Brazilian buffalo cows using information from genotyped and non-genotyped animals. We used the 90 K Axiom® Buffalo Genotyping array. A repeatability model was used. An iterative process was performed to calculate the weights of markers as a function of the squared effects of Single Nucleotide Polymorphism (SNP) and allele frequencies. The 10 SNPs with the largest effects for MY, %F and %P were studied and they explained 7·48, 9·94 and 6·56% of the genetic variance, respectively. These regions harbor genes with biological functions that could be related to the traits analyzed. The identification of such regions and genes will contribute to a better understanding of their influence on milk production and milk quality traits of buffaloes.

Study of the effect of genotype-environment interaction on age at first calving and production traits in Nellore cattle using multi-trait reaction norms and Bayesian inference.

This study investigated the effects of genotype-environment interaction on yearling weight, age at first calving and post-weaning weight gain in Nellore cattle using multi-trait reaction norm models. The environmental gradient was defined as a function of the mean yearling weight of the contemporary groups. A first-order random regression sire model with four classes of residual variance was used in the analyses and Bayesian methods were applied to estimate the (co)variance components. The heritability estimates ranged from 0.284 to 0.547, 0.222 to 0.316 and 0.256 to 0.522 for yearling weight, age at first calving and post-weaning weight gain, respectively. The lowest genetic correlations between environment groups for each trait were 0.38, 0.02 and 0.04 for yearling weight, age at first calving and post-weaning weight gain, respectively. Differences in the correlation estimates were observed between traits in the same environments, with the magnitude of the estimates tending toward zero as the environment improved. The results highlight the importance of including genotype-environment interactions in genetic evaluation programs considering the differences observed between environmental groups not only in terms of heritability, but also of genetic correlations.

Genomic selection in dairy cattle simulated populations.

Genomic selection is arguably the most promising tool for improving genetic gain in domestic animals to emerge in the last few decades, but is an expensive process. The aim of this study was to evaluate the economic impact related to the implementation of genomic selection in a simulated dairy cattle population. The software QMSim was used to simulate genomic and phenotypic data. The simulated genome contained 30 chromosomes with 100 cm each, 1666 SNPs markers equally spread and 266 QTLs randomly designated for each chromosome. The numbers of markers and QTLs were designated according to information available from Animal QTL (http://www.animalgenome.org/QTLdb) and Bovine QTL (http://bovineqtl.tamu.edu/). The allelic frequency changes were assigned in a gamma distribution with alpha parameters equal to 0·4. Recurrent mutation rates of 1·0e-4 were assumed to apply to markers and QTLs. A historic population of 1000 individuals was generated and the total number of animals was reduced gradually along 850 generations until we obtained a number of 200 animals in the last generation, characterizing a bottleneck effect. Progenies were created along generations from random mating of the male and female gametes, assuming the same proportion of both genders. Than the population was extended for another 150 generations until we obtained 17 000 animals, with only 320 male individuals in the last generation. After this period a 25 year of selection was simulated taking into account a trait limited by sex with heritability of 0·30 (i.e. milk yield), one progeny/cow/year and variance equal to 1·0. Annually, 320 bulls were mated with 16 000 dams, assuming a replacement rate of 60 and 40% for males and females, respectively. Selection and discard criteria were based in four strategies to obtain the EBVs assuming as breeding objective to maximize milk yield. The progeny replaced the discarded animals creating an overlapping generation structure. The selection strategies were: RS is selection based on random values; PS is selection based on phenotypic values; Blup is selection based on EBVs estimated by BLUP; and GEBV is selection based on genomic estimated breeding values in one step, using high (GBlup) and low (GBlupi) density panels. Results indicated that the genetic evaluation using the aid of genomic information could provide better genetic gain rates in dairy cattle breeding programs as well as reduce the average inbreeding coefficient in the population. The economic viability indicators showed that only Blup and GBlup/GBlupi strategies, the ones that used milk control and genetic evaluation were economic viable, considering a discount rate of 6·32% per year.

Random regression models to estimate genetic parameters for weights in Murrah buffaloes.

This article reports genetic analysis of the weight at different ages of Murrah water buffaloes, using random regression models (RRM). Models ranging from third to sixth order polynomial were used to describe direct genetic and animal permanent environmental effects. Contemporary group was included as a fixed effect, and a cubic polynomial was used to model the mean curve of the population. The residual was modeled considering a log-linear function. Two models were selected for study of genetic parameters. The first model included third and sixth order polynomials for direct genetic and animal permanent environmental effects (M36). The second model included sixth order polynomials for all random effects (M66). The estimates of heritability varied from 0.16 + 0.04 (44 days) to 0.38 + 0.04 (568 days) for model M36 and from 0.16 + 0.05 (33 days) to 0.42 + 0.05 (600 days) for model M66. Regarding estimates of the correlation for all effects, the magnitude tended to decline with the increase of the time span between measurements. These results indicate that the species has potential for genetic selection based on weight at different ages, since we found favorable genetic variability within the herd, with selection likely to be more efficient at ages near 600 days.

Genetic association between SNPs in the DGAT1 gene and milk production traits in Murrah buffaloes.

This study identified polymorphisms in the DGAT1 gene in Murrah buffaloes and investigated the associations to milk production and quality traits (milk, fat and protein yields and percentages, somatic cell count). Genomic DNA was extracted from hair follicles collected from the tail of 196 females. Three SNPs were identified in DGAT1 gene by sequencing. Statistical analyses were performed to verify the linkage and the association between polymorphisms and traits. The estimated value of r (2) between two SNPs in exon 17 (g.11,783G > A and g.11,785 T > C) was 0.029. SNP g.11,785 T > C was significantly associated (P < 0.05) to fat and protein percentage. Dominance effect was significant for milk and fat yields and protein percentage (P < 0.05). The additive effect of the SNP g.11,785 T > C was significant for protein production and somatic cell count (P < 0.05). It indicates that assisted marker selection might be done with considerations to balance production and udder health.

Polymorphisms in oxytocin and α1a adrenergic receptor genes and their effects on production traits in dairy buffaloes.

The use of molecular markers may auxiliary the buffalo breeding. The oxytocin (OXT) and the adrenergic receptor α1A (ADRA1A) may be involved in milk ejection in ruminants. The aim of this study was to verify the existence of polymorphisms in the OXT and ADRA1A genes and their associations with milk production traits. A total of 220 buffaloes were genotyped using PCR-RFLP for both genes. The SNP identified in the ADRA1A gene was associated with protein percentage in dairy buffaloes. This is the first report of such association in the literature, which has not been studied in other species.

Polymorphisms in the MTRN1A gene and their effects on the productive and reproductive traits in buffaloes.

The main environmental factor that affects the regulation of reproductive seasonality is photoperiod through its effects on melatonin secretion. The melatonin receptor MTRN1A appears to be involved in regulating the reproductive seasonality and milk production in the period. The aim of this study was to identify polymorphisms in the MTRN1A gene and their possible associations with milk, fat and protein productions, fat and protein percentages, age at first calving, and first calving interval in buffaloes. Three genotypes (CC, CT, and TT) were identified by PCR-RFLP, and there was a significant association with protein percentage (P < 0.0001). Further studies are necessary to better understand the influence of melatonin gene and their receptors in the productive functions of buffaloes.

Effects of a single nucleotide polymorphism in the leptin gene on the productive traits of dairy buffaloes (Bubalus bubalis).

The gene responsible for coding the leptin hormone has been associated with productive and reproductive traits in cattle. In dairy cattle, different polymorphisms found in the leptin gene have been associated with several traits of economic interest, such as energy balance, milk yield and composition, live weight, fertility and dry matter consumption. The aim of this study was to detect genetic variability in the leptin gene of buffaloes and to test possible associations with milk yield, fat and protein percentages, age at first calving and first calving interval. Three genotypes (AA, AG and GG) were identified by polymerase chain reaction-restriction fragment length polymorphism, which presented genotypic frequencies of 0.30, 0.54 and 0.16, respectively. The allele frequencies were 0.57 for the A allele and 0.43 for the G allele. No significant effects were found in the present study, but there is an indicative that leptin gene affects lipid metabolism.

First polymorphisms in JY-1 gene in cattle (Bos taurus indicus) and their association with sexual precocity and growth traits.

Protein JY-1 is a bovine oocyte-specific protein that regulates granulosa cell function and is involved in early embryonic development, influencing the chance of pregnancy. This study investigated molecular markers for the JY-1 gene. Seven SNPs were identified in exon 3 of the gene. The positions of the SNPs in the exon and the respective substitutions are: 163 (T/C), 281 (T/C), 321 (T/C), 532 (T/C), 652 (A/G), 679 (T/C), and 722 (G/C) (GenBank: JN592587 and JF262042.2). SNP 163 is located in a coding region and causes a proline-to-leucine substitution. The other SNPs are located in the 3'UTR region. SNPs 163, 281, 321, and 679 were genotyped in 297 Nellore heifers and the haplotypes were constructed. The haplotypes of JY-1 were not correlated with the traits studied at 5 %.

Genetic parameters for milk, fat and protein yields in Murrah buffaloes (Bubalus bubalis Artiodactyla, Bovidae).

The objective of the present study was to estimate genetic parameters for test-day milk, fat and protein yields and 305-day-yields in Murrah buffaloes. 4,757 complete lactations of Murrah buffaloes were analyzed. Co-variance components were estimated by the restricted maximum likelihood method. The models included additive direct genetic and permanent environmental effects as random effects, and the fixed effects of contemporary group, milking number and age of the cow at calving as linear and quadratic covariables. Contemporary groups were defined by herd-year-month of test for test-day yields and by herd-year-season of calving for 305-day yields. The heritability estimates obtained by two-trait analysis ranged from 0.15 to 0.24 for milk, 0.16 to 0.23 for protein and 0.13 to 0.22 for fat, yields. Genetic and phenotypic correlations were all positive. The observed population additive genetic variation indicated that selection might be an effective tool in changing population means in milk, fat and protein yields.

Genetic parameters for milk, fat and protein yields in Murrah buffaloes (Bubalus bubalis Artiodactyla, Bovidae)

The objective of the present study was to estimate genetic parameters for test-day milk, fat and protein yields and 305-day-yields in Murrah buffaloes. 4,757 complete lactations of Murrah buffaloes were analyzed. Co-variance components were estimated by the restricted maximum likelihood method. The models included additive direct genetic and permanent environmental effects as random effects, and the fixed effects of contemporary group, milking number and age of the cow at calving as linear and quadratic covariables. Contemporary groups were defined by herd-year-month of test for test-day yields and by herd-year-season of calving for 305-day yields. The heritability estimates obtained by two-trait analysis ranged from 0.15 to 0.24 for milk, 0.16 to 0.23 for protein and 0.13 to 0.22 for fat, yields. Genetic and phenotypic correlations were all positive. The observed population additive genetic variation indicated that selection might be an effective tool in changing population means in milk, fat and protein yields.