A genome‐wide association study for a proxy of intermuscular fat level in the Italian Large White breed identifies genomic regions affecting an important quality parameter for dry‐cured hams

Intermuscular fat content in protected designations of origin dry‐cured hams is a very important meat quality trait that affects the acceptability of the product by the consumers. An excess in intermuscular fat (defined as the level of fat deposition between leg muscles) is a defect that depreciates the final product. In this study we carried out a genome‐wide association study for visible intermuscular fat (VIF) of hams in the Italian Large White pig breed. This trait was evaluated on the exposed muscles of green legs in 1122 performance‐tested gilts by trained personnel, according to a classification scale useful for routine and cheap evaluation. All animals were genotyped with the Illumina PorcineSNP60 BeadChip. The genome‐wide association study identified three QTL regions on porcine chromosome 1 (SSC1; accounting for ~79% of the SNPs below the 5.0E?04 threshold) and SSC2, two on SSC7 and one each on SSC3, SSC6, SSC9, SSC11, SSC13, SSC15, SSC16 and SSC17. The most significant SNP (ALGA0004143 on SSC1 at 77.3 Mb; P_FDR < 0.05), included in the largest QTL region which spanned about 6.8 Mb on SSC1, is located within the glutamate ionotropic receptor kainate type subunit 2 (GRIK2) gene. Functional annotation of all genes included in QTL regions for VIF suggested that intermuscular fat in the Italian Large White breed is a complex trait apparently influenced by complex biological mechanisms also involving obesity‐related processes. These QTL target mainly chromosome regions different from those affecting subcutaneous and intramuscular fat deposition.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Genome‐wide association study for jumping performances in French sport horses

