Genome‐wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome‐wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F₂ progeny were analyzed for each blood‐related trait after filtering for quality control. Data were analyzed by the genome‐wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10^?6) on SSC3, 6, 8, 13 and 17 were identified for blood‐related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.     

Genome-wide association study of porcine hematological parameters in a Large White × Minzhu F2 resource population

Hematological traits, which are important indicators of immune function in animals, have been commonly examined as biomarkers of disease and disease severity in humans and animals. Genome-wide significant quantitative trait loci (QTLs) provide important information for use in breeding programs of animals such as pigs. QTLs for hematological parameters (hematological traits) have been detected in pig chromosomes, although these are often mapped by linkage analysis to large intervals making identification of the underlying mutation problematic. Single nucleotide polymorphisms (SNPs) are the common form of genetic variation among individuals and are thought to account for the majority of inherited traits. In this study, a genome-wide association study (GWAS) was performed to detect regions of association with hematological traits in a three-generation resource population produced by intercrossing Large White boars and Minzhu sows during the period from 2007 to 2011. Illumina PorcineSNP60 BeadChip technology was used to genotype each animal and seven hematological parameters were measured (hematocrit (HCT), hemoglobin (HGB), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV), red blood cell count (RBC) and red blood cell volume distribution width (RDW)). Data were analyzed in a three step Genome-wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) method. A total of 62 genome-wide significant and three chromosome-wide significant SNPs associated with hematological parameters were detected in this GWAS. Seven and five SNPs were associated with HCT and HGB, respectively. These SNPs were all located within the region of 34.6-36.5 Mb on SSC7. Four SNPs within the region of 43.7-47.0 Mb and fifty-five SNPs within the region of 42.2-73.8 Mb on SSC8 showed significant association with MCH and MCV, respectively. At chromosome-wide significant level, one SNP at 29.2 Mb on SSC1 and two SNPs within the region of 26.0-26.2 Mb were found to be significantly associated with RBC and RDW, respectively. Many of the SNPs were located within previously reported QTL regions and appeared to narrow down the regions compared with previously described QTL intervals. In current research, a total of seven significant SNPs were found within six candidate genes SCUBE3, KDR, TDO, IGFBP7, ADAMTS3 and AFP. In addition, the KIT gene, which has been previously reported to relate to hematological parameters, was located within the region significantly associated with MCH and MCV and could be a candidate gene. These results of this study may lead to a better understanding of the molecular mechanisms of hematological parameters in pigs.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Genetic association analysis highlights new loci that modulate hematological trait variation in Caucasians and African Americans

Red blood cell, white blood cell, and platelet measures, including their count, sub-type and volume, are important diagnostic and prognostic clinical parameters for several human diseases. To identify novel loci associated with hematological traits, and compare the architecture of these phenotypes between ethnic groups, the CARe Project genotyped 49,094 single nucleotide polymorphisms (SNPs) that capture variation in ~2,100 candidate genes in DNA of 23,439 Caucasians and 7,112 African Americans from five population-based cohorts. We found strong novel associations between erythrocyte phenotypes and the glucose-6 phosphate dehydrogenase (G6PD) A-allele in African Americans (rs1050828, P<2.0×10(-13), T-allele associated with lower red blood cell count, hemoglobin, and hematocrit, and higher mean corpuscular volume), and between platelet count and a SNP at the tropomyosin-4 (TPM4) locus (rs8109288, P=3.0×10(-7) in Caucasians; P=3.0×10(-7) in African Americans, T-allele associated with lower platelet count). We strongly replicated many genetic associations to blood cell phenotypes previously established in Caucasians. A common variant of the α-globin (HBA2-HBA1) locus was associated with red blood cell traits in African Americans, but not in Caucasians (rs1211375, P<7×10(-8), A-allele associated with lower hemoglobin, mean corpuscular hemoglobin, and mean corpuscular volume). Our results show similarities but also differences in the genetic regulation of hematological traits in European- and African-derived populations, and highlight the role of natural selection in shaping these differences.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Genome‐wide association study for reproductive traits in a Large White pig population

Using the PorcineSNP80 BeadChip, we performed a genome‐wide association study for seven reproductive traits, including total number born, number born alive, litter birth weight, average birth weight, gestation length, age at first service and age at first farrowing, in a population of 1207 Large White pigs. In total, we detected 12 genome‐wide significant and 41 suggestive significant SNPs associated with six reproductive traits. The proportion of phenotypic variance explained by all significant SNPs for each trait ranged from 4.46% (number born alive) to 11.49% (gestation length). Among them, 29 significant SNPs were located within known QTL regions for swine reproductive traits, such as corpus luteum number, stillborn number and litter size, of which one QTL region associated with litter size contained the ALGA0098819 SNP for total number born. Subsequently, we found that 376 functional genes contained or were near these significant SNPs. Of these, 14 genes—BHLHA15, OCM2, IL1B2, GCK, SMAD2, HABP2, PAQR5, GRB10, PRELID2, DMKN, GPI, GPIHBP1, ADCY2 and ACVR2B—were considered important candidates for swine reproductive traits based on their critical roles in embryonic development, energy metabolism and growth development. Our findings contribute to the understanding of the genetic mechanisms for reproductive traits and could have a positive effect on pig breeding programs.     

Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population

