Genome-wide association study and annotating candidate gene networks affecting age at first calving in Nellore cattle

Identification of selection signatures and genetic diversity in the sheep

In the study, data obtained from OvineSNP50K SNP chips using the Illumina® iScan platform for Eşme sheep were used. The integrated haplotype score (iHS) and runs of homozygosity (ROH) statistical approaches were used to identify selection signatures. Using the iHS analysis, it was discovered that there are 10 genomic regions and 51 genes on ovine chromosomes 1, 9, 11, and 12 that are under selection. Three genomic regions and 97 genes on ovine chromosomes 6 and 11 were found to be under selection using the ROH analysis. Candidate genes associated with economic and ecological traits were detected using both approaches. Among the genetic diversity parameters considered in this study, the minor allele frequency (MAF), the genetic distance between individuals (D), as well as observed (Ho) and expected heterozygosities (He) values were 0.300, 0.309, 0.388, and 0.390, respectively. The obtained Ho, He and D values indicate a moderate level of genetic diversity. The ratio of polymorphic SNPs (PN) was 0.947, and the average values of FROH and FHOM were 0.030 and 0.029, respectively. Considering the PN value obtained in the study, it is evident that the SNPs in the population exhibit a high level of polymorphism at 94.7%. While the FROH value obtained indicates high genetic diversity among the individuals in the present study, the FHOM value suggests that the population is predominantly composed of heterozygous individuals. As a result, evidence indicating genetic advancements have been made for target traits in breeding programs within the population. Additionally, candidate genes suitable for future molecular marker-supported breeding programs have been identified. In addition, a better understanding of the genetic structure and production potential of the population has been achieved. Findings have shown that Eşme sheep are a breed with high meat production potential and strong adaptation abilities.

Whole-genome resequencing deciphers patterns of genetic diversity, phylogeny, and evolutionary dynamics in Kashmir cattle

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

Genome-wide association study identifies genomic regions associated with key reproductive traits in Korean Hanwoo cows

BACKGROUND

Conducting genome-wide association studies (GWAS) for reproductive traits in Hanwoo cattle, including age at first calving (AFC), calving interval (CI), gestation length (GL), and number of artificial inseminations per conception (NAIPC), is of paramount significance. These analyses provided a thorough exploration of the genetic basis of these traits, facilitating the identification of key markers for targeted trait improvement. Breeders can optimize their selection strategies, leading to more efficient and sustainable breeding programs, by incorporating genetic insights. This impact extends beyond individual traits and contributes to the overall productivity and profitability of the Hanwoo beef cattle industry. Ultimately, GWAS is essential in ensuring the long-term genetic resilience and adaptability of Hanwoo cattle populations. The primary goal of this study was to identify significant single nucleotide polymorphisms (SNPs) or quantitative trait loci (QTLs) associated with the studied reproductive traits and subsequently map the underlying genes that hold promise for trait improvement.

RESULTS

A genome-wide association study of reproductive traits identified 68 significant single nucleotide polymorphisms (SNPs) distributed across 29 Bos taurus autosomes (BTA). Among them, BTA14 exhibited the highest number of identified SNPs (25), whereas BTA6, BTA7, BTA8, BTA10, BTA13, BTA17, and BTA20 exhibited 8, 5, 5, 3, 8, 2, and 12 significant SNPs, respectively. Annotation of candidate genes within a 500 kb region surrounding the significant SNPs led to the identification of ten candidate genes relevant to age at first calving. These genes were: FANCG, UNC13B, TESK1, TLN1, and CREB3 on BTA8; FAM110B, UBXN2B, SDCBP, and TOX on BTA14; and MAP3K1 on BTA20. Additionally, APBA3, TCF12, and ZFR2, located on BTA7 and BTA10, were associated with the calving interval; PAX1, SGCD, and HAND1, located on BTA7 and BTA13, were linked to gestation length; and RBM47, UBE2K, and GPX8, located on BTA6 and BTA20, were linked to the number of artificial inseminations per conception in Hanwoo cows.

CONCLUSIONS

The findings of this study enhance our knowledge of the genetic factors that influence reproductive traits in Hanwoo cattle populations and provide a foundation for future breeding strategies focused on improving desirable traits in beef cattle. This research offers new evidence and insights into the genetic variants and genome regions associated with reproductive traits and contributes valuable information to guide future efforts in cattle breeding.

Understanding the underlying genetic mechanisms for age at first calving, inter-calving period and scrotal circumference in Bonsmara cattle

BACKGROUND

Reproduction is a key feature of the sustainability of a species and thus represents an important component in livestock genetic improvement programs. Most reproductive traits are lowly heritable. In order to gain a better understanding of the underlying genetic basis of these traits, a genome-wide association was conducted for age at first calving (AFC), first inter-calving period (ICP) and scrotal circumference (SC) within the South African Bonsmara breed. Phenotypes and genotypes (120,692 single nucleotide polymorphisms (SNPs) post editing) were available on 7,128 South African Bonsmara cattle; the association analyses were undertaken using linear mixed models.

RESULTS

