A genome scan for quantitative trait loci affecting average daily gain and Kleiber ratio in Baluchi Sheep

A Genomic and Phenotypic Investigation of Feed Efficiency and Growth Traits in Targhee and Rambouillet Sheep

The U.S. range sheep industry uses estimated breeding values (EBVs) as part of their breeding objectives to increase post-weaning weight. The study objective was to quantify the relationship between lamb growth EBVs, feed intake, and feed efficiency. Eighty-one range ewe lambs were enrolled in the study to measure residual feed intake (RFI) over two 42-d periods at both the weaning and yearling stages. The ewe lambs' post-weaning weight EBVs (PWWT EBVs) were linearly associated with their phenotypic traits. Preliminary genome wide associations (GWAs) were also performed with Dry Matter Intake (DMI), RFI, mid-test body size, and average daily gain (ADG) and Ovine 50K SNP genotypes. Post-weaning weight EBVs were associated with dry matter intake (DMI) (p < 0.05) but had no association with residual feed intake (RFI) (p > 0.05) in both experimental periods. However, PWWT EBV was predictive of mid-test body weight in both periods (p < 0.05). A single SNP at Oar2:68,812,505, located within DMRT2, was associated with DMI and RFI in the second experimental period (Bonferroni corrected p <0.05). While selecting for higher post-weaning weight range ewes may increase feed consumed due to a larger body size, it was not associated with feed efficiency.

Genome-Wide Association Analysis of Growth Traits in Hu Sheep

(1) Background: The Hu sheep is a renowned breed characterized by high reproduction, year-round estrus, and resistance to high humidity and temperature conditions. However, the breed exhibits lower growth rates and meat yields, which necessitate improvements through selective breeding. The integration of molecular markers in sheep breeding programs has the potential to enhance growth performance, reduce breeding cycles, and increase meat production. Currently, the applications of SNP chips for genotyping in conjunction with genome-wide association studies (GWAS) have become a prevalent approach for identifying candidate genes associated with economically significant traits in livestock. (2) Methods: To pinpoint candidate genes influencing growth traits in Hu sheep, we recorded the birth weight, weaning weight, and weights at 3, 4, 5, 6, and 7 months for a total of 567 Hu sheep, and genotyping was performed using the Ovine 40K SNP chip. (3) Results: Through GWAS analysis and KEGG pathway enrichment, we identified three candidate genes associated with birth weight (CAMK2B, CACNA2D1, and CACNA1C). Additionally, we found two candidate genes linked to weaning weight (FGF9 and BMPR1B), with CACNA2D1 also serving as a shared gene between birth weight and weaning weight traits. Furthermore, we identified eight candidate genes related to monthly weight (FIGF, WT1, KCNIP4, JAK2, WWP1, PLCL1, GPRIN3, and CCSER1). (4) Conclusion: Our findings revealed a total of 13 candidate genetic markers that can be utilized for molecular marker-assisted selection, aiming to improve meat production in sheep breeding programs.

Investigation of growth traits in Turkish Merino lambs using multi-locus GWAS approaches: Middle Anatolian Merino

This study explored the genetic basis of growth traits in Middle Anatolian Merino lambs using multi-locus genome-wide association study (GWAS) analyses. Descriptive statistics indicated moderate heritability (h² = 0.363) for birth weight (BW) and (h² = 0.309) for weaning weight (WW), both statistically significant (p < 0.001). Strong genetic correlations were observed between WW and BW (rG = 0.922) and WW and Kleiber ratio (KR, rG = 0.896), implying that simultaneous improvements may be possible through targeted selection. Five multi-locus methods (mrMLM, FASTmrMLM, pLARmEB, FASTmrEMMA, and ISIS EM-BLASSO) were used to identify the polygenic basis of the traits. For BW and WW, 20 and 18 significant SNPs (LOD ≥ 5) were detected, respectively, with some SNPs co-detected by multiple methods. In contrast, only 10 significant SNPs were identified for KR, all exclusively by the ISIS EM-BLASSO approach. Pathway analyses within ± 100 Kb of associated SNPs revealed genes and pathways influencing these traits, which could be leveraged in future breeding programs for enhanced growth performance. The identified SNPs, particularly those associated with BW and WW traits, could facilitate genomic selection if validated in larger and more diverse populations.

A comprehensive review of livestock development: insights into domestication, phylogenetics, diversity, and genomic advances

