Genome-Wide Association Study Dissects the Genetic Architecture of Seed Weight and Seed Quality in Rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium^® SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using ‘pseudomolecules’ representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.     

Genetic variation of BnaA3.NIP5;1 expressing in the lateral root cap contributes to boron deficiency tolerance in Brassica napus

Boron (B) is essential for vascular plants. Rapeseed (Brassica napus) is the second leading crop source for vegetable oil worldwide, but its production is critically dependent on B supplies. BnaA3.NIP5;1 was identified as a B-efficient candidate gene in B. napus in our previous QTL fine mapping. However, the molecular mechanism through which this gene improves low-B tolerance remains elusive. Here, we report genetic variation in BnaA3.NIP5;1 gene, which encodes a boric acid channel, is a key determinant of low-B tolerance in B. napus. Transgenic lines with increased BnaA3.NIP5;1 expression exhibited improved low-B tolerance in both the seedling and maturity stages. BnaA3.NIP5;1 is preferentially polar-localized in the distal plasma membrane of lateral root cap (LRC) cells and transports B into the root tips to promote root growth under B-deficiency conditions. Further analysis revealed that a CTTTC tandem repeat in the 5’UTR of BnaA3.NIP5;1 altered the expression level of the gene, which is tightly associated with plant growth and seed yield. Field tests with natural populations and near-isogenic lines (NILs) confirmed that the varieties carried BnaA3.NIP5;1^Q allele significantly improved seed yield. Taken together, our results provide novel insights into the low-B tolerance of B. napus, and the elite allele of BnaA3.NIP5;1 could serve as a direct target for breeding low-B-tolerant cultivars.     

Genome-wide association study for candidate genes controlling seed yield and its components in rapeseed (Brassica napus subsp. napus)

Genetic improvement of seed yield per plant (SY) is one of the major objectives in Brassica napus breeding programme. SY, being a complex quantitative trait is directly and indirectly influenced by yield-component traits such as siliqua length (SL), number of seeds per siliqua (NSS), and thousand seed weight (TSW). Therefore, concurrent improvement in SL, NSS and TSW can lead to higher SY in B. napus. This study was conducted to identify significant SNPs and putative candidate genes governing SY and its component traits (SL, NSS, TSW). All these traits were evaluated in a diverse set of 200 genotypes representing diversity from wide geographical locations. Of these, a set of 125 genotypes were chosen based on pedigree diversity and multi-location trait variation for genotyping by sequencing (GBS). Best linear unbiased predictors (BLUPs) of all the traits were used for genome-wide association study (GWAS) with 85,126 SNPs obtained from GBS. A total of 16, 18, 27 and 18 SNPs were found to be significantly associated for SL, NSS, TSW and SY respectively. Based on linkage disequilibrium decay analysis, 150 kb genomic region flanking the SNP was used for the identification of underlying candidate genes for each test trait. Important candidate genes involved in phytohormone signaling (WAT1, OSR1, ARR8, CKX1, REM7, REM9, BG1) and seed storage proteins (Cruciferin) were found to have significant influence on seed weight and yield. Genes involved in sexual reproduction and fertilization (PERK7, PERK13, PRK3, GATA15, NFD6) were found to determine the number of seeds per siliqua. Several genes found in this study namely ATS3A, CKX1, SPL2, SPL6, SPL9, WAT1 showed pleiotropic effect with yield component traits. Significant SNPs and putative candidate genes identified for SL, NSS, TSW and SY could be used in marker-assisted breeding for improvement of crop yield in B. napus. Genotypes identified with high SL, NSS, TSW and SY could serve as donors in crop improvement programs in B. napus.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-01060-9.     

Genome-wide association study dissects the genetic control of plant height and branch number in response to low-phosphorus stress in Brassica napus