Genomic restricted maximum likelihood analysis of the 7,128 SA Bonsmara cattle yielded genomic heritability's of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC. A total of 16, 23 and 51 suggestive (P ≤ 4 × 10-6) SNPs were associated with AFC, ICP and SC, while 11, 11 and 44 significant (P ≤ 4 × 10-7) SNPs were associated with AFC, ICP and SC respectively. A total of 11 quantitative trait loci (QTL) and 11 candidate genes were co-located with these associated SNPs for AFC, with 10 QTL harbouring 11 candidate genes for ICP and 41 QTL containing 40 candidate genes for SC. The QTL identified were close to genes previously associated with carcass, fertility, growth and milk-related traits. The biological pathways influenced by these genes include carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation.

CONCLUSIONS

This was the first attempt to study the underlying polymorphisms associated with reproduction in South African beef cattle. Genes previously reported in cattle breeds for numerous traits bar AFC, ICP or SC were detected in this study. Over 20 different genes have not been previously reported in beef cattle populations and may have been associated due to the unique genetic composite background of the SA Bonsmara breed.

Genetic improvement of economic traits in Murrah buffalo using significant SNPs from genome-wide association study

GWAS helps to identify QTL and candidate genes of specific traits. Buffalo breeding has primarily focused on milk production, but its negative correlation with reproduction traits resulted in unfavorable decline of reproductive performance among buffaloes. A genome wide scan was performed on a total of 120 Murrah buffaloes genotyped by ddRAD sequencing for 13 traits related to female fertility, production, and growth. The identified 25 significant single nucleotide polymorphisms (SNPs) (P <1×10⁶) are associated with age at first calving (AFC), age at first service (AFS), period from calving to 1^st Artifical Insemination (AI), service period (SP) and 6 month body weight (6M). Fifteen genetic variants overlapped with different QTL regions of reported studies. Among the associated loci, outstanding candidate genes for fertility, including AQP1, TRNAE-CUC, NRIP1, CPNE4, and VOPP1, have effect in different fertility traits. AQP1 gene is expressed in ovulatory phase and various stages of pregnancy. TRNAE-CUC gene is associated with AFC and number . of calvings after 4 years of age. Glycogen content-associated gene CPNE4 regulates muscle glycogen and is upregulated during early pregnancy. NRIP1 generegulates ovulation, corpus luteum at pregnancy, and mammary gland development. The objective is to identify potential genomic regions and genetic variants associated with economic traits and to select the most significant SNP which have positive effect on all the traits.

Analysis of genetic diversity and selection signals in Chaling cattle of southern China using whole-genome scan

China has diversified resources of indigenous cattle, which are classified into Northern, Central, and Southern groups according to their geographical distribution. Chaling cattle belong to Southern group. This breed is famous for the production of good quality meat with elite meat grades. To analyze the genetic diversity of Chaling cattle, 20 samples were sequenced using whole-genome resequencing technology, along with 138 published whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle as control. It was found that Chaling cattle originated from Chinese indicine cattle. The genetic diversity of Chaling cattle is higher than that of Indian indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle, but lower than that of Chinese indicine cattle and Xiangxi cattle. Annotating the selection signals obtained by composite likelihood ratio, θπ, F_ST , π-ratio, and XP-EHH methods, several genes associated with immunity, heat tolerance, reproduction, growth, and meat quality showed strong selection signals. In general, this study provides a theoretical basis for analyzing the genetic mechanism of Chaling cattle with excellent adaptability, rough feeding tolerance, good immune performance, and good meat quality. This work lays a foundation for genetic breeding of Chaling cattle in future.

Progress and opportunities through use of genomics in animal production

The rearing of farmed animals is a vital component of global food production systems, but its impact on the environment, human health, animal welfare, and biodiversity is being increasingly challenged. Developments in genetic and genomic technologies have had a key role in improving the productivity of farmed animals for decades. Advances in genome sequencing, annotation, and editing offer a means not only to continue that trend, but also, when combined with advanced data collection, analytics, cloud computing, appropriate infrastructure, and regulation, to take precision livestock farming (PLF) and conservation to an advanced level. Such an approach could generate substantial additional benefits in terms of reducing use of resources, health treatments, and environmental impact, while also improving animal health and welfare.

Genome-Wide Association Study of Parasite Resistance to Gastrointestinal Nematodes in Corriedale Sheep

