Genome-wide selection signatures detection in Shanghai Holstein cattle population identified genes related to adaption, health and reproduction traits

Abstract

Background Over several decades, a wide range of natural and artificial selection events in response to subtropical environments, intensive pasture and intensive feedlot systems have greatly changed the customary behaviour, appearance, and important economic traits of Shanghai Holstein cattle. In particular, the longevity of the Shanghai Holstein cattle population is generally short, approximately the 2nd to 3rd lactation. In this study, two complementary approaches, integrated haplotype score (iHS) and runs of homozygosity (ROH), were applied for the detection of selection signatures within the genome using genotyping by genome-reduced sequence data from 1092 cows. Results In total, 101 significant iHS genomic regions containing selection signatures encompassing a total of 256 candidate genes were detected. There were 27 significant |iHS| genomic regions with a mean |iHS| score > 2. The average number of ROH per individual was 42.15 ± 25.47, with an average size of 2.95 Mb. The length of 78 % of the detected ROH was within the range of 1–2 MB and 2–4 MB, and 99 % were shorter than 8 Mb. A total of 168 genes were detected in 18 ROH islands (top 1 %) across 16 autosomes, in which each SNP showed a percentage of occurrence > 30 %. There were 160 and 167 genes associated with the 52 candidate regions within health-related QTL intervals and 59 candidate regions within reproduction-related QTL intervals, respectively. Annotation of the regions harbouring clustered |iHS| signals and candidate regions for ROH revealed a panel of interesting candidate genes associated with adaptation and economic traits, such as IL22RA1, CALHM3, ITGA9, NDUFB3, RGS3, SOD2, SNRPA1, ST3GAL4, ALAD, EXOSC10, and MASP2. In a further step, a total of 1472 SNPs in 256 genes were matched with 352 cis-eQTLs in 21 tissues and 27 trans-eQTLs in 6 tissues. For SNPs located in candidate regions for ROH, a total of 108 cis-eQTLs in 13 tissues and 4 trans-eQTLs were found for 1092 SNPs. Eighty-one eGenes were significantly expressed in at least one tissue relevant to a trait (P value < 0.05) and matched the 256 genes detected by iHS. For the 168 significant genes detected by ROH, 47 gene-tissue pairs were significantly associated with at least one of the 37 traits. Conclusions We provide a comprehensive overview of selection signatures in Shanghai Holstein cattle genomes by combining iHS and ROH. Our study provides a list of genes associated with immunity, reproduction and adaptation. For functional annotation, the cGTEx resource was used to interpret SNP-trait associations. The results may facilitate the identification of genes relevant to important economic traits and can help us better understand the biological processes and mechanisms affected by strong ongoing natural or artificial selection in livestock populations.

Genome-Wide Association Study and Selective Sweep Analysis Reveal the Genetic Architecture of Body Weights in a Chicken F2 Resource Population

“…As an importantly economic trait, body weight in chickens has undergone long-term artificial selection in the past decades, which is expected to left selective signatures on chicken genomics. The detection of selection signatures can expedite the identification of genes responsible for important economic traits and better understanding the biological mechanisms affected by strong ongoing natural or artificial selection in livestock populations ( 36 , 37 ). Accordingly, in this study selective signature analysis was utilized to confirm overlapping genomic regions detected by GWAS to screen important candidate genes for chicken BW.…”

Transcriptomic Analysis of Circulating Leukocytes Obtained during the Recovery from Clinical Mastitis Caused by Escherichia coli in Holstein Dairy Cows

“…The genes containing 12 of these variants were then found to be differentially expressed in leucocytes between the E. coli (+) EARLY and CONT cows. Among these, three DEGs ( WIPI1 , ARSG and SLC16A6 ) located together on BTA 19 came from a genomic region with an extended run of homozygosity (ROH) which was detected in a population of Shanghai Holstein cattle [ 59 ]. Analyses of ROH allow the identification of genomic regions with possible selection signatures for the breed, as their size and frequency vary according to population diversity and selection pressure.…”

Identification of Runs of Homozygosity Islands and Genomic Estimated Inbreeding Values in Caqueteño Creole Cattle (Colombia)

“…The database was composed of the genotypes of 127 Creole animals genotyped with the BovineHD BeadChip 777,000 SNPs (Illumina, San Diego, CA, USA). Assuming an MAF < 0.01, the PLINK program was applied for quality control of the genotypes [ 8 , 12 , 16 ], a call rate for SNPs < 0.90, and a call rate for samples < 0.98.…”

“…Genomic analysis of the population structure indicates the proportion of the total genetic variability represented by the associated principal component (PC), obtained by the plink-pca command line of the PLINK program [17]. 0.01, the PLINK program was applied for quality control of the genotypes [8,12,16], a call rate for SNPs < 0.90, and a call rate for samples < 0.98.…”