Improvement of protein efficiency (PE) is a key factor for a sustainable pig production, as nitrogen excretion contributes substantially to environmental pollution. Protein efficiency has been shown to be heritable and genetically correlated with performance traits such as feed conversion ratio (FCR) and average daily feed intake (ADFI). This study aimed to identify genomic regions associated with these traits through single-variant genome-wide association studies (GWAS) and regional heritability mapping (RHM) using whole-genome sequence variants from low-pass sequencing of more than 1000 Swiss Large White pigs. Genomic heritability estimates using 15 million variants were moderate to high, ranging from 0.33 to 0.47. GWAS did not identify significant variants for PE and FCR, but identified 45 variants at suggestive significance levels for ADFI on chromosome 1 and one for ADG on chromosome 14. Similarly, RHM detected no significant regions for PE and FCR, but five suggestive regions for ADFI (chromosome 1) and one for ADG (chromosome 14). However, by combining leading signals from GWAS and RHM, i.e. overlapping leading variants and significant regions, we highlighted putative candidate genes for PE, including PHYKPL, COL23A1, PPFIBP2, GVIN1, SYT9, RBMXL2, ZNF215, and olfactory receptor genes. Combining GWAS and RHM allowed us to identify genomic regions that may influence PE and production traits. Our apparent difficulty in detecting significant regions for these traits probably reflects the relatively small sample size, differences in genetic architecture across study designs and experimental conditions, and that polymorphisms explaining large proportions of the trait variation may not segregate in this population. Nevertheless, we identified plausible functional candidate genes in the highlighted regions, including those involved in nutrient sensing, the urea cycle, and metabolic pathways, in particular IGF1-insulin, and that have previously been reported to be associated with nitrogen metabolism in cattle and with muscle and adipose tissue metabolism and feed intake in pigs. We also highlighted a range of noncoding RNAs. Their targets and roles in gene regulation should be further investigated in this context.

BACKGROUND:Genetic correlations between complex traits suggest that pleiotropic variants contribute to trait variation. Genome-wide association studies (GWAS) aim to uncover the genetic underpinnings of traits. Multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS enable detecting variants associated with multiple phenotypes. In this study, we used array-derived genotypes and phenotypes for 24 reproduction, production, and conformation traits to explore differences between the two methods and used imputed sequence variant genotypes to fine-map six quantitative trait loci (QTL).RESULTS:We considered genotypes at 44,733 SNPs for 5,753 pigs from the Swiss Large White breed that had deregressed breeding values for 24 traits. Single-trait association analyses revealed eleven QTL that affected 15 traits. Multi-trait association testing and the meta-analysis of the single-trait GWAS revealed between 3 and 6 QTL, respectively, in three groups of traits. The multi-trait methods revealed three loci that were not detected in the single-trait GWAS. Four QTL that were identified in the single-trait GWAS, remained undetected in the multi-trait analyses. To pinpoint candidate causal variants for the QTL, we imputed the array-derived genotypes to the sequence level using a sequenced reference panel consisting of 421 pigs. This approach provided genotypes at 16 million imputed sequence variants with a mean accuracy of imputation of 0.94. The fine-mapping of six QTL with imputed sequence variant genotypes revealed four previously proposed causal mutations among the top variants.CONCLUSIONS:Our findings in a medium-size cohort of pigs suggest that multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS provide very similar results. Although multi-trait association methods provide a useful overview of pleiotropic loci segregating in mapping populations, the investigation of single-trait association studies is still advised, as multi-trait methods may miss QTL that are uncovered in single-trait GWAS.

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

We performed single-trait and multi-trait genome-wide association studies (GWAS) to detect quantitative trait loci (QTL) for body length and feed intake in 3,638 Swiss Large White pigs that had genotypes at 43,994 SNPs. Single-trait GWAS (ST-GWAS) revealed two non-overlapping QTL for each trait. Multivariate GWAS implemented with GEMMA software yielded 75 significant SNPs (-log10(p)>6) at the same four QTL. A meta-analysis of summary statistics from the two ST-GWAS resulted in 83 significant SNPs at the same QTL. After including all samples that had at least one phenotype recorded (in total 5,109), the P-values of top variants decreased suggesting higher statistical power due to a larger sample size. Fine mapping using imputed sequence variant genotypes confirmed the presence of three previously reported candidate causal variants in MC4R, VRTN and BMP2 genes among the most significantly associated variants.

