Estimation of direct and social effects of feeding duration in growing pigs using records from automatic feeding stations.

Automatic feeding systems in pig production allow for the recording of individual feeding behavior traits, which might be influenced by the social interactions among individuals. This study fitted mixed models to estimate the direct and social effects on visit duration at the feeder of group-housed pigs. The dataset included 74,413 records of each visit duration time (min) event at the automatic feeder from 135 pigs housed in 14 pens. The sequence of visits at the feeder was employed as a proxy for the social interaction between individuals. To estimate animal effects, the direct effect was apportioned to the animal feeding (feeding pig), and the social effect was apportioned to the animal that entered the feeder immediately after the feeding pig left the feeding station (follower). The data were divided into two subsets: "non-immediate replacement" time (NIRT, N = 6,256), where the follower pig occupied the feeder at least 600 s after the feeding pig left the feeder, and "immediate replacement" time (IRT, N = 58,255), where the elapsed time between replacements was less than or equal to 60 s. The marginal posterior distribution of the parameters was obtained by Bayesian method. Using the IRT subset, the posterior mean of the proportion of variance explained by the direct effect (PrpσTemefos) was 18% for all models. The proportion of variance explained by the follower social effect (Prpσ^f2) was 2%, and the residual variance (σ^e2) decreased, suggesting an improved model fit by including the follower effect. Fitting the models with the NIRT subset, the estimate of PrpσTemefos was 20% but the Prpσ^f2 was almost zero and σ^e2 was identical for all models. For the IRT subset, the predicted best linear unbiased predictor (BLUP) of direct (Direct BLUP) and social (Follower BLUP) random effects on visit duration at the feeder of an animal was calculated. Feeder visit duration time was not correlated with traits, such as weight gain or average feed intake (P > 0.05), whereas for the daily feeder occupation time, the estimated correlation was positive with the Direct BLUP (r^ = 0.51, P < 0.05) and negative with the Follower BLUP (r^= -0.26, P < 0.05). The results suggest that the visit duration of an animal at the single-space feeder was influenced by both direct and social effects when the replacement time between visits was less than 1 min. Finally, animals that spent a longer time per day at the feeder seemed to do so by shortening the meal length of the preceding individual at the feeder.

Emerging issues in genomic selection.

Genomic selection (GS) is now practiced successfully across many species. However, many questions remain, such as long-term effects, estimations of genomic parameters, robustness of genome-wide association study (GWAS) with small and large datasets, and stability of genomic predictions. This study summarizes presentations from the authors at the 2020 American Society of Animal Science (ASAS) symposium. The focus of many studies until now is on linkage disequilibrium between two loci. Ignoring higher-level equilibrium may lead to phantom dominance and epistasis. The Bulmer effect leads to a reduction of the additive variance; however, the selection for increased recombination rate can release anew genetic variance. With genomic information, estimates of genetic parameters may be biased by genomic preselection, but costs of estimation can increase drastically due to the dense form of the genomic information. To make the computation of estimates feasible, genotypes could be retained only for the most important animals, and methods of estimation should use algorithms that can recognize dense blocks in sparse matrices. GWASs using small genomic datasets frequently find many marker-trait associations, whereas studies using much bigger datasets find only a few. Most of the current tools use very simple models for GWAS, possibly causing artifacts. These models are adequate for large datasets where pseudo-phenotypes such as deregressed proofs indirectly account for important effects for traits of interest. Artifacts arising in GWAS with small datasets can be minimized by using data from all animals (whether genotyped or not), realistic models, and methods that account for population structure. Recent developments permit the computation of P-values from genomic best linear unbiased prediction (GBLUP), where models can be arbitrarily complex but restricted to genotyped animals only, and single-step GBLUP that also uses phenotypes from ungenotyped animals. Stability was an important part of nongenomic evaluations, where genetic predictions were stable in the absence of new data even with low prediction accuracies. Unfortunately, genomic evaluations for such animals change because all animals with genotypes are connected. A top-ranked animal can easily drop in the next evaluation, causing a crisis of confidence in genomic evaluations. While correlations between consecutive genomic evaluations are high, outliers can have differences as high as 1 SD. A solution to fluctuating genomic evaluations is to base selection decisions on groups of animals. Although many issues in GS have been solved, many new issues that require additional research continue to surface.

Estimation of genome-wide and locus-specific breed composition in pigs.

