Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle.

Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.

Characterization of runs of homozygosity, heterozygosity-enriched regions, and population structure in cattle populations selected for different breeding goals.

BACKGROUND: A decline in the level of genetic diversity in livestock can result in reduced response to selection, greater incidence of genetic defects, and inbreeding depression. In this context, various metrics have been proposed to assess the level of genetic diversity in selected populations. Therefore, the main goals of this study were to: 1) investigate the population structure of 16 cattle populations from 15 different pure breeds or composite populations, which have been selected for different breeds goals; and, 2) identify and compare runs of homozygosity (ROH) and heterozygosity-enriched regions (HER) based on different single nucleotide polymorphism (SNP) panels and whole-genome sequence data (WGS), followed by functional genomic analyses. RESULTS: A total of 24,187 ROH were found across all cattle populations, with 55% classified in the 2-4 Mb size group. Fourteen homozygosity islands were found in five populations, where four ROH islands located on BTA1, BTA5, BTA16, and BTA19 overlapped between the Brahman (BRM) and Gyr (GIR) breeds. A functional analysis of the genes found in these islands revealed candidate genes known to play a role in the melanogenesis, prolactin signaling, and calcium signaling pathways. The correlations between inbreeding metrics ranged from 0.02 to 0.95, where the methods based on homozygous genotypes (FHOM), uniting of gametes (FUNI), and genotype additive variance (FGRM) showed strong correlations among them. All methods yielded low to moderate correlations with the inbreeding coefficients based on runs of homozygosity (FROH). For the HER, 3576 runs and 26 islands, distributed across all autosomal chromosomes, were found in regions containing genes mainly related to the immune system, indicating potential balancing selection. Although the analyses with WGS did not enable detection of the same island patterns, it unraveled novel regions not captured when using SNP panel data. CONCLUSIONS: The cattle populations that showed the largest amount of ROH and HER were Senepol (SEN) and Montana (MON), respectively. Overlapping ROH islands were identified between GIR and BRM breeds, indicating a possible historical connection between the populations. The distribution and pattern of ROH and HER are population specific, indicating that different breeds have experienced divergent selection processes or different genetic processes.

Identification of novel candidate genes for age at first calving in Nellore cows using a SNP chip specifically developed for Bos taurus indicus cattle.

The age at first calving has a great economic impact on the beef cattle system and calving at 24 months is an objective of selection for a more efficient herd. However, an age at first calving around 36 months has been observed for Nellore cattle in Brazil. Thus, a genome-wide association study (GWAS) was carried out with 8376 records of age at first calving and 3239 animals genotyped with the GGP-Indicus 35K, which has been developed specifically for Bos taurus indicus. The weighted single-step genomic best linear unbiased prediction method was used, with adjacent SNPs (single nucleotide polymorphisms) in genomic windows of 1.0 Mb. After quality control, 3239 (2161 males and 1078 females) animals genotyped for 30,519 SNPs were used in GWAS analysis. The average and standard deviation of age at first calving were 1041.7 and 140.6 days, respectively. The heritability estimate was 0.10 ± 0.02. The GWAS analysis found seven genomic regions in BTA1, 2, 5, 12, 18, 21, and 24, which explained a total of 11.24% of the additive genetic variance of age at first calving. In these regions were found 62 protein coding genes, and the genes HSD17B2, SERPINA14, SERPINA1, SERPINA5, STAT1, NFATC1, ATP9B, CTDP1, THPO, ECE2, PSMD2, EIF4G1, EIF2B2, DVL3, POLR2H, TMTC2, and GPC6 are possible candidates for age at first birth due their function. Moreover, two molecular functions ("serine-type endopeptidase inhibitor activity" and "negative regulation of endopeptidase activity") were significant, which depend on several serpin genes. The use of a SNP chip developed especially for Bos taurus indicus allowed to find genomic regions for age at first calving, which are close to QTLs previously reported for other reproduction-related traits. Future studies can reveal the causal variants and their effects on reproductive precocity of Nellore cows.