Selection of genetically resistant animals is one alternative to reduce the negative impact of gastrointestinal nematodes (GIN) on sheep production. The aim of this study was to identify genomic regions associated with GIN resistance in Corriedale sheep by single-step genome-wide association studies (ssGWAS) using 170, 507 and 50K single nucleotide polymorphisms (SNPs). Analysis included 19,547 lambs with faecal egg counts (FEC) records, a pedigree file of 40,056 animals and 454, 711 and 383 genotypes from 170, 507 and 50K SNPs, respectively. Genomic estimated breeding values (GEBV) were obtained with single-step genomic BLUP methodology (ssGBLUP), using a univariate animal model, which included contemporary group, type of birth and age of dam as class fixed effects and age at FEC recording as covariate. The SNP effects as wells as p-values were estimated with POSTGSF90 program. Significance level was defined by a chromosome-wise False Discovery Rate of 5%. Significant genomic regions were identified in chromosomes 1, 3, 12 and 19 with the 170 SNP set, in chromosomes 7, 12 and 24 using the 507 SNP chip and only in chromosome 7 with the 50K SNP chip. Candidate genes located in these regions, using Oar_v4.0 as reference genome, were TIMP3, TLR5, LEPR and TLR9 (170 SNPs), SYNDIG1L and MGRN1 (507 SNP chip) and INO80, TLN2, TSHR and EIF2AK4 (50K SNP chip). These results validate genomic regions associated with FEC previously identified in Corriedale and other breeds and report new candidate regions for further investigation.

A Random Forest-Based Genome-Wide Scan Reveals Fertility-Related Candidate Genes and Potential Inter-Chromosomal Epistatic Regions Associated With Age at First Calving in Nellore Cattle

This study aimed to perform a genome-wide association analysis (GWAS) using the Random Forest (RF) approach for scanning candidate genes for age at first calving (AFC) in Nellore cattle. Additionally, potential epistatic effects were investigated using linear mixed models with pairwise interactions between all markers with high importance scores within the tree ensemble non-linear structure. Data from Nellore cattle were used, including records of animals born between 1984 and 2015 and raised in commercial herds located in different regions of Brazil. The estimated breeding values (EBV) were computed and used as the response variable in the genomic analyses. After quality control, the remaining number of animals and SNPs considered were 3,174 and 360,130, respectively. Five independent RF analyses were carried out, considering different initialization seeds. The importance score of each SNP was averaged across the independent RF analyses to rank the markers according to their predictive relevance. A total of 117 SNPs associated with AFC were identified, which spanned 10 autosomes (2, 3, 5, 10, 11, 17, 18, 21, 24, and 25). In total, 23 non-overlapping genomic regions embedded 262 candidate genes for AFC. Enrichment analysis and previous evidence in the literature revealed that many candidate genes annotated close to the lead SNPs have key roles in fertility, including embryo pre-implantation and development, embryonic viability, male germinal cell maturation, and pheromone recognition. Furthermore, some genomic regions previously associated with fertility and growth traits in Nellore cattle were also detected in the present study, reinforcing the effectiveness of RF for pre-screening candidate regions associated with complex traits. Complementary analyses revealed that many SNPs top-ranked in the RF-based GWAS did not present a strong marginal linear effect but are potentially involved in epistatic hotspots between genomic regions in different autosomes, remarkably in the BTAs 3, 5, 11, and 21. The reported results are expected to enhance the understanding of genetic mechanisms involved in the biological regulation of AFC in this cattle breed.

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.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

Genomic legacy of migration in endangered caribou

Wide-ranging animals, including migratory species, are significantly threatened by the effects of habitat fragmentation and habitat loss. In the case of terrestrial mammals, this results in nearly a quarter of species being at risk of extinction. Caribou are one such example of a wide-ranging, migratory, terrestrial, and endangered mammal. In populations of caribou, the proportion of individuals considered as "migrants" can vary dramatically. There is therefore a possibility that, under the condition that migratory behavior is genetically determined, those individuals or populations that are migratory will be further impacted by humans, and this impact could result in the permanent loss of the migratory trait in some populations. However, genetic determination of migration has not previously been studied in an endangered terrestrial mammal. We examined migratory behavior of 139 GPS-collared endangered caribou in western North America and carried out genomic scans for the same individuals. Here we determine a genetic subdivision of caribou into a Northern and a Southern genetic cluster. We also detect >50 SNPs associated with migratory behavior, which are in genes with hypothesized roles in determining migration in other organisms. Furthermore, we determine that propensity to migrate depends upon the proportion of ancestry in individual caribou, and thus on the evolutionary history of its migratory and sedentary subspecies. If, as we report, migratory behavior is influenced by genes, caribou could be further impacted by the loss of the migratory trait in some isolated populations already at low numbers. Our results indicating an ancestral genetic component also suggest that the migratory trait and their associated genetic mutations could not be easily re-established when lost in a population.

Sustainable Intensification of Beef Production in the Tropics: The Role of Genetically Improving Sexual Precocity of Heifers

Simple Summary

Tropical pasture-based beef production systems play a vital role in global food security. The importance of promoting sustainable intensification of such systems has been debated worldwide. Demand for beef is growing together with concerns over the impact of its production on the environment. Implementing sustainable livestock intensification programs relies on animal genetic improvement. In tropical areas, the lack of sexual precocity is a bottleneck for cattle efficiency, directly impacting the sustainability of production systems. In the present review we present and discuss the state of the art of genetic evaluation for sexual precocity in Bos indicus beef cattle, covering the definition of measurable traits, genetic parameter estimates, genomic analyses, and a case study of selection for sexual precocity in Nellore breeding programs.

Abstract

