Genome-wide association study for female fertility in Nordic Red cattle

Genomic insights into key genes and QTLs involved in cattle reproduction

From an economic standpoint, reproductive characteristics are fundamental for sustainable production, particularly for monotocous livestock like cattle. A longer inter-calving interval is indicative of low reproductive capacity. This issue changes the dynamics of current and future lactations since it necessitates more inseminations, veterinary care, and hormone interventions. Various reproductive phenotypes, including ovulation, mating, fertility, pregnancy, embryonic growth, and calving-related traits, are observed in dairy cattle, and these traits have been associated with several QTLs. Calving ease, age at puberty, scrotal circumference, and inseminations per conception have been associated with 4437, 10623, 10498, and 2476 Quantitative Trait Loci (QTLs), respectively. This data offers valuable insights into enhancing and comprehending reproductive traits in livestock breeding. Studying QTLs associated with reproductive traits has far-reaching implications across various fields, from agriculture and animal husbandry to human health, evolutionary biology, and conservation. It provides the foundation for informed breeding practices, advances in biotechnology, and a deeper understanding of the genetic underpinnings of reproduction.

Identifying Regions of the Genome Associated with Conception Rate to the First Service in Holstein Heifers Bred by Artificial Insemination and as Embryo Transfer Recipients

Heifer conception rate to the first service (HCR1) is defined as the number of heifers that become pregnant to the first breeding service compared to the heifers bred. This study aimed to identify loci associated and gene sets enriched for HCR1 for heifers that were bred by artificial insemination (AI, n = 2829) or were embryo transfer (ET, n = 2086) recipients, by completing a genome-wide association analysis and gene set enrichment analysis using SNP data (GSEA-SNP). Three unique loci, containing four positional candidate genes, were associated (p < 1 × 10-5) with HCR1 for ET recipients, while the GSEA-SNP identified four gene sets (NES ≥ 3) and sixty-two leading edge genes (LEGs) enriched for HCR1. While no loci were associated with HCR1 bred by AI, one gene set and twelve LEGs were enriched (NES ≥ 3) for HCR1 with the GSEA-SNP. This included one gene (PKD2) shared between HCR1 AI and ET services. Identifying loci associated or enriched for HCR1 provides an opportunity to use them as genomic selection tools to facilitate the selection of cattle with higher reproductive efficiency, and to better understand embryonic loss.

Effect of genotyping density on the detection of runs of homozygosity and heterozygosity in cattle

Runs of homozygosity (ROHom) are contiguous stretches of homozygous regions of the genome. In contrast, runs of heterozygosity (ROHet) are heterozygosity-rich regions. The detection of these two types of genomic regions (ROHom and ROHet) is influenced by the parameters involved in their identification and the number of available single-nucleotide polymorphisms (SNPs). The present study aimed to test the effect of chip density in detecting ROHom and ROHet in the Italian Simmental cattle breed. A sample of 897 animals were genotyped at low density (50k SNP; 397 individuals), medium density (140k SNP; 348 individuals), or high density (800k SNP; 152 individuals). The number of ROHom and ROHet per animal (nROHom and nROHet, respectively) and their average length were calculated. ROHom or ROHet shared by more than one animal and the number of times a particular SNP was inside a run were also computed (SNPROHom and SNPROHet). As the chip density increased, the nROHom increased, whereas their average length decreased. In contrast, the nROHet decreased and the average length increased as the chip density increased. The most repeated ROHom harbored no genes, whereas in the most repeated ROHet four genes (SNRPN, SNURF, UBE3A, and ATP10A) previously associated with reproductive traits were found. Across the 3 datasets, 31 SNP, located on Bos taurus autosome (BTA) 6, and 37 SNP (located on BTA21) exceeded the 99th percentile in the distribution of the SNPROHom and SNPROHet, respectively. The genomic region on BTA6 mapped the SLIT2, PACRGL, and KCNIP4 genes, whereas 19 and 18 genes were mapped on BTA16 and BTA21, respectively. Interestingly, most of genes found through the ROHet analysis were previously reported to be related to health, reproduction, and fitness traits. The results of the present study confirm that the detection of ROHom is more reliable when the chip density increases, whereas the ROHet trend seems to be the opposite. Genes and quantitative trait loci (QTL) mapped in the highlighted regions confirm that ROHet can be due to balancing selection, thus related to fitness traits, health, and reproduction, whereas ROHom are mainly involved in production traits. The results of the present study strengthened the usefulness of these parameters in analyzing the genomes of livestock and their biological meaning.