Advances in pig genomic technologies enable implementation of new methods to estimate breed composition, allowing innovative and efficient ways to evaluate and ensure breed and line background. Existing methods to test for homozygosity at key loci involve test mating the animal in question and observing phenotypic patterns among offspring, requiring extensive resources. In this study, whole-genome pig DNA microarray data from over 8,000 SNP was used to profile the composition of U.S. registered purebred pigs using a refined linear regression method that enhances the interpretation of coefficients. In a simulation analysis, a strong correlation between true and estimated breed composition was observed (R2 = 0.94). Applying these methods to 930 Yorkshire animals registered with the National Swine Registry, 95% were estimated to have a "genome-wide" Yorkshire breed composition of at least 0.825 or 82.5%, with similar performance for evaluating datasets of registered Duroc (n = 88) Landrace (n = 129), and Hampshire (n = 17) breeds. We also developed new methods to evaluate locus-based breed probabilities. Such methods have been applied to multi-locus SNP genotypes flanking the KIT gene known to predominantly control coat color, thereby inferring the probability that an animal has haplotypes in the KIT region that are predominant in white breeds. These methods have been adopted by the National Swine Registry as a means to identify purebred Yorkshire animals.

Differential Effects of Sex Pheromone Compounds on Adult Female Sea Lamprey (Petromyzon marinus) Locomotor Patterns.

Synchronization of male and female locomotor activity plays a critical role in ensuring reproductive success, especially in semelparous species. The goal of this study was to elucidate the effects of individual chemical signals, or pheromones, on the locomotor activity in the sea lamprey (Petromyzon marinus). In their native habitat, adult preovulated females (POF) and ovulated females (OF) are exposed to sex pheromone compounds that are released from spermiated males and attract females to nests during their migration and spawning periods. In this study, locomotor activity of individual POF and OF was measured hourly in controlled laboratory conditions using an automated video-tracking system. Differences in the activity between a baseline day (no treatment exposure) and a treatment day (sex pheromone compound or control exposure) were examined for daytime and nighttime periods. Results showed that different pheromone compound treatments affected both POF and OF sea lamprey (p < 0.05) but in different ways. Spermiated male washings (SMW) and one of its main components, 7α,12α,24-trihydroxy-5α-cholan-3-one 24 sulfate (3kPZS), decreased activity of POF during the nighttime. SMW also reduced activity in POF during the daytime. In contrast, SMW increased activity of OF during the daytime, and an additional compound found in SMW, petromyzonol sulfate (PZS), decreased the activity during the nighttime. In addition, we examined factors that allowed us to infer the overall locomotor patterns. SMW increased the maximum hourly activity during the daytime, decreased the maximum hourly activity during the nighttime, and reduced the percentage of nocturnal activity in OF. Our findings suggest that adult females have evolved to respond to different male compounds in regards to their locomotor activity before and after final maturation. This is a rare example of how species-wide chemosensory stimuli can affect not only the amounts of activity but also the overall locomotor pattern in a vertebrate species.

Effects of sex pheromones and sexual maturation on locomotor activity in female sea lamprey (Petromyzon marinus).

Synchronization of male and female locomotor rhythmicity can play a vital role in ensuring reproductive success. Several physiological and environmental factors alter these locomotor rhythms. As sea lamprey, Petromyzon marinus, progress through their life cycle, their locomotor activity rhythm changes multiple times. The goal of this study was to elucidate the activity patterns of adult female sea lamprey during the sexual maturation process and discern the interactions of these patterns with exposure to male pheromones. During these stages, preovulated and ovulated adult females are exposed to sex pheromone compounds, which are released by spermiated males and attract ovulated females to the nest for spawning. The locomotor behavior of adult females was monitored in a natural stream with a passive integrated tag responder system as they matured, and they were exposed to a sex pheromone treatment (spermiated male washings) or a control (prespermiated male washings). Results showed that, dependent on the hour of day, male sex pheromone compounds reduce total activity (p < 0.05) and cause increases in activity during several daytime hours in preovulated and ovulated females. These results are one of the first examples of how sex pheromones modulate a locomotor rhythm in a vertebrate, and they suggest that the interaction between maturity stage and sex pheromone exposure contributes to the differential locomotor rhythms found in adult female sea lamprey. This phenomenon may contribute to the reproductive synchrony of mature adults, thus increasing reproductive success in this species.

Probing genetic control of swine responses to PRRSV infection: current progress of the PRRS host genetics consortium.