Pork quality is an economically important trait and one of the main selection criteria for breeding in the swine industry. In this genome-wide association study (GWAS), 455 pigs from a porcine Large White × Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for intramuscular fat content (IMF), marbling, moisture, color L*, color a*, color b* and color score in the longissimus muscle (LM). Association tests between each trait and the SNPs were performed via the Genome Wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approach. From the Ensembl porcine database, SNP annotation was implemented using Sus scrofa Build 9. A total of 45 SNPs showed significant association with one or multiple meat quality traits. Of the 45 SNPs, 36 were located on SSC12. These significantly associated SNPs aligned to or were in close approximation to previously reported quantitative trait loci (QTL) and some were located within introns of previously reported candidate genes. Two haplotype blocks ASGA0100525-ASGA0055225-ALGA0067099-MARC0004712-DIAS0000861, and ASGA0085522-H3GA0056170 were detected in the significant region. The first block contained the genes MYH1, MYH2 and MYH4. A SNP (ASGA0094812) within an intron of the USP43 gene was significantly associated with five meat quality traits. The present results effectively narrowed down the associated regions compared to previous QTL studies and revealed haplotypes and candidate genes on SSC12 for meat quality traits in pigs.     

QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F₂ intercross between Landrace and Korean native pigs. Eight blood‐related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F₂ progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome‐wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood‐related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker‐assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.     

A genome‐wide association study suggests several novel candidate genes for carcass traits in Chinese Simmental beef cattle

A genome‐wide association study (GWAS) was conducted for two carcass traits in Chinese Simmental beef cattle. The experimental population consisted of 1301 individuals genotyped with the Illumina BovineHD SNP BeadChip (770K). After quality control, 671 990 SNPs and 1217 individuals were retained for the GWAS. The phenotypic traits included carcass weight and bone weight, which were measured after the cattle were slaughtered at 16 to 18 months of age. Three statistical models—a fixed polygene model, a random polygene model and a composite interval mapping polygene model—were used for the GWAS. The genome‐wide significance threshold after Bonferroni correction was 7.44E‐08 (= 0.05/671 990). In this study, we detected eight and seven SNPs significantly associated with carcass weight and bone weight respectively. In total, 11 candidate genes were identified within or close to these significant SNPs. Of these, we found several novel candidate genes, including PBX1, GCNT4, ALDH1A2, LCORL and WDFY3, to be associated with carcass weight and bone weight in Chinese Simmental beef cattle, and their functional roles need to be verified in further studies.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

猪2、7和8号染色体上影响血常规指标的数量性状基因座(QTL)检测

以3个品种（长白猪、大白猪、松辽黑猪）16个公猪家系共计368头仔猪组成资源群体,在猪2、7和8号染色体上共选取35个微卫星标记,采用基于线性混合模型的方差组分分析方法,对影响与猪白细胞、红细胞和血小板相关的共计18项血常规指标的数量性状基因座（quantitative trait loci,QTL）进行了检测．通过似然比检验,并以自由度为2的卡方分布作为检验统计量的分布,共发现22个在P〈5％水平下显著的QTL,其中在2号染色体上有9个,分别影响白细胞总数、中性粒细胞数、平均红细胞体积、血红蛋白含量、平均红细胞血红蛋白浓度、血小板总数、平均血小板体积、血小板分布宽度和血小板压积,在7号染色体有7个,分别影响白细胞总数、中性粒细胞数、平均红细胞血红蛋白浓度、血红蛋白含量、血小板总数、平均红细胞体积和红细胞分布宽度变异,在8号染色体上有6个,分别影响中性粒细胞百分比、淋巴细胞百分比、平均红细胞血红蛋白浓度、血小板总数、血小板压积和平均红细胞体积．为尽可能地避免由于多重检验所造成的假阳性率的升高,我们采用了控制假检出率（false discovery rate,FDR）的方法来对这22个QTL进行进一步检验,发现有14个达到FDR〈5％显著水平,其中又有9个达到FDR〈1％显著水平．     

Association between the ABO locus and hematological traits in Korean