Genomic Regions Associated with Milk Composition and Fertility Traits in Spring-Calved Dairy Cows in New Zealand

The objective of this study was to identify genomic regions and genes that are associated with the milk composition and fertility traits of spring-calved dairy cows in New Zealand. Phenotypic data from the 2014–2015 and 2021–2022 calving seasons in two Massey University dairy herds were used. We identified 73 SNPs that were significantly associated with 58 potential candidate genes for milk composition and fertility traits. Four SNPs on chromosome 14 were highly significant for both fat and protein percentages, and the associated genes were DGAT1, SLC52A2, CPSF1, and MROH1. For fertility traits, significant associations were detected for intervals from the start of mating to first service, the start of mating to conception, first service to conception, calving to first service, and 6-wk submission, 6-wk in-calf, conception to first service in the first 3 weeks of the breeding season, and not in calf and 6-wk calving rates. Gene Ontology revealed 10 candidate genes (KCNH5, HS6ST3, GLS, ENSBTAG00000051479, STAT1, STAT4, GPD2, SH3PXD2A, EVA1C, and ARMH3) that were significantly associated with fertility traits. The biological functions of these genes are related to reducing the metabolic stress of cows and increasing insulin secretion during the mating period, early embryonic development, foetal growth, and maternal lipid metabolism during the pregnancy period.

The study of candidate genes in the improvement of egg production in ducks – a review

Duck is the second-largest poultry species aside from chicken. The rate of egg production is a major determinant of the economic income of poultry farmers. Among the reproductive organs, the ovary is a major part of the female reproductive system which is highly important for egg production. Based on the importance of this organ, several studies have been carried out to identify candidate genes at the transcriptome level, and also the expression level of these genes at different tissues or egg-laying conditions, and single nucleotide polymorphism (SNPs) of genes associated with egg production in duck. In this review, expression profile and association study analyses at SNPs level of different candidate genes with egg production traits of duck were highlighted. Furthermore, different studies on transcriptome analysis, Quantitative Trait Loci (QTL) mapping, and Genome Wide Association Study (GWAS) approach used to identify potential candidate genes for egg production in ducks were reported. This review would widen our knowledge on molecular markers that are associated or have a positive correlation to improving egg production in ducks, for the increasing world populace.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Genome-wide association study with imputed whole-genome sequence variants including large deletions for female fertility in 3 Nordic dairy cattle breeds

Fertility is an economically important trait in livestock. Poor fertility in dairy cattle can be due to loss-of-function variants affecting any essential gene that causes early embryonic mortality in homozygotes. To identify fertility-associated quantitative trait loci, we performed single-marker association analyses for 8 fertility traits in Holstein, Jersey, and Nordic Red Dairy cattle using imputed whole-genome sequence variants including SNPs, indels, and large deletion. We then performed stepwise selection of independent markers from GWAS loci using conditional and joint association analyses. From single-marker analyses for fertility traits, we reported genome-wide significant associations of 30,384 SNPs, 178 indels, and 3 deletions in Holstein; 23,481 SNPs, 189 indels, and 13 deletions in Nordic Red; and 17 SNPs in Jersey cattle. Conditional and joint association analyses identified 37 and 23 independent associations in Holstein and Nordic Red Dairy cattle, respectively. Fertility-associated GWAS loci were enriched for developmental and cellular processes (Gene Ontology enrichment, false discovery rate < 0.05). For these quantitative trait loci regions (top marker and 500 kb of surrounding regions), we proposed several candidate genes with functional annotations corresponding to embryonic lethality and various fertility-related phenotypes in mouse and cattle. The inclusion of these top markers in future releases of the custom SNP chip used for genomic evaluations will enable their validation in independent populations and improve the accuracy of genomic predictions.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Fertility-Associated Polymorphism within Bovine ITGβ5 and Its Significant Correlations with Ovarian and Luteal Traits

Simple Summary

The fertility of bovines is essential for cattle husbandry. ITGβ5, which is suggested to be closely related to fertility, is known to mediate cell adhesion and affect a variety of cellular activities. To investigate the relationship between the ITGβ5 gene and the fertility of bovines, 696 ovarian samples were collected and six potential indel (insertion/deletion) within ITGβ5 were analysed, from which a deletion mutation was found to be polymorphic. The genotype frequency and allele frequency of this locus in the investigated population were calculated and the population genetic parameters analyzed. In addition, this locus was found to be significantly correlated with ovarian width and corpus luteum diameter. Considering the importance of ovary and corpus luteum in reproduction, it is tempting to speculate the crucial effects of ITGβ5 on bovine fertility, which still need further validation. The results of our study might provide a theoretical basis for future breeding to enhance bovine reproduction.

