Strong selection signatures for Aleutian disease tolerance acting on novel candidate genes linked to immune and cellular responses in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by Aleutian mink disease virus (AMDV). This study aimed to identify candidate regions and genes underlying selection for response against AMDV using whole-genome sequence (WGS) data. Three case-control selection signatures studies were conducted between animals (N = 85) producing high versus low antibody levels against AMDV, grouped by counter immunoelectrophoresis (CIEP) test and two enzyme-linked immunosorbent assays (ELISA). Within each study, selection signals were detected using fixation index (FST) and nucleotide diversity (θπ ratios), and validated by cross-population extended haplotype homozygosity (XP-EHH) test. Within- and between-studies overlapping results were then evaluated. Within-studies overlapping results indicated novel candidate genes related to immune and cellular responses (e.g., TAP2, RAB32), respiratory system function (e.g., SPEF2, R3HCC1L), and reproduction system function (e.g., HSF2, CFAP206) in other species. Between-studies overlapping results identified three large segments under strong selection pressure, including two on chromosome 1 (chr1:88,770-98,281 kb and chr1:114,133-120,473) and one on chromosome 6 (chr6:37,953-44,279 kb). Within regions with strong signals, we found novel candidate genes involved in immune and cellular responses (e.g., homologous MHC class II genes, ITPR3, VPS52) in other species. Our study brings new insights into candidate regions and genes controlling AD response.

Epidemiology, pathogenesis, and diagnosis of Aleutian disease caused by Aleutian mink disease virus: A literature review with a perspective of genomic breeding for disease control in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by the Aleutian mink disease virus (AMDV). Commonly referred to as mink plasmacytosis, AD is an economically significant disease in mink-breeding countries. Aleutian disease mainly induces weight loss, lower fertility, and dropped pelt quality in adults and can result in acute interstitial pneumonia with high mortality rates in kits. In this review, we employed the scientific literature on AD over the last 70 years to discuss the historical and contemporary status of AD outbreaks and seroprevalence in mink farming countries. We also explained different forms of AD and the differences between the pathogenicity of the virus in kits and adults. The application of the available AD serological tests in AD control strategies was argued. We explained how selection programs could help AD control and proposed different approaches to selecting animals for building AD-tolerant herds. The advantages of genomic selection for AD tolerance over traditional breeding strategies were discussed in detail. We also explained how genomic selection could help AD control by selecting tolerant animals for the next generation based on genome-wide single nucleotide polymorphisms (SNP) data and the challenges of implementing genomic selection for AD tolerance in the mink industry. This review collected the information required for designing successful breeding programs for AD tolerance. Examples of the application of information are presented, and data gaps are highlighted. We showed that AD tolerance is necessary to be among the traits that animals are selected for in the mink industry.

Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls.

BACKGROUND: The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., FROH > 1 Mb) and excess of SNP homozygosity (FSNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants. RESULTS: Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in FROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting FROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility. CONCLUSIONS: Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

Population structure and genomic footprints of selection in five major Iranian horse breeds.

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log10 rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log10 rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle.

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The fdM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Detection of Common Copy Number of Variants Underlying Selection Pressure in Middle Eastern Horse Breeds Using Whole-Genome Sequence Data.

Dareshouri, Arabian, and Akhal-Teke are 3 Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole genome of 4 Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n = 3) and Akhal-Teke (n = 3) horses. Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size = 100 kb, step size = 50 kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8 837 950 single-nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and 3 quantitative trait loci, including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted P-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs that were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Weighted Single-Step GWAS for Body Mass Index and Scans for Recent Signatures of Selection in Yorkshire Pigs.

