Next Article in Journal
Integrative Analysis of Genomic and Clinical Data Reveals Intrinsic Characteristics of Bladder Urothelial Carcinoma Progression
Previous Article in Journal
Ankrd45 Is a Novel Ankyrin Repeat Protein Required for Cell Proliferation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 6
10.3390/genes10060463
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Hematological Traits in Polled Yak by Genome-Wide Association Studies Using Individual SNPs and Haplotypes

by
Xiaoming Ma
1,2,
Congjun Jia
1,2,
Donghai Fu
1,2,
Min Chu
1,2,
Xuezhi Ding
1,2,
Xiaoyun Wu
1,2,
Xian Guo
1,2,
Jie Pei
1,2,
Pengjia Bao
1,2,
Chunnian Liang
1,2,* and
Ping Yan
1,2,*
1
Animal Science Department, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
2
Key Laboratory for Yak Genetics, Breeding, and Reproduction Engineering of Gansu Province, Lanzhou 730050, China
*
Authors to whom correspondence should be addressed.
Genes 2019, 10(6), 463; https://doi.org/10.3390/genes10060463
Submission received: 10 April 2019 / Revised: 11 June 2019 / Accepted: 12 June 2019 / Published: 17 June 2019
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Yak (Bos grunniens) is an important domestic animal living in high-altitude plateaus. Due to inadequate disease prevention, each year, the yak industry suffers significant economic losses. The identification of causal genes that affect blood- and immunity-related cells could provide preliminary reference guidelines for the prevention of diseases in the population of yaks. The genome-wide association studies (GWASs) utilizing a single-marker or haplotype method were employed to analyze 15 hematological traits in the genome of 315 unrelated yaks. Single-marker GWASs identified a total of 43 significant SNPs, including 35 suggestive and eight genome-wide significant SNPs, associated with nine traits. Haplotype analysis detected nine significant haplotype blocks, including two genome-wide and seven suggestive blocks, associated with seven traits. The study provides data on the genetic variability of hematological traits in the yak. Five essential genes (GPLD1, EDNRA, APOB, HIST1H1E, and HIST1H2BI) were identified, which affect the HCT, HGB, RBC, PDW, PLT, and RDWSD traits and can serve as candidate genes for regulating hematological traits. The results provide a valuable reference to be used in the analysis of blood properties and immune diseases in the yak.

1. Introduction

Yak (Bos grunniens) is an iconic symbol of the Qinghai–Tibet Plateau and adjacent high-altitude areas [1]. This region is well known for its high-altitude, pristine natural environment, and extreme seasonal changes [2]. During the early Holocene period, the ancient Qiang people began to domesticate wild yaks [3,4]. While only a few animal species can survive in the Qinghai–Tibet Plateau, yaks thrive in the extreme conditions present in this area. They provide food such as meat and milk, a mode of transportation, shelter, and fuel for the residents of this high-altitude region [5,6]. As typical representatives of indigenous animals, yaks are characterized by a high degree of adaptability and low sensitivity to harsh ecological pressures [7,8]. However, due to relatively primitive breeding technology used by the herdsmen, and the difficulty or even lack of implementing immunizations and other measures preventing epidemics, yaks are very susceptible to parasites, which has caused significant economic losses to the yak-breeding industry.
Blood cells have a critical function in immune response and resistance to diseases [9,10]. They include three major cell types: leukocytes (white blood cells, WBCs), erythrocytes (red blood cells, RBCs), and platelets. WBCs are involved in the phagocytosis of bacteria and disease prevention. RBCs are essential for oxygen transport within the organism and the removal of the byproducts of respiration. The platelets play an important role in blood coagulation, wound healing, inflammatory response, and organ transplant rejection [11]. Blood homeostasis is the property responsible for the maintenance of routine indicators, like the number, volume, and biological activity of peripheral blood cells, within a relatively narrow physiological range. In pathologic conditions, these indicators are frequently outside of the normal physiological range [12]. In humans, these traits are to a large extent inherited but also influenced by environmental factors, such as age, gender, obesity, smoking, and alcohol consumption. Genetic and environmental factors can produce significant differences among individuals [13,14,15,16]. Therefore, the hematological indicators are often used as clinical indicators reflecting the overall health of patients.
Genome-wide association studies (GWASs) represent a powerful tool for mapping genetic loci associated with common disease and quantitative traits. Hematologic parameters are important indicators reflecting the immune status and health of the human or animal organism and are used as biomarkers in GWASs [17]. In recent years, GWASs have been used in human studies to locate more than 60 quantitative trait loci (QTLs) for blood traits among people of European ancestry, the Japanese and African Americans [18,19,20,21,22,23,24,25]. In the field of animal husbandry, pigs have been bred for disease resistance earlier than other species. Many QTLs affecting blood parameters traits in this species were also identified. A search of the animal QTL database [26] yielded 2482 QTLs affecting pig blood traits. Using whole genome-wide scanning in the White Duroc*Erhualian F2 resource population, Yang and coworkers located QTLs that affect the characteristics of WBCs and platelets on the Sus scrofa chromosomes (SSCs) 2, 7, 8, 12, 13, 15, and X [27]. Zou and collaborators identified QTLs affecting RBC traits at three growth stages in SSCs 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 13 [28]. Wang and collaborators used the porcine SNP60 gene chip for GWASs of 18 hematological traits (seven leukocyte traits, seven red blood cell traits, and four platelet traits) affecting piglet vaccine immunization. After the alignment based on multiple tests to adjust for each property characteristic of leukocytes, erythrocytes, and platelets, chromosomes 10, 24 and 77 on the SNP were found to be significant [29]. To the best of our knowledge, the GWAS approach was not used for bovine blood traits. The only available information is related to validation tests on haplotypes of related Leptin genes [30].
The current study employs the bovine HD SNP chip to locate the QTLs of blood traits in Polled yaks by GWAS analysis. The objective of this work is to provide crucial preliminary information necessary for the identification of genes affecting blood traits and causal mutations in yak.

2. Materials and Methods

All yaks involved in this study were in compliance with the guidelines for the care and use of laboratory animals issued by the State Council of the People’s Republic of China. Additionally, the current investigation was approved by the Animal Management and Ethics Committee of the Lanzhou Institute of Animal Science and Veterinary Medicine, Chinese Academy of Agricultural Sciences (Permit No. SYXK-2016-0039)

