奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状,随着平衡育种理念的提出和发展,繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展,全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状, 随着平衡育种理念的提出和发展, 繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展, 全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

Keeping profitable cows in the herd: A lifetime cost-benefit assessment to support culling decisions

Increasing the productive lifespan of dairy cows is important to achieve a sustainable dairy industry, but making strategic culling decisions based on cow profitability is challenging for farmers. The objective of this study was to carry out a lifetime cost-benefit analysis based on production and health records and to explore different culling decisions among farmers. The cost-benefit analysis was conducted for 22 747 dairy cows across 114 herds in Quebec, Canada for which feed costs and the occurrence of diseases were reported. Costs and revenues related to productive lifespan were compared among cohorts of cows that left their respective herd at the end of their last completed lactation or stayed for a complete additional lactation. Hierarchical clustering analysis was carried out based on costs and revenues to explore different culling decisions among farmers. Our results showed that the knowledge of lifetime cumulative costs and revenues was of great importance to identify low-profitable cows at an earlier lactation, while only focusing on current lactation costs and revenues can lead to an erroneous assessment of profitability. While culling decisions were mostly based on current lactation costs and revenues and disregarded the occurrence of costly events on previous lactations, there was variation among farmers as we identified three different culling decision clusters. Monitoring cumulative costs and revenues would help farmers to identify low-profitable cows at an earlier lactation and make the decision to increase herd productive lifespan and farm profitability by keeping the most profitable cows.     

A high recombination rate in eusocial Hymenoptera: evidence from the common wasp Vespula vulgaris

Background

Longevity expressed as the number of days between birth and death is a trait of great importance for both human and animal populations. In our analysis we use dairy cattle to demonstrate how the association of Single Nucleotide Polymorphisms (SNPs) located within selected genes with longevity can be modeled. Such an approach can be extended to any genotyped population with time to endpoint information available. Our study is focused on selected genes in order to answer the question whether genes, known to be involved into the physiological determination of milk production, also influence individual's survival.

Results

Generally, the highest risk differences among animals with different genotypes are observed for polymorphisms located within the leptin gene. The polymorphism with a highest effect on functional longevity is LEP-R25C, for which the relative risk of culling for cows with genotype CC is 3.14 times higher than for the heterozygous animals. Apart from LEP-R25C, also FF homozygotes at the LEP-Y7F substitution attribute 3.64 times higher risk of culling than the YY homozygotes and VV homozygotes at LEP-A80V have 1.83 times higher risk of culling than AA homozygotes. Differences in risks between genotypes of polymorphisms within the other genes (the butyrophilin subfamily 1 member A1 gene, BTN1A1; the acyl-CoA:diacylglycerol acyltransferase 1 gene, DGAT1; the leptin receptor gene, LEPR; the ATP-binding cassette sub-family G member 2, ABCG2) are much smaller.

Conclusions

Our results indicate association between LEP and longevity and are very well supported by results of other studies related to dairy cattle. In view of the growing importance of functional traits in dairy cattle, LEP polymorphisms should be considered as markers supporting selection decisions. Furthermore, since the relationship between both LEP polymorphism and its protein product with longevity in humans is well documented, with our result we were able to demonstrate that livestock with its detailed records of family structure, genetic, and environmental factors as well as extensive trait recording can be a good model organism for research aspects related to humans.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Significant differences in fertility between dairy cows selected for one QTL located on bovine chromosome 3 are not attributable to energy balance,although eating behaviour is affected

