Genome-wide association for the outcome of fixed-time artificial insemination of Brahman heifers in northern Australia

Analysis of Runs of Homozygosity in Aberdeen Angus Cattle

A large number of cattle breeds have marked phenotypic differences. They are valuable models for studying genome evolution. ROH analysis can facilitate the discovery of genomic regions that may explain phenotypic differences between breeds affecting traits of economic importance. This paper investigates genome-wide ROH of 189 Aberdeen Angus bulls using the Illumina Bovine GGP HD Beadchip150K to structurally and functionally annotate genes located within or in close ROH of the Aberdeen Angus cattle genome. The method of sequential SNP detection was used to determine the ROH. Based on this parameter, two ROH classes were allocated. The total length of all ROH islands was 11,493 Mb. As a result of studying the genomic architecture of the experimental population of Aberdeen Angus bulls, nine ROH islands and 255 SNPs were identified. Thirteen of these overlapped with regions bearing 'selection imprints' previously identified in other breeds of cattle, and five of these regions were identified in other Aberdeen Angus populations. The total length of the ROH islands was 11,493 Mb. The size of individual islands ranged from 0.038 to 1.812 Mb. Structural annotation showed the presence of 87 genes within the identified ROH islets.

Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production traits.

Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers

BACKGROUND

In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production.

RESULTS

Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified.

CONCLUSIONS

In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.

Pregnancy rate per cycle is heritable and reduces with cycle in naturally mated tropically adapted beef cows

Pregnancy in cattle is the outcome of the complex process of initiation of cycling, fertilization, maternal recognition of pregnancy and foeto-placental development. Though much is known about initiation of cycling and associated risk factors, there are virtually no data on pregnancy rate per cycle for naturally mated cattle, especially for extensively managed, tropically adapted genotypes, which this study aimed to determine. Tropical composite (Bos indicus and African Sanga crosses with Bos taurus) and Brahman cattle (n = 2,181) of known pedigree in four-year groups at four sites were mated annually for 84 days. Body condition, ovarian function, pregnancies, calving and lactation were monitored through six full reproductive cycles using 4-8 weekly ultrasound of the reproductive tract outside the calving period and daily monitoring during calving. From this, dates of commencement of cycling and conception in each year were estimated for each animal, enabling calculation of established pregnancy for consecutive 21-day periods while cycling and of pregnancies within four months of calving while lactating (P4M). Pregnancy per 21-day period (cycle) during mating for cycling animals averaged 63%, 71%, 41% and 28% in four consecutive cycles. Pregnant per cycle was 2%-11% higher in tropical composites than in Brahmans. The only other consistently significant risk to becoming pregnant was if cycling commenced later than three weeks before mating commenced. P4M averaged 62% and was lower for cows in sub-optimal body condition and in first-parity and later-calving cows. Pregnant per cycle was moderately heritable (~20%), while heritability was moderate to high (33%) for P4M. Selection for pregnant per cycle could be achieved indirectly by selection for P4M, a trait that is readily measured.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Genome-Wide Association Study Provides Insights into Important Genes for Reproductive Traits in Nelore Cattle

Simple Summary

In this study, we investigated the association between single nucleotide polymorphisms (SNPs) and reproductive traits in order to identify candidate genes and biological pathways associated with these traits in Nelore beef cattle. The genome-wide association analysis revealed genomic regions that could explain part of the genetic variance of the studied traits. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, being important for reproductive efficiency. The identification of candidate genes that were associated with the studied traits as well as genes enriched in the functional terms and pathways may be useful for exploring the genetic architecture underlying reproductive traits and may be used in Nelore breeding programs.

Abstract

The identification of genomic regions associated with reproductive traits as well as their biological processes allows a better understanding of the phenotypic variability of these traits. This information could be applied to animal breeding programs to accelerate genetic gain. The aim of this study was to evaluate the association between single nucleotide polymorphisms (SNP) with a scrotal circumference at 365 days of age (SC365) and at 450 days of age (SC450), gestation length (GL) as a calf trait, age at first calving (AFC), accumulated productivity (ACP), heifer early calving until 30 months (HC30), and stayability (STAY) traits, in order to identify candidate genes and biological pathways associated with reproductive traits in Nelore cattle. The data set consisted of pedigree, phenotypes, and genotypes of Nelore cattle from the “Associação Nacional de Criadores e Pesquisadores” (ANCP). The association analyses were performed using the Weighted Single-Step Genome-Wide Association method; the regions, consisting of 10 consecutive SNP, which explained more than 0.5% of additive genetic variance, were considered as a significant association. A total of 3, 6, 7, 5, 10, 25, and 12 windows were associated with SC355, SC450, GL, AFC, ACP, HC30, and STAY, respectively. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, among them CAMK1D, TASP1, ACOXL, RAB11FIP5, and SFXN5. Moreover, the genes were enriched in functional terms, like negative regulation of fat cell differentiation, fatty acid alpha-oxidation, and sphingolipids signaling pathway. The identification of the genes associated with the traits, as well as genes enriched in the terms and pathway mentioned above, should contribute to future biological validation studies and may be used as candidate genes in Nelore breeding programs.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

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.