Abstract

There is an urgent need to improve bovine fertility, and molecular marker-assisted selection (MAS) can accelerate this process. Genome-wide association studies suggest that Integrin β5 (ITGβ5) might affect fertility in bovines. As a member of the integrins family, ITGβ5 can bind to the extracellular matrix and mediate various cellular processes. In our study, primers spanning six potential insertion/deletion (indel) polymorphisms within the ITGβ5 gene were designed and 696 ovary samples from different individuals, the vast majority not in oestrum were collected for genetic variation detection. A deletion locus, rs522759246, namely P1-D_13-bp, was found to be polymorphic. The allele D frequency was 0.152 and the polymorphism information content (PIC) value was 0.224, indicating a low-degree PIC. This locus did not follow the Hardy–Weinberg equilibrium (p = 1.200E-23). Importantly, associations between P1-D_13-bp and ovarian morphological traits were established. Polymorphisms of this locus had significant correlations with ovarian width (p = 0.015). The corpus luteum is also linked to fertility and P1-D_13-bp was significantly correlated with corpus luteum diameter (p = 0.005). In conclusion, an indel mutation within the bovine ITGβ5 gene was identified, which was significantly associated with several ovarian and luteal traits.

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.

Identification of Major Loci and Candidate Genes for Meat Production-Related Traits in Broilers

Background

Carcass traits are crucial characteristics of broilers. However, the underlying genetic mechanisms are not well understood. In the current study, significant loci and major-effect candidate genes affecting nine carcass traits related to meat production were analyzed in 873 purebred broilers using an imputation-based genome-wide association study.

Results

The heritability estimates of nine carcass traits, including carcass weight, thigh muscle weight, and thigh muscle percentage, were moderate to high and ranged from 0.21 to 0.39. Twelve genome-wide significant SNPs and 118 suggestively significant SNPs of 546,656 autosomal variants were associated with carcass traits. All SNPs for six weight traits (body weight at 42 days of age, carcass weight, eviscerated weight, whole thigh weight, thigh weight, and thigh muscle weight) were clustered around the 24.08 Kb region (GGA24: 5.73–5.75 Mb) and contained only one candidate gene (DRD2). The most significant SNP, rs15226023, accounted for 4.85–7.71% of the estimated genetic variance of the six weight traits. The remaining SNPs for carcass composition traits (whole thigh percentage and thigh percentage) were clustered around the 42.52 Kb region (GGA3: 53.03–53.08 Mb) and contained only one candidate gene (ADGRG6). The most significant SNP in this region, rs13571431, accounted for 11.89–13.56% of the estimated genetic variance of two carcass composition traits. Some degree of genetic differentiation in ADGRG6 between large and small breeds was observed.

Conclusion

We identified one 24.08 Kb region for weight traits and one 42.52 Kb region for thigh-related carcass traits. DRD2 was the major-effect candidate gene for weight traits, and ADGRG6 was the major-effect candidate gene for carcass composition traits. Our results supply essential information for causative mutation identification of carcass traits in broilers.

Assessing the genetic background and genomic relatedness of red cattle populations originating from Northern Europe

Background

Local cattle breeds need special attention, as they are valuable reservoirs of genetic diversity. Appropriate breeding decisions and adequate genomic management of numerically smaller populations are required for their conservation. At this point, the analysis of dense genome-wide marker arrays provides encompassing insights into the genomic constitution of livestock populations. We have analyzed the genetic characterization of ten cattle breeds originating from Germany, The Netherlands and Denmark belonging to the group of red dairy breeds in Northern Europe. The results are intended to provide initial evidence on whether joint genomic breeding strategies of these populations will be successful.

Results

Traditional Danish Red and Groningen White-Headed were the most genetically differentiated breeds and their populations showed the highest levels of inbreeding. In contrast, close genetic relationships and shared ancestry were observed for the populations of German Red and White Dual-Purpose, Dutch Meuse-Rhine-Yssel, and Dutch Deep Red breeds, reflecting their common histories. A considerable amount of gene flow from Red Holstein to German Angler and to German Red and White Dual-Purpose was revealed, which is consistent with frequent crossbreeding to improve productivity of these local breeds. In Red Holstein, marked genomic signatures of selection were reported on chromosome 18, suggesting directed selection for important breeding goal traits. Furthermore, tests for signatures of selection between Red Holstein, Red and White Dual-Purpose, and Meuse-Rhine-Yssel uncovered signals for all investigated pairs of populations. The corresponding genomic regions, which were putatively under different selection pressures, harboured various genes which are associated with traits such as milk and beef production, mastitis and female fertility.