Background and AimsOilseed rape (Brassica napus) is one of the most important oil crops worldwide. Phosphorus (P) deficiency severely decreases the plant height and branch number of B. napus. However, the genetic bases controlling plant height and branch number in B. napus under P deficiency remain largely unknown. This study aims to mine candidate genes for plant height and branch number by genome-wide association study (GWAS) and determine low-P-tolerance haplotypes.MethodsAn association panel of B. napus was grown in the field with a low P supply (P, 0 kg ha⁻¹) and a sufficient P supply (P, 40 kg ha⁻¹) across 2 years and plant height and branch number were investigated. More than five million single-nucleotide polymorphisms (SNPs) were used to conduct GWAS of plant height and branch number at two contrasting P supplies.Key ResultsA total of 2127 SNPs were strongly associated (P < 6·25 × 10⁻⁰⁷) with plant height and branch number at two P supplies. There was significant correlation between phenotypic variation and the number of favourable alleles of associated loci on chromosomes A10 (chrA10_821671) and C08 (chrC08_27999846), which will contribute to breeding improvement by aggregating these SNPs. BnaA10g09290D and BnaC08g26640D were identified to be associated with chrA10_821671 and chrC08_27999846, respectively. Candidate gene association analysis and haplotype analysis showed that the inbred lines carrying ATT at BnaA10g09290Hap1 and AAT at BnaC08g26640Hap1 had greater plant height than lines carrying other haplotype alleles at low P supply.ConclusionOur results demonstrate the power of GWAS in identifying genes of interest in B. napus and provided insights into the genetic basis of plant height and branch number at low P supply in B. napus. Candidate genes and favourable haplotypes may facilitate marker-based breeding efforts aimed at improving P use efficiency in B. napus.     

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.     

Lack of association between genetic polymorphisms in ROBO1, MRPL19/C2ORF3 and THEM2 with Developmental Dyslexia

Developmental Dyslexia (DD) is a heritable, complex genetic disorder characterized by specific impairment in reading and writing ability that is substantially below the expected reading ability given the person's chronological age, measured intelligence and age-appropriate education. More than ten susceptible genes have been identified for DD. A Single Nucleotide Polymorphism (SNP) of these genes was found to be associated with various phenotypes of DD. To identify the role of SNPs of four candidate genes namely, MRPL19/C2ORF3, ROBO1 and THEM2 in an Indian population, we genotyped eight SNPs of these genes in 157 children with DD and 212 normal readers using a MassARRAY technique with a MALDI-TOF MS analyzer. Power analysis of some of these SNPs showed > 80% of power. Chi-square test, Odds Ratios (ORs), 95% Confidence Intervals (CIs) and Bonferroni's correction were applied to identify the significance of the genotyped SNPs and haplotypes. Our study failed to show any association of SNPs and haplotypes of these genes with DD in an Indian population.     

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.     

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome‐wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue‐specific associated loci were related to total leaf and seed GSL content. Thirty‐six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double‐low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.     

植物种子植酸研究进展

磷是植物生长发育所必需的大量营养元素。在种子发育过程中,植酸是磷的贮存库,对维持植物体内磷平衡有重要的作用。在种子萌发过程中,植酸酶分解植酸盐,释放磷、矿质营养和肌醇供幼苗生长。本文综述了近年来植物(作物)种子中植酸的生物合成途径、种子植酸含量的遗传、低植酸作物的育种等研究进展。首先,植酸生物合成途径中最初的反应底物为葡萄糖-6-磷酸,形成肌醇后,以肌醇为底物合成植酸共有两条路径:依赖脂类与不依赖脂类,目前,已分离鉴定若干植酸合成所需的关键酶及其编码基因,包括肌醇-3-磷酸合成酶、肌醇激酶、肌醇多磷酸盐激酶,以及参与植酸运输的ATP结合盒转运子。其次,利用作图群体及关联分析群体,分别在水稻(Oryza sativa L.)、白菜(Brassica rapa L.)、菜豆(Phaseolus vulgaris L.)等植物中鉴定出多个与种子植酸磷含量相关的遗传位点。第三,筛选获得有价值的低植酸突变体是培育低植酸作物的主要途径。当把低植酸作为育种目标时,可能会忽略种子植酸含量的降低给植物带来的不利影响,如何消除低植酸造成的不利影响,成为科学家们亟需解决的问题。     

Imprints of Natural Selection Along Environmental Gradients in Phenology-Related Genes of Quercus petraea

We explored single nucleotide polymorphism (SNP) variation in candidate genes for bud burst from Quercus petraea populations sampled along gradients of latitude and altitude in Western Europe. SNP diversity was monitored for 106 candidate genes, in 758 individuals from 32 natural populations. We investigated whether SNP variation reflected the clinal pattern of bud burst observed in common garden experiments. We used different methods to detect imprints of natural selection (F_ST outlier, clinal variation at allelic frequencies, association tests) and compared the results obtained for the two gradients. F_ST outlier SNPs were found in 15 genes, 5 of which were common to both gradients. The type of selection differed between the two gradients (directional or balancing) for 3 of these 5. Clinal variations were observed for six SNPs, and one cline was conserved across both gradients. Association tests between the phenotypic or breeding values of trees and SNP genotypes identified 14 significant associations, involving 12 genes. The results of outlier detection on the basis of population differentiation or clinal variation were not very consistent with the results of association tests. The discrepancies between these approaches may reflect the different hierarchical levels of selection considered (inter- and intrapopulation selection). Finally, we obtained evidence for convergent selection (similar for gradients) and clinal variation for a few genes, suggesting that comparisons between parallel gradients could be used to screen for major candidate genes responding to natural selection in trees.     

<Emphasis Type="Italic">LMI1</Emphasis>-like genes involved in leaf margin development of <Emphasis Type="Italic">Brassica napus</Emphasis>

In rapeseed (Brassica napus L.), leaf margins are variable and can be entire, serrate, or lobed. In our previous study, the lobed-leaf gene (LOBED-LEAF 1, BnLL1) was mapped to a 32.1 kb section of B. napus A10. Two LMI1-like genes, BnaA10g26320D and BnaA10g26330D, were considered the potential genes that controlled the lobed-leaf trait in rapeseed. In the present study, these two genes and another homologous gene (BnaC04g00850D) were transformed into Arabidopsis thaliana (L.) Heynh. plants to identify their functions. All three LMI1-like genes of B. napus produced serrate leaf margins. The expression analysis indicated that the expression level of BnaA10g26320D determined the difference between lobed- and entire-leaved lines in rapeseed. Therefore, it is likely that BnaA10g26320D corresponds to BnLL1.     

Sequence diversity, haplotype analysis, association mapping and functional marker development in the waxy and starch synthase IIa genes for grain-yield-related traits in hexaploid wheat (Triticum aestivum L.)

Development of high-yielding cereal crops could meet increasing global demands for food, feed and bio-fuels. Wheat is one of the world??s most important cereal crops. The biosynthesis of starch is the major determinant of yield in wheat. Two starch biosynthesis genes, the waxy (Wx) genes and the starch synthase IIa (SSIIa) genes, were amplified and sequenced in 92 diverse wheat genotypes using genome-specific primers. Nucleotide diversity, haplotype analysis and association mapping were performed. The first exon (5??-UTR) and the first intron of the three homoeologous Wx genes were isolated using expressed sequence tag sequences. The Wx genes contained 12 exons separated by 11 introns. SNP (single nucleotide polymorphism) frequency ranged from 1 SNP/3,648?bp for Wx-D1 to 1 SNP/135?bp for SSIIa-A1, with an average of 1 SNP/230?bp. The average SNP frequencies in exon and intron regions were 1 SNP/322?bp and 1 SNP/228?bp, respectively. Thirty, 23 and 5 SNPs were identified and formed five, six and five haplotypes for SSIIa-A1, SSIIa-B1 and SSIIa-D1, respectively. However, no association was found between these SNPs and seven yield-related traits. Twenty-two, 15 and 1 SNPs were detected and formed nine, five and two haplotypes for Wx-A1, Wx-B1 and Wx-D1, respectively. Three unique nucleotides C+A+T at SNP5, SNP6 and SNP12 formed Wx-B1-H3, which was significantly associated with increased grain weight, thousand kernel weight, and total starch content in three spring wheat genotypes and five winter wheat genotypes. Cost-effective and co-dominant SNP markers were developed using temperature-switch (TS)-PCR and are being used for marker-assisted selection of doubled haploid lines with enhanced grain yield and starch content in winter wheat breeding programs.     

Effect of phytase supplementation on phosphorus digestibility in low-phytate barley fed to finishing pigs

Forty crossbred barrows (Camborough 15 Line female×Canabred sire) weighing an average of 79.6±8.0?kg were used in a factorial design experiment (5 barleys×2 enzyme levels) conducted to determine the effects of phytase supplementation on nutrient digestibility in low-phytate barleys fed to finishing pigs. The pigs were assigned to one of 10 dietary treatments comprised of a normal 2-rowed, hulled variety of barley (CDC Fleet, 0.26% phytate) or 2 low-phytate hulled genotypes designated as LP422 (0.14% phytate) and LP635 (0.09% phytate). A normal, hulless barley (CDC Dawn, 0.26% phytate) and a hulless genotype designated as LP422H (0.14% phytate) were also included. All barleys were fed with and without phytase (Natuphos 5000 FTU/kg). The diets fed contained 98% barley, 0.5% vitamin premix, 0.5% trace mineral premix, 0.5% NaCl and 0.5% chromic oxide but no supplemental phosphorus. The marked feed was provided for a 7-day acclimatization period, followed by a 3-day faecal collection. In the absence of phytase, phosphorus digestibility increased substantially (P<0.05) as the level of phytate in the barley declined. For the hulled varieties, phosphorus digestibility increased from 12.9% for the normal barley (0.26% phytate) to 35.3 and 39.8% for the two low-phytate genotypes (0.14 and 0.09% phytate respectively). For the hulless varieties, phosphorus digestibility increased from 9.2% for the normal barley (0.26% phytate) to 34.7% for the hulless variety with 54% of the normal level of phytate (0.14% phytate). In contrast, when phytase was added to the diet, there was little difference in phosphorus digestibility between pigs fed normal barley and those fed the low-phytate genotypes (significant barley×enzyme interaction, P=0.01). For the hulled varieties, phosphorus digestibility was 50.1% for the barley with the normal level of phytate (0.26% phytate) compared with 51.1 and 52.4% for the varieties with 54 and 35% of the normal level of phytate (0.14 and 0.09% phytate respectively). For the hulless varieties, phosphorus digestibility increased from 47.1% for the normal barley (0.26% phytate) to 54.4% for the hulless variety with 54% of the normal level of phytate (0.14% phytate). In conclusion, both supplementation with phytase and selection for low-phytate genotypes of barley were successful in increasing the digestibility of phosphorus for pigs. Unfortunately, the effects did not appear to be additive. Whether or not swine producers will choose low-phytate barley or supplementation with phytase as a means to improve phosphorus utilization, will likely depend on the yield potential of low-phytate barley and the additional costs associated with supplementation with phytase.     

A genome-wide association study uncovers novel genomic regions and candidate genes of yield-related traits in upland cotton

Key message

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield.

Abstract

Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

     

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

Haplotype distribution and association of candidate genes with salt tolerance in Indian wild rice germplasm

Key message

The association of natural genetic variations of salt-responsive candidate genes belonging to different gene families with salt-tolerance phenotype and their haplotype variation in different geographic regions.

Abstract

Soil salinity covers a large part of the arable land of the world and is a major factor for yield losses in salt-sensitive crops, such as rice. Different gene families that respond to salinity have been identified in rice, but limited success has been achieved in developing salt-tolerant cultivars. Therefore, 21 salt stress-responsive candidate genes belonging to different gene families were re-sequenced to analyse their genetic variation and association with salt tolerance. The average single nucleotide polymorphism (SNP) density was 16 SNPs per kbp amongst these genes. The identified nucleotide and haplotype diversity showed comparatively higher genetic variation in the transporter family genes. Linkage disequilibrium (LD) analysis showed significant associations of SNPs in BADH2, HsfC1B, MIPS1, MIPS2, MYB2, NHX1, NHX2, NHX3, P5CS1, P5CS2, PIP1, SIK1, SOS1, and SOS2 genes with the salt-tolerant phenotype. A combined analysis of SNPs in the 21 candidate genes and eight other HKT transporter genes produced two separate clusters of tolerant genotypes, carrying unique SNPs in the ion transporter and osmoticum-related genes. Haplotype network analysis showed all the major and few minor alleles distributed over distant geographic regions. Minor haplotypes may be recently evolved alleles which migrated to distant geographic regions and may represent recent expansion of Indian wild rice. The analysis of genetic variation in different gene families identified the relationship between adaptive variations and functional significance of the genes. Introgression of the identified alleles from wild relatives may enhance the salt tolerance and consequently rice production in the salinity-affected areas.

     

SNPs associated with body weight and backfat thickness in two pig breeds identified by a genome-wide association study

Growth and fat deposition are important economic traits due to the influence on production in pigs. In this study, a dataset of 1200 pigs with 345,570 SNPs genotyped by sequencing (GBS) was used to conduct a GWAS with single-marker regression method to identify SNPs associated with body weight and backfat thickness (BFT) and to search for candidate genes in Landrace and Yorkshire pigs. A total of 27 and 13 significant SNPs were associated with body weight and BFT, respectively. In the region of 149.85–149.89 Mb on SSC6, the SNP (SSC6: 149876737) for body weight and the SNP (SSC6: 149876507) for BFT were in the same locus region (a gap of 230 bp). Two SNPs were located in the DOCK7 gene, which is a protein-coding gene that plays an important role in pigmentation. Two SNPs located on SSC8: 54567459 and SSC11: 33043081 were found to overlap weight and BFT; however, no candidate gene was found in these regions. In addition, based on other significant SNPs, two positional candidate genes, NSRP1 and CADPS, were proposed to influence weight. In conclusion, this is the first study report using GBS data to identify the significant SNPs for weight and BFT. A total of four particularly interesting SNPs and one potential candidate genes (DOCK7) were found for these traits in domestic pigs. This study improves our knowledge to better understand the complex genetic architecture of weight and BFT, but further validation studies of these candidate loci and genes are recommended in pigs