Increasing productivity through continued animal genetic improvement is a crucial part of implementing sustainable livestock intensification programs. In Zebu cattle, the lack of sexual precocity is one of the main obstacles to improving beef production efficiency. Puberty-related traits are complex, but large-scale data sets from different “omics” have provided information on specific genes and biological processes with major effects on the expression of such traits, which can greatly increase animal genetic evaluation. In addition, genetic parameter estimates and genomic predictions involving sexual precocity indicator traits and productive, reproductive, and feed-efficiency related traits highlighted the feasibility and importance of direct selection for anticipating heifer reproductive life. Indeed, the case study of selection for sexual precocity in Nellore breeding programs presented here show that, in 12 years of selection for female early precocity and improved management practices, the phenotypic means of age at first calving showed a strong decreasing trend, changing from nearly 34 to less than 28 months, with a genetic trend of almost −2 days/year. In this period, the percentage of early pregnancy in the herds changed from around 10% to more than 60%, showing that the genetic improvement of heifer’s sexual precocity allows optimizing the productive cycle by reducing the number of unproductive animals in the herd. It has a direct impact on sustainability by better use of resources. Genomic selection breeding programs accounting for genotype by environment interaction represent promising tools for accelerating genetic progress for sexual precocity in tropical beef cattle.

Genome-wide association study using haplotype libraries and repeated measures model to identify candidate genomic regions for stillbirth in Holstein cattle

Reduced fertility is one of the main causes of economic losses on dairy farms, resulting in economic losses estimated at $938 per stillbirth case in Holstein herds. The identification of genomic regions associated with stillbirth could help to develop better management and breeding strategies aimed to reduce the frequency of undesirable gestation outcomes. Here, 10,570 cows and 50,541 birth records were used to perform a haplotype-based GWAS. A total of 41 significantly associated pseudo-SNPs (haplotypes within haplotype blocks converted to a binary classification) were identified after Bonferroni adjustment for multiple tests. A total of 117 positional candidate genes were annotated within or close (in a 200-kb interval) to significant pseudo-SNPs (haplotype blocks). The guilt-by-association functional prioritization identified 31 potential functional candidate genes for reproductive performance out of the 117 positional candidate genes annotated. These genes play crucial roles in biological processes associated with pregnancy persistence, fetus development, immune response, among others. These results helped us to better understand the genetic basis of stillbirth in dairy cattle and may be useful for the prediction of stillbirth in Holstein cattle, helping to reduce the related economic losses caused by this phenotype.

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.

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.

Genetic Features of Reproductive Traits in Bovine and Buffalo: Lessons From Bovine to Buffalo

Bovine and buffalo are important livestock species that have contributed to human lives for more than 1000 years. Improving fertility is very important to reduce the cost of production. In the current review, we classified reproductive traits into three categories: ovulation, breeding, and calving related traits. We systematically summarized the heritability estimates, molecular markers, and genomic selection (GS) for reproductive traits of bovine and buffalo. This review aimed to compile the heritability and genome-wide association studies (GWASs) related to reproductive traits in both bovine and buffalos and tried to highlight the possible disciplines which should benefit buffalo breeding. The estimates of heritability of reproductive traits ranged were from 0 to 0.57 and there were wide differences between the populations. For some specific traits, such as age of puberty (AOP) and calving difficulty (CD), the majority beef population presents relatively higher heritability than dairy cattle. Compared to bovine, genetic studies for buffalo reproductive traits are limited for age at first calving and calving interval traits. Several quantitative trait loci (QTLs), candidate genes, and SNPs associated with bovine reproductive traits were screened and identified by candidate gene methods and/or GWASs. The IGF1 and LEP pathways in addition to non-coding RNAs are highlighted due to their crucial relevance with reproductive traits. The distribution of QTLs related to various traits showed a great differences. Few GWAS have been performed so far on buffalo age at first calving, calving interval, and days open traits. In addition, we summarized the GS studies on bovine and buffalo reproductive traits and compared the accuracy between different reports. Taken together, GWAS and candidate gene approaches can help to understand the molecular genetic mechanisms of complex traits. Recently, GS has been used extensively and can be performed on multiple traits to improve the accuracy of prediction even for traits with low heritability, and can be combined with multi-omics for further analysis.

Genome-wide association study for age at puberty and resumption of cyclicity in a crossbred dairy cattle population

Fertility is of primary economic importance in dairy cattle and the most common reason for involuntary culling. However, standard fertility traits have very low heritability that renders genetic selection slow and difficult. In this study, we explored fertility from an endocrine standpoint. A total of 1,163 crossbred Holstein-Normande females in a 3-generation familial design were studied for progesterone level measured every 10 d to determine age at puberty (PUB) and commencement of postpartum luteal activity (CPLA). Genetic parameters were estimated using REML with WOMBAT software. The heritability estimates were 0.38 ± 0.10 and 0.16 ± 0.07 for PUB and CPLA, respectively. Moreover, the 2 traits were genetically correlated (0.45 ± 0.23), suggesting a partially common determinism. Because of the family structure, a linkage disequilibrium and linkage analysis approach was preferred over standard genome-wide association study to map genomic regions associated with these traits. Ten quantitative trait loci (QTL) were detected for PUB on chromosomes 1, 3, 11, 13, 14, 21, and 29, whereas 3 QTL were associated with CPLA on chromosomes 21 and 26. Only the QTL on chromosome 21 was common to both traits. Four functional candidate genes (NCOA2, GAS2, OVOL1, and FOSL1) were identified in the detected regions. These findings will contribute to a clearer understanding of fertility determinism and enhance the value of introducing endocrinological data in fertility studies.