Conclusions

This study provides comprehensive knowledge on the genetic constitution and genomic connectedness of divergent red cattle populations in Northern Europe. The results will help to design and optimize breeding strategies. A joint genomic evaluation including some of the breeds studied here seems feasible.

Genome-wide association study for mammary structure in Canadian Angus cows

Functional and enduring mammary structure is pivotal for producer profitability, and animal health and welfare in beef production. Genetic evaluations for teat and udder score in Canadian Angus cattle have previously been developed. The aim of this study was to identify genomic regions associated with teat and udder structure in Canadian Angus cows thereby enhancing knowledge of the biological architecture of these traits. Thus, we performed a weighted single-step genome wide association study (WssGWAS) to identify candidate genes for teat and udder score in 1,582 Canadian Angus cows typed with the GeneSeek^® Genomic Profiler Bovine 130K SNP array. Genomically enhanced estimated breeding values (GEBVs) were converted to SNP marker effects using unequal variances for markers to calculate weights for each SNP over three iterations. At the genome wide level, we detected windows of 20 consecutive SNPs that explained more than 0.5% of the variance observed in these traits. A total of 35 and 28 windows were identified for teat and udder score, respectively, with two SNP windows in common for both traits. Using Ensembl, the SNP windows were used to search for candidate genes and quantitative trait loci (QTL). A total of 94 and 71 characterized genes were identified in the regions for teat and udder score, respectively. Of these, 7 genes were common for both traits. Gene network and enrichment analysis, using Ingenuity Pathway Analysis (IPA), signified key pathways unique to each trait. Genes of interest were associated with immune response and wound healing, adipose tissue development and morphology, and epithelial and vascular development and morphology. Genetic architecture from this GWAS confirms that teat and udder score are distinct, polygenic traits involving varying and complex biological pathways, and that genetic selection for improved teat and udder score is possible.

Association analysis of loci implied in "buffering" epistasis

The existence of buffering mechanisms is an emerging property of biological networks, and this results in the buildup of robustness through evolution. So far, there are no explicit methods to find loci implied in buffering mechanisms. However, buffering can be seen as interaction with genetic background. Here we develop this idea into a tractable model for quantitative genetics, in which the buffering effect of one locus with many other loci is condensed into a single statistical effect, multiplicative on the total additive genetic effect. This allows easier interpretation of the results and simplifies the problem of detecting epistasis from quadratic to linear in the number of loci. Using this formulation, we construct a linear model for genome-wide association studies that estimates and declares the significance of multiplicative epistatic effects at single loci. The model has the form of a variance components, norm reaction model and likelihood ratio tests are used for significance. This model is a generalization and explanation of previous ones. We test our model using bovine data: Brahman and Tropical Composite animals, phenotyped for body weight at yearling and genotyped at high density. After association analysis, we find a number of loci with buffering action in one, the other, or both breeds; these loci do not have a significant statistical additive effect. Most of these loci have been reported in previous studies, either with an additive effect or as footprints of selection. We identify buffering epistatic SNPs present in or near genes reported in the context of signatures of selection in multi-breed cattle population studies. Prominent among these genes are those associated with fertility (INHBA, TSHR, ESRRG, PRLR, and PPARG), growth (MSTN, GHR), coat characteristics (KIT, MITF, PRLR), and heat resistance (HSPA6 and HSPA1A). In these populations, we found loci that have a nonsignificant statistical additive effect but a significant epistatic effect. We argue that the discovery and study of loci associated with buffering effects allow attacking the difficult problems, among others, of the release of maintenance variance in artificial and natural selection, of quick adaptation to the environment, and of opposite signs of marker effects in different backgrounds. We conclude that our method and our results generate promising new perspectives for research in evolutionary and quantitative genetics based on the study of loci that buffer effect of other loci.

Genome-Wide Association Study and Pathway Analysis for Female Fertility Traits in Iranian Holstein Cattle

