Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep

Whole genome sequencing reveals that five genes are related to BW trait in sheep

BW is an important economic trait in sheep that influences growth and development. Currently, most studies have used a single approach to screen genes associated with BW traits in sheep. To address this limitation, we conducted a genome-wide association study (GWAS) covering four different BW periods: birth, weaning, 6 months, and 12 months. Five new candidate genes: MAP3K1, ANKRD55, ABCB1, MEF2C and TRNAW-CCA-87 were screened using a combination of GWAS and quantitative trait loci analysis in sheep. Additionally, five genes were subjected to Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. These genes were primarily enriched in pathways related to growth hormone and energy metabolism. The results demonstrated that the above genes potentially influenced the growth and development of sheep. The five new candidate genes are closely related to the BW trait in sheep, which will be valuable for understanding the genetic mechanisms underlying BW traits and for guiding sheep breeding.

Characterization of the Hepatic Transcriptome for Divergent Immune-Responding Sheep Following Natural Exposure to Gastrointestinal Nematodes

Infections with gastrointestinal nematodes (GINs) reduce the economic efficiency of sheep operations and compromise animal welfare. Understanding the host's response to GIN infection can help producers identify animals that are naturally resistant to infection. The objective of this study was to characterize the hepatic transcriptome of sheep that had been naturally exposed to GIN parasites. The hepatic transcriptome was studied using RNA-Sequencing technology in animals characterized as high (n = 5) or medium (n = 6) based on their innate immune acute-phase (AP) response phenotype compared with uninfected controls (n = 4), and with biased antibody-mediated (AbMR, n = 5) or cell-mediated (CMR, n = 5) adaptive immune responsiveness compared to uninfected controls (n = 3). Following the assessment of sheep selected for innate responses, 0, 136, and 167 genes were differentially expressed (DE) between high- and medium-responding animals, high-responding and uninfected control animals, and medium-responding and uninfected control animals, respectively (false discovery rate (FDR) < 0.05, and fold change |FC| > 2). When adaptive immune responses were assessed, 0, 53, and 57 genes were DE between antibody- and cell-biased animals, antibody-biased and uninfected control animals, and cell-biased and uninfected control animals, respectively (FDR < 0.05, |FC| > 2). Functional analyses identified enriched gene ontology (GO) terms and metabolic pathways related to the innate immune response and energy metabolism. Six functional candidate genes were identified for further functional and validation studies to better understand the underlying biological mechanisms of host responses to GINs. These, in turn, can potentially help improve decision making and management practices to increase the overall host immune response to GIN infection.

Genomic Regions Associated with Resistance to Gastrointestinal Nematode Parasites in Sheep-A Review

Gastrointestinal nematodes (GINs) can be a major constraint and global challenge to the sheep industry. These nematodes infect the small intestine and abomasum of grazing sheep, causing symptoms such as weight loss, diarrhea, hypoproteinemia, and anemia, which can lead to death. The use of anthelmintics to treat infected animals has led to GIN resistance, and excessive use of these drugs has resulted in residue traced in food and the environment. Resistance to GINs can be measured using multiple traits, including fecal egg count (FEC), Faffa Malan Chart scores, hematocrit, packed cell volume, eosinophilia, immunoglobulin (Ig), and dagginess scores. Genetic variation among animals exists, and understanding these differences can help identify genomic regions associated with resistance to GINs in sheep. Genes playing important roles in the immune system were identified in several studies in this review, such as the CFI and MUC15 genes. Results from several studies showed overlapping quantitative trait loci (QTLs) associated with multiple traits measuring resistance to GINs, mainly FEC. The discovery of genomic regions, positional candidate genes, and QTLs associated with resistance to GINs can help increase and accelerate genetic gains in sheep breeding programs and reveal the genetic basis and biological mechanisms underlying this trait.

Genome-wide association studies for economically important traits in mink using copy number variation

Copy number variations (CNVs) are structural variants consisting of duplications and deletions of DNA segments, which are known to play important roles in the genetics of complex traits in livestock species. However, CNV-based genome-wide association studies (GWAS) have remained unexplored in American mink. Therefore, the purpose of the current study was to investigate the association between CNVs and complex traits in American mink. A CNV-based GWAS was performed with the ParseCNV2 software program using deregressed estimated breeding values of 27 traits as pseudophenotypes, categorized into traits of growth and feed efficiency, reproduction, pelt quality, and Aleutian disease tests. The study identified a total of 10,137 CNVs (6968 duplications and 3169 deletions) using the Affymetrix Mink 70K single nucleotide polymorphism (SNP) array in 2986 American mink. The association analyses identified 250 CNV regions (CNVRs) associated with at least one of the studied traits. These CNVRs overlapped with a total of 320 potential candidate genes, and among them, several genes have been known to be related to the traits such as ARID1B, APPL1, TOX, and GPC5 (growth and feed efficiency traits); GRM1, RNASE10, WNT3, WNT3A, and WNT9B (reproduction traits); MYO10, and LIMS1 (pelt quality traits); and IFNGR2, APEX1, UBE3A, and STX11 (Aleutian disease tests). Overall, the results of the study provide potential candidate genes that may regulate economically important traits and therefore may be used as genetic markers in mink genomic breeding programs.

Weighted single-step genome-wide association study and functional enrichment analyses for gastrointestinal nematode resistance traits in Santa Ines sheep

This study aimed to identify genomic regions, pathways, and putative candidate genes associated with resistance to gastrointestinal nematode in Santa Ines sheep. The phenotypic information comprised 5529 records from 1703 naturally infected animals. After genomic data quality control, 37,511 SNPs from 589 animals were available. The weighted single-step approach for genome-wide association study was performed to estimate the SNP effects and variances accounted by 10-SNP sliding windows. Confirming the polygenic nature of the studied traits, 20, 22, 21, and 19 genomic windows that explained more than 0.5% of the additive genetic variance were identified for fecal egg counts (FEC), Famacha© (FAM), packed cell volume (PCV), and total plasma protein (TPP), respectively. A total of 81, 122, 106, and 101 protein-coding genes were found in windows associated with FEC, FAM, PCV, and TPP, respectively. Several protein-coding genes related to the immune system and inflammatory response functions were identified within those genomic regions, such as ADCY9, ADRB2, BRAF, CADM1, CCL20, CD70, CREBBP, FNBP1, HTR4, IL16, IL22, IL26, MAPK8, NDFIP1, NLRC3, PAK5, PLCB1, PLCB4, ROCK1, TEK, TNFRSF12A, and VAV1. Functional enrichment analysis by DAVID tool also revealed many significant (P < 0.05) pathways and Gene Ontology terms that could be related to resistance to gastrointestinal nematode in Santa Ines sheep, such as chemokine signaling pathway (oas04062), cAMP signaling pathway (oas04024), cGMP-PKG signaling pathway (Oas04022), platelet activation (Oas04611), Rap1 signaling pathway (oas04015), and oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705). These results contribute to improving the knowledge of the genetic architecture of resistance to gastrointestinal nematode in Santa Ines sheep.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.