Controlling extra fat deposition is economically favorable in modern swine industry. Understanding the genetic architecture of fat deposition traits such as body mass index (BMI) can help in improving genomic selection for such traits. We utilized a weighted single-step genome-wide association study (WssGWAS) to detect genetic regions and candidate genes associated with BMI in a Yorkshire pig population. Three extended haplotype homozygosity (EHH)-related statistics were also incorporated within a de-correlated composite of multiple signals (DCMS) framework to detect recent selection signatures signals. Overall, the full pedigree consisted of 7016 pigs, of which 5561 had BMI records and 598 pigs were genotyped with an 80 K single nucleotide polymorphism (SNP) array. Results showed that the most significant windows (top 15) explained 9.35% of BMI genetic variance. Several genes were detected in regions previously associated with pig fat deposition traits and treated as potential candidate genes for BMI in Yorkshire pigs: FTMT, SRFBP1, KHDRBS3, FOXG1, SOD3, LRRC32, TSKU, ACER3, B3GNT6, CCDC201, ADCY1, RAMP3, TBRG4, CCM2. Signature of selection analysis revealed multiple candidate genes previously associated with various economic traits. However, BMI genetic variance explained by regions under selection pressure was minimal (1.31%). In conclusion, candidate genes associated with Yorkshire pigs' BMI trait were identified using WssGWAS. Gene enrichment analysis indicated that the identified candidate genes were enriched in the insulin secretion pathway. We anticipate that these results further advance our understanding of the genetic architecture of BMI in Yorkshire pigs and provide information for genomic selection for fat deposition in this breed.

Genomic Prediction of Average Daily Gain, Back-Fat Thickness, and Loin Muscle Depth Using Different Genomic Tools in Canadian Swine Populations.

Improvement of prediction accuracy of estimated breeding values (EBVs) can lead to increased profitability for swine breeding companies. This study was performed to compare the accuracy of different popular genomic prediction methods and traditional best linear unbiased prediction (BLUP) for future performance of back-fat thickness (BFT), average daily gain (ADG), and loin muscle depth (LMD) in Canadian Duroc, Landrace, and Yorkshire swine breeds. In this study, 17,019 pigs were genotyped using Illumina 60K and Affymetrix 50K panels. After quality control and imputation steps, a total of 41,304, 48,580, and 49,102 single-nucleotide polymorphisms remained for Duroc (n = 6,649), Landrace (n = 5,362), and Yorkshire (n = 5,008) breeds, respectively. The breeding values of animals in the validation groups (n = 392-774) were predicted before performance test using BLUP, BayesC, BayesCπ, genomic BLUP (GBLUP), and single-step GBLUP (ssGBLUP) methods. The prediction accuracies were obtained using the correlation between the predicted breeding values and their deregressed EBVs (dEBVs) after performance test. The genomic prediction methods showed higher prediction accuracies than traditional BLUP for all scenarios. Although the accuracies of genomic prediction methods were not significantly (P > 0.05) different, ssGBLUP was the most accurate method for Duroc-ADG, Duroc-LMD, Landrace-BFT, Landrace-ADG, and Yorkshire-BFT scenarios, and BayesCπ was the most accurate method for Duroc-BFT, Landrace-LMD, and Yorkshire-ADG scenarios. Furthermore, BayesCπ method was the least biased method for Duroc-LMD, Landrace-BFT, Landrace-ADG, Yorkshire-BFT, and Yorkshire-ADG scenarios. Our findings can be beneficial for accelerating the genetic progress of BFT, ADG, and LMD in Canadian swine populations by selecting more accurate and unbiased genomic prediction methods.

Prediction of biological age and evaluation of genome-wide dynamic methylomic changes throughout human aging.

