A survey of single nucleotide polymorphisms identified from whole-genome sequencing and their functional effect in the porcine genome.

Genetic variants detected from sequence have been used to successfully identify causal variants and map complex traits in several organisms. High and moderate impact variants, those expected to alter or disrupt the protein coded by a gene and those that regulate protein production, likely have a more significant effect on phenotypic variation than do other types of genetic variants. Hence, a comprehensive list of these functional variants would be of considerable interest in swine genomic studies, particularly those targeting fertility and production traits. Whole-genome sequence was obtained from 72 of the founders of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the founding boars (12 Duroc and 12 Landrace) and 48 Yorkshire-Landrace composite sows. Sequence reads were mapped to the Sscrofa10.2 genome build, resulting in a mean of 6.1 fold (×) coverage per genome. A total of 22 342 915 high confidence SNPs were identified from the sequenced genomes. These included 21 million previously reported SNPs and 79% of the 62 163 SNPs on the PorcineSNP60 BeadChip assay. Variation was detected in the coding sequence or untranslated regions (UTRs) of 87.8% of the genes in the porcine genome: loss-of-function variants were predicted in 504 genes, 10 202 genes contained nonsynonymous variants, 10 773 had variation in UTRs and 13 010 genes contained synonymous variants. Approximately 139 000 SNPs were classified as loss-of-function, nonsynonymous or regulatory, which suggests that over 99% of the variation detected in our pigs could potentially be ignored, allowing us to focus on a much smaller number of functional SNPs during future analyses.

Genome-wide copy number variation in the bovine genome detected using low coverage sequence of popular beef breeds.

Copy number variations (CNVs) are large insertions, deletions or duplications in the genome that vary between members of a species and are known to affect a wide variety of phenotypic traits. In this study, we identified CNVs in a population of bulls using low coverage next-generation sequence data. First, in order to determine a suitable strategy for CNV detection in our data, we compared the performance of three distinct CNV detection algorithms on benchmark CNV datasets and concluded that using the multiple sample read depth approach was the best method for identifying CNVs in our sequences. Using this technique, we identified a total of 1341 copy number variable regions (CNVRs) from genome sequences of 154 purebred sires used in Cycle VII of the USMARC Germplasm Evaluation Project. These bulls represented the seven most popular beef breeds in the United States: Hereford, Charolais, Angus, Red Angus, Simmental, Gelbvieh and Limousin. The CNVRs covered 6.7% of the bovine genome and spanned 2465 protein-coding genes and many known quantitative trait loci (QTL). Genes harbored in the CNVRs were further analyzed to determine their function as well as to find any breed-specific differences that may shed light on breed differences in adaptation, health and production.

Linkage disequilibrium among commonly genotyped SNP variants detected from bull sequence.

Genomic prediction utilizing causal variants could increase selection accuracy above that achieved with SNPs genotyped by currently available arrays used for genomic selection. A number of variants detected from sequencing influential sires are likely to be causal, but noticeable improvements in prediction accuracy using imputed sequence variant genotypes have not been reported. Improvement in accuracy of predicted breeding values may be limited by the accuracy of imputed sequence variants. Using genotypes of SNPs on a high-density array and non-synonymous SNPs detected in sequence from influential sires of a multibreed population, results of this examination suggest that linkage disequilibrium between non-synonymous and array SNPs may be insufficient for accurate imputation from the array to sequence. In contrast to 75% of array SNPs being strongly correlated to another SNP on the array, less than 25% of the non-synonymous SNPs were strongly correlated to an array SNP. When correlations between non-synonymous and array SNPs were strong, distances between the SNPs were greater than separation that might be expected based on linkage disequilibrium decay. Consistently near-perfect whole-genome linkage disequilibrium between the full array and each non-synonymous SNP within the sequenced bulls suggests that whole-genome approaches to infer sequence variants might be more accurate than imputation based on local haplotypes. Opportunity for strong linkage disequilibrium between sequence and array SNPs may be limited by discrepancies in allele frequency distributions, so investigating alternate genotyping approaches and panels providing greater chances of frequency-matched SNPs strongly correlated to sequence variants is also warranted. Genotypes used for this study are available from https://www.animalgenome.org/repository/pub/;USDA2017.0519/.

Porcine single nucleotide polymorphisms and their functional effect: an update.

OBJECTIVE: To aid in the development of a comprehensive list of functional variants in the swine genome, single nucleotide polymorphisms (SNP) were identified from whole genome sequence of 240 pigs. Interim data from 72 animals in this study was published in 2017. This communication extends our previous work not only by utilizing genomic sequence from additional animals, but also by the use of the newly released Sscrofa 11.1 reference genome. RESULTS: A total of 26,850,263 high confidence SNP were identified, including 19,015,267 reported in our previously published results. Variation was detected in the coding sequence or untranslated regions (UTR) of 78% of the genes in the porcine genome: 1729 loss-of-function variants were predicted in 1162 genes, 12,686 genes contained 64,232 nonsynonymous variants, 250,403 variants were present in UTR of 15,739 genes, and 15,284 genes contained 90,939 synonymous variants. In total, approximately 316,000 SNP were classified as being of high to moderate impact (i.e. loss-of-function, nonsynonymous, or regulatory). These high to moderate impact SNP will be the focus of future genome-wide association studies.

Beef steers with average dry matter intake and divergent average daily gain have altered gene expression in the jejunum.