Effect of indel variants within the sorting nexin 29 (SNX29) gene on growth traits of goats

The sorting nexin 29 gene (SNX29) is a well-known regulator of myocyte differentiation and proliferation. In this work, two indels (17-bp and 21-bp) were identified in the goat SNX29 gene, and their effects on the growth traits of 1,759 Shaanbei white cashmere (SBWC) goats were analyzed. Both indels had three genotypes [homozygote wild type (II), heterozygote (ID), and homozygote mutation (DD)] and displayed medium genetic diversity (0.25 < polymorphism information content (PIC) < 0.50) in the population. The 17-bp indel was significantly associated with chest width (p = 0.009), body weight (p = 0.021), and chest depth (p = 0.032), with the II genotype dominant. The 21-bp indel was significantly associated with chest width (p = 0.001), chest depth (p = 4.8E-5), heart girth (p = 0.007), and hip width (p = 0.002). Because the two indels were in the upstream (17-bp) and intron (21-bp) regions of the SNX29 gene, transcription factor binding sites were predicted. The IRF5 and MYC could bind with the 17-bp indel and 21-bp indel sequences, respectively. This study indicates that SNX29 is a promising candidate gene that can be used to improve meat production in goat breeding.

Identification of genomic regions related to age at first calving and first calving interval in water buffalo using single-step GBLUP

In Brazil, water buffaloes have been used to produce milk for mozzarella cheese production. Consequently, the main selection criterion applied for the buffalo genetic improvement is the estimated mozzarella yield as a function of milk, fat and protein production. However, given the importance of reproductive traits in production systems, this study aimed to use techniques for identifying genomic regions that affect the age at first calving (AFC) and first calving interval (FCI) in buffalo cows and to select candidate genes for the identification of QTL and gene expression studies. The single-step GBLUP method was used for the identification of genomic regions. Windows of 1 Mb containing single-nucleotide polymorphisms were constructed and the 10 windows that explained the greatest proportion of genetic variance were considered candidate regions for each trait. Genes present into the selected windows were identified using the UOA_WB_1 assembly as the reference, and their ontology was defined with the Panther tool. Candidate regions for both traits were identified on BBU 3, 12, 21 and 22; for AFC, candidates were detected on BBU 6, 7, 8, 9 and 15 and for first calving interval on BBU 4, 14 and 19. This study identified regions with great contribution to the additive genetic variance of age at first calving and first calving interval in the population of buffalo cows studied. The ROCK2, PMVK, ADCY2, MAP2K6, BMP10 and GFPT1 genes are main candidates for reproductive traits in water dairy buffaloes, and these results may have future applications in animal breeding programs or in gene expression studies of the species.

Determination of genetic variation within the DYRK2 gene and its associations with milk traits in cattle

To speed up the progress of marker-assisted selection (MAS) in cattle breeding, the dual-specificity tyrosine phosphorylation-regulated kinase 2 (DYRK2), cadherin 2 (CDH2), and kinesin family member 1A (KIF1A) genes were chosen based on our pervious genome-wide association study (GWAS) analysis results. DYRK2 is a kinase that may participate in cell growth and/or development; it shows phosphorylation activity toward serine, threonine, and tyrosine fragments of proteins, and it is different from other protein kinases. The CDH2 gene encodes a classic cadherin, which is a member of the cadherin superfamily. The protein encoded by KIF1A is a member of the kinesin family and plays a role in the transportation of membrane organelles along axon microtubules. We detected insertion/deletion (InDel) variation in these three candidate genes in 438 individual cattle (Xinjiang Brown cattle and Wagyu × Luxi crossbreed cattle). Only DYRK2-P3-11 bp was polymorphic and genotyped. The polymorphism information content of DYRK2-P3-11 bp was 0.336. Correlation analyses showed that InDel polymorphism was significantly associated with six different milk traits. These findings may aid future analyses of InDel genotypes in cattle breeds, and speed up the progress of MAS in cattle breeding.

Genome-Wide Assessment of Runs of Homozygosity in Chinese Wagyu Beef Cattle

