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.

Gamevar.f90: a software package for calculating individual gametic diversity.

BACKGROUND: Traditional selection in livestock and crops focuses on additive genetic values or breeding values of the individuals. While traditional selection utilizes variation between individuals, differences between gametes within individuals have been less frequently exploited in selection programs. With the successful implementation of genomic selection in livestock and crops, estimation and selection for gametic variation is becoming possible. RESULTS: The gamevar.f90 software is designed to estimate individual-level variance of genetic values of gametes for complex traits in large populations. The software estimates the (co)variances of gametic diversity as well as other diversity parameters that are useful for selection programs and mating designs. The calculation is carried out chromosome by chromosome and can be easily parallelized. The gamevar.f90 program is written in Fortran with efficient computing algorithms in a user-friendly software package with easily-handled input and output files. Finally, we applied the program to estimate gametic variance for hundreds of bulls for lifetime net merit, productive life, and livability. The RPTA (relative predicted transmitting ability), assuming a future selection intensity (if) of 1.5, showed larger variance than GEBV/2, indicating that greater future genetic gains can be obtained using an index that includes gametic variances. We also used the relative coefficient of variation to estimate with 95% confidence the sample sizes required to observe 90% variability of the progeny for lifetime net merit (or to allow at maximum 10% of change in the EBV predicted from progeny data). CONCLUSIONS: Collectively, we develop an efficient computer program package, gamevar.f90, for estimating gametic variance for large numbers of individuals. The novel information on gametic variation will be useful in future animal and crop breeding programs.

Genome-wide association study applied to type traits related to milk yield in water buffaloes (Bubalus bubalis).

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

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.

Bos taurus-indicus hybridization correlates with intralocus sexual-conflict effects of PRDM9 on male and female fertility in Holstein cattle.

BACKGROUND: Crossover localization during meiotic recombination is mediated by the fast-evolving zinc-finger (ZnF) domain of gene PRDM9. To study its impact on dairy cattle performance, we compared its genetic variation between the relatively small Israeli (IL) Holsteins and the North American (US) Holsteins that count millions. RESULTS: Initially, we analyzed the major BTA1 haplotypes present in IL Holsteins based on the 10 most telomeric SNPs of the BovineSNP50 BeadChip. Sequencing of representative haplotype carriers indicated that for all frequent haplotypes (> 6%), the variable PRDM9 ZnF array consisted of seven tandem ZnF repeats. Two rare haplotypes (frequency < 4%) carried an indicine PRDM9, whereas all others were variants of the taurine type. These two haplotypes included the minor SNP allele, which was perfectly linked with a previously described PRDM9 allele known to induce unique localization of recombination hotspots. One of them had a significant (p = 0.03) negative effect on IL sire fertility. This haplotype combined the rare minor alleles of the only SNPs with significant (p < 0.05) negative substitution effects on US sire fertility (SCR). Analysis of telomeric SNPs indicated general agreement of allele frequencies (R = 0.95) and of the substitution effects on sire fertility (SCR, R = 0.6) between the US and IL samples. Surprisingly, the alleles that had a negative impact on male fertility had the most positive substitution effects on female fertility traits (DPR, CCR and HCR). CONCLUSIONS: A negative genetic correlation between male and female fertility is encoded within the BTA1 telomere. Cloning the taurine PRDM9 gene, which is the common form carried by Holsteins, we encountered the infiltration of an indicine PRDM9 variant into this population. During meiosis, in heterozygous males, the indicine PRDM9 variant may induce incompatibility of recombination hotspots and male infertility. However, this variant is associated with favorable female fertility, which would explain its survival and the general negative correlation (R = - 0.3) observed between male and female fertility in US Holsteins. Further research is needed to explain the mechanism underlying this positive effect and to devise a methodology to unlink it from the negative effect on male fertility during breeding.

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.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle.

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

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.