Improvement of reproduction in dairy cows has become a major challenge in dairy production. We have recently shown that dairy cows carrying the ‘fertil−’ haplotype for one quantitative trait locus (QTL), affecting female fertility and located on the bovine chromosome 3, had a significantly lower conception rate after the first artificial insemination than cows carrying the ‘fertil+’ haplotype. The objective of this paper was to study other phenotypic modifications linked to this QTL. In the present study, 23 ‘fertil+’ and 18 ‘fertil−’ cows were characterized for live weight, milk production, food intake, eating behaviour and plasma metabolites. These parameters were measured during the first lactation, from calving to 40 weeks postpartum (wkpp). In the first 7 weeks of lactation, ‘fertil+’ primiparous cows had a significantly higher live BW and milk production than ‘fertil−’ cows. Dry matter intake tended to be slightly higher for ‘fertil+’ than for ‘fertil−’ primiparous cows in this period. However, energy balance was similar for the two haplotypes in the whole lactation, except in the first wkpp, and consequently, could not explain their different fertility. The major observation concerned the eating behaviour. ‘Fertil+’ primiparous cows had a significantly lower eating rate than ‘fertil−’ cows during the 40 weeks of lactation. In parallel, ‘fertil+’ cows spent significantly more time at the feeder for a similar number of visits than ‘fertil−’ cows. Furthermore, no differences in plasma concentrations of non-esterified fatty acids and insulin were observed between the two haplotypes. Plasma glucose was significantly lower in ‘fertil+’ than in ‘fertil−’ cows in the second wkpp. Taken together, our results show that ‘fertil+’ and ‘fertil−’ dairy cows, with different fertility, have also different eating behaviour without any variation in energy balance, except in the first week of lactation.     

Genome-wide association studies for agronomical traits in winter rice accessions of Assam

The 297 winter rice accessions of Assam, North East India were genotyped by sequencing (GBS). The 50,985 high-quality SNPs were filtered and assigned to 12 rice chromosomes. The population structure analysis revealed three major subgroups SG1, SG2, and SG3 consisting of 30, 8, and 143 accessions respectively. The remaining 116 accessions were grouped as admixture population. Phenotypic data were recorded on13 agronomical traits for genome-wide association studies (GWAS). The 60 significant marker-trait associations (MTAs) were identified for 11 agronomical traits, which explained 0 to 15% of phenotypic variance (PV). A QTL ‘hot spot’ was detected near the centromeric region on chromosome 6. The identified QTLs may be validated and utilized in ‘genomics assisted breeding’ for improvement of existing rice cultivars of Assam and North East India.     

Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample population size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein–Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with traditional measures of fertility, as well as a fertility phenotype derived from milk progesterone profiles. Traditional fertility measures investigated were days to first heat, days to first service, pregnancy rate to first service, number of services and calving interval (CI); post-partum interval to the commencement of luteal activity (CLA) was derived using routine milk progesterone assays. Phenotypic and genotypic data on 37 590 single nucleotide polymorphisms (SNPs) were available for up to 1570 primiparous cows. Genetic parameters were estimated using linear animal models, and univariate and bivariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. Heritability estimates of the traditional fertility traits varied from 0.03 to 0.16; the heritability for CLA was 0.13. The posterior quantitative trait locus (QTL) probabilities, across the genome, for the traditional fertility measures were all <0.021. Posterior QTL probabilities of 0.060 and 0.045 were observed for CLA on SNPs each on chromosome 2 and chromosome 21, respectively, in the univariate analyses; these probabilities increased when CLA was included in the bivariate analyses with the traditional fertility traits. For example, in the bivariate analysis with CI, the posterior QTL probability of the two aforementioned SNPs were 0.662 and 0.123. Candidate genes in the vicinity of these SNPs are discussed. The results from this study suggest that the power of genome-wide association studies in cattle may be increased by sharing of data and also possibly by using physiological measures of the trait under investigation.     

Associations between the breeding values of Holstein-Friesian bulls and longevity and culling reasons of their daughters