2.1. Animal Blood Sample Collection

The cohort of yaks used in this study was from the Ashdan Mountain area of Qinghai Province, China. Blood samples, 5 mL each, were collected from the jugular vein of 315 yaks. The animals had no genetic relationship to avoid inclusion of yaks with the same pedigree. DNA was extracted by the phenol–chloroform method [31], and the isolated DNA was dissolved in a Tris-EDTA (TE) buffer solution at pH = 8.0. The concentration and quality of DNA were measured by the NANODROP 1000 Nucleic Acid Protein Analyzer (Thermo Fisher Scientific Inc., Waltham, MA USA). The ratio of A 260/280 between 1.8 and 2.0 was considered to reflect acceptable quality. Qualified DNA samples were diluted to yield the concentration of 50 ng/μL and stored at −20 °C for subsequent SNP chip typing.
Standard sets of hematological data were recorded using the Sysmex pocH-100iV Diff whole blood analyzer (Sysmex Corporation, Japan) within 24 h of specimen collection at the Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences. The remaining portion of the blood sample was stored at −20 °C until DNA extraction (see above). To perform single-marker GWASs, 15 blood traits were used. They included 8 erythroid traits: hematocrit (HCT), hemoglobin (HGB), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV), red blood cell count (RBC), red blood cell volume distribution width-SD (RDW-SD), and red blood cell volume distribution width (RDW). Three leukocyte traits were also analyzed: lymphocyte count (LYM), white blood cell count (WBC), and medium white blood cell count (OTHR). The remaining four traits were related to the platelets and included platelet distribution width (PDW), platelet count (PLT), platelet-large cell ratio (P-LCR), and mean platelet volume (MPV). The presence of correlations between the 15 hematological parameters analyzed was established using the R psych package (http://personality-project.org/r/psych.manual.pdf).

2.2. SNP Array Genotyping and Quality Control

Genomic DNA was processed according to Illumina’s chip protocol and the chip was scanned using the iScan instrument (Illumina Inc., San Diego, CA, USA). The results of the chip scan were analyzed with GenomeStudio (version 2011.1) software. Quality control was performed using PLINK (version 1.07) [32] and results with a parameter of detection rates >95.0% and minor allele frequencies >0.01, HWE (p-value < 1×10−5), were excluded. Moreover, yaks with more than 10% of missing genotypes or with Mendelian errors greater than >5% were excluded. Only SNP markers that were consistent with family separation rules have been included in subsequent statistical analyses.

2.3. Statistical Analyses

The Bonferroni calibrated multiple test [33] was used to determine the significance threshold, with a genomic significant level threshold of 0.05 per effective SNP locus and a chromosomal significant horizontal threshold of 1 per effective SNP locus. For the purpose of this study, the corresponding thresholds were set as 9.5 × 10−6 (0.05/104,679) and 4.7 × 10−7 (1/104,679).

2.4. Single-Marker GWAS

For each SNP, the linear trend of alleles and phenotypic traits was assessed by the general linear mixed model [34,35]. The linear model contains polygenic random effects, and a variance–covariance matrix to form a structural Q matrix. The Q matrix was computed using admixture [36]. Since the CV value (Corss validation) was found to be the smallest when K = 3, therefore K = 3 was chosen for the Q matrix.
The model was described as follows:
y = Xbi + Qv + e
where y is the phenotype, bi is the regression coefficient of phenotype on SNP genotypes, X represents the vector of the SNP genotype indicators and assumes values of 0, 1, or 2 corresponding, respectively, to the three genotypes 11, 12 and 22, v is a variable dependent on population structure, Q is the corresponding principal components matrix, and e is the vector of residual errors with e, N (0, Iσ e2). The single-marker GWASs were conducted using the Tassel 5.2.43 software package [35].

2.5. Haplotype Analysis

Individual haplotype blocks and haplotype frequency were detected using the standard expectation maximization (EM) [37] algorithm. Reference thresholds and partitioning criteria for SNPs within closely linked haplotypes were established following the recommendations of the user manual of the PLINK software [32]. The recessive multi-alleles were used as a random effect in the association study model. The haplotype analysis model used in the current research is similar to that employed in the single-marker model, except for the replacement of the SNP effect in the single-marker GWAS by the haplotype effect.

3. Results

3.1. SNP Characteristics and Phenotype Statistics

The mean, standard deviation and coefficient of variation (CV) for 15 phenotypic observations of blood traits in the experimental population are presented in Table 1. The CV values ranged from a minimum of 2.92 to a maximum of 57.71 for MCHC and P-LCR, respectively.
After quality control, none of the animals had a genotyping call rate lower than 95%; all 315 yaks were included in the association analyses. A total of 42,276 SNPs failed the missingness test (GENO > 0.05), 653,185 SNPs had a minor allele frequency (MAF) of less than 0.05, and 8710 SNPs were severely departing from the Hardy Weinberg Equilibrium (HWE) (p-value less than 10−5). These SNPs were excluded from further analysis so that a total of 104,864 SNPs remained. Additional 185 SNPs, including unmapped SNPs, were also removed, yielding a total of 104,679 SNPs used for subsequent GWASs.

3.2. Loci for Erythrocyte Traits

Sixteen significant SNPs were found by using the single-marker GWAS method. They included nine suggested and seven genome-wide significant SNPs. The SNPs were identified for three erythrocyte traits, of which six were for HCT, five for HGB, and five for RBC (Table 2). Together, 16 single markers were detected which were located within five annotated genes. They include nine markers in the region of genes region, and one marker was located in a region 69,046 bp away from the nearest annotated gene. All identified markers were distributed on BTA (Bos taurus chromosome) 5, 16, 14, and 23. The most significant association of SNP was with the HCT trait (p = 6.21 × 10−8) (Figure 1).
In haplotype analysis, a total of six significant haplotype blocks, including two genome-wide and four suggestive blocks, were identified for four erythrocyte traits. One block each was associated with HCT, HGB, and RDWSD, while three blocks were associated with RBCs (Table 3). These significant SNPs were located on chromosome 23 (Figure 2). The most significant haplotype block was located in the GPLD1 gene and was associated with the HGB trait.

3.3. White Blood Cell Counts

By single-marker GWAS, seven single suggestive SNPs were found for two white blood cell counts traits. Three of them were identified for LYM and four for OTHR. The most significant (p = 7.71 × 10−7) marker, bovine HD0100040184 was associated with OTHR and was located in the ATP2C1 gene on chromosome 1 (Figure 1).
Haplotype analysis revealed only 1 significant haplotype block, which was associated with OTHR traits. The most significant annotation blocks located in the LOC785941 gene on chromosome 8 was also associated with OTHR traits (Figure 2).

3.4. Platelet Traits

Single-marker GWAS detected 14 significant suggestive SNPs. All of them were associated with PLT. Expected PLT traits markers were included within genes or within the nearest regions of genes; other trait makers were outside the gene region. The most significant (p = 3.87 × 10−7) PLT trait marker was bovine HD2100007261, located in the SH3GL3 gene on chromosome 4 (Figure 1).
Haplotype analysis detected two haplotype blocks for PLT. Both blocks were located on chromosome 23 at position from 31,565,453 to 31,686,075 bp position (Figure 2).

4. Discussion

The present investigation employed a single-marker and haplotype GWAS method based on the BovineHD BeadChip to analyze 15 hematological traits in polled yaks. The accumulated data resulted in the identification of 43 significant SNPs by single-marker GWAS, including 35 suggestive and eight genome-wide significant SNPs, as well as nine significant haplotype blocks, including seven suggestive and two genome-wide blocks, by haplotype analysis. The combination of single-marker and haplotype analyses revealed significant associations of 15 hematological quantitative traits in comparable genomic regions. However, the distribution of genomic P values obtained with these two methods were slightly different. This apparent lack of consistency points to the fact that the efficacy of the analysis is data dependent. Additionally, the results obtained with the single-marker method were characterized by more significant p-values than those acquired with the haplotype method. The reason for this discrepancy may be due to the close relationship between continuous SNPs located in the same region; clustering of SNPs may reduce the number of significant correlations. It is often difficult to weigh between increasing LD and increasing freedom which results in reduced detection efficiency. Furthermore, the LDs distributed in the genome are not uniform, i.e., in some places, the LD is high, whereas in others places the LD is low, so we compared the two methods to capture more associated SNPs.
The current investigation represents the first GWAS of hematological traits in the yak. Previous studies using this methodology in livestock were focused on hematological traits in pigs. By employing a similar method, Zhang and coworkers have determined that four SNPs on SSCs 2 have pleiotropic effects when single-marker analysis was used and are located at 24,777,963 bp on chromosome 7 when haplotype analysis was employed [38]. Seven genes have been selected as potential candidates after SNP annotations and expression variation were analyzed. Zhang studied selective White Duroc 6 Erhualian F2 intercross in three growth stages and found 185 significant SNPs for 18 hematological traits. In the GWAS analysis of the dairy cow, hematologic traits were not addressed, but the association of leptin with immunological [39] and hematological traits was identified. The study of Orrù and colleagues [30] revealed that haplotype polymorphism affects the hematologic variables of the dairy cow during the perinatal period. Our study identified several novel SNPs affecting blood traits. However, they did not overlap with significant loci of the same traits in other species. The possible reasons for this discrepancy include differences among species, the use of cattle chips that yielded only 10,486 SNPs after the quality control caused insufficient marker density, and the genetic complexity of blood cells.
The results of this GWAS showed a significant association pattern of erythrocyte phenotypes. The Manhattan plot pattern was similar for HGB, HCT, and RBC and the three erythrocyte traits which shared common regions at BTA23, from 32,982 Mb to 32,983 Mb, containing two SNPs (BovineHD2300009627 and BovineHD2300009628). Both of these parameters are directly related to Erythrocytes, respectively correlation analysis of hematological traits by comparison (Additional file 1: Table S1). Moreover, the results indicated a high correlation between the three erythrocyte traits (r1 = 0.96, r2 = 0.93, r3 = 0.94; p < 0.001). A significant degree of correlation implies that the QTL on BTA23 might be pleiotropic. However, this association was not found in leukocyte and platelet traits.
By applying a single-marker analysis, the current study identified 43 significant SNPs on 12 different chromosomes associated with 15 hematological traits (Additional file 1: Table S1). Together, 15 genes occupying regions from 63 to 77,309 bp were annotated. The subsequent verification of the function of these genes resulted in the selection of three genes as potential candidates. They included endothelin receptor type A (EDNRA), apolipoprotein B (APOB), and glycosylphosphatidylinositol-specific phospholipase D1 (GPLD1)—all functionally associated with blood-related cells or immune function.
The SNP BovineHD1700003099 located in the gene explained 13.35% of phenotypic variants of PDW (Table 2). EDNRA, a homologous isomer of the endothelin receptor, belongs to the rhodopsin family of the G protein-coupled receptor superfamily [40]. EDNRA is an ET-1 selective receptor, expressed mostly in vascular smooth muscle cells. Activation of EDNRA results in prolonged vasoconstriction. The endothelin receptors are closely related to the onset and development of cardiovascular and kidney diseases, diabetes mellitus, autoimmune diseases, and cancer [41,42,43]. The EDNRA gene is also involved in the vasoconstriction mechanism, which is closely related to pulmonary hypertension [44]. In addition, the results indicated that 14 markers were significantly correlated with the PLT traits, and the association of the BovineHD 1100022374 marker was most significant (p =1.68 × 10−6). Moreover, all these markers were located within the APOB gene. The expression of the APOB gene is positively correlated with the degree of atherosclerosis, and the composition and function of platelets in patients with hyperlipidemia is altered. These changes may be related to the effect of circulating lipoproteins on platelet phenotype [45]. Since yaks are chronically exposed to low atmospheric pressure and hypoxia, the compensatory increase in the number of erythrocytes leads to a higher blood viscosity, which results in the disorder of lipid metabolism and an elevation in blood lipid levels. In this regard, hypoxia promotes the formation of atherosclerotic plaque through the increase of permeability of vascular walls for lipids [46]. There is a coordination relationship between the platelet activation mechanism and the cholesterol-rich microdomain [47]. In the species of cattle with close relatives of yak, the base mutation and insertion mutation of the APOB gene can be found to cause cholesterol deficiency in the individual [48,49]. Furthermore, the data obtained in the present study demonstrate that two significant markers (BovineHD4100016160, BovineHD2300009636) overlapped and fell into the GPLD1 gene, which is also associated with erythrocyte traits. The product of the GPLD1 gene, also known as GPI-PLD and PLD, belongs to the glycosyl phosphatidyl inositol (GPI) family of proteolytic enzymes. The concentration of GPLD1 in human serum can be as high as 5–10 μg/mL [50]. GPI has a variety of anchoring proteins widely expressed in the cardiovascular system. During the onset and development of atherosclerosis, GPI acts not only as a physiological receptor, but also activates signaling pathways leading to growth or apoptosis of vascular cells. The deficiency of the homologous complement increases the sensitivity of erythrocytes, becoming the most important factor in the pathogenesis of paroxysmal nocturnal hemoglobinuria. It has been documented that the overexpression of GPLD1 increases its ability to anchor CD55 and CD59 by nearly 50%, resulting in the inhibition of cell proliferation, and increase in the rate of apoptosis and complement-mediated cell death. These effects are attributed to the blockade of the CD55 and CD59 signaling pathway and the decrease or inhibition of the resistance to complement activity, constituting an important factor in the immune escape of cancer cells [51]. Together, these findings point to the involvement of the GPLD1 gene in the pathways affecting the properties of red blood cells, making GPLD1 a candidate gene controlling the characteristics of red blood cells.
In the present work, nine significant haplotype blocks located mostly in BTA8, BTA23, and BTA24 were identified by haplotype analysis. Among them, three significant blocks for RBC were found within two annotated genes, GPLD1 and TINAg, located on BTA23. Also, RBC, HCT, and HGB were shown to have common significant blocks located on the GPLD1 gene. Another erythrocyte trait, RDW-SD, was associated with the block at BTA23, from 3150 to 3156 Mb containing seven SNPs (BovineHD2300009102, BovineHD2300009103, BovineHD2300009105, BovineHD2300009110, BovineHD2300009112, BovineHD2300009113 and BovineHD2300009117), a region that also includes the HIST1H1E and HIST1H2BI genes. These genes are associated with the function of blood and immune cells. Kordbacheh and coworkers have demonstrated that histones promote erythrocyte aggregation and sedimentation in a concentration-dependent manner, increasing their osmotic fragility and lysis [52]. In addition, histones negatively interfere with the denaturation of red blood cells by reducing their ability to pass through an artificial spleen. In vivo, histones induce anemia, increase hemoglobin content in the spleen, and decrease the number of platelets and white blood cells within minutes.

5. Conclusions

In summary, our research provides important information about the genetic determinants of hematological traits in yaks and points out five genes, namely, GPLD1, EDNRA, APOB, HIST1H1E, and HIST1H2BI, affecting the HCT, HGB, RBC, PDW, PLT, and RDWSD traits, and serving as candidate genes for the future dissection of molecular mechanisms regulating hematological traits. The present investigation provides a valuable reference for the analysis of blood properties in the yak and immune diseases impacting this species.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/6/463/s1, Table S1: Description of correlation and p-value among the hematological traits, Table S2: Description of all identified SNPs showing significant association with hematological traits by single-marker GWAS.

Author Contributions

X.M. and C.J. conceived and designed the experiments; D.F. performed the experiments; M.C. and J.P. analyzed the data; X.G. and X.D. contributed reagents/materials/analysis tools; X.W. and P.B. prepared the figures and/or tables; C.L. and P.Y. approved the final draft of the manuscript submitted for review and publication.

Funding

This study, including experiments design, samples collection, data analysis and manuscript writing were funded by the Agricultural Science and Technology Innovation Program (CAAS-ASTIP-2014-LIHPS-01), the National Beef Cattle Industry Technology & System (CARS-37), and the Herbivorous Livestock Industry and Technical System of Gansu Province (GARS08).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qiu, Q.; Zhang, G.; Ma, T.; Qian, W.; Wang, J.; Ye, Z.; Cao, C.; Hu, Q.; Kim, J.; Larkin, D.M. The yak genome and adaptation to life at high altitude. Nat. Genet. 2012, 44, 946. [Google Scholar] [CrossRef]
  2. Mizuno, S.; Ishizaki, T.; Toga, H.; Sakai, A.; Isakova, J.; Taalaibekova, E.; Baiserkeev, Z.; Kojonazarov, B.; Aldashev, A. Endogenous Asymmetric Dimethylarginine Pathway in High Altitude Adapted Yaks. Biomed. Res. Int. 2015, 2015, 196904. [Google Scholar] [CrossRef]
  3. Wang, Z.; Shen, X.; Liu, B.; Su, J.; Yonezawa, T.; Yu, Y.; Guo, S.; Ho, S.Y.; Vila, C.; Hasegawa, M. Phylogeographical analyses of domestic and wild yaks based on mitochondrial DNA: New data and reappraisal. J. Biogeogr. 2010, 37, 2332–2344. [Google Scholar] [CrossRef]
  4. Guo, S.; Savolainen, P.; Su, J.; Qian, Z.; Qi, D.; Jie, Z.; Yang, Z.; Zhao, X.; Liu, J. Origin of mitochondrial DNA diversity of domestic yaks. BMC Evol. Biol. 2006, 6, 73. [Google Scholar] [CrossRef]
  5. Wang, L.L.; Yu, Q.L.; Han, L.; Ma, X.L.; Song, R.D.; Zhao, S.N.; Zhang, W.H. Study on the effect of reactive oxygen species-mediated oxidative stress on the activation of mitochondrial apoptosis and the tenderness of yak meat. Food Chem. 2018, 244, 394. [Google Scholar] [CrossRef]
  6. Zi, X.D. Reproduction in female yaks (Bos grunniens) and opportunities for improvement. Theriogenology 2003, 59, 1303–1312. [Google Scholar] [CrossRef]
  7. Recavarren, S.; Arias-Stella, J. Right Ventricular Hypertrophy in People Born and Living at High Altitudes. Br. Heart J. 1964, 26, 806. [Google Scholar] [CrossRef]
  8. Alexander, A.F.; Jensen, R. Gross cardiac changes in cattle with high mountain (brisket) disease and in experimental cattle maintained at high altitudes. Am. J. Vet. Res. 1959, 20, 680–689. [Google Scholar]
  9. Greenough, P.R. Breeding for Disease Resistance in Farm Animals. Can. Vet. J. 2010, 33, 763. [Google Scholar] [CrossRef]
  10. Müller, M.; Brem, G. Disease resistance in farm animals. Experientia 1991, 47, 923. [Google Scholar] [CrossRef]
  11. Hoffman, R.; Benz Jr, E.J.; Silberstein, L.E.; Heslop, H.; Anastasi, J.; Weitz, J. Hematology: Basic Principles and Practice; Elsevier Health Sciences: Amsterdam, The Netherlands, 2013. [Google Scholar]
  12. Okada, Y.; Kamatani, Y. Common genetic factors for hematological traits in Humans. J. Hum. Genet. 2012, 57, 161–169. [Google Scholar] [CrossRef] [Green Version]
  13. Evans, D.M.; Frazer, I.H.; Martin, N.G. Genetic and environmental causes of variation in basal levels of blood cells. Twin Res. 1999, 2, 250–257. [Google Scholar] [CrossRef]
  14. Fisch, I.R.; Freedman, S.H. Smoking, oral contraceptives, and obesity: Effects on white blood cell count. JAMA 1975, 234, 500–506. [Google Scholar] [CrossRef]
  15. Hall, M.A.; Ahmadi, K.R.; Norman, P.; Snieder, H.; MacGregor, A.J.; Vaughan, R.W.; Spector, T.D.; Lanchbury, J.S. Genetic influence on peripheral blood T lymphocyte levels. Genes Immun. 2000, 1, 423. [Google Scholar] [CrossRef]
  16. Whitfield, J.B.; Martin, N.G. Genetic and environmental influences on the size and number of cells in the blood. Genet. Epidemiol. 2010, 2. [Google Scholar] [CrossRef]
  17. Kamatani, Y.; Matsuda, K.; Okada, Y.; Kubo, M.; Hosono, N.; Daigo, Y.; Nakamura, Y.; Kamatani, N. Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat. Genet. 2010, 42, 210–215. [Google Scholar] [CrossRef]
  18. Ganesh, S.K.; Zakai, N.A.; Van Rooij, F.J.; Soranzo, N.; Smith, A.V.; Nalls, M.A.; Chen, M.H.; Kottgen, A.; Glazer, N.L.; Dehghan, A.; et al. Multiple loci influence erythrocyte phenotypes in the CHARGE Consortium. Nat. Genet. 2009, 41, 1191–1198. [Google Scholar] [CrossRef] [Green Version]
  19. Jin, L.; Glessner, J.T.; Haitao, Z.; Cuiping, H.; Zhi, W.; Bradfield, J.P.; Mentch, F.D.; Yiran, G.; Cecilia, K.; Qianghua, X. GWAS of blood cell traits identifies novel associated loci and epistatic interactions in Caucasian and African-American children. Hum. Mol. Genet. 2013, 22, 1457–1464. [Google Scholar] [CrossRef]
  20. Zhao, Y.H.; Liu, Z.H.; Wang, L.; Xiao, H.M.; Du, C.G.; Zhang, Y.J.; Su, R.; Li, J.Q. Expression of the RORα gene in Inner Mongolian cashmere goat hair follicles. Genet. Mol. Res. 2015, 14, 380–388. [Google Scholar] [CrossRef]
  21. Soranzo, N.; Rendon, A.; Gieger, C.; Jones, C.; Watkins, N.; Menzel, S.; Doring, A.; Stephens, J.; Prokisch, H.; Erber, W.; et al. A novel variant on chromosome 7q22.3 associated with mean platelet volume, counts, and function. Blood 2009, 113, 3831. [Google Scholar] [CrossRef]
  22. Nalls, M.A.; Wilson, J.G.; Patterson, N.J.; Tandon, A.; Zmuda, J.M.; Huntsman, S.; Garcia, M.; Hu, D.; Li, R.; Beamer, B.A.; et al. Admixture mapping of white cell count: Genetic locus responsible for lower white blood cell count in the Health ABC and Jackson Heart studies. Am. J. Hum. Genet. 2008, 82, 81–87. [Google Scholar] [CrossRef]
  23. Meisinger, C.; Prokisch, H.; Gieger, C.; Soranzo, N.; Mehta, D.; Rosskopf, D.; Lichtner, P.; Klopp, N.; Stephens, J.; Watkins, N.A. A Genome-wide Association Study Identifies Three Loci Associated with Mean Platelet Volume. Am. J. Hum. Genet. 2009, 84, 66–71. [Google Scholar] [CrossRef] [Green Version]
  24. Yoichiro, K.; Sukanya, W.; Hidenori, O.; Takahisa, K.; Atsushi, T.; Naoya, H.; Michiaki, K.; Tatsuhiko, T.; Naoyuki, K.; Hiromitsu, K. A genome-wide association study identifies variants in the HLA-DP locus associated with chronic hepatitis B in Asians. Nat. Genet. 2009, 41, 591–595. [Google Scholar] [CrossRef]
  25. Soranzo, N.; Spector, T.D.; Mangino, M.; Kühnel, B.; Rendon, A.; Teumer, A.; Willenborg, C.; Wright, B.; Chen, L.; Li, M.; et al. A genome-wide meta-analysis identifies 22 loci associated with eight hematological parameters in the HaemGen consortium. Nat. Genet. 2009, 41, 1182–1190. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, Z.; Fritz-Waters, E.; Reecy, J. AnimalQTLdb: A livestock QTL database tool set for positional QTL information mining and beyond. Nucleic Acids Res. 2007, 35, 604–609. [Google Scholar] [CrossRef]
  27. Yang, S.; Ren, J.; Yan, X.; Huang, X.; Zou, Z.; Zhang, Z.; Yang, B.; Huang, L. Quantitative trait loci for porcine white blood cells and platelet-related traits in a White Duroc×Erhualian F2 resource population. Anim. Genet. 2009, 40, 273–278. [Google Scholar] [CrossRef]
  28. Zou, Z.; Ren, J.; Yan, X.; Huang, X.; Yang, S.; Zhang, Z.; Yang, B.; Li, W.; Huang, L. Quantitative trait loci for porcine baseline erythroid traits at three growth ages in a White Duroc × Erhualian F 2 resource population. Mamm. Genome 2008, 19, 640–646. [Google Scholar] [CrossRef]
  29. Wang, J.Y.; Luo, Y.R.; Fu, W.X.; Lu, X.; Zhou, J.P.; Ding, X.D.; Liu, J.F.; Zhang, Q. Genome-wide association studies for hematological traits in swine. Anim. Genet. 2013, 44, 34–43. [Google Scholar] [CrossRef]
  30. Orrù, L.; Abeni, F.; Catillo, G.; Grandoni, F.; Crisà, A.; De Matteis, G.; Scatà, M.; Napolitano, F.; Moioli, B. Leptin gene haplotypes are associated with change in immunological and hematological variables in dairy cow during the peripartum period. J. Anim. Sci. 2012, 90, 16–26. [Google Scholar] [CrossRef]
  31. Moore, D.D.; Dowhan, D. Preparation and analysis of DNA. In Current Protocols in Molecular Biology; Ausubel, F.M., Brent, R., Kingston, R.E., Moore, D.D., Seidman, J.G., Smith, J.A., Struhl, R., Eds.; John Wiley & Sons: New York, NY, USA, 1997; Volume 1, pp. 241–245. ISBN 0-471-50338-X. [Google Scholar]
  32. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, Q.; Cui, J.; Chazaro, I.; Cupples, L.A.; Demissie, S. Power and type I error rate of false discovery rate approaches in genome-wide association studies. BMC Genet. 2005, 6, 1–4. [Google Scholar] [CrossRef]
  34. Breslow, N.E.; Clayton, D.G. Approximate Inference in Generalized Linear Mixed Models. J. Am. Stat. Assoc. 1993, 88, 9–25. [Google Scholar] [CrossRef]
  35. Bradbury, P.; Zhang, Z.; Kroon, D.; Casstevens, T.Y.; Buckler, E. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  36. Alexander, D.H.; Novembre, J.K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef] [Green Version]
  37. Long, J.; Williams, R.; Urbanek, M. An E-M algorithm and testing strategy for multiple-locus haplotypes. Am. J. Hum. Genet. 1995, 56, 799–810. [Google Scholar] [CrossRef]
  38. Zhang, F.; Zhang, Z.; Yan, X.; Chen, H.; Zhang, W.; Yuan, H.; Huang, L. Genome-wide association studies for hematological traits in Chinese Sutai pigs. BMC Genet. 2014, 15, 41. [Google Scholar] [CrossRef]
  39. Zhang, Z.; Hong, Y.; Gao, J.; Xiao, S.; Ma, J.; Zhang, W.; Ren, J.; Huang, L. Genome-wide association study reveals constant and specific loci for hematological traits at three time stages in a White Duroc × Erhualian F2 resource population. PLoS ONE 2013, 8, e63665. [Google Scholar] [CrossRef]
  40. Hoeltke, C.; Kopka, W.K.; Heindel, W.; Schober, O.; Schaefers, M.; Bremer, C. Biodistribution of A Nonpeptidic Fluorescent Endothelin A Receptor Imaging Probe. Mol. Imaging 2009, 8, 27–34. [Google Scholar] [CrossRef]
  41. Yanagisawa, M.; Kurihara, H.; Kimura, S.; Tomobe, Y.; Kobayashi, M.; Mitsui, Y.; Yazaki, Y.; Goto, K.; Masaki, T. A novel potent vasoconstrictor peptide produced by vascular endothelial cells. Nature 1988, 6163. [Google Scholar] [CrossRef]
  42. Maclean, M.R.; Mcculloch, K.M.; Baird, M. Endothelin ETA- and ETB-receptor-mediated vasoconstriction in rat pulmonary arteries and arterioles. J. Cardiovasc. Pharm. 1994, 23, 838–845. [Google Scholar] [CrossRef]
  43. Watts, S.W. Endothelin receptors: what’s new and what do we need to know? Am. J. Physiol. Regul. Integr. Comp. Physiol. 2010, 298, 254–260. [Google Scholar] [CrossRef]
  44. Nagalakshmi, S.; Derk, C.T. Endothelin receptor antagonists as disease modifiers in systemic sclerosis. Inflamm. Allergy Drug Targets 2011, 10. [Google Scholar] [CrossRef]
  45. Cleland, J.G.; Daubert, J.C.; Erdmann, E.; Freemantle, N.; Gras, D.; Kappenberger, L.; Tavazzi, L. The effect of cardiac resynchronization on morbidity and mortality in heart failure. Dig. World Core Med. J. 2005, 114, 216–218. [Google Scholar] [CrossRef]
  46. Liu, J.Q. Variations of levels of lipoprotein, platelet and zinc in 95 acute myocardial infarction patients in high plateau. China Trop. Med. 2010, 1001, 99–109. [Google Scholar] [CrossRef]
  47. Bodin, S.; Tronchère, H.; Payrastre, B. Lipid rafts are critical membrane domains in blood platelet activation processes. Biochim. Et Biophys. Acta 2003, 1610, 247–257. [Google Scholar] [CrossRef] [Green Version]
  48. Menzi, F.; Besuchet-Schmutz, N.; Fragnière, M.; Hofstetter, S.; Jagannathan, V.; Mock, T.; Raemy, A.; Studer, E.; Mehinagic, K.; Regenscheit, N.; et al. A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle. Anim. Genet. 2016, 47. [Google Scholar] [CrossRef]
  49. Gross, J.; Schwinn, A.C.; Schmitz-Hsu, F.; Menzi, F.; Drögemüller, C.; Albrecht, C.; Bruckmaier, R. Rapid Communication: Cholesterol deficiency–associated APOB mutation impacts lipid metabolism in Holstein calves and breeding bulls. J. Anim. Sci. 2016, 94. [Google Scholar] [CrossRef]
  50. Gondi, C.S.; Lakka, S.N.; Siddique, K.; Dinh, D.H.; Olivero, W.C.; Gujrati, M.; Rao, J.S. Expression of antisense uPAR and antisense uPA from a bicistronic adenoviral construct inhibits glioma cell invasion, tumor growth, and angiogenesis. Oncogene 2003, 22, 5967–5975. [Google Scholar] [CrossRef] [Green Version]
  51. Tang, J.H.; He, W.J.; Huang, H.; Tan, C.C.; Duan, Q.; Wang, K.J.; Yuan, X.Y.; Zhu, X.J. Important roles of glycosylphosphatidylinositol (GPI)-specific phospholipase D and some GPI-anchored proteins in the pathogenesis of hepatocellular carcinoma. Clin. Biochem. 2009, 42, 400–407. [Google Scholar] [CrossRef]
  52. Kordbacheh, F.; O’Meara, C.H.; Coupland, L.A.; Lelliott, P.M.; Parish, C.R. Extracellular histones induce erythrocyte fragility and anemia. Blood 2017, 130. [Google Scholar] [CrossRef]