Female fertility is an important trait that contributes to cow’s profitability and it can be improved by genomic information. The objective of this study was to detect genomic regions and variants affecting fertility traits in Iranian Holstein cattle. A data set comprised of female fertility records and 3,452,730 pedigree information from Iranian Holstein cattle were used to predict the breeding values, which were then employed to estimate the de-regressed proofs (DRP) of genotyped animals. A total of 878 animals with DRP records and 54k SNP markers were utilized in the genome-wide association study (GWAS). The GWAS was performed using a linear regression model with SNP genotype as a linear covariate. The results showed that an SNP on BTA19, ARS-BFGL-NGS-33473, was the most significant SNP associated with days from calving to first service. In total, [69] significant SNPs were located within 27 candidate genes. Novel potential candidate genes include OSTN, DPP6, EphA5, CADPS2, Rfc1, ADGRB3, Myo3a, C10H14orf93, KIAA1217, RBPJL, SLC18A2, GARNL3, NCALD, ASPH, ASIC2, OR3A1, CHRNB4, CACNA2D2, DLGAP1, GRIN2A and ME3. These genes are involved in different pathways relevant to female fertility and other characteristics in mammals. Gene set enrichment analysis showed that thirteen GO terms had significant overrepresentation of genes statistically associated with female fertility traits. The results of network analysis identified CCNB1 gene as a hub gene in the progesterone-mediated oocyte maturation pathway, significantly associated with age at first calving. The candidate genes identified in this study can be utilized in genomic tests to improve reproductive performance in Holstein cattle.

Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle

The present study was carried out in Tharparkar cattle for identification of genome-wide SNPs and microsatellites, and then annotate the identified high-quality SNPs to milk production, fertility, carcass, adaptability and immune response of economically important traits. A total of 146,011 SNPs were identified with respect to Bos taurus reference genome which are indicus specific, out of which 10,519 SNPs were found to be novel. Similarly, a total of 87,047 SNPs were identified with respect to Bos indicus reference genome. After final annotation of SNPs identified with respect to Bos indicus reference genome, 2871 SNPs were found to be associated in 383 candidate genes having to do with milk production, fertility, carcass, immune response and adaptability traits. Following that, 2571 microsatellites were identified. The information mined from the data might be of importance for the future breed improvement programs, conservation efforts and for enhancing the SNPs density of the existing bovine SNP chips.

New insights into the associations among feed efficiency, metabolizable efficiency traits and related QTL regions in broiler chickens

Background

Improving the feed efficiency would increase profitability for producers while also reducing the environmental footprint of livestock production. This study was conducted to investigate the relationships among feed efficiency traits and metabolizable efficiency traits in 180 male broilers. Significant loci and genes affecting the metabolizable efficiency traits were explored with an imputation-based genome-wide association study. The traits measured or calculated comprised three growth traits, five feed efficiency related traits, and nine metabolizable efficiency traits.

Results

The residual feed intake (RFI) showed moderate to high and positive phenotypic correlations with eight other traits measured, including average daily feed intake (ADFI), dry excreta weight (DEW), gross energy excretion (GEE), crude protein excretion (CPE), metabolizable dry matter (MDM), nitrogen corrected apparent metabolizable energy (AMEn), abdominal fat weight (AbF), and percentage of abdominal fat (AbP). Greater correlations were observed between growth traits and the feed conversion ratio (FCR) than RFI. In addition, the RFI, FCR, ADFI, DEW, GEE, CPE, MDM, AMEn, AbF, and AbP were lower in low-RFI birds than high-RFI birds (P < 0.01 or P < 0.05), whereas the coefficients of MDM and MCP of low-RFI birds were greater than those of high-RFI birds (P < 0.01). Five narrow QTLs for metabolizable efficiency traits were detected, including one 82.46-kb region for DEW and GEE on Gallus gallus chromosome (GGA) 26, one 120.13-kb region for MDM and AMEn on GGA1, one 691.25-kb region for the coefficients of MDM and AMEn on GGA5, one region for the coefficients of MDM and MCP on GGA2 (103.45–103.53 Mb), and one 690.50-kb region for the coefficient of MCP on GGA14. Linkage disequilibrium (LD) analysis indicated that the five regions contained high LD blocks, as well as the genes chromosome 26 C6orf106 homolog (C26H6orf106), LOC396098, SH3 and multiple ankyrin repeat domains 2 (SHANK2), ETS homologous factor (EHF), and histamine receptor H3-like (HRH3L), which are known to be involved in the regulation of neurodevelopment, cell proliferation and differentiation, and food intake.

Conclusions

Selection for low RFI significantly decreased chicken feed intake, excreta output, and abdominal fat deposition, and increased nutrient digestibility without changing the weight gain. Five novel QTL regions involved in the control of metabolizable efficiency in chickens were identified. These results, combined through nutritional and genetic approaches, should facilitate novel insights into improving feed efficiency in poultry and other species.

Genome-Wide Association Analysis Reveals Key Genes Responsible for Egg Production of Lion Head Goose

