Xiao-Lin Wu

Evaluation of candidate gene effects for beef backfat via Bayesian model selection

Candidate gene approaches provide tools for exploring and localizing causative genes affecting quantitative traits and the underlying variation may be better understood by determining the relative magnitudes of effects of their polymorphisms. Diacyglycerol O-acyltransferase 1 (DGAT1), fatty acid binding protein (heart) 3 (FABP3), growth hormone 1 (GH1), leptin (LEP) and thyroglobulin (TG) have been previously identified as genes contributing to genetic control of subcutaneous fat thickness (SFT) in beef cattle. In the present research, Bayesian model selection was used to evaluate effects of these five candidate genes by comparing competing non-nested models and treating candidate gene effects as either random or fixed. The analyses were implemented in SAS to simplify the programming and computation. Phenotypic data were gathered from a F2 population of Wagyu × Limousin cattle. The five candidate genes had significant but varied effects on SFT in this population. Bayesian model selection identified the DGAT1 model as the one with the greatest model probability, whether candidate gene effects were considered random or fixed, and DGAT1 had the greatest additive effect on SFT. The SAS codes developed in the study are freely available and can be downloaded at: http://www.ansci.wsu.edu/programs/.

Genome-Wide Genetic Diversity and Differentially Selected Regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee Sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.

Reactomes of porcine alveolar macrophages infected with porcine reproductive and respiratory syndrome virus.

Porcine reproductive and respiratory syndrome (PRRS) has devastated pig industries worldwide for many years. It is caused by a small RNA virus (PRRSV), which targets almost exclusively pig monocytes or macrophages. In the present study, five SAGE (serial analysis of gene expression) libraries derived from 0 hour mock-infected and 6, 12, 16 and 24 hours PRRSV-infected porcine alveolar macrophages (PAMs) produced a total 643,255 sequenced tags with 91,807 unique tags. Differentially expressed (DE) tags were then detected using the Bayesian framework followed by gene/mRNA assignment, arbitrary selection and manual annotation, which determined 699 DE genes for reactome analysis. The DAVID, KEGG and REACTOME databases assigned 573 of the DE genes into six biological systems, 60 functional categories and 504 pathways. The six systems are: cellular processes, genetic information processing, environmental information processing, metabolism, organismal systems and human diseases as defined by KEGG with modification. Self-organizing map (SOM) analysis further grouped these 699 DE genes into ten clusters, reflecting their expression trends along these five time points. Based on the number one functional category in each system, cell growth and death, transcription processes, signal transductions, energy metabolism, immune system and infectious diseases formed the major reactomes of PAMs responding to PRRSV infection. Our investigation also focused on dominant pathways that had at least 20 DE genes identified, multi-pathway genes that were involved in 10 or more pathways and exclusively-expressed genes that were included in one system. Overall, our present study reported a large set of DE genes, compiled a comprehensive coverage of pathways, and revealed system-based reactomes of PAMs infected with PRRSV. We believe that our reactome data provides new insight into molecular mechanisms involved in host genetic complexity of antiviral activities against PRRSV and lays a strong foundation for vaccine development to control PRRS incidence in pigs.

The Reverse Cholesterol Transport Pathway Improves Understanding of Genetic Networks for Fat Deposition and Muscle Growth in Beef Cattle

In the present study, thirteen genes involved in the reverse cholesterol transport (RCT) pathway were investigated for their associations with three fat depositions, eight fatty acid compositions and two growth-related phenotypes in a Wagyu x Limousin reference population, including 6 F1 bulls, 113 F1 dams, and 246 F2 progeny. A total of 37 amplicons were used to screen single nucleotide polymorphisms (SNPs) on 6 F1 bulls. Among 36 SNPs detected in 11 of these 13 genes, 19 were selected for genotyping by the Sequenom assay design on all F2 progeny. Single-marker analysis revealed seven SNPs in ATP binding cassette A1, apolipoproteins A1, B and E, phospholipid transfer protein and paraoxinase 1 genes significantly associated with nine phenotypes (P<0.05). Previously, we reported genetic networks associated with 19 complex phenotypes based on a total of 138 genetic polymorphisms derived from 71 known functional genes. Therefore, after Bonferroni correction, these significant (adjusted P<0.05) and suggestive (adjusted P<0.10) associations were then used to identify genetic networks related to the RCT pathway. Multiple-marker analysis suggested possible genetic networks involving the RCT pathway for kidney-pelvic-heart fat percentage, rib-eye area, and subcutaneous fat depth phenotypes with markers derived from paraoxinase 1, apolipoproteins A1 and E, respectively. The present study confirmed that genes involved in cholesterol homeostasis are useful targets for investigating obesity in humans as well as for improving meat quality phenotypes in a livestock production.

The Complementary Neighborhood Patterns and Methylation-to-Mutation Likelihood Structures of 15,110 Single-Nucleotide Polymorphisms in the Bovine Genome