Genes 10 00463 g001 550
Figure 1. A portion of the GWAS results of hematological traits at a significant level in Polled yak. X axis indicates chromosome, Y axis indicates −lg (p-value), the solid line indicates genome-wide significant level and the dashed line indicates the chromosome-wide significance level, the bars below each chromosome are the density of the markers. RBC: red blood cell count. HCT: hematocrit; HGB: hemoglobin; OTHR: medium white blood cell count; PLT: platelet count; MPV: mean platelet volume; PLCR: platelet count; LYM: lymphocyte count; PDW: platelet distribution width.(A–I) Manhattan plot representing HCT, LYM, PDW, HGB, OTHR, PLCR, RBC, MPV and PLT traits, respectively.
Figure 1. A portion of the GWAS results of hematological traits at a significant level in Polled yak. X axis indicates chromosome, Y axis indicates −lg (p-value), the solid line indicates genome-wide significant level and the dashed line indicates the chromosome-wide significance level, the bars below each chromosome are the density of the markers. RBC: red blood cell count. HCT: hematocrit; HGB: hemoglobin; OTHR: medium white blood cell count; PLT: platelet count; MPV: mean platelet volume; PLCR: platelet count; LYM: lymphocyte count; PDW: platelet distribution width.(A–I) Manhattan plot representing HCT, LYM, PDW, HGB, OTHR, PLCR, RBC, MPV and PLT traits, respectively.
Genes 10 00463 g001
Genes 10 00463 g002 550
Figure 2. Manhattan plots for the haplotype analysis of hematological traits surpass the genome-wide significant threshold. Log 10 (1/p-value) values are shown for all SNPs that passed quality control. The numbers indicate the chromosomes in the genome. The solid line and dotted line denote the Bonferroni-corrected genome-wide and suggestive significant threshold, respectively. SNPs surpassing the genome-wide threshold are highlighted in red and SNPs reaching the suggestive threshold in green. The bars below each chromosome are the density of the markers. E.g., Manhattan of HCT means Manhattan map obtained by haplotype analysis of HCT traits. HCT: hematocrit; RBC: red blood cell count; HGB: hemoglobin; OTHR: medium white blood cell count; PLT: platelet count; RDW-SD: red blood cell volume distribution width-SD; PLCR: platelet count.(A–G) Manhattan plot representing HCT, HGB, OTHR, RBC, PLT, RDWSD and PLCR traits, respectively.
Figure 2. Manhattan plots for the haplotype analysis of hematological traits surpass the genome-wide significant threshold. Log 10 (1/p-value) values are shown for all SNPs that passed quality control. The numbers indicate the chromosomes in the genome. The solid line and dotted line denote the Bonferroni-corrected genome-wide and suggestive significant threshold, respectively. SNPs surpassing the genome-wide threshold are highlighted in red and SNPs reaching the suggestive threshold in green. The bars below each chromosome are the density of the markers. E.g., Manhattan of HCT means Manhattan map obtained by haplotype analysis of HCT traits. HCT: hematocrit; RBC: red blood cell count; HGB: hemoglobin; OTHR: medium white blood cell count; PLT: platelet count; RDW-SD: red blood cell volume distribution width-SD; PLCR: platelet count.(A–G) Manhattan plot representing HCT, HGB, OTHR, RBC, PLT, RDWSD and PLCR traits, respectively.
Genes 10 00463 g002
Table
Table 1. Descriptive statistics of 15 hematological traits in the yak population.
Table 1. Descriptive statistics of 15 hematological traits in the yak population.
TraitsAbbreviationValue 1 (No. 2)CV 3
Mean corpuscular hemoglobin concentration (g/L)MCHC318.061 ± 9.281 (313)2.92
Red blood cell volume distribution width-SD (fL)RDWSD31.427 ± 2.362 (315)7.52
Mean corpuscular volume (fl)MCV42.725 ± 3.25 (315)7.61
Red cell distribution coefficient of variationRDWCV0.187 ± 0.015 (315)7.75
Mean platelet volume (fl)MPV7.146 ± 0.627 (171)8.77
Mean corpuscular hemoglobin (pg)MCH13.611 ± 1.312 (315)9.64
Red blood cell count (1012/L)RBC10.031 ± 1.266 (314)12.62
Hematocrit (%)HCT0.429 ± 0.056 (309)12.83
Hemoglobin (g/L)HGB136.217 ± 18.577 (314)13.64
Platelet distribution width (%)PDW8.777 ± 1.314 (169)14.97
White blood cell count (109/L)WBC9.04 ± 2.197 (278)24.31
Medium white blood cell count (109/L)OTHR4.603 ± 1.505 (270)32.69
Lymphocyte count (109/L)LYM4.538 ± 1.588 (276)34.99
Platelet count (109/L)PLT326.193 ± 170.149 (311)52.17
Platelet-large cell ratio (%)PLCR0.068 ± 0.039 (111)57.71
1 Values are shown in mean ± standard deviation. 2 The numbers of recorded individuals are given in parentheses. 3 CV: coefficient of variation.
Table
Table 2. Associated SNPs and nearby candidate genes for hematological traits.
Table 2. Associated SNPs and nearby candidate genes for hematological traits.
Traits 1SNP 2CHR 3Position 4p ValueVAR (%) 5Nearest Gene 6
NameDistance (bp)
HCT *BOVINEHD05000239785847849898.51× 10−6 7.44IFLTD169,046
HCT **BOVINEHD410001616023329946556.21× 10−8 10.48GPLD1WITHIN
HCT **BOVINEHD230000963623329965216.55 × 10−8 10.38GPLD1WITHIN
HGB **BOVINEHD410001616023329946552.63× 10−7 9.48GPLD1WITHIN
HGB **BOVINEHD230000963623329965212.76 × 10−7 9.38GPLD1WITHIN
RBC *BOVINEHD140001675314603951189.25 × 10−7 7.52LOC782927WITHIN
RBC *BOVINEHD230000963623329965215.40 × 10−6 7.59GPLD1WITHIN
RBC *BOVINEHD410001616023329946555.84 × 10−6 7.59GPLD1WITHIN
RBC *BOVINEHD23000015452363752806.80 × 10−6 7.45C23H6ORF14242
RBC *BOVINEHD23000015502363831727.65 × 10−6 7.38C23H6ORF14242
MPV *BOVINEHD100002855110986904343.00 × 10−6 14.31LOC78848113,887
MPV *BOVINEHD200000700920233057435.68 × 10−6 11.73IL31RA63
PDW *BOVINEHD170000309917107781561.98 × 10−6 13.35EDNRAWITHIN
PLT *BOVINEHD110002236211778976261.68 × 10−6 7.29APOBWITHIN
PLT *BOVINEHD210000726121250449563.87 × 10−79.3SH3GL311,156
PLT *BOVINEHD24000027332496635135.51 × 10−6 7.69LOC53895877,309
LYM *BOVINEHD200000412020128404342.75 × 10−6 9.13LOC529061WITHIN
LYM *BOVINEHD210001767121610247854.74 × 10−6 8.76LOC781495WITHIN
LYM *BOVINEHD230000578223218704231.11 × 10−6 8.48MUT68,056
OTHR *BOVINEHD010004018411405432317.71 × 10−710.17ATP2C120,601
OTHR *BOVINEHD0300001040334615676.51 × 10−6 8.7LOC1002951702804
OTHR *BOVINEHD03000213953735432624.99 × 10−6 8.89NEGR1WITHIN
OTHR *BOVINEHD08000057808184171079.50 × 10−6 8.44TUSC1151,827
1 The abbreviations of hematological traits are given in Table 1, e.g., MCV is mean corpuscular volume. 2 The number of significant SNPs for each hematological trait. 3,4 Chromosomal locations and positions of the most significant SNP associated with hematological traits in Bos taurus UMD_3.1 assembly. 5 Interpretation rate of phenotypic variation of the marker. 6 The nearest annotated gene to the significant SNP. The annotated gene database is from http://asia.ensembl.org/index.html. SNP designated as in a gene or distance (bp) from a gene region in Bos taurus UMD_3.1 assembly. ** Genome-wide significant. * Chromosome-wide significant.
Table
Table 3. Associated haplotype blocks and nearby candidate genes for hematological traits.
Table 3. Associated haplotype blocks and nearby candidate genes for hematological traits.
Traits 1NSNP 2CHR3START 4END 5SNP1 6SNP2 7VAR (%) 8p ValueNearest Gene 9
NameDistance (bp)
RBC *22360451416051410BovineHD2300001442BovineHD23000014447.417.16 × 10−6TINAGWITHIN
RBC *22360451416051410BovineHD2300001442BovineHD23000014447.397.45 × 10−6TINAGWITHIN
RBC *2233298219932983832BovineHD2300009627BovineHD23000096287.585.40 × 10−6GPLD1WITHIN
HGB **2233298219932983832BovineHD2300009627BovineHD23000096289.382.76 × 10−7GPLD1WITHIN
HCT **2233298219932983832BovineHD2300009627BovineHD230000962810.376.55 × 10−8GPLD1WITHIN
PLT *22477519397753085BovineHD2400002173ARS-BFGL-NGS-1125397.695.51 × 10−6LOC100336967WITHIN
OTHR *581771684517825120BovineHD4100006566ARS-BFGL-NGS-100898.439.50 × 10−6LOC7859418868
RDWSD *7233150445231558355BovineHD2300009102BovineHD23000091036.438.17 × 10−6HIST1H1E,HIST1H2BIWITHIN
PLCR *2233293576932941851BovineHD2300009596BovineHD230000959919.082.53 × 10−6LOC520603WITHIN
1 The abbreviations of hematological traits are given in Table 1, e.g., RBC is red blood cell count. 2 Number of SNPs included in haplotype block. 3,4,5 Chromosomal locations and positions of the most significant SNP associated with hematological traits in Bos taurus UMD_3.1 assembly, and starting position and ending position on the chromosome. 6,7 The first SNP and the second SNP in the haplotype block. 8 Interpretation rate of phenotypic variation of the marker. 9 The nearest annotated gene to the significant SNP. The annotated gene database is from http://asia.ensembl.org/index.html. SNP designated as in a gene or distance (bp) from a gene region in Bos taurus UMD_3.1 assembly. ** Genome-wide significant. * Chromosome-wide significant.

Share and Cite

MDPI and ACS Style

Ma, X.; Jia, C.; Fu, D.; Chu, M.; Ding, X.; Wu, X.; Guo, X.; Pei, J.; Bao, P.; Liang, C.; et al. Analysis of Hematological Traits in Polled Yak by Genome-Wide Association Studies Using Individual SNPs and Haplotypes. Genes 2019, 10, 463. https://doi.org/10.3390/genes10060463

AMA Style

Ma X, Jia C, Fu D, Chu M, Ding X, Wu X, Guo X, Pei J, Bao P, Liang C, et al. Analysis of Hematological Traits in Polled Yak by Genome-Wide Association Studies Using Individual SNPs and Haplotypes. Genes. 2019; 10(6):463. https://doi.org/10.3390/genes10060463

Chicago/Turabian Style

Ma, Xiaoming, Congjun Jia, Donghai Fu, Min Chu, Xuezhi Ding, Xiaoyun Wu, Xian Guo, Jie Pei, Pengjia Bao, Chunnian Liang, and et al. 2019. "Analysis of Hematological Traits in Polled Yak by Genome-Wide Association Studies Using Individual SNPs and Haplotypes" Genes 10, no. 6: 463. https://doi.org/10.3390/genes10060463

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop