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.

Estimation of genetic parameters for the tick and hemoparasite burden in Angus cattle.

The study was conducted with the objective of estimating genetic and phenotypic parameters for tick (CRM) and Babesia bigemina (IBBi), Babesia bovis (IBBo), and Anaplasma marginale (IAM) burden in Angus female breed in Brazil. The sample group was composed of Angus females raised in herds located in a region of endemic instability for cattle tick fever in the state of Rio Grande Sul (RS), Brazil. The variance components were estimated using Bayesian inference and Gibbs sampling algorithm, considering a multi-trait animal model. Heritability estimates showed values of low magnitude, ranging from 0.03 (IBBo) to 0.16 (CRM), while repeatability estimates ranged between 0.07 (IBBo) and 0.21 (CRM). Regarding the genetic correlation estimates, the values showed low (-0.01 for IBBo × IAM) to moderate (0.55 between IBBi × IAM) magnitudes. The results indicate that it is possible to use tick count and hemoparasite infection levels as selection criteria, with small genetic gains.

Selection effect for growth traits on energy requirements in beef Nelore steers.

Growth data of 77,372 Nelore steers were used to estimate the selection effect on energy requirements considering two beef production systems: cow-calf and slaughter cycles. All the animals had measures from 120 days to 7 years old. The parameters necessary to evaluate the selection effect on energy requirements were obtained by random regression analysis using Legendre polynomials. The models included additive direct and maternal effects, and animal and maternal permanent environmental effects as random terms. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of animal age (cubic regression) were considered random covariables. The coefficients from the model M3353_5 were used to calculate the genetic gains necessary to predict the increase in phenotypes. The selection was simulated for body weight (BW) and weight gain (WG) at different ages and energy requirements were calculated using NRC equations. The cost of feed was calculated for a cow-calf and slaughter cycle of production considering a system of Brachiaria decumbens pasture without supplementation. In slaughter system, the selection for weight of 365 days of age is the best option. In cow-calf systems, the selection W120 is the best choice.

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.

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.

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.

Population structure of Simmental beef cattle using pedigree analysis.

The purpose of this study was to estimate parameters related to the population structure and genetic diversity in the Simental breed based on the pedigree information of 77,553 animals. The individual coefficients of inbreeding and average relatedness, number of complete generations, coefficient of change of inbreeding, effective size, effective number of founders, number of ancestors and generation interval were calculated. Using the Simmental cattle information, the mean inbreeding and average relatedness coefficients were 1.49% and 0.99%, respectively. The effective number of founders and ancestors was 163 and 132, respectively, and the effective population size was 48.03. Despite the relatively small inbreeding coefficient, some of the estimated population parameters indicated the need to adopt measures to maintain the genetic variability of the population.

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.

Genomic study of the resilience of buffalo cows to a negative energy balance.

The research article was carried out with the objective of studying the genetic variation on the resilience of buffaloes to negative energy balance-NEB (measured by changes in body weight in early lactation)-as well as investigating genomic regions of interest for this trait. A model of reaction norms was used, considering milk production as the trait to be analyzed and solutions of the contemporary groups to weight changes as environmental gradient. In this methodology, the genetic value of the slope represents the measure of resilience of the animals. After the estimation step, a genome-wide association analysis was performed for the slope of the reaction norms model, to obtain a list of windows and associated genes. The heritability estimates for milk production over the resilience gradient ranged from 0.13 to 0.28, with lower values in the intermediate environmental groups. Regarding the productive resilience of dairy buffalo cows to NEB, the genomic windows with the highest contribution to the genetic variance were detected on chromosomes BBU 1, 2, 3, 4, 9, 12, 19, and 21. A functional analysis of the genes described in the selected windows indicated association with metabolic routes related to growth and immunity of the animals, with an emphasis on the STAT6 gene. The results presented indicate that there is for this trait genetic variation to be used as selection criteria, in addition to genomic regions that can increase the precision of the selection.

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.

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.

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.

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.