Livestock plays an essential role in sustaining human livelihoods, offering a diverse range of species integral to food security, economic stability, and cultural traditions. The domestication of livestock, which began over 10,000 years ago, has driven significant genetic changes in species such as cattle, buffaloes, sheep, goats, and pigs. Recent advancements in genomic technologies, including next-generation sequencing (NGS), genome-wide association studies (GWAS), and genomic selection, have dramatically enhanced our understanding of these genetic developments. This review brings together key research on the domestication process, phylogenetics, genetic diversity, and selection signatures within major livestock species. It emphasizes the importance of admixture studies and evolutionary forces like natural selection, genetic drift, and gene flow in shaping livestock populations. Additionally, the integration of machine learning with genomic data offers new perspectives on the functional roles of genes in adaptation and evolution. By exploring these genomic advancements, this review provides insights into genetic variation and evolutionary processes that could inform future approaches to improving livestock management and adaptation to environmental challenges, including climate change.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Disease-driven top predator decline affects mesopredator population genomic structure

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

A SNP in the ovine cathepsin K (CTSK) gene is associated with yearling growth performance in a crossbred sheep population

Cathepsin K (CTSK) is a lysosomal protease existent in the skeletal muscles which is involved in biochemical processes related to obesity. Several studies have reported the effects of CTSK gene on body weight and fat deposition in human, mice and pigs. However, information about its structure and functions in sheep is very limited. Thus, this study was performed to evaluate the association between CTSK gene variants and yearling growth performance in Afshari × Booroola-Merino crossbred sheep. A fragment of 500 bp in exon 6 and partial of intron 5 of CTSK gene was amplified with polymerase chain reaction (PCR). All animals were genotyped by single-stranded conformation polymorphism (SSCP) and further confirmed by sequencing. Association analysis using a fixed linear model indicated that g.106510225G > A SNP was significantly related to average daily weight gain (ADWG) per year, fat-tail weight to carcass weight ratio (FW/CW), muscle thickness (MT) and muscle cross-sectional area (MCSA) of animals (p < 0.05). Due to the low polymorphic information content (PIC <0.25) for targeted locus in studied population, more association studies are needed to confirm the CTSK gene effects on growth traits in sheep.

New genomic insights into the conformation of Lipizzan horses

Conformation traits are important selection criteria in equine breeding, as they describe the exterior aspects of the horse (height, joint angles, shape). However, the genetic architecture of conformation is not well understood, as data of these traits mainly consist of subjective evaluation scores. Here, we performed genome-wide association studies on two-dimensional shape data of Lipizzan horses. Based on this data, we identified significant quantitative trait loci (QTL) associated with cresty neck on equine chromosome (ECA)16 within the MAGI1 gene, and with type, hereby differentiating heavy from light horses on ECA5 within the POU2F1 gene. Both genes were previously described to affect growth, muscling and fatty deposits in sheep, cattle and pigs. Furthermore, we pin-pointed another suggestive QTL on ECA21, near the PTGER4 gene, associated with human ankylosing spondylitis, for shape differences in the back and pelvis (roach back vs sway back). Further differences in the shape of the back and abdomen were suggestively associated with the RYR1 gene, involved in core muscle weakness in humans. Therefore, we demonstrated that horse shape space data enhance the genomic investigations of horse conformation.

Variation in the ovine FOXP3 gene and its effect on growth traits in Egyptian Barki sheep

The aim of the present study was to detect the FOXP3 gene polymorphisms in Barki sheep at a variable region covering exon 13, intron 13 and the coding sequence in exon 14 and to test the association of these polymorphisms with growth traits. 122 Barki lambs were phenotyped for various growth traits, viz., birth weight (BW), weaning weight (WW), pre-weaning daily gain in weight (ADG1), post-weaning daily gain in weight (ADG2) and marketing bodyweight (MW). The polymerase chain reaction - single-strand conformational polymorphisms (PCR-SSCP) and DNA sequencing methods were used to identify the genetic variants in the FOXP3 gene. The associations between the variation in FOXP3 gene and growth traits were tested using a general linear model. Two variants (F1 and F2 with gene frequencies of 0.64 and 0.36, respectively), and three genotypes (F1F1, F1F2 and F2F2 with frequencies of 0.37, 0.53 and 0.10, respectively) were detected. The association of FOXP3 genotype was significant (p < 0.05) with ADG2 and MW. It is concluded that FOXP3 genotype might be helpful for sheep breeders to produce fast-growing lambs. However, further studies are needed in a large population to confirm the association we found.

Genome-Wide Association Studies for Growth Curves in Meat Rabbits Through the Single-Step Nonlinear Mixed Model