BACKGROUND: Understanding the role of host genetics in resistance to porcine reproductive and respiratory syndrome virus (PRRSV) infection, and the effects of PRRS on pig health and related growth, are goals of the PRRS Host Genetics Consortium (PHGC). METHODS: The project uses a nursery pig model to assess pig resistance/susceptibility to primary PRRSV infection. To date, 6 groups of 200 crossbred pigs from high health farms were donated by commercial sources. After acclimation, the pigs were infected with PRRSV in a biosecure facility and followed for 42 days post infection (dpi). Blood samples were collected at 0, 4, 7, 10, 14, 21, 28, 35 and 42 dpi for serum and whole blood RNA gene expression analyses; weekly weights were recorded for growth traits. All data have been entered into the PHGC relational database. Genomic DNAs from all PHGC1-6 pigs were prepared and genotyped with the Porcine SNP60 SNPchip. RESULTS: Results have affirmed that all challenged pigs become PRRSV infected with peak viremia being observed between 4-21 dpi. Multivariate statistical analyses of viral load and weight data have identified PHGC pigs in different virus/weight categories. Sera are now being compared for factors involved in recovery from infection, including speed of response and levels of immune cytokines. Genome-wide association studies (GWAS) are underway to identify genes and chromosomal locations that identify PRRS resistant/susceptible pigs and pigs able to maintain growth while infected with PRRSV. CONCLUSIONS: Overall, the PHGC project will enable researchers to discover and verify important genotypes and phenotypes that predict resistance/susceptibility to PRRSV infection. The availability of PHGC samples provides a unique opportunity to continue to develop deeper phenotypes on every PRRSV infected pig.

Deciphering the luteal transcriptome: potential mechanisms mediating stage-specific luteolytic response of the corpus luteum to prostaglandin F2α.

The objective of this study was to identify prostaglandin F(2α) (PG)-induced changes in the transcriptome of bovine corpora lutea (CL) that are specific to mature, PG-responsive (day 11) CL vs. developing (day 4) CL, which do not undergo luteolysis in response to PG administration. CL were collected at 0, 4, and 24 h after PG injection on days 4 and 11 of the estrous cycle (n = 5 per day and time point), and microarray analysis was performed with GeneChip Bovine Genome Arrays. Data normalization was performed with affy package and significance testing with maanova from Bioconductor. Significance (relative to 0 h time point) was declared at fold change >2.0 or <0.5 and false discovery rate of <5%. At 4 and 24 h after PG, 221 (day 4) and 661 (day 11) and 248 (day 4) and 1,421 (day 11) regulated genes, respectively, were identified. The accentuated gene expression response in day 11 CL was accompanied by specific enrichment of PG-regulated genes in distinctive gene ontology categories (immune related and other), particularly at 24 h after injection. Specificity in putative transcription factor binding sites was observed among PG-regulated genes on day 11 vs. day 4, including a potential association of ETS transcription factors with acute PG-induced gene expression specific to day 11 CL. Temporal and PG-induced regulation of abundance of mRNA for ETS transcription factor family members linked to the stage-specific response to PG was not observed. Increased abundance of protein and/or mRNA for six PG-regulated putative ETS-responsive genes was noted in day 11 but not day 4 CL. Results reveal insight into stage-specific gene expression in bovine CL in response to PG and potential transcriptional mediators of luteolysis.

Genome-wide linkage analysis of global gene expression in loin muscle tissue identifies candidate genes in pigs.

BACKGROUND: Nearly 6,000 QTL have been reported for 588 different traits in pigs, more than in any other livestock species. However, this effort has translated into only a few confirmed causative variants. A powerful strategy for revealing candidate genes involves expression QTL (eQTL) mapping, where the mRNA abundance of a set of transcripts is used as the response variable for a QTL scan. METHODOLOGY/PRINCIPAL FINDINGS: We utilized a whole genome expression microarray and an F(2) pig resource population to conduct a global eQTL analysis in loin muscle tissue, and compared results to previously inferred phenotypic QTL (pQTL) from the same experimental cross. We found 62 unique eQTL (FDR <10%) and identified 3 gene networks enriched with genes subject to genetic control involved in lipid metabolism, DNA replication, and cell cycle regulation. We observed strong evidence of local regulation (40 out of 59 eQTL with known genomic position) and compared these eQTL to pQTL to help identify potential candidate genes. Among the interesting associations, we found aldo-keto reductase 7A2 (AKR7A2) and thioredoxin domain containing 12 (TXNDC12) eQTL that are part of a network associated with lipid metabolism and in turn overlap with pQTL regions for marbling, % intramuscular fat (% fat) and loin muscle area on Sus scrofa (SSC) chromosome 6. Additionally, we report 13 genomic regions with overlapping eQTL and pQTL involving 14 local eQTL. CONCLUSIONS/SIGNIFICANCE: Results of this analysis provide novel candidate genes for important complex pig phenotypes.