A genome‐wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo‐Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54 602 SNPs available on this chip, 44 424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single‐marker mixed model and a haplotype‐based mixed model (significant: P < 1E‐05; suggestive: P < 1E‐04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.     

Analysis of non‐synonymous SNPs of the porcine SERPINA6 gene as potential causal variants for a QTL affecting plasma cortisol levels on SSC7

Recently, the SERPINA6 gene encoding corticosteroid‐binding globulin (CBG) has been proposed as a candidate gene for a quantitative trait locus (QTL) affecting cortisol level on pig chromosome 7. The QTL was repeatedly detected in different lines, including a Piétrain × (German Landrace × German Large White) cross (PiF1) and purebred German Landrace (LR). In this study, we investigated whether the known non‐synonymous polymorphisms c.44G>T, c.622C>T, c.770C>T, c.793G>A, c.832G>A and c.919G>A of SERPINA6 are sufficient to explain the QTL in these two populations. Our investigations revealed that SNPs c.44G>T, c.622C>T, c.793G>A and c.919G>A are associated with cortisol level in PiF1 (P < 0.01). Haplotype analysis showed that these associations are largely attributable to differences between a major haplotype carrying SNPs c.793G>A and c.919G>A and a haplotype carrying SNPs c.44G>T and c.622C>T. Furthermore, some SNPs, particularly c.44G>T and c.622C>T and the carrier haplotype, showed association with meat quality traits including pH and conductivity (P < 0.05). In LR, the non‐synonymous SNPs segregate at very low frequency (<5%) and/or show only weak association with cortisol level (SNPs c.832G>A and c.919G>A; P < 0.05). These findings suggest that the non‐synonymous SNPs are not sufficient to explain the QTL across different breeds. Therefore, we examined whether the expression of SERPINA6 is affected by cis‐regulatory polymorphisms in liver, the major organ for CBG production. We found allelic expression imbalance of SERPINA6, which suggests that its expression is indeed affected by genetic variation in cis‐acting elements. This represents candidate causal variation for future studies of the molecular background of the QTL.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

Genomic analysis of the differential response to experimental infection with porcine circovirus 2b

Porcine circovirus type 2 (PCV2) is the etiological agent of a group of associated diseases (PCVAD) that affect production efficiency and can lead to mortality. Using different crossbred lines of pigs, we analyzed host genetic variation of viral load, immune response and weight change following experimental infection with a PCV2b strain (n = 386). Pigs expressed variation in the magnitude and initiation of viremia and immune response recorded weekly until 28 days post‐infection. A higher viral load was correlated with weight gain (r = ?0.26, P < 0.0001) and presence of PCV2‐specific antibodies (IgM, r = 0.26–0.34, P < 0.0001; IgG, r = 0.17–0.20, P < 0.01). In genome‐wide association analyses of the responses at different time points, the proportions of phenotypic variation explained by combined effects of 56 433 SNPs were 34.8–59.4% for viremia, 10.1–59.5% for antibody response and 5.6–14.9% for weight change. Relationships between genomic prediction of overall viral load and weight gain during the first weeks of challenge were negative (?0.21 and ?0.24 respectively, P < 0.0001). Individuals that carried more favorable alleles across three SNPs on SSC9 (0.60 Mb) and SSC12 (6.8 and 18.2 Mb) partially explained this relationship, having lower viral load (P < 0.0001); lower viremia at day 14 (P < 0.0001), day 21 (P < 0.01) and day 28 (P < 0.05) and greater overall average daily gain during infection (ADG_i; P < 0.01), ADG_i at week 3 (P < 0.001) and week 4 (P < 0.01). These additive genetic relationships could lead to molecular solutions to improve animal health and reduce production costs.     

Exploiting whole genome sequence data to fine map and characterize candidate genes within a quantitative trait loci region affecting androstenone on porcine chromosome 5

Male piglets are routinely castrated to eliminate boar taint. However, this treatment is undesirable, and alternative approaches, including genetic strategies to reduce boar taint, are demanded. Androstenone is one of the causative agents of boar taint, and a QTL region affecting this pheromone has previously been reported on SSC5: 22.6–24.8 Mb in Duroc. The QTL region is one of the few reported for androstenone that does not simultaneously affect levels of other sex steroids. The main objective of this study was to fine map this QTL. Whole genome sequence data from 23 Norwegian Duroc boars were analyzed to detect new polymorphisms within the QTL region. A subset of 161 SNPs was genotyped in 834 Duroc sires and analyzed for association with androstenone in adipose tissue and testosterone, estrone sulphate and 17β‐estradiol in blood plasma. Our results revealed 100 SNPs significantly associated with androstenone levels in fat (P < 0.001) with 94 of the SNPs being in strong linkage disequilibrium in the region 23.03–24.27 Mb. This haplotype block contains at least four positional candidate genes (HSD17B6, SDR9C7, RDH16 and STAT6) involved in androstenone biosynthesis. No significant associations were found between any of the SNPs and levels of testosterone and estrogens, confirming previous findings. The amount of phenotypic variance explained by single SNPs within the haplotype block was as high as 5.4%. As the SNPs in this region significantly affect levels of androstenone without affecting levels of other sex steroids, they are especially interesting as genetic markers for selection against boar taint.     

Genome‐wide prediction of age at puberty and reproductive longevity in sows

Traditional selection for sow reproductive longevity is ineffective due to low heritability and late expression of the trait. Incorporation of DNA markers into selection programs is potentially a more practical approach for improving sow lifetime productivity. Using a resource population of crossbred gilts, we explored pleiotropic sources of variation that influence age at puberty and reproductive longevity. Of the traits recorded before breeding, only age at puberty significantly affected the probability that females would produce a first parity litter. The genetic variance explained by 1‐Mb windows of the sow genome, compared across traits, uncovered regions that influence both age at puberty and lifetime number of parities. Allelic variants of SNPs located on SSC5 (27–28 Mb), SSC8 (36–37 Mb) and SSC12 (1.2–2 Mb) exhibited additive effects and were associated with both early expression of puberty and a greater than average number of lifetime parities. Combined analysis of these SNPs showed that an increase in the number of favorable alleles had positive impact on reproductive longevity, increasing number of parities by up to 1.36. The region located on SSC5 harbors non‐synonymous alleles in the arginine vasopressin receptor 1A (AVPR1A) gene, a G‐protein‐coupled receptor associated with social and reproductive behaviors in voles and humans and a candidate for the observed effects. This region is characterized by high levels of linkage disequilibrium in different lines and could be exploited in marker‐assisted selection programs across populations to increase sow reproductive longevity.     

Single nucleotide polymorphisms in the insulin‐like growth factor 1 (IGF1) gene are associated with growth‐related traits in farmed Atlantic salmon

Understanding the genetic basis of variation in traits related to growth and fillet quality in Atlantic salmon is of importance to the aquaculture industry. Several growth‐related QTL have been identified via the application of genetic markers. The IGF1 gene is considered a highly conserved and crucial growth‐regulating gene in salmonid species. However, the association between polymorphisms in the IGF1 gene and growth‐related traits in Atlantic salmon is unknown. Therefore, in this study, regions of the Atlantic salmon IGF1 gene were sequenced, aligned and compared across individuals. Three SNPs were identified in the putative promoter (SNP1, g.5763G>T; GenBank no. AGKD01012745 ), intron 1 (SNP2, g.7292C>T; GenBank no. AGKD01012745 ) and intron 3 (SNP3, g.4671A>C; GenBank no. AGKD01133398 ) regions respectively. These SNPs were genotyped in a population of 4800 commercial Atlantic salmon with data on several weight and fillet traits measured at harvest (at approximately 3 years of age). In a mixed model, association analysis of individual SNPs, SNP1 and SNP3 were both significantly associated with several weight traits (P < 0.05). The estimated additive effect on overall harvest weight was approximately 35 and 110 g for SNPs 1 and 3 respectively. A haplotype analysis confirmed the association between genetic variation in the IGF1 gene with overall body weight (P < 0.05) and fillet component traits (P < 0.05). Our findings suggest the identified nucleotide polymorphisms of the IGF1 gene may either affect farmed Atlantic salmon growth directly or be in population‐wide linkage disequilibrium with causal variation, highlighting their possible utility as candidates for marker‐assisted selection in the aquaculture industry.     

Genomic regions affecting backfat thickness and cannon bone circumference identified by genome‐wide association study in a Duroc pig population

We performed a genome‐wide association study using the porcine 60K SNP array to detect QTL regions for nine traits in a three‐generational Duroc samples (n = 651), viz. generations 1, 2 and 3 from a population selected over five generations using a closed nucleus breeding scheme. We applied a linear mixed model for association mapping to detect SNP effects, adjusting for fixed effects (sex and season) and random polygenic effects (reflecting genetic relatedness), and derived a likelihood ratio statistic for each SNP using the efficient mixed‐model association method. We detected a region on SSC6 for backfat thickness (BFT) and on SSC7 for cannon bone circumference (CANNON), with a genome‐wide significance of P < 0.01 after Bonferroni correction. These regions had been detected previously in other pig populations. Six genes are located in the BFT‐associated region, while the CANNON‐associated region includes 66 genes. In the future, significantly associated SNPs, derived by sequencing the coding regions of the six genes in the BFT region, can be used in marker‐assisted selection of BFT, whereas haplotypes constructed from the SSC7 region with strong LD can be used to select for the CANNON trait in our resource family.     

Three novel quantitative trait loci for skin thickness in swine identified by linkage and genome‐wide association studies

Skin is the largest organ in the pig body and plays a key role in protecting the body against pathogens and excessive water loss. Deciphering the genetic basis of swine skin thickness would enrich our knowledge about the skin. To identify the loci for porcine skin thickness, we first performed a genome scan with 194 microsatellite markers in a White Duroc × Erhualian F₂ intercross. We identified three genome‐wide significant QTL on pig chromosomes (SSC) 4, 7 and 15 using linkage analysis. The most significant QTL was found on SSC7 with a small confidence interval of ~5 cM, explaining 23.9 percent of phenotypic variance. Further, we conducted a genome‐wide association study (GWAS) using Illumina PorcineSNP60 Beadchips for the F₂ pedigree and a population of Chinese Sutai pigs. We confirmed significant QTL in the F₂ pedigree and replicated QTL on SSC15 in Chinese Sutai pigs. A meta‐analysis of GWASs on both populations detected a genomic region associated with skin thickness on SSC4. GWAS results were generally consistent with QTL mapping. Identical‐by‐descent analysis defined QTL on SSC7 in a 683‐kb region harboring an interesting candidate gene: HMGA1. On SSC15, the linkage disequilibrium analysis showed a haplotype block of 2.20 Mb that likely harbors the gene responsible for skin thickness. Our findings provide novel insights into the genetic basis of swine skin thickness, which would benefit further understanding of porcine skin function.     

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.     

Genome‐wide association study reveals regions associated with gestation length in two pig populations

Reproduction traits, such as gestation length (GLE), play an important role in dam line breeding in pigs. The objective of our study was to identify single nucleotide polymorphisms (SNPs) that are associated with GLE in two pig populations. Genotypes and deregressed breeding values were available for 2081 Dutch Landrace‐based (DL) and 2301 Large White‐based (LW) pigs. We identified two QTL regions for GLE, one in each population. For DL, three associated SNPs were detected in one QTL region spanning 0.52 Mbp on Sus scrofa chromosome (SSC) 2. For LW, four associated SNPs were detected in one region of 0.14 Mbp on SSC5. The region on SSC2 contains the heparin‐binding EGF‐like growth factor (HBEGF) gene, which promotes embryo implantation and has been described to be involved in embryo survival throughout gestation. The associated SNP can be used for marker‐assisted selection in the studied populations, and further studies of the HBEGF gene are warranted to investigate its role in GLE.     

A genome‐wide association study to detect genetic variation for postpartum dysgalactia syndrome in five commercial pig breeding lines

Postpartum dysgalactia syndrome (PDS) in sows is an important disease after parturition with a relevant economic impact, affecting the health and welfare of both sows and piglets. The genetic background of this disease has been discussed and its heritability estimated, but further genetic analyses are lacking in detail. The aim of the current study was to detect loci affecting the susceptibility to PDS through a genome‐wide association approach. The study was designed as a family‐based association study with matched sampling of affected sows and healthy half‐ or full‐sib control sows on six farms. For the study, 597 sows (322 affected vs. 275 healthy control sows) were genotyped on 62 163 single nucleotide polymorphisms (SNPs) using the Illumina PorcineSNP60 BeadChip. After quality control, 585 sows (314 affected vs. 271 healthy control sows) and 49 740 SNPs remained for further analysis. Statistics were performed mainly with the r package genabel and included a principal component analysis. A statistically significant genome‐wide associated SNP was identified on porcine chromosome (SSC) 17. Further promising results with moderate significance were detected on SSC 13 and on an unplaced scaffold with an older annotation on SSC 15. The PRICKLE2 and NRP2 genes were identified as candidate genes near associated SNPs. Several quantitative trait loci (QTL) have been previously described in these genomic regions, including QTL for mammary gland condition, as teat number and non‐functional nipples QTL, as well as QTL for body temperature and gestation length.     

Genome‐wide detection of CNVs associated with beak deformity in chickens using high‐density 600K SNP arrays

Beak deformity (crossed beaks) is found in several indigenous chicken breeds including Beijing‐You studied here. Birds with deformed beaks have reduced feed intake and poor production performance. Recently, copy number variation (CNV) has been examined in many species and is recognized as a source of genetic variation, especially for disease phenotypes. In this study, to unravel the genetic mechanisms underlying beak deformity, we performed genome‐wide CNV detection using Affymetrix chicken high‐density 600K data on 48 deformed‐beak and 48 normal birds using penncnv . As a result, two and eight CNV regions (CNVRs) covering 0.32 and 2.45 Mb respectively on autosomes were identified in deformed‐beak and normal birds respectively. Further RT‐qPCR studies validated nine of the 10 CNVRs. The ratios of six CNVRs were significantly different between deformed‐beak and normal birds (P < 0.01). Within these six regions, three and 21 known genes were identified in deformed‐beak and normal birds respectively. Bioinformatics analysis showed that these genes were enriched in six GO terms and one KEGG pathway. Five candidate genes in the CNVRs were further validated using RT‐qPCR. The expression of LRIG2 (leucine rich repeats and immunoglobulin like domains 2) was lower in birds with deformed beaks (P < 0.01). Therefore, the LRIG2 gene could be considered a key factor in view of its known functions and its potential roles in beak deformity. Overall, our results will be helpful for future investigations of the genomic structural variations underlying beak deformity in chickens.     

Genome‐wide association study for the level of serum electrolytes in Italian Large White pigs

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome‐wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance‐tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome‐wide efficient mixed‐model association. For the level of Ca²⁺, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (P_{nominal value} < 5.00E‐05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.     

Whole‐genome resequencing‐based QTL‐seq identified candidate genes and molecular markers for fresh seed dormancy in groundnut

The subspecies fastigiata of cultivated groundnut lost fresh seed dormancy (FSD) during domestication and human‐made selection. Groundnut varieties lacking FSD experience precocious seed germination during harvest imposing severe losses. Development of easy‐to‐use genetic markers enables early‐generation selection in different molecular breeding approaches. In this context, one recombinant inbred lines (RIL) population (ICGV 00350 × ICGV 97045) segregating for FSD was used for deploying QTL‐seq approach for identification of key genomic regions and candidate genes. Whole‐genome sequencing (WGS) data (87.93 Gbp) were generated and analysed for the dormant parent (ICGV 97045) and two DNA pools (dormant and nondormant). After analysis of resequenced data from the pooled samples with dormant parent (reference genome), we calculated delta‐SNP index and identified a total of 10,759 genomewide high‐confidence SNPs. Two candidate genomic regions spanning 2.4 Mb and 0.74 Mb on the B05 and A09 pseudomolecules, respectively, were identified controlling FSD. Two candidate genes—RING‐H2 finger protein and zeaxanthin epoxidase—were identified in these two regions, which significantly express during seed development and control abscisic acid (ABA) accumulation. QTL‐seq study presented here laid out development of a marker, GMFSD1, which was validated on a diverse panel and could be used in molecular breeding to improve dormancy in groundnut.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.     

A whole genome sequence association study for puberty in a large Duroc × Erhualian F2 population

A large proportion of gilts and sows are culled from reproduction populations because of anestrus and pubertal reproductive failure. Selecting early onset of puberty gilts has a favorable effect on sows’ reproductivity. However, age at puberty is hard to be routinely measured in commercial herds. With molecular genetic predictors, identifying individuals that have a propensity for early onset of puberty can be simplified. We previously performed genome scanning and a genome‐wide association study for puberty in an F2 resource population using 183 microsatellites and 62 125 SNPs respectively. The detection power and resolution of identified quantitative trait loci were very low. Herein, we re‐sequenced 19 founders of the F2 resource population in high coverage, and whole genome sequences of F2 individuals were imputed to perform an association study for reproductive traits. A total of 2339 SNPs associated with pubertal reproductive failure were identified in the region of 30.94–40.74 Mb on SSC7, with the top one, positioned at 33.36 Mb, explaining 16% of the phenotypic variances. We improved the magnitude of the P‐value by 10E+5 to 10E+7 using the whole genome sequence rather than using low/middle density markers as in previous studies, and we narrowed down the QTL confidence interval to 5.25 Mb. Combining the annotation of gene function, RAB23 and BAK1 were perceived as the most compelling candidate genes. The identified loci may be useful in culling sows failing to show estrus by marker‐assisted selection to increase reproductive efficiency of swine herds.