Genome-Wide Association Study Based on Random Regression Model Reveals Candidate Genes Associated with Longitudinal Data in Chinese Simmental Beef Cattle

Genome-Wide Association Study Reveals the Genetic Architecture of Growth and Meat Production Traits in a Chicken F2 Resource Population

BACKGROUND/OBJECTIVES

For genomic selection to enhance the efficiency of broiler production, finding SNPs and candidate genes that define the manifestation of main selected traits is essential. We conducted a genome-wide association study (GWAS) for growth and meat productivity traits of roosters from a chicken F2 resource population (n = 152).

METHODS

The population was obtained by crossing two breeds with contrasting phenotypes for performance indicators, i.e., Russian White (slow-growing) and Cornish White (fast-growing). The birds were genotyped using the Illumina Chicken 60K SNP iSelect BeadChip. After LD filtering of the data, 54,188 SNPs were employed for the GWAS analysis that allowed us to reveal significant specific associations for phenotypic traits of interest and economic importance.

RESULTS

At the threshold value of p < 9.2 × 10-7, 83 SNPs associated with body weight at the age of 28, 42, and 63 days were identified, as well as 171 SNPs associated with meat qualities (average daily gain, slaughter yield, and dressed carcass weight and its components). Moreover, 34 SNPs were associated with a group of three or more traits, including 15 SNPs significant for a group of growth traits and 5 SNPs for a group of meat productivity indicators. Relevant to these detected SNPs, nine prioritized candidate genes associated with the studied traits were revealed, including WNT2, DEPTOR, PPA2, UNC80, DDX51, PAPPA, SSC4D, PTPRU, and TLK2.

CONCLUSIONS

The found SNPs and candidate genes can serve as genetic markers for growth and meat performance characteristics in chicken breeding in order to achieve genetic improvement in broiler production.

Population structure and identification of genomic regions associated with productive traits in five Italian beef cattle breeds

Italy has a long history in beef production, with local breeds such as Marchigiana, Chianina, Romagnola, Maremmana, and Podolica which produce high-quality meat. Selection has improved meat production, precocity, growth ability and muscle development, but the genetic determinism of such traits is mostly unknown. Using 33K SNPs-data from young bulls (N = 4064) belonging to these five Italian breeds, we demonstrated that the Maremmana and Podolica rustic breeds are closely related, while the specialised Marchigiana, Chianina, and Romagnola breeds are more differentiated. A genome-wide association study for growth and muscle development traits (average daily gain during the performance test, weight at 1 year old, muscularity) was conducted in the five Italian breeds. Results indicated a region on chromosome 2, containing the myostatin gene (MSTN), which displayed significant genome-wide associations with muscularity in Marchigiana cattle, a breed in which the muscle hypertrophy phenotype is segregating. Moreover, a significant SNP on chromosome 14 was associated, in the Chianina breed, to muscularity. The identification of diverse genomic regions associated with conformation traits might increase our knowledge about the genomic basis of such traits in Italian beef cattle and, eventually, such information could be used to implement marker-assisted selection of young bulls tested in the performance test.

Including microbiome information in a multi-trait genomic evaluation: a case study on longitudinal growth performance in beef cattle

BACKGROUND

Growth rate is an important component of feed conversion efficiency in cattle and varies across the different stages of the finishing period. The metabolic effect of the rumen microbiome is essential for cattle growth, and investigating the genomic and microbial factors that underlie this temporal variation can help maximize feed conversion efficiency at each growth stage.

RESULTS

By analysing longitudinal body weights during the finishing period and genomic and metagenomic data from 359 beef cattle, our study demonstrates that the influence of the host genome on the functional rumen microbiome contributes to the temporal variation in average daily gain (ADG) in different months (ADG1, ADG2, ADG3, ADG4). Five hundred and thirty-three additive log-ratio transformed microbial genes (alr-MG) had non-zero genomic correlations (rg) with at least one ADG-trait (ranging from |0.21| to |0.42|). Only a few alr-MG correlated with more than one ADG-trait, which suggests that a differential host-microbiome determinism underlies ADG at different stages. These alr-MG were involved in ribosomal biosynthesis, energy processes, sulphur and aminoacid metabolism and transport, or lipopolysaccharide signalling, among others. We selected two alternative subsets of 32 alr-MG that had a non-uniform or a uniform rg sign with all the ADG-traits, regardless of the rg magnitude, and used them to develop a microbiome-driven breeding strategy based on alr-MG only, or combined with ADG-traits, which was aimed at shaping the rumen microbiome towards increased ADG at all finishing stages. Combining alr-MG information with ADG records increased prediction accuracy of genomic estimated breeding values (GEBV) by 11 to 22% relative to the direct breeding strategy (using ADG-traits only), whereas using microbiome information, only, achieved lower accuracies (from 7 to 41%). Predicted selection responses varied consistently with accuracies. Restricting alr-MG based on their rg sign (uniform subset) did not yield a gain in the predicted response compared to the non-uniform subset, which is explained by the absence of alr-MG showing non-zero rg at least with more than one of the ADG-traits.

CONCLUSIONS

Our work sheds light on the role of the microbial metabolism in the growth trajectory of beef cattle at the genomic level and provides insights into the potential benefits of using microbiome information in future genomic breeding programs to accurately estimate GEBV and increase ADG at each finishing stage in beef cattle.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Genome-wide identification and evolutionary analysis of the FGF gene family in buffalo