Bayesian analysis was performed to examine the single-nucleotide polymorphism (SNPs) neighborhood patterns in cattle using 15,110 SNPs, each with a flanking sequence of 500 bp. Our analysis confirmed three well-known features reported in plants and/or other animals: (1) the transition is the most abundant type of SNPs, accounting for 69.8% in cattle; (2) the transversion occurs most frequently (38.56%) in cattle when the A + T content equals two at their immediate adjacent sites; and (3) C ↔ T and A ↔ G transitions have reverse complementary neighborhood patterns and so do A ↔ C and G ↔ T transversions. Our study also revealed several novel SNP neighborhood patterns that have not been reported previously. First, cattle and humans share an overall SNP pattern, indicating a common mutation system in mammals. Second, unlike C ↔ T/A ↔ G and A ↔ C/G ↔ T, the true neighborhood patterns for A ↔ T and C ↔ G might remain mysterious because the sense and antisense sequences flanking these mutations are not actually recognizable. Third, among the reclassified four types of SNPs, the neighborhood ratio between A + T and G + C was quite different. The ratio was lowest for C ↔ G, but increased for C ↔ T/A ↔ G, further for A ↔ C/G ↔ T, and the most for A ↔ T. Fourth, when two immediate adjacent sites provide structures for CpG, it significantly increased transitions compared to the structures without the CpG. Finally, unequal occurrence between A ↔ G and C ↔ T in five paired neighboring structures indicates that the methylation-induced deamination reactions were responsible for ∼20% of total transitions. In addition, conversion can occur at both CpG sites and non-CpG sites. Our study provides new insights into understanding molecular mechanisms of mutations and genome evolution.

A Novel Type of Sequence Variation: Multiple-Nucleotide Length Polymorphisms Discovered in the Bovine Genome

Three types of sequence variations—single-nucleotide polymorphisms (SNPs), insertions and deletions (indels), and short tandem repeats (STRs)—have been extensively reported in mammalian genomes. In this study, we discovered a novel type of sequence variation, i.e., multiple-nucleotide length polymorphisms (MNLPs) in bovine UCN3 (Urocortin 3) and its receptor CRHR2 (corticotropin-releasing hormone receptor 2) genes. Both MNLPs featured involvement of multiple-nucleotide length polymorphisms (5–18 bases), low sequence identity, and 1.7- to 11-fold changes in promoter activity between two alleles. Therefore, this novel genetic complexity would contribute significantly to the evolutionary, functional, and phenotypic complexity of genomes within or among species.

Accurate Profiling of Gene Expression and Alternative Polyadenylation with Whole Transcriptome Termini Site Sequencing (WTTS-Seq).

Construction of next-generation sequencing (NGS) libraries involves RNA manipulation, which often creates noisy, biased, and artifactual data that contribute to errors in transcriptome analysis. In this study, a total of 19 whole transcriptome termini site sequencing (WTTS-seq) and seven RNA sequencing (RNA-seq) libraries were prepared from Xenopus tropicalis adult and embryo samples to determine the most effective library preparation method to maximize transcriptomics investigation. We strongly suggest that appropriate primers/adaptors are designed to inhibit amplification detours and that PCR overamplification is minimized to maximize transcriptome coverage. Furthermore, genome annotation must be improved so that missing data can be recovered. In addition, a complete understanding of sequencing platforms is critical to limit the formation of false-positive results. Technically, the WTTS-seq method enriches both poly(A)+ RNA and complementary DNA, adds 5′- and 3′-adaptors in one step, pursues strand sequencing and mapping, and profiles both gene expression and alternative polyadenylation (APA). Although RNA-seq is cost prohibitive, tends to produce false-positive results, and fails to detect APA diversity and dynamics, its combination with WTTS-seq is necessary to validate transcriptome-wide APA.

Improving accuracy of genomic prediction in Brangus cattle by adding animals with imputed low-density SNP genotypes

Reliable genomic prediction of breeding values for quantitative traits requires the availability of sufficient number of animals with genotypes and phenotypes in the training set. As of 31 October 2016, there were 3,797 Brangus animals with genotypes and phenotypes. These Brangus animals were genotyped using different commercial SNP chips. Of them, the largest group consisted of 1,535 animals genotyped by the GGP-LDV4 SNP chip. The remaining 2,262 genotypes were imputed to the SNP content of the GGP-LDV4 chip, so that the number of animals available for training the genomic prediction models was more than doubled. The present study showed that the pooling of animals with both original or imputed 40K SNP genotypes substantially increased genomic prediction accuracies on the ten traits. By supplementing imputed genotypes, the relative gains in genomic prediction accuracies on estimated breeding values (EBV) were from 12.60% to 31.27%, and the relative gain in genomic prediction accuracies on de-regressed EBV was slightly small (i.e. 0.87%-18.75%). The present study also compared the performance of five genomic prediction models and two cross-validation methods. The five genomic models predicted EBV and de-regressed EBV of the ten traits similarly well. Of the two cross-validation methods, leave-one-out cross-validation maximized the number of animals at the stage of training for genomic prediction. Genomic prediction accuracy (GPA) on the ten quantitative traits was validated in 1,106 newly genotyped Brangus animals based on the SNP effects estimated in the previous set of 3,797 Brangus animals, and they were slightly lower than GPA in the original data. The present study was the first to leverage currently available genotype and phenotype resources in order to harness genomic prediction in Brangus beef cattle.