The lion head goose is one of the most important agricultural resources in China; however, its breeding process is relatively slow. In the present study, a genome-wide association study was performed for the genetic selection of egg production characters in lion head geese. We detected 30 single-nucleotide polymorphisms located in or near 30 genes that might be associated with egg production character, and quantitative real-time polymerase chain reaction was used to verify their expression level in lion head geese. The results showed that the expression levels of CRTC1 (encoding CREB-regulated transcription coactivator 1), FAAH2 (encoding fatty acid amide hydrolase 2), GPC3 (encoding glypican 3), and SERPINC1 (encoding serpin family C member 1) in high egg production population were significantly lower than those in the low egg production populations (*P < 0.05). The expression levels of CLPB (encoding caseinolytic peptidase B protein homolog), GNA12 (encoding guanine nucleotide-binding protein subunit alpha-12), and ZMAT5 (encoding zinc finger, matrin type 5) in the high egg production population were significantly higher than those in the low egg production populations (*P < 0.05). The expression of BMP4 (encoding bone morphogenetic protein 4), FRMPD3 (encoding FERM and PDZ domain containing 3), LIF (encoding leukemia inhibitory factor), and NFYC (encoding nuclear transcription factor Y subunit gamma) in the high egg production population were very significantly lower than those in the low egg production population (**P < 0.01). Our findings provide an insight into the economic traits of lion head goose. These candidate genes might be valuable for future breeding improvement.

Genome-wide association studies for the number of animals born alive and dead in duroc pigs

Litter size is one of the most important economic traits for pig production as it is directly related to the production efficiency. As an important litter size trait in pigs, the number of piglets born alive at birth (NBA) receives widespread interests in the pig industry. However, traits of piglets born dead, including the number of stillborn piglets (NS) and the piglets mummified at birth (NM) should be noted to explain the loss of reproduction. Herein, in the present study, a total of 803 producing sows were sampled and 2807 farrowing records for NBA, NM, and NS traits were collected in a Duroc swine population. Subsequently, a genome-wide association study (GWAS) was performed for NBA, NS and NM in parity groups 1 to 5. In total, 10 putative regions were found associated with these traits. After stepwise conditional analyses around the putative regions, eight independent signals were ultimately identified for NBA, NS, and NM, and there were seven promising candidate genes related to these traits, including ARID1A, RXRG, NFATC4, ABTB2, GRAMD1B, NDRG1, and APC. Our findings contribute to the understanding of the significant genetic causes of piglets born alive and dead, and could have a positive effect on pig production efficiency and economic profits.

The microRNA gene <i>bta-mir-2313</i> in cattle: an atlas of regulatory elements and an association analysis with growth and carcass traits in the Slovenian Simental cattle breed

Abstract. MicroRNAs (miRNA) are a class of non-coding RNAs important in posttranscriptional regulation of target genes. Regulation requires complementarity between the target mRNA and the miRNA region responsible for their recognition and binding, also called the seed region. Previous studies have proven that expression profiles and genetic variations of miRNA genes (miR-SNP; SNP – single nucleotide polymorphism) and their target sites (miR-TS-SNPs) have an impact on phenotypic variation and disease susceptibility in human, animal models, and livestock. MicroRNA-associated polymorphisms therefore represent biomarker potential for phenotypic traits in livestock. Effects of miRNA gene polymorphisms on phenotypic traits have been studied in several animal species but much less in cattle. The aim of the present study was therefore to analyze the genetic variability in the bta-mir-2313 gene and test associations with growth and carcass traits of the Slovenian Simmental cattle breed. Additionally, validated and predicted genomic information related to the miRNA gene bta-mir-2313 has been obtained and presented as an atlas of miRNA regulatory elements. Sanger sequencing has been used for biomarker development and genotyping of 145 animals of Slovenian dual-purpose Simmental cattle. Out of nine known polymorphisms located within pre-miRNA regions, one mature miRNA seed SNP was polymorphic in the Slovenian Simmental cattle breed. An additional three polymorphisms were identified within the flanking pri-miRNA regions. There was no significant effect of polymorphisms on 18 tested fattening and carcass traits; however, validated polymorphisms could now be tested in association with other traits in other cattle populations. The microRNA gene bta-mir-2313 warrants further genetic and functional analyses since it overlaps with a large number of quantitative trait loci (QTL), has over 3100 predicted targets and highly polymorphic mature seed regions, and is located within protein-coding gene GRAMD1B, previously associated with production traits in cattle. Mature miRNA seed SNPs present important genomic loci for functional studies because they could affect the gain/loss of downstream targets and should be systematically studied in cattle.