ABSTRACT: BACKGROUND: Recently, genomewide association studies identified a pleiotropic gene locus, ABO, as being significantly associated with hematological traits. To confirm the effects of ABO on hematological traits, we examined the link between the ABO locus and hematological traits in Korean population-based cohorts. RESULTS: Six tagging SNPs for ABO were analyzed with regard to their effects on hematological traits [white blood cell count (WBC), red blood cell count (RBC), platelet (Plat), mean corpuscular volume (MCV), and mean corpuscular haemoglobin concentration (MCHC)]. Linear regression analyses were performed, controlling for recruitment center, sex, and age as covariates. Of the 6 tagging SNPs, 3 (rs2073823, rs8176720, and rs495828) and 3 (rs2073823, rs8176717, and rs687289) were significantly associated with RBC and MCV, respectively (Bonferroni correction p-value criteria < 0.05/6 = 0.008). rs2073823 and a reported SNP (rs8176746), as well as rs495828 and a reported SNP (rs651007), showed perfect linkage disequilibrium status (r2s = 0.99). Of the remaining 3 SNPs (rs8176720, rs8176717 and rs687289), rs8176717 generated an independent signal with moderate p-value (= 0.045) when it was adjusted for by rs2073823 (the most significant SNP). We also identified a copy number variation (CNV) that was tagged by the SNP rs8176717, the minor allele of which correlated with the deletion allele of CNV. Our haplotype analysis indicated that the haplotype that contained the CNV deletion was significantly associated with MCV (beta +/- se = 0.363 +/- 0.118, p =2.09 x 10-3). CONCLUSIONS: Our findings confirm that ABO is one of the genetic factors that are associated with hematological traits in the Korean population. This result is notable, because GWASs fail to evaluate the link between a CNV and phenotype traits.     

Association analysis of polymorphisms of porcine LMP2 and LMP7 genes with haematological traits

Low molecular weight polypeptides 2 (LMP2) and low molecular weight polypeptides 7 (LMP7) are located within the major histocompatibility complex and have been associated with autoimmune disease. In this study, polymorphisms of porcine LMP2 and LMP7 genes were analyzed by PCR-SSCP and DNA sequencing methods. Four SNPs (DQ659151:g.2115T>C; DQ659151:g.4343A>G; DQ872631:g.1232C>G; DQ872631:g.2847C>T) were identified. Four SNPs of genes were analyzed for association with 22 haematological traits in Large White (n = 195), Landrace (n = 84) and Songliao Black (n = 86) pig population. Of all the 22 traits, seven were significant associated with the SNPs of LMP2/LMP7 gene (P < 0.05). They included white blood cell count (WBC) (P = 0.028), neutrophilic granulocyte count (GRAN) (P = 0.037), monocytes percentage (MO%) (P = 0.015), red blood cell (RBC) (P = 0.004), red blood cell volume distribution width (RDW) (P = 0.004), mean platelet volume (MPV) (P = 0.016) and CD4(+)CD8(+)% (P = 0.045). These results suggest LMP2/LMP7 gene should be regarded as molecular marker to estimate animal's immune status for their effects on hematological traits.     

Genome-wide association study-based identification genes influencing agronomic traits in rice (Oryza sativa L.)

Rice is one of the most important cereal crops, providing the daily dietary intake for approximately 50% of the global human population. Here, we re-sequenced 259 rice accessions, generating 1371.65 Gb of raw data. Furthermore, we performed genome-wide association studies (GWAS) on 13 agronomic traits using 2.8 million single nucleotide polymorphisms (SNPs) characterized in 259 rice accessions. Phenotypic data and best linear unbiased prediction (BLUP) values of each of the 13 traits over two years of each trait were used for the GWAS. The results showed that 816 SNP signals were significantly associated with the 13 agronomic traits. Then we detected candidate genes related to target traits within 200 kb upstream and downstream of the associated SNP loci, based on linkage disequilibrium (LD) blocks in the whole rice genome. These candidate genes were further identified through haplotype block constructions. This comprehensive study provides a timely and important genomic resource for breeding high yielding rice cultivars.     

Whole‐genome association study for energy balance and fat/protein ratio in German Holstein bull dams

The metabolic status of cows is important to health and fertility, especially in early lactation, and energy balance (EB) and fat/protein ratio (FPR) are considered as appropriate indicators for metabolic disorders. The aim of this study was to detect SNPs (single nucleotide polymorphisms) associated with EB and FPR in German Holstein bull dams belonging to the research herd Karkendamm. Bull dams were genotyped using the Illumina Bovine SNP 50K Bead chip^® comprising 54 001 SNPs. A total of 43 593 SNPs and 586 (EB) and 668 (FPR) bull dams passed the quality control criteria. Phenotypes were deregressed breeding values estimated via random regression animal models for lactation days 11, 20, 30, and 42 for EB and FPR. Whole‐genome association analyses were carried out fitting principal components as covariates to adjust for genetic substructure. Permutation tests were applied to estimate genome‐wise significance. Across all observed lactation days, 19 SNPs located in four different intervals on chromosomes 1, 14, 16, and 27 were detected. For EB, seven markers across four chromosomes were identified. There was no overlap between markers associated with FPR and EB. SNPs associated with FPR were mostly located in QTL regions for milk production traits, especially in the region of DGAT1, whereas SNPs associated with EB mainly showed positional relationships to previously described QTL regions affecting functional traits, especially fertility.     

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.     

Next‐generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan)

To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing‐based bulked segregant analysis (Seq‐BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R‐ and S‐bulks with the help of draft genome sequence and reference‐guided assembly of ICPL 20096 (resistant parent). Seq‐BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re‐sequenced and their combined analysis with R‐ and S‐bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2‐Mb flanking regions of seven candidate SNPs identified through Seq‐BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re‐sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics‐assisted breeding in pigeonpea.