The objective of this study was to determine the association of differentially expressed genes (DEG) in the jejunum of steers with average DMI and high or low ADG. Feed intake and growth were measured in a cohort of 144 commercial Angus steers consuming a finishing diet containing (on a DM basis) 67.8% dry-rolled corn, 20% wet distillers grains with solubles, 8% alfalfa hay, and 4.2% vitamin/mineral supplement. From the cohort, a subset of steers with DMI within ±0.32 SD of the mean for DMI and the greatest (high) and least (low) ADG were chosen for slaughter and jejunum mucosa collection ( = 8 for each group). Dry matter intake (10.1 ± 0.05 kg/d) was not different ( = 0.41) but ADG was greater in the high-gain group (2.17 and 1.72 ± 0.02 kg/d for the high- and low-ADG groups, respectively; < 0.01). A total of 13,747 genes were found to be expressed in the jejunum, of which 64 genes were differentially expressed between the 2 groups (corrected < 0.05). Ten of the DEG were upregulated in the low-ADG group and 54 were upregulated in the high-ADG group. Gene ontology analysis determined that 24 biological process terms were overrepresented ( < 0.05), including digestion, drug and xenobiotic metabolism, and carbohydrate metabolism. Additionally, 89 molecular function terms were enriched ( < 0.05), including metallopeptidase activity, transporter activity, steroid hydrolase activity, glutathione transferase activity, and chemokine receptor binding. Metabolic pathways (28 pathways) impacted by the DEG ( < 0.05) included drug and xenobiotic metabolism by cytochrome P450, carbohydrate digestion and absorption, vitamin digestion and absorption, galactose metabolism, and linoleic acid metabolism. Results from this experiment indicate that cattle with average DMI and greater ADG likely have a greater capacity to handle foreign substances (xenobiotics). It is also possible that cattle with a greater ADG have a greater potential to digest and absorb nutrients in the small intestine.

GENOMICS SYMPOSIUM: Using genomic approaches to uncover sources of variation in age at puberty and reproductive longevity in sows.

Genetic variants associated with traits such as age at puberty and litter size could provide insight into the underlying genetic sources of variation impacting sow reproductive longevity and productivity. Genomewide characterization and gene expression profiling were used using gilts from the University of Nebraska-Lincoln swine resource population ( = 1,644) to identify genetic variants associated with age at puberty and litter size traits. From all reproductive traits studied, the largest fraction of phenotypic variation explained by the Porcine SNP60 BeadArray was for age at puberty (27.3%). In an evaluation data set, the predictive ability of all SNP from high-ranked 1-Mb windows (1 to 50%), based on genetic variance explained in training, was greater (12.3 to 36.8%) compared with the most informative SNP from these windows (6.5 to 23.7%). In the integrated data set ( = 1,644), the top 1% of the 1-Mb windows explained 6.7% of the genetic variation of age at puberty. One of the high-ranked windows detected (SSC2, 12-12.9 Mb) showed pleiotropic features, affecting both age at puberty and litter size traits. The RNA sequencing of the hypothalami arcuate nucleus uncovered 17 differentially expressed genes (adjusted < 0.05) between gilts that became pubertal early (<155 d of age) and late (>180 d of age). Twelve of the differentially expressed genes are upregulated in the late pubertal gilts. One of these genes is involved in energy homeostasis (), a function in which the arcuate nucleus plays an important contribution, linking nutrition with reproductive development. Energy restriction during the gilt development period delayed age at puberty by 7 d but increased the probability of a sow to produce up to 3 parities ( < 0.05). Identification of pleotropic functional polymorphisms may improve accuracy of genomic prediction while facilitating a reduction in sow replacement rates and addressing welfare concerns.

Genomewide association study of liver abscess in beef cattle.

Fourteen percent of U.S. cattle slaughtered in 2011 had liver abscesses, resulting in reduced carcass weight, quality, and value. Liver abscesses can result from a common bacterial cause, , which inhabits rumen lesions caused by acidosis and subsequently escapes into the blood stream, is filtered by the liver, and causes abscesses in the liver. Our aim was to identify SNP associated with liver abscesses in beef cattle. We used lung samples as a DNA source because they have low economic value, they have abundant DNA, and we had unrestricted access to sample them. We collected 2,304 lung samples from a beef processing plant: 1,152 from animals with liver abscess and 1,152 from animals without liver abscess. Lung tissue from pairs of animals, 1 with abscesses and another without, were collected from near one another on the viscera table to ensure that pairs of phenotypically extreme animals came from the same lot. Within each phenotype (abscess or no abscess), cattle were pooled by slaughter sequence into 12 pools of 96 cattle for each phenotype for a total of 24 pools. The pools were constructed by equal volume of frozen lung tissue from each animal. The DNA needed to allelotype each pool was then extracted from pooled lung tissue and the BovineHD Bead Array (777,962 SNP) was run on all 24 pools. Total intensity (TI), an indicator of copy number variants, was the sum of intensities from red and green dyes. Pooling allele frequency (PAF) was red dye intensity divided TI. Total intensity and PAF were weighted by the inverse of their respective genomic covariance matrices computed over all SNP across the genome. A false discovery rate ≤ 5% was achieved for 15 SNP for PAF and 20 SNP for TI. Genes within 50 kbp from significant SNP were in diverse pathways including maintenance of pH homeostasis in the gastrointestinal tract, maintain immune defenses in the liver, migration of leukocytes from the blood into infected tissues, transport of glutamine into the kidney in response to acidosis to facilitate production of bicarbonate to increase pH, aggregate platelets to liver injury to facilitate liver repair, and facilitate axon guidance. Evidence from the 35 detected SNP associations combined with evidence of polygenic variation indicate that there is adequate genetic variation in incidence rate of liver abscesses, which could be exploited to select sires for reduced susceptibility to subacute acidosis and associated liver abscess.