The use of DNA methylation signatures to predict chronological age and aging rate is of interest in many fields, including disease prevention and treatment, forensics, and anti-aging medicine. Although a large number of methylation markers are significantly associated with age, most age-prediction methods use a few markers selected based on either previously published studies or datasets containing methylation information. Here, we implemented reproducing kernel Hilbert spaces (RKHS) regression and a ridge regression model in a Bayesian framework that utilized phenotypic and methylation profiles simultaneously to predict chronological age. We used over 450,000 CpG sites from the whole blood of a large cohort of 4409 human individuals with a range of 10-101 years of age. Models were fitted using adjusted and un-adjusted methylation measurements for cell heterogeneity. Un-adjusted methylation scores delivered a significantly higher prediction accuracy than adjusted methylation data, with a correlation between age and predicted age of 0.98 and a root mean square error (RMSE) of 3.54 years in un-adjusted data, and 0.90 (correlation) and 7.16 (RMSE) years in adjusted data. Reducing the number of predictors (CpG sites) through subset selection improved predictive power with a correlation of 0.98 and an RMSE of 2.98 years in the RKHS model. We found distinct global methylation patterns, with a significant increase in the proportion of methylated cytosines in CpG islands and a decreased proportion in other CpG types, including CpG shore, shelf, and open sea (P < 5e-06). Epigenetic drift seemed to be a widespread phenomenon as more than 97% of the age-associated methylation sites had heteroscedasticity. Apparent methylomic aging rate (AMAR) had a sex-specific pattern, with an increase in AMAR in females with age related to males.

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep.

Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation FST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the FST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.

Genetic and phenotypic parameters for body weights, harvest length, and growth curve parameters in American mink.

Understanding the genetics underlying growth curve is important for selection of animals with better growth potential, but little is known about the genetics of growth curve parameters in mink. This study estimated the genetic parameters for body weights (BWs), harvest length (HL), and growth parameters derived from the Richards model. For this purpose, individual BW of 1,088 mink measured seven times in 3-wk intervals (weeks 13 to 31 of life) were used for growth curve modeling using the Richards model. The BW traits included BW at week 13 (BW13), 16 (BW16), 19 (BW19), 22 (BW22), 25 (BW25), 28 (BW28), and 31 (BW31). Univariate analyses indicated that sex and birth-year had significant effects (P < 0.05) on BW, HL, asymptotic weight (α), growth rate at mature (k), shape parameter (m), weight at the inflection point (WIP), and age at the inflection point (AIP). In contrast, the color type had only significant effect (P < 0.05) on BW31 and HL. Estimated heritabilities (±SE) were ranged from 0.36 ± 0.13 (BW13) to 0.46 ± 0.10 (BW22) for BW and were 0.51 ± 0.09, 0.29 ± 0.09, 0.30 ± 0.09, 0.33 ± 0.1, 0.44 ± 0.10, and 0.47 ± 0.10 for HL, α, k, m, WIP, and AIP, respectively. The parameter α had non-significant (P > 0.05) genetic correlations (±SE) with k (-0.21 ± 0.23) and m (-0.10 ± 0.22), suggesting that changing shape parameters (k and m) will not influence asymptotic weight (α). Strong significant (P < 0.05) phenotypic (from 0.46 ± 0.03 to 0.60 ± 0.03) and genetic (0.70±0.13 to 0.88±0.09) correlations were observed between HL and different BW measures. The α, AIP, and WIP parameters had significant (P < 0.05) genetic correlations with HL indicated that selection for higher α, AIP, and WIP values would increase HL. Parameters k and m had nonsignificant (P > 0.05) genetic correlations with HL, indicating the change of the curve shape could not influence HL. Overall, the results suggest that growth curve parameters are heritable and can respond to genetic or genomic selection for optimizing the performance in mink.

Signatures of selection analysis using whole-genome sequence data reveals novel candidate genes for pony and light horse types.

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33, CRB1, CPN1, FAM13A, and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.

Whole-Genome Signatures of Selection in Sport Horses Revealed Selection Footprints Related to Musculoskeletal System Development Processes.

Selective breeding has led to gradual changes at the genome level of horses. Deciphering selective pressure patterns is progressive to understand how breeding strategies have shaped the sport horse genome; although, little is known about the genomic regions under selective pressures in sport horse breeds. The major goal of this study was to shed light on genomic regions and biological pathways under selective pressures in sport horses. In this study, whole-genome sequences of 16 modern sport and 35 non-sport horses were used to investigate the genomic selective signals of sport performance, by employing fixation index, nucleotide diversity, and Tajima's D approaches. A total number of 49 shared genes were identified using these approaches. The functional enrichment analysis for candidate genes revealed novel significant biological processes related to musculoskeletal system development, such as limb development and morphogenesis, having been targeted by selection in sport breeds.