Selection signature reveals genes associated with susceptibility loci affecting respiratory disease due to pleiotropic and hitchhiking effect in Chinese indigenous pigs

Identifying selection signatures for immune response and resilience to Aleutian disease in mink using genotype data

Aleutian disease (AD) brings tremendous financial losses to the mink industry. Selecting AD-resilient mink has been conducted to control AD. Such selections could have altered the patterns of genetic variation responding to selection pressures. This study aimed to identify selection signatures for immune response (IRE) and resilience to AD. A total of 1,411 mink from an AD-positive facility were used. For IRE, 264 animals were categorized according to the combined results of enzyme-linked immunosorbent assay (ELISA) and counterimmunoelectrophoresis (CIEP). For resilience, two grouping methods were used: 1) general resilience performance (GRP, n = 30) was evaluated based on the feed conversion ratio, Kleiber ratio, and pelt quality; and 2) female reproductive performance (FRP, n = 36) was measured based on the number of kits alive 24 h after birth. Detection methods were the pairwise fixation index, nucleotide diversity, and cross-population extended haplotype homozygosity. A total of 619, 569, and 526 SNPs were identified as candidates for IRE, GRP, and FRP, respectively. The annotated genes were involved in immune system process, growth, reproduction, and pigmentation. Two olfactory-related Gene Ontology (GO) terms were significant (q < 0.05) for all traits, suggesting the impact of AD on the sense of smell of infected mink. Differences in detected genes and GO terms among different color types for IRE indicated variations in immune response to AD among color types. The mitogen-activated protein kinase (MAPK) signaling pathway was significant (q < 0.05) for FRP, suggesting that AD may disrupt MAPK signaling and affect FRP. The findings of this research contribute to our knowledge of the genomic architecture and biological mechanisms underlying AD resilience in mink.

Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Identifying selection signatures for litter size in Guangxi Bama Xiang pigs

Litter size is an important economic trait in pig production. However, the genetic mechanisms underlying varying litter size in Guangxi Bama Xiang pigs remain unknown. To identify selection signatures for litter size in Guangxi Bama Xiang pigs, we obtained 297 Illumina PorcineSNP50 BeadChip array data and the average born number (ABN) from parity one to nine in Guangxi Bama Xiang pigs. Fixation index (Fst) methods were used to identify the selection signature of the litter size, and three phenotypic gradient differential population pairs (according to the ABN) in individuals were used to reduce the false positives of signature selections. Single nucleotide polymorphisms (SNPs) were identified in the VEGFA promoter and exons. The general linear model was used to analyse the differences in distinct genotypes after they were typed using three-round multiplex PCR technology. Finally, the transcriptome factor and CpG island in the VEGFA promoter were predicted. A total of 328, 328 and 317 significant loci were identified in the 1st, 2nd and 3rd population pairs, respectively. After removing the false positives, 25 SNPs were defined as the selection signatures in relation to litter size. Ten (VEGFA, USP49, USP25, SRPK1, SLC26A8, RPL10A, PPARD, MAPK14, HMGA1 and CHRDL2) out of 52 genes in the selection regions were annotated as the candidate genes of litter size, respectively, VEGFA. There were no SNPs in the VEGFA exon region, but we obtained three SNPs (rs786889605, rs343769603 and rs323942424) in the VEGFA promoter regions. The ABN in CC was significantly higher than that in TT in rs786889605, and the ABN in TT was significantly lower than that in GG in rs323942424. Meanwhile, the mutation of the VEGFA promoter result in the loss of Sp1 and NF-1 and the formation of Oct-1. In summary, we obtained ten candidate genes, and two mutations in the VEGFA promoter that could be important potential molecular biomarkers for litter size in Bama Xiang pigs.

Hitchhiking Mapping of Candidate Regions Associated with Fat Deposition in Iranian Thin and Fat Tail Sheep Breeds Suggests New Insights into Molecular Aspects of Fat Tail Selection

Simple Summary

Fatness-related traits are economically very important in sheep production and are associated with serious diseases in humans. Using a denser set of SNP markers and a variety of statistical approaches, our results were able to refine the regions associated with fat deposition and to suggest new insights into molecular aspects of fat tail selection. These results may provide a strong foundation for studying the regulation of fat deposition in sheep and do offer hope that the causal mutations and the mode of inheritance of this trait will soon be discovered by further investigation.

Abstract

The fat tail is a phenotype that divides indigenous Iranian sheep genetic resources into two major groups. The objective of the present study is to refine the map location of candidate regions associated with fat deposition, obtained via two separate whole genome scans contrasting thin and fat tail breeds, and to determine the nature of the selection occurring in these regions using a hitchhiking approach. Zel (thin tail) and Lori-Bakhtiari (fat tail) breed samples that had previously been run on the Illumina Ovine 50 k BeadChip, were genotyped with a denser set of SNPs in the three candidate regions using a Sequenom Mass ARRAY platform. Statistical tests were then performed using different and complementary methods based on either site frequency (F_ST and Median homozygosity) or haplotype (iHS and XP-EHH). The results from candidate regions on chromosome 5 and X revealed clear evidence of selection with the derived haplotypes that was consistent with selection to near fixation for the haplotypes affecting fat tail size in the fat tail breed. An analysis of the candidate region on chromosome 7 indicated that selection differentiated the beneficial alleles between breeds and homozygosity has increased in the thin tail breed which also had the ancestral haplotype. These results enabled us to confirm the signature of selection in these regions and refine the critical intervals from 113 kb, 201 kb, and 2831 kb to 28 kb, 142 kb, and 1006 kb on chromosome 5, 7, and X respectively. These regions contain several genes associated with fat metabolism or developmental processes consisting of TCF7 and PPP2CA (OAR5), PTGDR and NID2 (OAR7), AR, EBP, CACNA1F, HSD17B10,SLC35A2, BMP15, WDR13, and RBM3 (OAR X), and each of which could potentially be the actual target of selection. The study of core haplotypes alleles in our regions of interest also supported the hypothesis that the first domesticated sheep were thin tailed, and that fat tail animals were developed later. Overall, our results provide a comprehensive assessment of how and where selection has affected the patterns of variation in candidate regions associated with fat deposition in thin and fat tail sheep breeds.