Runs of homozygosity (ROH) are continuous homozygous regions that generally exist in the DNA sequence of diploid organisms. Identifications of ROH leading to reduction in performance can provide valuable insight into the genetic architecture of complex traits. Here, we evaluated genome-wide patterns of homozygosity and their association with important traits in Chinese Wagyu beef cattle. We identified a total of 29,271 ROH segments from 462 animals. Within each animal, an average number of ROH was 63.36 while an average length was 62.19 Mb. To evaluate the enrichment of ROH across genomes, we initially identified 280 ROH regions by merging ROH events across all individuals. Of these, nine regions containing 154 candidate genes, were significantly associated with six traits (body height, chest circumference, fat coverage, backfat thickness, ribeye area, and carcass length; p < 0.01). Moreover, we found 26 consensus ROH regions with frequencies exceeding 10%, and several regions overlapped with QTLs, which are associated with body weight, calving ease, and stillbirth. Among them, we observed 41 candidate genes, including BCKDHB, MAB21L1, SLC2A13, FGFR3, FGFRL1, CPLX1, CTNNA1, CORT, CTNNBIP1, and NMNAT1, which have been previously reported to be related to body conformation, meat quality, susceptibility, and reproductive traits. In summary, we assessed genome-wide autozygosity patterns and inbreeding levels in Chinese Wagyu beef cattle. Our study identified many candidate regions and genes overlapped with ROH for several important traits, which could be unitized to assist the design of a selection mating strategy in beef cattle.

Effects of eight InDel variants in FHIT on milk traits in Xinjiang brown cattle

In our previous genome-wide association study (GWAS), we identified the fragile histidine triad diadenosine triphosphatase (FHIT) gene in Xinjiang brown cattle (XJBC) as a candidate gene associated with cattle productive traits, with potential application in mark-assisted selection (MAS) in cattle breeding. FHIT is a prototype of a class of tumor suppressor genes that contain genomic loci mapped to common fragile loci. Here, 388 healthy and unrelated XJBC were selected to identify insertion/deletion (InDel) variants in the bovine FHIT and assess their effects on milk traits. Eight of the thirteen InDel loci were found to be polymorphic in FHIT. The polymorphism information content of the eight loci ranged from 0.061 to 0.375. The correlation analysis showed that all the new InDel variants were significantly related to six different milk traits (p < 0.05). The following variants presented a significant relationship with productive traits: P2-23bp with the 305 milk yield (p = 0.005) in the sixth parity; P3-24bp with the milk fat yield (p = 0.009) in the third parity; P5-21bp with the somatic cell score (p = 0.001) in the first parity and with the milk protein percentage (p = 0.002) in the sixth parity; and P7-26bp with the somatic cell score (p = 0.003) in the sixth parity. These findings will help evaluate InDel genotypes, within and between cattle breeds and identify potential target loci to accelerate progress in MAS in cattle breeding.

Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers

Age at first calving (AFC) plays an important role in the economic efficiency of beef cattle production. This trait can be affected by a combination of genetic and environmental factors, leading to physiological changes in response to heifers’ adaptation to a wide range of environments. Genome-wide association studies through the reaction norm model were carried out to identify genomic regions associated with AFC in Nellore heifers, raised under different environmental conditions (EC). The SNP effects for AFC were estimated in three EC levels (Low, Medium, and High, corresponding to average contemporary group effects on yearling body weight equal to 159.40, 228.6 and 297.6 kg, respectively), which unraveled shared and unique genomic regions for AFC in Low, Medium, and High EC levels, that varied according to the genetic correlation between AFC in different EC levels. The significant genomic regions harbored key genes that might play an important biological role in controlling hormone signaling and metabolism. Shared genomic regions among EC levels were identified on BTA 2 and 14, harboring candidate genes associated with energy metabolism (IGFBP2, IGFBP5, SHOX, SMARCAL1, LYN, RPS20, MOS, PLAG1, CHCD7, and SDR16C6). Gene set enrichment analyses identified important biological functions related to growth, hormone levels affecting female fertility, physiological processes involved in female pregnancy, gamete generation, ovulation cycle, and age at puberty. The genomic regions highlighted differences in the physiological processes linked to AFC in different EC levels and metabolic processes that support complex interactions between the gonadotropic axes and sexual precocity in Nellore heifers.

Transcriptome profile and association study revealed STAT3 gene as a potential quality marker of bovine gametes

The present study was aimed to investigate differences in molecular signatures in oocytes derived from Holstein-Friesian heifers with different genetic merit for fertility, euthanized during day 0 or day 12 of the estrous cycle. Moreover, association between single nucleotide polymorphisms (SNPs) of ODC1 and STAT3 genes and bull fertility traits was investigated. The gene expression patterns were analyzed using cDNA array and validated with quantitative real-time polymerase chain reaction (PCR). The result revealed that several genes have shown not only to be regulated by fertility merit but also by the day of oocyte recovery during the estrous cycle. The STAT3 gene was found to be upregulated in oocytes recovered from animals with high fertility merit at both day 0 and day 12. Some other genes like PTTG1, ODC1 and TUBA1C were downregulated at day 0 and upregulated at day 12 in high, compared with low, fertility merit recovered oocytes. In contrast, the transcript abundance of TPM3 was upregulated at day 0 and downregulated at day 12 in high, compared with low, fertility merit recovered oocytes. In addition, ODC1 and STAT3 were found to be associated (P < 0.05) with sperm quality traits as well as flow cytometry parameters. Therefore, the expression of several candidate genes including ODC1 and STAT3 was related to the genetic merit of the cow. In addition polymorphisms in these two genes were found to be associated with bull semen quality.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin

Background

Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance.

Results

Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively.

Conclusion

This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Genetic parameters for fertility traits in Nellore bulls

Nellore is the main cattle breed used in Brazil, being the largest commercial herd in the world. Beyond the importance of male reproductive efficiency for farm profit, the use of reproductive techniques, mainly artificial insemination, turns the evaluation of male reproductive traits even more important. Estimation of genetic parameters increases the knowledge on traits variances and allows envisaging the possibility of the inclusion of new traits as selection criterion. Genetic parameters for fifteen traits that can be classified as testicular biometry or physical and morphological semen traits were estimated for a Nellore bull population ranging from 18 to 36 months. Single-trait and bi-trait animal models were used for (co)variance components estimation. The contemporary group was considered as fixed effect and age at measurement as covariable. Scrotal circumference presented heritability of 0.47 ± 0.12. This value is similar to the heritabilities found for all testicular biometry traits (0.34-0.48). Sperm progressive motility, which has a direct effect on bull fertility, presented low heritability (0.07 ± 0.08). Major and total sperm defects presented moderate to high heritabilities (0.49 ± 0.18 and 0.39 ± 0.15, respectively), indicating that great genetic gain can be obtained through selection against sperm defects. High and positive genetic correlations were observed among testicular biometry traits, which also presented favourable genetic correlations with physical and morphological traits of the semen with magnitude ranging from high to low. Scrotal circumference presented moderate to high and favourable genetic correlations with sperm progressive motility, sperm turbulence, major sperm defects and total sperm defects. Thus, the selection for scrotal circumference results in favourable correlated genetic response for semen quality. The results show that the use of scrotal circumference as reference trait for bull fertility is appropriate, since it presents high heritability and favourable genetic correlation with semen quality.

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.

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.

Identifying Loci Under Positive Selection in Yellow Korean Cattle (Hanwoo)

Jeju Black cattle is one of the aboriginal Korean cattle breeds that has been isolated in Jeju island for a long time, while Yellow Hanwoo cattle has been extensively selected for beef production traits for the last several decades. Aiming to investigate broader patterns of selection, we genotyped 352 Yellow Hanwoo and 169 Jeju Black cattle using a customized 150 K bovine chip. Our composite selection signals’ analysis to identify selection signatures (cross-population extended haplotype homozygosity [XP-EHH], ΔSAF, and F_ST) identified recent and strong signature of selection near many loci with mutations affecting the traits under strong selection as outlier in Yellow Hanwoo, including SCP2 (P = 8.41 × 10⁻¹⁰) that may be involved in the meat quality. We found nine candidate regions with significant clusters of selection signals, and further bioinformatics analyses of the genes located within these regions revealed mainly genes involved in G-protein coupled receptor signaling pathway (GO:0007186) or olfactory transduction (bta04740), which may be due to adaptation to natural environments in Jeju island. Based on the stronger correlation of Ne₁₀/Ne₁₀₀ ratio between Yellow Hanwoo (0.61) and Jeju Black (0.66) cattle, our results suggest that the difference of chromosomal regions of selection signature between the 2 cattle breeds was due to a consequence of selection processes to adapt to environmental differences between Jeju island and the main inland, Korean peninsula.

Genotype-by-environment interaction of fertility traits in Danish Holstein cattle using a single-step genomic reaction norm model

Genotype-by-environment (G × E) interactions could play an important role in cattle populations, and it should be considered in breeding programmes to select the best sires for different environments. The objectives of this study were to study G × E interactions for female fertility traits in the Danish Holstein dairy cattle population using a reaction norm model (RNM), and to detect the particular genomic regions contributing to the performance of these traits and the G × E interactions. In total 4534 bulls were genotyped by an Illumina BovineSNP50 BeadChip. An RNM with a pedigree-based relationship matrix and a pedigree-genomic combined relationship matrix was used to explore the existence of G × E interactions. In the RNM, the environmental gradient (EG) was defined as herd effect. Further, the genomic regions affecting interval from calving to first insemination (ICF) and interval from first to last insemination (IFL) were detected using single-step genome-wide association study (ssGWAS). The genetic correlations between extreme EGs indicated that G × E interactions were sizable for ICF and IFL. The genomic RNM (pedigree-genomic combined relationship matrix) had higher prediction accuracy than the conventional RNM (pedigree-based relationship matrix). The top genomic regions affecting the slope of the reaction norm included immunity-related genes (IL17, IL17F and LIF), and growth-related genes (MC4R and LEP), while the top regions influencing the intercept of the reaction norm included fertility-related genes such as EREG, AREG and SMAD4. In conclusion, our findings validated the G × E interactions for fertility traits across different herds and were helpful in understanding the genetic background of G × E interactions for these traits.

Genome-wide scan reveals divergent selection among taurine and zebu cattle populations from different regions