Growth is a complex trait with moderate to high heritability in livestock and must be described by the longitudinal data measured over multiple time points. Therefore, the used phenotype in genome-wide association studies (GWAS) of growth traits could be either the measures at the preselected time point or the fitted parameters of whole growth trajectory. A promising alternative approach was recently proposed that combined the fitting of growth curves and estimation of single-nucleotide polymorphism (SNP) effects into single-step nonlinear mixed model (NMM). In this study, we collected the body weights at 35, 42, 49, 56, 63, 70, and 84 days of age for 401 animals in a crossbred population of meat rabbits and compared five fitting models of growth curves (Logistic, Gompertz, Brody, Von Bertalanffy, and Richards). The logistic model was preferably selected and subjected to GWAS using the approach of single-step NMM, which was based on 87,704 genome-wide SNPs. A total of 45 significant SNPs distributed on five chromosomes were found to simultaneously affect the two growth parameters of mature weight (A) and maturity rate (K). However, no SNP was found to be independently associated with either A or K. Seven positional genes, including KCNIP4, GBA3, PPARGC1A, LDB2, SHISA3, GNA13, and FGF10, were suggested to be candidates affecting growth performances in meat rabbits. To the best of our knowledge, this is the first report of GWAS based on single-step NMM for longitudinal traits in rabbits, which also revealed the genetic architecture of growth traits that are helpful in implementing genome selection.

Genome-wide association studies of preweaning growth and in vivo carcass composition traits in Esme sheep

Sheep are considered as a major contributor of global food security. Moreover, sheep preweaning growth traits as well as in vivo carcass composition traits such as ultrasonic measurements of Longissimus dorsi muscle depth (UMD) and back-fat thickness (UFD) are crucially important indicators of meat yield and hot carcass composition. Despite their relative importance for productivity and profitability of a sheep production system, detected QTL for these traits are quite scarce. Therefore, we implemented GWAS for these traits using animal mixed model-based association approach provided by GenABEL in Esme sheep. Three genome-wide and 14 individual chromosome-wide associated SNPs were discovered. As a result, ESRP1, LOC105613082, ZNF641, DUSP5, TEAD1, SMOX, PTPRT, RALYL, POM121C, PHIP, LOC101106051, ZIM3, PEG3, TRPC7, FBXL4, LOC105610397, LOC105616489 and DNAAF2 were suggested as candidates. Some of the discovered genes and involved pathways were already annotated to contribute growth and development in various species including human, mice and cattle. All in all, the results of this study are expected to strongly contribute to shed a light on the underlying molecular mechanisms behind growth and carcass composition traits, with potential implications on studies aiming faster genetic improvement, targeted low-resolution SNP panel designs and genome-editing studies.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits

Background

Variants that regulate transcription, such as expression quantitative trait loci (eQTL), have shown enrichment in genome-wide association studies (GWAS) for mammalian complex traits. However, no study has reported eQTL in sheep, although it is an important agricultural species for which many GWAS of complex meat traits have been conducted. Using RNA sequence data produced from liver and muscle from 149 sheep and imputed whole-genome single nucleotide polymorphisms (SNPs), our aim was to dissect the genetic architecture of the transcriptome by associating sheep genotypes with three major molecular phenotypes including gene expression (geQTL), exon expression (eeQTL) and RNA splicing (sQTL). We also examined these three types of eQTL for their enrichment in GWAS of multi-meat traits and fatty acid profiles.

Results

Whereas a relatively small number of molecular phenotypes were significantly heritable (h² > 0, P < 0.05), their mean heritability ranged from 0.67 to 0.73 for liver and from 0.71 to 0.77 for muscle. Association analysis between molecular phenotypes and SNPs within ± 1 Mb identified many significant cis-eQTL (false discovery rate, FDR < 0.01). The median distance between the eQTL and transcription start sites (TSS) ranged from 68 to 153 kb across the three eQTL types. The number of common variants between geQTL, eeQTL and sQTL within each tissue, and the number of common variants between liver and muscle within each eQTL type were all significantly (P < 0.05) larger than expected by chance. The identified eQTL were significantly (P < 0.05) enriched in GWAS hits associated with 56 carcass traits and fatty acid profiles. For example, several geQTL in muscle mapped to the FAM184B gene, hundreds of sQTL in liver and muscle mapped to the CAST gene, and hundreds of sQTL in liver mapped to the C6 gene. These three genes are associated with body composition or fatty acid profiles.

Conclusions

We detected a large number of significant eQTL and found that the overlap of variants between eQTL types and tissues was prevalent. Many eQTL were also QTL for meat traits. Our study fills a gap in the knowledge on the regulatory variants and their role in complex traits for the sheep model.

Multi-trait GWAS using imputed high-density genotypes from whole-genome sequencing identifies genes associated with body traits in Nile tilapia

Background

Body traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits.

Results

A multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively.

Conclusions

The high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07341-z.