Identifying selection signatures and runs of homozygosity for spine curvature in Chinese indigenous pigs

In China there are approximately 100 pig breeds, which show great diversity in their appearance. However, information on genome selection signatures, such as spine curvature, is scarce. Therefore, we used the fixation index (F_ST ) and cross-population extended haplotype homozygosity (XPEHH) methods to explore the genome selection signatures of spine curvature in six breeds of Chinese indigenous pig. We identified 396 and 389 single nucleotide polymorphisms using the F_ST and XPEHH methods, respectively. We detected 19 selection signatures and 28 genes located in the selected regions. Five candidate genes (MAP3K7, CUX1, GRIN2B, ALPL and MACF1) were identified in the selection signatures. Additionally, 719 high-frequency runs of homozygosity regions, 17 unique runs of homozygosity regions, 78 genes and 27 pathways were identified in the runs of homozygosity analysis. The TGF-beta signaling pathway and eight genes related to the spine formation, spine defects and intervertebral disk degeneration were identified, comprising ACVR1, FMOD, ITGA4, MAPK8, PDGF, RPL3, SULF1 and UBE2D1. In summary, we identified 13 candidate genes related to spine curvature in Chinese indigenous pigs.

Establishment of a model of Mycoplasma hyopneumoniae infection using Bama miniature pigs

Mycoplasma hyopneumoniae (M. hyopneumoniae), is the primary aetiological agent of enzootic pneumonia leading to chronic respiratory disease prevalent worldwide. Conventional pigs are the only animals used for pathogenicity studies and vaccine evaluations of M. hyopneumoniae. Considering that the challenge animals have better genetic stability and a smaller body size to operate with, an alternative experimental animal model of M. hyopneumoniae infection with Bama miniature pigs was established. Nine seven-week-old snatch-farrowed, porcine colostrum-deprived (SF-pCD) Bama miniature pigs and nine conventional pigs were randomly divided into two infected groups (Bama miniature-infected (BI) and conventional-infected groups (CI), BI and CI, n = 6) and two control groups (Bama miniature control (BC) and conventional control (CC) groups, BC and CC, n = 3). Every piglet was tracheally inoculated with 5 × 108 CCU/mL containing 10% suspension of a stock of frozen lung homogenate from SF-pCD pigs infected with virulent strain JS or sterilized KM2 medium. Typical lung lesions appeared in all infected pigs after necropsy, and the mean gross lung lesions was 17.3 and 13.7 in groups of BI and CI. Serum IgG and nasal sIgA antibody titres were increased significantly. Cilia shedding and mucus staining increased greatly in JS-infected bronchi. Obvious reddish gross lesions and M. hyopneumoniae antigen were detected, especially apparently observed in group of BI. Moreover, DNA copies of M. hyopneumoniae from bronchoalveolar lavage fluid (BALF) of each JS-infected piglet reached more than 108, and M. hyopneumoniae could be re-isolated from each infected BALF. These results indicate that Bama miniature pigs could be used as an alternative and more maneuverable experimental infection model for M. hyopneumoniae and display typical clinical and pathological features consistent with those in conventional pigs.

Signature of genome wide gene expression in classical swine fever virus infected macrophages and PBMCs of indigenous vis-a-vis crossbred pigs

The genetic basis of differential host immune response vis-à-vis transcriptome profile was explored in PBMCs of indigenous (Ghurrah) and crossbred pigs after classical swine fever vaccination and in monocyte derived macrophages (MDMs) challenged with virulent classical swine fever (CSF) virus. The humoral immune response (E2 antibody) was higher (74.87%) in crossbred than indigenous pigs (58.20%) at 21st days post vaccination (21dpv). The rate of reduction of ratio of CD4⁺/CD8⁺ was higher in crossbred pigs than indigenous pigs at 7th days post vaccination (7dpv). The immune genes IFIT1, IFIT5, RELA, NFKB2, TNF and LAT2 were up regulated at 7dpv in RNA seq data set and was in concordance during qRT-PCR validation. The Laminin Subunit Beta 1 (LAMB1) was significantly (p ≤ 0.05) down-regulated in MDMs of indigenous pigs and consequently a significantly (p ≤ 0.01) higher copy number of virulent CSF virus was evidenced in macrophages of crossbred pigs than indigenous pigs. Activation of LXR:RXR pathway at 60 h post infection (60hpi) in MDMs of indigenous versus crossbred pigs inhibited nuclear translocation of NF-κB, resulted into transrepression of proinflammatory genes. But it helped in maintenance of HDL level by lowering down cholesterol/LDL level in MDMs of indigenous pigs. The key immune genes (TLR2, TLR4, IL10, IL8, CD86, CD54, CASP1) of TREM1 signaling pathway were upregulated at 7dpv in PBMCs but those genes were downregulated at 60hpi in MDMs indigenous pigs. Using qRT-PCR, the validation of differentially expressed, immunologically important genes (LAMB1, OAS1, TLR 4, TLR8 and CD86) in MDMs revealed that expression of these genes were in concordance with RNA-seq data.