A powerful and flexible linear mixed model framework for the analysis of relative quantification RT-PCR data.

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is currently viewed as the most precise technique to quantify levels of messenger RNA. Relative quantification compares the expression of a target gene under two or more experimental conditions normalized to the measured expression of a control gene. The statistical methods and software currently available for the analysis of relative quantification of RT-PCR data lack the flexibility and statistical properties to produce valid inferences in a wide range of experimental situations. In this paper we present a novel method for the analysis of relative quantification of qRT-PCR data, which consists of the analysis of cycles to threshold values (C(T)) for a target and a control gene using a general linear mixed model methodology. Our method allows testing of a broader class of hypotheses than traditional analyses such as the classical comparative C(T). Moreover, a simulation study using plasmode datasets indicated that the estimated fold-change in pairwise comparisons was the same using either linear mixed models or a comparative C(T) method, but the linear mixed model approach was more powerful. In summary, the method presented in this paper is more accurate, powerful and flexible than the traditional methods for analysis of qRT-PCR data. This new method is especially useful for studies involving multiple experimental factors and complex designs.

Characterization and in vitro control of MPF activity in zebrafish eggs.

We describe the characterization of maturation-promoting factor (MPF) in zebrafish eggs and used different defined conditions to maintain its activity in vitro. MPF activity levels are high in freshly ovulated mature eggs and decline rapidly within 5 min after either fertilization or parthenogenetic activation. The MPF activity of eggs matured in vitro declines faster when the eggs are incubated in Hank's culture medium supplemented with 0.5% BSA (H-BSA) than when incubated in Chinook salmon ovarian fluid (CSOF). MPF activity in nonactivated, aged eggs remains high in H-BSA supplemented with 75 microM MG132 or 10 mM caffeine, but neither MG132 nor caffeine can sustain high MPF activity in activated eggs. MG132-treated eggs showed delayed completion of metaphase and extrusion of the second polar body. Nuclear staining of the activated eggs confirmed the correlation between their cell cycle stage and MPF activity at each time point. An embryotoxic effect was found when matured eggs were held in 100 microM of MG132 or 20 mM caffeine for 1 h. Calcium-depleted medium and 1,2-bis(o-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid also showed detrimental effects on the embryos. Conversely, nonactivated, aged matured eggs maintained high MPF activity and developmental potential when CSOF was used as a holding medium.

Characterization of a canine model of autosomal recessive retinitis pigmentosa due to a PDE6A mutation.

PURPOSE: To characterize a canine model of autosomal recessive RP due to a PDE6A gene mutation. METHODS: Affected and breed- and age-matched control puppies were studied by electroretinography (ERG), light and electron microscopy, immunohistochemistry, and assay for retinal PDE6 levels and enzymatic activity. RESULTS: The mutant puppies failed to develop normal rod-mediated ERG responses and had reduced light-adapted a-wave amplitudes from an early age. The residual ERG waveforms originated primarily from cone-driven responses. Development of photoreceptor outer segments stopped, and rod cells were lost by apoptosis. Immunohistochemistry demonstrated a marked reduction in rod opsin immunostaining outer segments and relative preservation of cones early in the disease process. With exception of rod bipolar cells, which appeared to be reduced in number relatively early in the disease process, other inner retinal cells were preserved in the early stages of the disease, although there was marked and early activation of Müller glia. Western blot analysis showed that the PDE6A mutation not only resulted in a lack of PDE6A protein but the affected retinas also lacked the other PDE6 subunits, suggesting expression of PDE6A is essential for normal expression of PDE6B and PDE6G. Affected retinas lacked PDE6 enzymatic activity. CONCLUSIONS: This represents the first characterization of a PDE6A model of autosomal recessive retinitis pigmentosa, and the PDE6A mutant dog shows promise as a large animal model for investigation of therapies to rescue mutant rod photoreceptors and to preserve cone photoreceptors in the face of a rapid loss of rod cells.

Full conjugate analysis of normal multiple traits with missing records using a generalized inverted Wishart distribution.

A Markov chain Monte Carlo (MCMC) algorithm to sample an exchangeable covariance matrix, such as the one of the error terms (R0) in a multiple trait animal model with missing records under normal-inverted Wishart priors is presented. The algorithm (FCG) is based on a conjugate form of the inverted Wishart density that avoids sampling the missing error terms. Normal prior densities are assumed for the 'fixed' effects and breeding values, whereas the covariance matrices are assumed to follow inverted Wishart distributions. The inverted Wishart prior for the environmental covariance matrix is a product density of all patterns of missing data. The resulting MCMC scheme eliminates the correlation between the sampled missing residuals and the sampled R0, which in turn has the effect of decreasing the total amount of samples needed to reach convergence. The use of the FCG algorithm in a multiple trait data set with an extreme pattern of missing records produced a dramatic reduction in the size of the autocorrelations among samples for all lags from 1 to 50, and this increased the effective sample size from 2.5 to 7 times and reduced the number of samples needed to attain convergence, when compared with the 'data augmentation' algorithm.