Comparative transcriptomic analysis of high and low egg-producing duck ovaries

The egg-laying rate is an important indicator of egg production of laying ducks. Egg production directly impacts the economic benefits of the duck industry. In order to obtain better insight into the molecular mechanisms associated with the process of egg production, comparative transcriptomic analysis of the ovaries of Jinding ducks with high and low egg production was performed using the Illumina HiSeq 2500 system. A total of 843 differentially expressed genes (DEGs) was identified, 367 that were down-regulated and 476 that were up-regulated in high egg production (HEP) ovaries, as compared with low egg production (LEP) ovaries. Some genes, such as MC5R, APOD, ORAI1, and DYRK4, were more active in HEP ovaries, indicating that these genes may play important roles in regulation of egg production. Among these 843 DEGs, 685 were assigned to gene ontology (GO) categories. Of these, 25 genes were related to reproduction, and 30 were related to the reproductive process, including some associated with ovarian follicle development, circadian regulation of gene expression, circadian rhythm, and estrogen receptor binding. Furthermore, some important functional pathways were revealed, such as the steroid biosynthesis pathway, the endocrine and other factor-regulated calcium reabsorption pathways, circadian rhythm, the neuroactive ligand-receptor interaction pathway, fatty acid biosynthesis, and the calcium-signaling pathway, which appear to be much more active in the HEP group, as compared to those of the LEP group. The results of this study provide very useful information that may contribute to future functional studies of genes involved in bird reproduction.

Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle

Cattle internal organs as accessible raw materials have a long history of being widely used in beef processing, feed and pharmaceutical industry. These traits not only are of economic interest to breeders, but they are intrinsically linked to many valuable traits, such as growth, health, and productivity. Using the Illumina Bovine HD 770K SNP array, we performed a genome-wide association study for heart weight, liver weight, spleen weight, lung weight, and kidney weight in 1,217 Simmental cattle. In our research, 38 significant single nucleotide polymorphisms (SNPs) ( P < 1.49 × 10−6) were identified for five internal organ weight traits. These SNPs are within or near 13 genes, and some of them have been reported previously, including NDUFAF4, LCORL, BT.94996, SLIT2, FAM184B, LAP3, BBS12, MECOM, CD300LF, HSD17B3, TLR4, MXI1, and MB21D2. In addition, we detected four haplotype blocks on BTA6 containing 18 significant SNPs associated with spleen weight. Our results offer worthy insights into understanding the genetic mechanisms of internal organs' development, with potential application in breeding programs of Simmental beef cattle.

Human-Mediated Introgression of Haplotypes in a Modern Dairy Cattle Breed

Domestic animals can serve as model systems of adaptive introgression and their genomic signatures. In part, their usefulness as model systems is due to their well-known histories. Different breeding strategies such as introgression and artificial selection have generated numerous desirable phenotypes and superior performance in domestic animals. The modern Danish Red Dairy Cattle is studied as an example of an introgressed population. It originates from crossing the traditional Danish Red Dairy Cattle with the Holstein and Brown Swiss breeds, both known for high milk production. This crossing happened, among other things due to changes in the production system, to raise milk production and overall performance. The genomes of modern Danish Red Dairy Cattle are heavily influenced by regions introgressed from the Holstein and Brown Swiss breeds and under subsequent selection in the admixed population. The introgressed proportion of the genome was found to be highly variable across the genome. Haplotypes introgressed from Holstein and Brown Swiss contained or overlapped known genes affecting milk production, as well as protein and fat content (CD14, ZNF215, BCL2L12, and THRSP for Holstein origin and ITPR2, BCAT1, LAP3, and MED28 for Brown Swiss origin). Genomic regions with high introgression signals also contained genes and enriched QTL associated with calving traits, body confirmation, feed efficiency, carcass, and fertility traits. These introgressed signals with relative identity-by-descent scores larger than the median showing Holstein or Brown Swiss introgression are mostly significantly correlated with the corresponding test statistics from signatures of selection analyses in modern Danish Red Dairy Cattle. Meanwhile, the putative significant introgressed signals have a significant dependency with the putative significant signals from signatures of selection analyses. Artificial selection has played an important role in the genomic footprints of introgression in the genome of modern Danish Red Dairy Cattle. Our study on a modern cattle breed contributes to an understanding of genomic consequences of selective introgression by demonstrating the extent to which adaptive effects contribute to shape the specific genomic consequences of introgression.