In this study, to identify genomic signatures of divergent selection, we genotyped 10 cattle breeds/populations (n = 275), representing eight Ethiopian cattle populations (n = 229) and two zebu populations (n = 46) adapted to tropical and sub-tropical environments, using the high-density single-nucleotide polymorphisms (SNPs) derived mainly from Bos indicus breeds, and using five reference taurine breeds (n = 212). Population genetic differentiation (F_ST ) values across sliding windows were estimated between zebu and reference combined taurine breeds. The most differentiated regions (F_ST  ≥ 0.53), representing the top 1% smoothed F_ST values, were considered to represent regions under diversifying selection. In total, 285 and 317 genes were identified in the comparisons of Ethiopian cattle with taurine and Asian zebu with taurine respectively. Some of these genes are involved in stress responses/thermo-tolerance and DNA damage repair (HSPA4, HSF1, CMPK1 and EIF2AK4), pigmentation (ERBB3 and MYO1A), reproduction/fertility (UBE2D3, ID3 and PSPC1), immune response (PIK3CD and AKIRIN2) and body stature and size (MBP2, LYN and NPM1). Additionally, the candidate genes were associated with functional terms (e.g. cellular response to stress, DNA repair, inflammatory response) important for physiological adaptation to environmental stresses. The results of our study may shed light on the influence of artificial and natural selection in shaping the genomic diversity of modern cattle breeds and also may serve as a basis for further genetic investigation of traits of tropical adaptation in cattle.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits

Background

The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (F_ROH), genomic relationship matrix (F_GRM), and pedigree-based coefficient (F_PED).

Results

The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of F_PED and F_ROH, and their correlations (− 0.05 to 0.26) were low. Estimates of correlation between F_GRM-F_PED was zero, while the correlation (− 0.01 to − 0.07) between F_GRM-F_ROH decreased as a function of ROH length, except for F_{ROH > 8Mb} (− 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands.

Conclusions

Low F_PED-F_ROH correlation estimates indicate that F_PED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and F_ROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Electronic supplementary material

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

Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates.

Genome-wide association study using haplotype alleles for the evaluation of reproductive traits in Nelore cattle

Zebu cattle (Bos taurus indicus) are highly adapted to tropical regions. However, females reach puberty after taurine heifers, which affects the economic efficiency of beef cattle breeding in the tropical regions. The aims of this study were to establish associations between the haplotype alleles of the bovine genome and age at first calving (AFC) in the Nelore cattle, and to identify the genes and quantitative trait loci (QTL) related to this phenotype. A total of 2,273 Nelore cattle (995 males and 1,278 females) genotyped using the Illumina BovineHD BeadChip were used in the current study. The association analysis included females with valid first calving records as well as open heifers. Linkage disequilibrium (LD) analysis among the markers was performed using blocks of 5, 10, and 15 markers, which were determined by sliding windows shifting one marker at a time. Then, the haplotype block size to be used in the association study was chosen based on the highest r² average among the SNPs in the block. The five HapAlleles most strongly associated with the trait (top five) were considered as significant associations. The results of the analysis revealed four genomic regions related to AFC, which overlapped with 20 QTL of the reproductive traits reported previously. Furthermore, there were 19 genes related to reproduction in those regions. In conclusion, the use of haplotypes allowed the detection of chromosomal regions associated with AFC in Nelore cattle, and provided the basis for elucidating the mechanisms underlying this trait.

Insertion/Deletion Within the KDM6A Gene Is Significantly Associated With Litter Size in Goat

A previous whole-genome association analysis identified lysine demethylase 6A (KDM6A), which encodes a type of histone demethylase, as a candidate gene associated to goat fecundity. KDM6A gene knockout mouse disrupts gametophyte development, suggesting that it has a critical role in reproduction. In this study, goat KDM6A mRNA expression profiles were determined, insertion/deletion (indel) variants in the gene identified, indel variants effect on KDM6A gene expression assessed, and their association with first-born litter size analyzed in 2326 healthy female Shaanbei white cashmere goats. KDM6A mRNA was expressed in all tissues tested (heart, liver, spleen, lung, kidney, muscle, brain, skin and testis); the expression levels in testes at different developmental stages [1-week-old (wk), 2, 3 wk, 1-month-old (mo), 1.5 and 2 mo] indicated a potential association with the mitosis-to-meiosis transition, implying that KDM6A may have an essential role in goat fertility. Meanwhile, two novel intronic indels of 16 bp and 5 bp were identified. Statistical analysis revealed that only the 16 bp indel was associated with first-born litter size (P < 0.01), and the average first-born litter size of individuals with an insertion/insertion genotype higher than that of those with the deletion/deletion genotype (P < 0.05). There was also a significant difference in genotype distributions of the 16 bp indel between mothers of single-lamb and multi-lamb litters in the studied goat population (P = 0.001). Consistently, the 16 bp indel also had a significant effect on KDM6A gene expression. Additionally, there was no significant linkage disequilibrium (LD) between these two indel loci, consistent with the association analysis results. Together, these findings suggest that the 16 bp indel in KDM6A may be useful for marker-assisted selection (MAS) of goats.