Taking into account functional traits in the breeding practice should lead to a longer productive life of cows. However, despite the increased contribution of these traits in bull selection indices, their daughters are frequently culled as early as the 2nd or 3rd lactation. The problem is whether and to what extent the genetic potential of animals is realized in the production practice. Therefore, the purpose of this study was to determine the associations between the breeding value (BV) of bulls and their daughters for cow longevity and culling reasons in the Holstein-Friesian cattle population in Poland. Data for 532 062 cows culled in 2012, 2015, and 2018 were analyzed. A majority of 5 045 cow sires originated from Poland, Germany, France, the Netherlands, and the United States. The highest variation in the contribution of culling reasons was for the cows culled at the age of 2–4 years. The contribution of the culling reasons, analyzed in relation to the cow culling age, remained similar and the only exception was culling because of old age, for which a significant increase was observed only for the culling age of at least 9 years (13.8%), which was reached by only 7.3% of the cows. The sires were characterized by generally high BV for conformation and reproductive traits. However, they had, at most, the average genetic potential for functional longevity. There were a number of beneficial associations found between the BV of bulls and the distribution of culling reasons in their daughters. For example, it concerns relations between the somatic cell score in milk and culling due to udder diseases and low milk yield, between the interval from calving to first insemination and low milk yield, between the protein yield and old age, or between the BV for certain conformation traits (size, udder) and cow culling due to age. In these cases, as the BV increased for a given trait, the contribution of the corresponding cow culling reason tended to decrease. Our study showed that it seems reasonable to consider Holstein-Friesian cows aged at least 9 years at culling to be long-living animals. This is primarily evidenced by the rapid increase in the culling due to old age in relation to younger cows. Nowadays the above age limit can be suggested as a criterion of longevity for Holstein-Friesian cows but the criterion should be updated to the relation genotype-environment-economy that tends to change over time.     

Genome-wide association study of milk production traits in a crossbred dairy sheep population using three statistical models

Milk production is one of the most important characteristics of dairy sheep, and the identification of genes affecting milk production traits is critical to understanding the genetics and improve milk production in future generations. Three statistical techniques, namely GWAS, ridge-regression BLUP and BayesC , were used to identify SNPs in significant association with three milk production traits (milk yield, fat yield and protein yield) in a crossbred dairy sheep population. The results suggested that chromosomes 1, 3, 4, 5, 7 and 11 were likely to harbor genes important to milk production because these chromosomes had the greatest top-100-SNP variance contributions on the three milk production traits. The GWAS analysis identified between 74 and 288 genome-wide significant SNP (P < 0.05) whereas the BayesCπ model revealed between six and 63 SNPs, each with >95% posterior probability of inclusion as having a non-zero association effect on at least one of the three milk production traits. Positional candidate genes for milk production in sheep were searched, based on the sheep genomic assembly OAR version 3.1, such as those which map position coincided with or was located within 0.1 Mbp of a genome-wide suggestive or significant SNP. These identified SNPs and candidate genes supported some previous findings and also added new information about genetic markers for genetic improvement of lactation in dairy sheep, but keeping in mind that the majority of these positional candidate genes are not necessarily true causative loci for these traits and future validations are thus necessary.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Genome Wide Association Study Identifies 20 Novel Promising Genes Associated with Milk Fatty Acid Traits in Chinese Holstein

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 Chinese Holstein cows with the PLINK software. Genotypes were obtained with the Illumina BovineSNP50 Bead chip and a total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant associated with multiple SNPs, especially C18:0 (105 SNPs), C18 index (93 SNPs), and C14 index (84 SNPs); Several SNPs are close to or within the DGAT1, SCD1 and FASN genes which are well-known to affect milk composition traits of dairy cattle. Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising candidates for C10:0, C12:0, C14:0, C14:1, C14 index, C18:0, C18:1n9c, C18 index, SFA, UFA and SFA/UFA were found, which composed of HTR1B, CPM, PRKG1, MINPP1, LIPJ, LIPK, EHHADH, MOGAT1, ECHS1, STAT1, SORBS1, NFKB2, AGPAT3, CHUK, OSBPL8, PRLR, IGF1R, ACSL3, GHR and OXCT1. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting milk fatty acid traits in dairy cattle.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1595-0) contains supplementary material, which is available to authorized users.     