Background The key-ancestor approach has been frequently applied to prioritize individuals for whole-genome sequencing based on their marginal genetic contribution to current populations. Using this approach, we selected 70 key ancestors from two lines of the Swiss Large White breed that have been selected divergently for fertility and fattening traits and sequenced their genomes with short paired-end reads. Results Using pedigree records, we estimated the effective population size of the dam and sire line to 72 and 44, respectively. In order to assess sequence variation in both lines, we sequenced the genomes of 70 boars at an average coverage of 16.69-fold. The boars explained 87.95 and 95.35% of the genetic diversity of the breeding populations of the dam and sire line, respectively. Reference-guided variant discovery using the GATK revealed 26,862,369 polymorphic sites. Principal component, admixture and fixation index (F ST ) analyses indicated considerable genetic differentiation between the lines. Genomic inbreeding quantified using runs of homozygosity was higher in the sire than dam line (0.28 vs 0.26). Using two complementary approaches, we detected 51 signatures of selection. However, only six signatures of selection overlapped between both lines. We used the sequenced haplotypes of the 70 key ancestors as a reference panel to call 22,618,811 genotypes in 175 pigs that had been sequenced at very low coverage (1.11-fold) using the GLIMPSE software. The genotype concordance, non-reference sensitivity and non-reference discrepancy between thus inferred and Illumina PorcineSNP60 BeadChip-called genotypes was 97.60, 98.73 and 3.24%, respectively. The low-pass sequencing-derived genomic relationship coefficients were highly correlated ( r > 0.99) with those obtained from microarray genotyping. Conclusions We assessed genetic diversity within and between two lines of the Swiss Large White pig breed. Our analyses revealed considerable differentiation, even though the split into two populations occurred only few generations ago. The sequenced haplotypes of the key ancestor animals enabled us to implement genotyping by low-pass sequencing which offers an intriguing cost-effective approach to increase the variant density over current array-based genotyping by more than 350-fold.

SummaryA recessive form of arthrogryposis multiplex congenita (AMC) was detected 20 years ago in the Swiss Large White (SLW) pig population. A diagnostic marker test enabled the identification of carrier animals, but the underlying causal mutation remains unknown. To identify the mutation underlying AMC, we collected SNP chip genotyping data for 11 affected piglets and 23 healthy pigs. Association testing using 47 829 SNPs confirmed that AMC maps to SSC5 (P = 9.4 × 10−13). Subsequent autozygosity mapping revealed a common 6.06 Mb region (from 66 757 970 to 72 815 151 bp) of extended homozygosity in 11 piglets affected by AMC. Using WGS data, we detected a 63‐bp insertion compatible with the recessive inheritance of AMC in the second exon of KIF21A gene encoding Kinesin Family Member 21A. The 63‐bp insertion is predicted to introduce a premature stop codon in KIF21A gene (p.Val41_Phe42insTer) that truncates 1614 amino acids (~97%) from the protein. We found that this deleterious allele still segregates at a frequency of 0.1% in the SLW pig population. Carrier animals can now be detected unambiguously and excluded from breeding.

A recessive sperm defect of Yorkshire boars was detected more than a decade ago. Affected boars produce ejaculates that contain spermatozoa with defective acrosomes, resulting in low fertility. The acrosome defect was mapped to porcine chromosome 15 but the causal mutation has not been identified. We re-analyzed microarray-derived genotypes of affected boars and confirmed that the acrosome defect maps to a 12.24 Mb segment of porcine chromosome 15. To detect the mutation causing defective acrosomes, we sequenced the genomes of two affected and three unaffected boars to an average coverage of 11-fold. Read depth analysis revealed a 55 kb deletion that is associated with the acrosome defect. The deletion encompasses the BOLL gene encoding the boule homolog, an RNA binding protein which is an evolutionarily conserved member of the DAZ (Deleted in AZoospermia) gene family. Lack of BOLL expression causes spermatogenic arrest and sperm maturation failure in many species. Boars that carry the deletion in the homozygous state produce sperm but their acrosomes are defective, suggesting that lack of porcine BOLL compromises acrosome formation. Our findings warrant further research to investigate the role of BOLL during spermatogenesis and sperm maturation in pigs.