Meta-analysis of sequence-based association studies across three cattle breeds reveals 25 QTL for fat and protein percentages in milk at nucleotide resolution

Background

Genotyping and whole-genome sequencing data have been generated for hundreds of thousands of cattle. International consortia used these data to compile imputation reference panels that facilitate the imputation of sequence variant genotypes for animals that have been genotyped using dense microarrays. Association studies with imputed sequence variant genotypes allow for the characterization of quantitative trait loci (QTL) at nucleotide resolution particularly when individuals from several breeds are included in the mapping populations.

Results

We imputed genotypes for 28 million sequence variants in 17,229 cattle of the Braunvieh, Fleckvieh and Holstein breeds in order to compile large mapping populations that provide high power to identify QTL for milk production traits. Association tests between imputed sequence variant genotypes and fat and protein percentages in milk uncovered between six and thirteen QTL (P < 1e-8) per breed. Eight of the detected QTL were significant in more than one breed. We combined the results across breeds using meta-analysis and identified a total of 25 QTL including six that were not significant in the within-breed association studies. Two missense mutations in the ABCG2 (p.Y581S, rs43702337, P = 4.3e-34) and GHR (p.F279Y, rs385640152, P = 1.6e-74) genes were the top variants at QTL on chromosomes 6 and 20. Another known causal missense mutation in the DGAT1 gene (p.A232K, rs109326954, P = 8.4e-1436) was the second top variant at a QTL on chromosome 14 but its allelic substitution effects were inconsistent across breeds. It turned out that the conflicting allelic substitution effects resulted from flaws in the imputed genotypes due to the use of a multi-breed reference population for genotype imputation.

Conclusions

Many QTL for milk production traits segregate across breeds and across-breed meta-analysis has greater power to detect such QTL than within-breed association testing. Association testing between imputed sequence variant genotypes and phenotypes of interest facilitates identifying causal mutations provided the accuracy of imputation is high. However, true causal mutations may remain undetected when the imputed sequence variant genotypes contain flaws. It is highly recommended to validate the effect of known causal variants in order to assess the ability to detect true causal mutations in association studies with imputed sequence variants.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4263-8) contains supplementary material, which is available to authorized users.

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.

Synergies between assisted reproduction technologies and functional genomics

This review, is a synopsis of advanced reproductive technologies in farm animals, including the discussion of their limiting factors as revealed by the study of offspring derived from embryos produced in vitro and through cloning. These studies show that the problems of epigenetic mis-programming, which were reported in the initial stages of assisted reproduction, still persist. The importance of whole-genome analyses, including the methylome and transcriptome, in improving embryo biotechnologies in farm animals, are discussed. Genome editing approaches for the improvement of economically-relevant traits in farm animals are also described. Efficient farm animal embryo biotechnologies, including cloning and the most recent technologies such as genome editing, will effectively complement the latest strategies to accelerate genetic improvement of farm animals.

Genome-wide association analysis of milk yield traits in Nordic Red Cattle using imputed whole genome sequence variants

BACKGROUND

The Nordic Red Cattle consisting of three different populations from Finland, Sweden and Denmark are under a joint breeding value estimation system. The long history of recording of production and health traits offers a great opportunity to study production traits and identify causal variants behind them. In this study, we used whole genome sequence level data from 4280 progeny tested Nordic Red Cattle bulls to scan the genome for loci affecting milk, fat and protein yields.

RESULTS

Using a genome-wise significance threshold, regions on Bos taurus chromosomes 5, 14, 23, 25 and 26 were associated with fat yield. Regions on chromosomes 5, 14, 16, 19, 20 and 25 were associated with milk yield and chromosomes 5, 14 and 25 had regions associated with protein yield. Significantly associated variations were found in 227 genes for fat yield, 72 genes for milk yield and 30 genes for protein yield. Ingenuity Pathway Analysis was used to identify networks connecting these genes displaying significant hits. When compared to previously mapped genomic regions associated with fertility, significantly associated variations were found in 5 genes common for fat yield and fertility, thus linking these two traits via biological networks.

CONCLUSION

This is the first time when whole genome sequence data is utilized to study genomic regions affecting milk production in the Nordic Red Cattle population. Sequence level data offers the possibility to study quantitative traits in detail but still cannot unambiguously reveal which of the associated variations is causative. Linkage disequilibrium creates difficulties to pinpoint the causative genes and variations. One solution to overcome these difficulties is the identification of the functional gene networks and pathways to reveal important interacting genes as candidates for the observed effects. This information on target genomic regions may be exploited to improve genomic prediction.