Associations between strain,herd size,age at first calving,culling reason and lifetime performance characteristics in Holstein-Friesian cows

Cow longevity and lifetime performance traits are good indicators of breeding effectiveness and animal welfare. They are also interrelated with the economics of dairy herd. Unfortunately, a high milk yield is often associated with deteriorated cow health and fertility and, consequently, with an increased culling rate. This situation, observed also in the Polish population of Holstein-Friesian cattle, inspired us to undertake a study on the associations between some factors and lifetime performance characteristics. The data set consisted of the records on 135 496 cows, including 131 526 of the Black and White strain (BW), and 3970 of the Red and White strain (RW) covered by performance recording and culled in 2012. It was found that cows of the BW strain and those from the largest herds (>100 cows) reached higher lifetime and mean daily energy-corrected milk (ECM) yields than cows of the RW strain and those from smaller herds culled at a similar age. Cows youngest at first calving (<2.0 years) were characterised by the highest lifetime ECM yield. It indicates that heifers can be bred even when they are younger than 15 to 16 months with no significant negative effect on their later performance. Infertility and reproduction problems (39.6%) and udder diseases (15.5%) constituted the most frequent reasons for cow culling. Cow longevity and lifetime productivity were considerably affected by the interactions between the studied factors.     

Analyzes of genome-wide association follow-up study for calving traits in dairy cattle

ABSTRACT: BACKGROUND: There is often a pronounced disagreement between results obtained from different genome-wide association studies in cattle. There are multiple reasons for this disagreement. Particularly the presence of false positives leads to a need to validate detected QTL before they are optimally incorporated or weighted in selection decisions or further studied for causal gene. In dairy cattle progeny testing scheme new data is routinely accumulated which can be used to validate previously discovered associations. However, the data is not an independent sample and the sample size may not be sufficient to have enough power to validate previous discoveries. Here we compared two strategies to validate previously detected QTL when new data is added from the same study population. We compare analyzing a combined dataset (COMB) including all data presently available to only analyzing a validation dataset (VAL) i.e. a new dataset not previously analyzed as an independent replication. Secondly, we confirm SNP detected in the Reference population (REF) (i.e. previously analyzed dataset consists of older bulls) in the VAL dataset. RESULTS: Clearly the results from the combined (COMB) dataset which had nearly twice the sample size of other two subsets allowed the detection of far more significant associations than the two smaller subsets. The number of significant SNPs in REF (older bulls) was about four times higher compare to VAL (younger bulls) though both had similar sample sizes, 2,219 and 2,039 respectively. A total of 424 SNP-trait combinations on 22 chromosomes showed genome-wide significant association involving 284 unique SNPs in the COMB dataset. In the REF data set 101 associations (73 unique SNPs) and in the VAL 24 associations (18 unique SNPs) were found genome-wide significant.Sixty-eight percent of the SNPs in the REF dataset could be confirmed in the VAL dataset. Out of 469 unique SNPs showing chromosome-wide significant association with calving traits in the REF dataset 321 could be confirmed in the VAL dataset at P < 0.05. CONCLUSIONS: The follow-up study for GWAS in cattle will depend on the aim of the study. If the aim is to discover novel QTL, analyses of the COMB dataset is recommended, while in case of identification of the causal mutation underlying a QTL, confirmation of the discovered SNPs are necessary to avoid following a false positive.     

Associations between milk performance traits in Holstein cows and 16 candidate SNPs identified by arrayed primer extension (APEX) microarray

An oligonucleotide microarray-which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis-was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples. The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.     

The role of bovine causal genes underlying dairy traits in Spanish Churra sheep

In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.     

Associations Between Milk Performance Traits in Holstein Cows and 16 Candidate SNPs Identified by Arrayed Primer Extension (APEX) Microarray

An oligonucleotide microarray—which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis—was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples.

The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.     

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6–14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.