Association of lipid-related genes implicated in conceptus elongation with female fertility traits in dairy cattle

Elongation of the preimplantation conceptus is a requirement for pregnancy success in ruminants, and failures in this process are highly associated with subfertility in dairy cattle. Identifying genetic markers that are related to early conceptus development and survival and utilizing these markers in selective breeding can improve the reproductive efficiency of dairy herds. Here, we evaluated the association of 1,679 SNP markers within or close to 183 candidate genes involved in lipid metabolism of the elongating conceptus with different fertility traits in US Holstein cattle. A total of 27,371 bulls with predicted transmitting ability records for daughter pregnancy rate, cow conception rate, and heifer conception rate were used as the discovery population. The associations found in the discovery population were validated using 2 female populations (1,122 heifers and 2,138 lactating cows) each with 4 fertility traits, including success to first insemination, number of services per conception, age at first conception for heifers, or days open for cows. Marker effects were estimated using a linear mixed model with SNP genotype as a linear covariate and a random polygenic effect. After multiple testing correction, 39 SNP flagging 27 candidate genes were associated with at least one fertility trait in the discovery population. Of these 39 markers, 3 SNP were validated in the heifer population and 4 SNP were validated in the cow population. The 3 SNP validated in heifers are located within or near genes CAT, MYOF, and RBP4, and the 4 SNP validated in lactating cows are located within or close to genes CHKA, GNAI1, and HMOX2. These validated genes seem to be relevant for reducing pregnancy losses, and the SNP within these genes are excellent candidates for inclusion in genomic tests to improve reproductive performance in dairy cattle.

Across-breed validation study confirms and identifies new loci associated with sexual precocity in Brahman and Nellore cattle

The aim of this study was to identify candidate regions associated with sexual precocity in Bos indicus. Nellore and Brahman were set as validation and discovery populations, respectively. SNP selected in Brahman to validate in Nellore were from gene regions affecting reproductive traits (G1) and significant SNP (p ≤ 10^-3 ) from a meta-analysis (G2). In the validation population, early pregnancy (EP) and scrotal circumference (SC) were evaluated. To perform GWAS in validation population, we used regression and Bayes C. SNP with p ≤ 10^-3 in regression and Bayes factor ≥3 in Bayes C were deemed significant. Significant SNP (for EP or SC) or SNP in their ±250 Kb vicinity region, which were in at least one discovery set (G1 or G2), were considered validated. SNP identified in both G1 and G2 were considered candidate. For EP, 145 SNP were validated in G1 and 41 in G2, and for SC, these numbers were 14 and 2. For EP, 21 candidate SNP were detected (G1 and G2). For SC, no candidate SNP were identified. Validated SNP and their vicinity region were located close to quantitative trait loci or genes related to reproductive traits and were enriched in gene ontology terms related to reproductive success. These are therefore strong candidate regions for sexual precocity in Nellore and Brahman.

Genomic analyses of tropical beef cattle fertility based on genotyping pools of Brahman cows with unknown pedigree

We introduce an innovative approach to lowering the overall cost of obtaining genomic EBV (GEBV) and encourage their use in commercial extensive herds of Brahman beef cattle. In our approach, the DNA genotyping of cow herds from 2 independent properties was performed using a high-density bovine SNP chip on DNA from pooled blood samples, grouped according to the result of a pregnancy test following their first and second joining opportunities. For the DNA pooling strategy, 15 to 28 blood samples from the same phenotype and contemporary group were allocated to pools. Across the 2 properties, a total of 183 pools were created representing 4,164 cows. In addition, blood samples from 309 bulls from the same properties were also taken. After genotyping and quality control, 74,584 remaining SNP were used for analyses. Pools and individual DNA samples were related by means of a "hybrid" genomic relationship matrix. The pooled genotyping analysis of 2 large and independent commercial populations of tropical beef cattle was able to recover significant and plausible associations between SNP and pregnancy test outcome. We discuss 24 SNP with significant association ( < 1.0 × 10) and mapped within 40 kb of an annotated gene. We have established a method to estimate the GEBV in young herd bulls for a trait that is currently unable to be predicted at all. In summary, our novel approach allowed us to conduct genomic analyses of fertility in 2 large commercial Brahman herds managed under extensive pastoral conditions.

Polymorphisms and genes associated with puberty in heifers

Puberty onset is a multifactorial process influenced by genetic determinants and environmental conditions, especially nutritional status. Genes, genetic variations, and regulatory networks compose the molecular basis of achieving puberty. In this article, we reviewed the discovery of multiple polymorphisms and genes associated with heifer puberty phenotypes and discuss the opportunities to use this evolving knowledge of genetic determinants for breeding early pubertal Bos indicus-influenced cattle. The discovery of polymorphisms and genes was mainly achieved through candidate gene studies, quantitative trait loci analyses, genome-wide association studies, and recently, global gene expression studies (transcriptome). These studies are recapitulated and summarized in the current review.