Artificial insemination in pig (Sus scrofa domesticus) breeding involves the evaluation of the semen quality of breeding boars. Ejaculates that fulfill predefined quality requirements are processed, diluted and used for inseminations. Within short time, eight Swiss Large White boars producing immotile sperm that had multiple morphological abnormalities of the sperm flagella were noticed. The eight boars were inbred on a common ancestor suggesting that the novel sperm flagella defect is a recessive trait. Transmission electron microscopy cross-sections revealed that the immotile sperm had disorganized flagellar axonemes. Haplotype-based association testing involving microarray-derived genotypes at 41,094 SNPs of six affected and 100 fertile boars yielded strong association (P=4.22 x 10−15) at chromosome 12. Autozygosity mapping enabled us to pinpoint the causal mutation on a 1.11 Mb haplotype located between 3,473,632 and 4,587,759 bp. The haplotype carries an intronic 13-bp deletion (Chr12:3,556,401-3,556,414 bp) that is compatible with recessive inheritance. The 13-bp deletion excises the polypyrimidine tract upstream exon 56 of DNAH17 (XM_021066525.1:c.8510-17_8510-5del) encoding dynein axonemal heavy chain 17. Transcriptome analysis of the testis of two affected boars revealed that the loss of the polypyrimidine tract causes exon skipping which results in the in-frame loss of 89 amino acids from DNAH17. Disruption of DNAH17 impairs the assembly of the flagellar axoneme and manifests in multiple morphological abnormalities of the sperm flagella. Direct gene testing may now be implemented to monitor the defective allele in the Swiss Large White population and prevent the frequent manifestation of a sterilizing sperm tail disorder in breeding boars.

Artificial insemination in pig (Sus scrofa domesticus) breeding involves the evaluation of the semen quality of breeding boars. Ejaculates that fulfill predefined quality requirements are processed, diluted and used for inseminations. Within short time, eight Swiss Large White boars producing immotile sperm that had multiple morphological abnormalities of the sperm flagella were noticed at a semen collection center. The eight boars were inbred on a common ancestor suggesting that the novel sperm flagella defect is a recessive trait. Transmission electron microscopy cross-sections revealed that the immotile sperm had disorganized flagellar axonemes. Haplotype-based association testing involving microarray-derived genotypes at 41,094 SNPs of six affected and 100 fertile boars yielded strong association (P = 4.22 × 10-15) at chromosome 12. Autozygosity mapping enabled us to pinpoint the causal mutation on a 1.11 Mb haplotype located between 3,473,632 and 4,587,759 bp. The haplotype carries an intronic 13-bp deletion (Chr12:3,556,401-3,556,414 bp) that is compatible with recessive inheritance. The 13-bp deletion excises the polypyrimidine tract upstream exon 56 of DNAH17 (XM_021066525.1: c.8510-17_8510-5del) encoding dynein axonemal heavy chain 17. Transcriptome analysis of the testis of two affected boars revealed that the loss of the polypyrimidine tract causes exon skipping which results in the in-frame loss of 89 amino acids from DNAH17. Disruption of DNAH17 impairs the assembly of the flagellar axoneme and manifests in multiple morphological abnormalities of the sperm flagella. Direct gene testing may now be implemented to monitor the defective allele in the Swiss Large White population and prevent the frequent manifestation of a sterilizing sperm tail disorder in breeding boars.

The phenotypic relationships between type traits and functional traits were analysed in Czech Holstein dairy cows born between years 2002 and 2015, with minimum proportion of Holstein genes 88%. Two slightly different models were used to evaluate the effects of 1 measured trait (in cm), 20 linear type traits (9 classes each), 4 composite traits, and final score (both with scales 0 – 100 points) on milk production traits and longevity. Sample for analysis of effect on milk yield included observations from 247 790 cows within one to four lactations (570 671 rows in total), second sample for longevity study included data from 228 161 cows with sums of one to maximum six lactations. Cows were required to obtain type classification scores between 30th and 210th day of the first lactation in age between 650 and 1206 days. Milk yield (in kg) records below 5079 kg and above 16 622 kg were set to absent.