Fibroblast growth factors (FGFs) are important polypeptide growth factors that play a critical role in many developmental processes, including differentiation, cell proliferation, and migration in mammals. This study employs in silico analyses to characterize the FGF gene family in buffalo, investigating their genome-wide identification, physicochemical properties, and evolutionary patterns. For this purpose, genomic and proteomic sequences of buffalo, cattle, goat, and sheep were retrieved from NCBI database. We identified a total of 22 FGF genes in buffalo. Physicochemical properties observed through ProtParam tool showed notable features of these proteins including in-vitro instability, thermostability, hydrophilicity, and basic nature. Phylogenetic analysis grouped 22 identified genes into nine sub-families based on evolutionary relationships. Additionally, analysis of gene structure, motif patterns, and conserved domains using TBtools revealed the remarkable conservation of this gene family across selected species throughout the course of evolution. Comparative amino acid analysis performed through ClustalW demonstrated significant conservation between buffalo and cattle FGF proteins. Mutational analysis showed three non-synonymous mutations at positions R103 > G, P7 > L, and E98 > Q in FGF4, FGF6, and FGF19, respectively in buffalo. Duplication events revealed only one segmental duplication (FGF10/FGF22) in buffalo and two in cattle (FGF10/FGF22 and FGF13/FGF13-like) with Ka/Ks values <1 indicating purifying selection pressure for these duplications. Comparison of protein structures of buffalo, goat, and sheep exhibited more similarities in respective structures. In conclusion, our study highlights the conservation of the FGF gene family in buffalo during evolution. Furthermore, the identified non-synonymous mutations may have implications for the selection of animals with better performance.Communicated by Ramaswamy H. Sarma.

Longitudinal genome-wide association analysis using a single-step random regression model for height in Japanese Holstein cattle

Growth traits, such as body weight and height, are essential in the design of genetic improvement programs of dairy cattle due to their relationship with feeding efficiency, longevity, and health. We investigated genomic regions influencing height across growth stages in Japanese Holstein cattle using a single-step random regression model. We used 72,921 records from birth to 60 mo of age with 4,111 animals born between 2000 and 2016. The analysis included 1,410 genotyped animals with 35,319 single nucleotide polymorphisms, consisting of 883 females with records and 527 bulls, and 30,745 animals with pedigree information. A single genomic region at the 58.4 megabase pair on chromosome 18 was consistently identified across 6 age points of 10, 20, 30, 40, 50, and 60 mo after multiple testing corrections for the significance threshold. Twelve candidate genes, previously reported for longevity and gestation length, were found near the identified genomic region. Another location near the identified region was also previously associated with body conformation, fertility, and calving difficulty. Functional Gene Ontology enrichment analysis suggested that the candidate genes regulate dephosphorylation and phosphatase activity. Our findings show that further study of the identified candidate genes will contribute to a better understanding of the genetic basis of height in Japanese Holstein cattle.

Estimation of the Genetic Components of (Co)variance and Preliminary Genome-Wide Association Study for Reproductive Efficiency in Retinta Beef Cattle

In this study, we analyzed the variation of reproductive efficiency, estimated as the deviation between the optimal and real parity number of females at each stage of the cow's life, in 12,554 cows belonging to the Retinta Spanish cattle breed, using classical repeatability and random regression models. The results of the analyses using repeatability model and the random regression model suggest that reproductive efficiency is not homogeneous throughout the cow's life. The h² estimate for this model was 0.30, while for the random regression model it increased across the parities, from 0.24 at the first calving to 0.51 at calving number 9. Additionally, we performed a preliminary genome-wide association study for this trait in a population of 252 Retinta cows genotyped using the Axiom Bovine Genotyping v3 Array. The results showed 5 SNPs significantly associated with reproductive efficiency, located in two genomic regions (BTA4 and BTA28). The functional analysis revealed the presence of 5 candidate genes located within these regions, which were previously involved in different aspects related to fertility in cattle and mice models. This new information could give us a better understanding of the genetic architecture of reproductive traits in this species, as well as allow us to accurately select more fertile cows.

Genome-Wide Association Study of Body Weight Trait in Yaks

The yak is the largest meat-producing mammal around the Tibetan Plateau, and it plays an important role in the economic development and maintenance of the ecological environment throughout much of the Asian highlands. Understanding the genetic components of body weight is key for future improvement in yak breeding; therefore, genome-wide association studies (GWAS) were performed, and the results were used to mine plant and animal genetic resources. We conducted whole genome sequencing on 406 Maiwa yaks at 10 × coverage. Using a multiple loci mixed linear model (MLMM), fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK), we found that a total of 25,000 single-nucleotide polymorphisms (SNPs) were distributed across chromosomes, and seven markers were identified as significantly (p-values < 3.91 × 10−7) associated with the body weight trait,. Several candidate genes, including MFSD4, LRRC37B, and NCAM2, were identified. This research will help us achieve a better understanding of the genotype−phenotype relationship for body weight.

Identification of Candidate Genes Regulating Carcass Depth and Hind Leg Circumference in Simmental Beef Cattle Using Illumina Bovine Beadchip and Next-Generation Sequencing Analyses

Genome-wide association studies are a robust means of identifying candidate genes that regulate economically important traits in farm animals. The aim of this study is to identify single-nucleotide polymorphisms (SNPs) and candidate genes potentially related to carcass depth and hind leg circumference in Simmental beef cattle. We performed Illumina Bovine HD Beadchip (~670 k SNPs) and next-generation sequencing (~12 million imputed SNPs) analyses of data from 1252 beef cattle, to which we applied a linear mixed model. Using a statistical threshold (p = 0.05/number of SNPs identified) and adopting a false discovery rate (FDR), we identified many putative SNPs on different bovine chromosomes. We identified 12 candidate genes potentially annotated with the markers identified, including CDKAL1 and E2F3, related to myogenesis and skeletal muscle development. The identification of such genes in Simmental beef cattle will help breeders to understand and improve related traits, such as meat yield.