Genetic analysis of the blood transcriptome of young healthy pigs to improve disease resilience.

BACKGROUND: Disease resilience is the ability of an animal to maintain productive performance under disease conditions and is an important selection target. In pig breeding programs, disease resilience must be evaluated on selection candidates without exposing them to disease. To identify potential genetic indicators for disease resilience that can be measured on selection candidates, we focused on the blood transcriptome of 1594 young healthy pigs with subsequent records on disease resilience. Transcriptome data were obtained by 3'mRNA sequencing and genotype data were from a 650 K genotyping array. RESULTS: Heritabilities of the expression of 16,545 genes were estimated, of which 5665 genes showed significant estimates of heritability (p < 0.05), ranging from 0.05 to 0.90, with or without accounting for white blood cell composition. Genes with heritable expression levels were spread across chromosomes, but were enriched in the swine leukocyte antigen region (average estimate > 0.2). The correlation of heritability estimates with the corresponding estimates obtained for genes expressed in human blood was weak but a sizable number of genes with heritable expression levels overlapped. Genes with heritable expression levels were significantly enriched for biological processes such as cell activation, immune system process, stress response, and leukocyte activation, and were involved in various disease annotations such as RNA virus infection, including SARS-Cov2, as well as liver disease, and inflammation. To estimate genetic correlations with disease resilience, 3205 genotyped pigs, including the 1594 pigs with transcriptome data, were evaluated for disease resilience following their exposure to a natural polymicrobial disease challenge. Significant genetic correlations (p < 0.05) were observed with all resilience phenotypes, although few exceeded expected false discovery rates. Enrichment analysis of genes ranked by estimates of genetic correlations with resilience phenotypes revealed significance for biological processes such as regulation of cytokines, including interleukins and interferons, and chaperone mediated protein folding. CONCLUSIONS: These results suggest that expression levels in the blood of young healthy pigs for genes in biological pathways related to immunity and endoplasmic reticulum stress have potential to be used as genetic indicator traits to select for disease resilience.

Adjusting for gene-specific covariates to improve RNA-seq analysis.

SUMMARY: This article suggests a novel positive false discovery rate (pFDR) controlling method for testing gene-specific hypotheses using a gene-specific covariate variable, such as gene length. We suppose the null probability depends on the covariate variable. In this context, we propose a rejection rule that accounts for heterogeneity among tests by using two distinct types of null probabilities. We establish a pFDR estimator for a given rejection rule by following Storey's q-value framework. A condition on a type 1 error posterior probability is provided that equivalently characterizes our rejection rule. We also present a suitable procedure for selecting a tuning parameter through cross-validation that maximizes the expected number of hypotheses declared significant. A simulation study demonstrates that our method is comparable to or better than existing methods across realistic scenarios. In data analysis, we find support for our method's premise that the null probability varies with a gene-specific covariate variable. AVAILABILITY AND IMPLEMENTATION: The source code repository is publicly available at https://github.com/hsjeon1217/conditional_method.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils.

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Plasma protein levels of young healthy pigs as indicators of disease resilience.

Selection for disease resilience, which refers to the ability of an animal to maintain performance when exposed to disease, can reduce the impact of infectious diseases. However, direct selection for disease resilience is challenging because nucleus herds must maintain a high health status. A possible solution is indirect selection of indicators of disease resilience. To search for such indicators, we conducted phenotypic and genetic quantitative analyses of the abundances of 377 proteins in plasma samples from 912 young and visually healthy pigs and their relationships with performance and subsequent disease resilience after natural exposure to a polymicrobial disease challenge. Abundances of 100 proteins were significantly heritable (false discovery rate (FDR) <0.10). The abundance of some proteins was or tended to be genetically correlated (rg) with disease resilience, including complement system proteins (rg = -0.24, FDR = 0.001) and IgG heavy chain proteins (rg = -0.68, FDR = 0.22). Gene set enrichment analyses (FDR < 0.2) based on phenotypic and genetic associations of protein abundances with subsequent disease resilience revealed many pathways related to the immune system that were unfavorably associated with subsequent disease resilience, especially the innate immune system. It was not possible to determine whether the observed levels of these proteins reflected baseline levels in these young and visually healthy pigs or were the result of a response to environmental disturbances that the pigs were exposed to before sample collection. Nevertheless, results show that, under these conditions, the abundance of proteins in some immune-related pathways can be used as phenotypic and genetic predictors of disease resilience and have the potential for use in pig breeding and management.

A challenge of selection for disease resilience is that it is difficult to directly select pigs that have greater resilience to multiple diseases in the healthy nucleus herd environment which is essential for breeding programs. A possible alternative is to select an indicator trait or marker that can be measured in a healthy setting, is heritable, and is associated with the genetics of disease resilience. In this study, we investigated plasma protein levels measured on young healthy pigs as indicator traits to select for disease resilience. For this purpose, we used plasma proteome data collected prior to the natural exposure of nursery pigs to multiple diseases, performed phenotypic and genetic quantitative analyses, and investigated their relationships with disease resilience. Our results suggest that plasma protein levels of young healthy pigs have the potential as biomarkers to select for disease resistance.

Integrative methylome and transcriptome analysis of porcine abdominal fat indicates changes in fat metabolism and immune responses during different development.

Fat is involved in synthesizing fatty acids (FAs), FA circulation, and lipid metabolism. Various genetic studies have been conducted on porcine fat but understanding the growth and specific adipose tissue is insufficient. The purpose of this study is to investigate the epigenetic difference in abdominal fat according to the growth of porcine. The samples were collected from the porcine abdominal fat of different developmental stages (10 and 26 weeks of age). Then, the samples were sequenced using MBD-seq and RNA-seq for profiling DNA methylation and RNA expression. In 26 weeks of age pigs, differentially methylated genes (DMGs) and differentially expressed genes (DEGs) were identified as 2,251 and 5,768, compared with 10 weeks of age pigs, respectively. Gene functional analysis was performed using GO and KEGG databases. In functional analysis results of DMGs and DEGs, immune responses such as chemokine signaling pathways, B cell receptor signaling pathways, and lipid metabolism terms such as PPAR signaling pathways and fatty acid degradation were identified. It is thought that there is an influence between DNA methylation and gene expression through changes in genes with similar functions. The effects of DNA methylation on gene expression were investigated using cis-regulation and trans-regulation analysis to integrate and interpret different molecular layers. In the cis-regulation analysis using 629 overlapping genes between DEGs and DMGs, immune response functions were identified, while in trans-regulation analysis through the TF-target gene network, the co-expression network of lipid metabolism-related functions was distinguished. Our research provides an understanding of the underlying mechanisms for epigenetic regulation in porcine abdominal fat with aging.

Fat is involved in the synthesis of new fatty acids (FAs), FA circulation, and lipid metabolism. Various genetic studies have been conducted on porcine fat but understanding the growth and specific adipose tissue is insufficient. The purpose of this study is to investigate the epigenetic difference in abdominal fat according to the growth of porcine. Modifications in DNA methylation and expression values were confirmed epigenetically with growth. Changed genes in each DNA and RNA showed identical trends in the function of immune response and lipid metabolism. The effects of DNA methylation on gene expression were investigated using cis-regulation (functional enrichment analysis of overlapping genes) and trans-regulation (transcription factor and target gene networking) analysis to integrate and interpret different molecular layers. Our research provides an understanding of the underlying mechanisms for epigenetic regulation in porcine abdominal fat with aging.

Regulation of Swine Growth by Backfat Tissue during Growing and Finishing Stages.

Recently, interest in the function of pig backfat (BF) has increased in the field of livestock animals, and many transcriptome-based studies using commercial pig breeds have been conducted. However, there is a lack of comprehensive studies regarding the biological mechanisms of Korean native pigs (KNPs) and Yorkshire pig crossbreeds. In this study, therefore, BF samples of F1 crossbreeds of KNPs and Yorkshire pigs were investigated to identify differentially expressed genes (DEGs) and their related terms using RNA-sequencing analysis. DEG analysis identified 611 DEGs, of which 182 were up-regulated and 429 were down-regulated. Lipid metabolism was identified in the up-regulated genes, whereas growth and maturation-related terminologies were identified in the down-regulated genes. LEP and ACTC1 were identified as highly connected core genes during functional gene network analysis. Fat tissue was observed to affect lipid metabolism and organ development due to hormonal changes driven by transcriptional alteration. This study provides a comprehensive understanding of BF contribution to crossbreeds of KNPs and Yorkshire pigs during growth periods.

Genetic association of polymorphisms in porcine RGS16 with porcine circovirus viral load in naturally infected Yorkshire pigs.

Regulator of G protein signaling 16 (RGS16) is known to be associated with porcine circovirus type 2 (PCV2). PCV2 associated disease (PCVAD) is a serious problem in the swine industry. The representative symptoms of PCVAD are high viral titer proliferation and decreased average daily gain. In this study, we identified single nucleotide polymorphisms (SNPs) in the RGS16 region, including the upstream region. Of the 22 identified SNPs, rs332913874, rs326071195, and rs318298586 were genotyped in 142 Yorkshire pigs. These SNPs were significantly associated with the PCV2 viral load. Moreover, the haplotype combination was also related to the PCV2 viral load. The haplotype and diplotype analysis also had a significant difference with the PCV2 viral load. Taken together, our results suggest that RGS16 SNPs considerably affect the PCV2 viral load.

Associations of natural variation in the CD163 and other candidate genes on host response of nursery pigs to porcine reproductive and respiratory syndrome virus infection.

Pigs with complete resistance to porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) have been produced by genetically knocking out the CD163 gene that encodes a receptor of the PRRSV for entry into macrophages. The objectives of this study were to evaluate associations of naturally occurring single nucleotide polymorphisms (SNPs) in the CD163 gene and in three other candidate genes (CD169, RGS16, and TRAF1) with host response to PRRSV-only infection and to PRRS vaccination and PRRSV/porcine circovirus 2b (PCV2b) coinfection. SNPs in the CD163 gene were not included on SNP genotyping panels that were used for previous genome-wide association analyses of these data. An additional objective was to identify the potential genetic interaction of variants at these four candidate genes with a mutation in the GBP5 gene that was previously identified to be associated with host response to PRRSV infection. Finally, the association of SNPs with expression level of the nearby gene was tested. Several SNPs in the CD163, CD169, and RGS16 genes were significantly associated with host response under PRRSV-only and/or PRRSV/PCV2b coinfection. The effects of all SNPs that were significant in the PRRSV-only infection trials depend on genetic background. The effects of some SNPs in the CD163, CD169, and RGS16 genes depend on genotype at the putative causative mutation in the GBP5 gene, which indicates a potential biological interaction of these genes with GBP5. In addition, genome-wide association results for the PRRSV-only infection trials revealed that SNPs located in the CDK5RAP2 or MEGF9 genes, near the TRAF1 gene, had suggestive effects on PRRS viral load, which indicates that these SNPs might contribute to PRRSV neuropathogenesis. In conclusion, natural genetic variants in the CD163, CD169, and RGS16 genes are associated with resistance to PRRSV and/or PCV2b infection and appear to interact with the resistance quantitative trait locus in the GBP5 gene. The identified SNPs can be used to select for increased natural resistance to PRRSV and/or PRRSV-PCV2b coinfection.

Quantitative analysis of the blood transcriptome of young healthy pigs and its relationship with subsequent disease resilience.

BACKGROUND: Disease resilience, which is the ability of an animal to maintain performance under disease, is important for pigs in commercial herds, where they are exposed to various pathogens. Our objective was to investigate population-level gene expression profiles in the blood of 912 healthy F1 barrows at ~ 27 days of age for associations with performance and health before and after their exposure to a natural polymicrobial disease challenge at ~ 43 days of age. RESULTS: Most significant (q < 0.20) associations of the level of expression of individual genes in blood of young healthy pigs were identified for concurrent growth rate and subjective health scores prior to the challenge, and for mortality, a combined mortality-treatment trait, and feed conversion rate after the challenge. Gene set enrichment analyses revealed three groups of gene ontology biological process terms that were related to disease resilience: 1) immune and stress response-related terms were enriched among genes whose increased expression was unfavorably associated with both pre- and post-challenge traits, 2) heme-related terms were enriched among genes that had favorable associations with both pre- and post-challenge traits, and 3) terms related to protein localization and viral gene expression were enriched among genes that were associated with reduced performance and health traits after but not before the challenge. CONCLUSIONS: Gene expression profiles in blood from young healthy piglets provide insight into their performance when exposed to disease and other stressors. The expression of genes involved in stress response, heme metabolism, and baseline expression of host genes related to virus propagation were found to be associated with host response to disease.

Gene expression in tonsils in swine following infection with porcine reproductive and respiratory syndrome virus.

BACKGROUND: Porcine reproductive and respiratory syndrome (PRRS) is a threat to pig production worldwide. Our objective was to understand mechanisms of persistence of PRRS virus (PRRSV) in tonsil. Transcriptome data from tonsil samples collected at 42 days post infection (dpi) were generated by RNA-seq and NanoString on 51 pigs that were selected to contrast the two PRRSV isolates used, NVSL and KS06, high and low tonsil viral level at 42 dpi, and the favorable and unfavorable genotypes at a genetic marker (WUR) for the putative PRRSV resistance gene GBP5. RESULTS: The number of differentially expressed genes (DEGs) differed markedly between models with and without accounting for cell-type enrichments (CE) in the samples that were predicted from the RNA-seq data. This indicates that differences in cell composition in tissues that consist of multiple cell types, such as tonsil, can have a large impact on observed differences in gene expression. Based on both the NanoString and the RNA-seq data, KS06-infected pigs showed greater activation, or less inhibition, of immune response in tonsils at 42 dpi than NVSL-infected pigs, with and without accounting for CE. This suggests that the NVSL virus may be better than the KS06 virus at evading host immune response and persists in tonsils by weakening, or preventing, host immune responses. Pigs with high viral levels showed larger CE of immune cells than low viral level pigs, potentially to trigger stronger immune responses. Presence of high tonsil virus was associated with a stronger immune response, especially innate immune response through interferon signaling, but these differences were not significant when accounting for CE. Genotype at WUR was associated with different effects on immune response in tonsils of pigs during the persistence stage, depending on viral isolate and tonsil viral level. CONCLUSIONS: Results of this study provide insights into the effects of PRRSV isolate, tonsil viral level, and WUR genotype on host immune response and into potential mechanisms of PRRSV persistence in tonsils that could be targeted to improve strategies to reduce viral rebreaks. Finally, to understand transcriptome responses in tissues that consist of multiple cell types, it is important to consider differences in cell composition.

Identification of Monoallelically Expressed Genes Associated with Economic Traits in Hanwoo (Korean Native Cattle).

Hanwoo, an indigenous Korean cattle breed, has been genetically improved by selecting superior sires called Korean-proven bulls. However, cows still contribute half of the genetic stock of their offspring, and allelic-specific expressed genes have potential, as selective targets of cows, to enhance genetic gain. The aim of this study is to identify genes that have MAEs based on both the genome and transcriptome and to estimate their effects on breeding values (BVs) for economically important traits in Hanwoo. We generated resequencing data for the parents and RNA-sequencing data for the muscle, fat, and brain tissues of the offspring. A total of 3801 heterozygous single nucleotide polymorphisms (SNPs) in offspring were identified and they were located in 1569 genes. Only 14 genes showed MAE (seven expressing maternal alleles and seven expressing paternal alleles). Tissue-specific MAE was observed, and LANCL1 showed maternal allele expression across all tissues. MAE genes were enriched for the biological process of cell death and angiogenesis, which included ACKR3 and PDCL3 genes, whose SNPs were significantly associated with BVs of lean meat production-related traits, such as weight at 12 months of age, carcass weight, and loin eye area. In the current study, monoallelically expressed genes were identified in various adult tissues and these genes were associated with genetic capacity in Hanwoo.

Integrated time-serial transcriptome networks reveal common innate and tissue-specific adaptive immune responses to PRRSV infection.

Porcine reproductive and respiratory syndrome virus (PRRSV) infection is the most important viral disease causing severe economic losses in the swine industry. However, mechanisms underlying gene expression control in immunity-responsible tissues at different time points during PRRSV infection are poorly understood. We constructed an integrated gene co-expression network and identified tissue- and time-dependent biological mechanisms of PRRSV infection through bioinformatics analysis using three tissues (lungs, bronchial lymph nodes [BLNs], and tonsils) via RNA-Seq. Three groups with specific expression patterns (i.e., the 3-dpi, lung, and BLN groups) were discovered. The 3 dpi-specific group showed antiviral and innate-immune signalling similar to the case for influenza A infection. Moreover, we observed adaptive immune responses in the lung-specific group based on various cytokines, while the BLN-specific group showed down-regulated AMPK signalling related to viral replication. Our study may provide comprehensive insights into PRRSV infection, as well as useful information for vaccine development.

The genetic basis of natural antibody titers of young healthy pigs and relationships with disease resilience.

BACKGROUND: Disease resilience is the ability to maintain performance under pathogen exposure but is difficult to select for because breeding populations are raised under high health. Selection for resilience requires a trait that is heritable, easy to measure on healthy animals, and genetically correlated with resilience. Natural antibodies (NAb) are important parts of the innate immune system and are found to be heritable and associated with disease susceptibility in dairy cattle and poultry. Our objective was to investigate NAb and total IgG in blood of healthy, young pigs as potential indicator traits for disease resilience. RESULTS: Data were from Yorkshire x Landrace pigs, with IgG and IgM NAb (four antigens) and total IgG measured by ELISA in blood plasma collected ~ 1 week after weaning, prior to their exposure to a natural polymicrobial challenge. Heritability estimates were lower for IgG NAb (0.12 to 0.24, + 0.05) and for total IgG (0.19 + 0.05) than for IgM NAb (0.33 to 0.53, + 0.07) but maternal effects were larger for IgG NAb (0.41 to 0.52, + 0.03) and for total IgG (0.19 + 0.05) than for IgM NAb (0.00 to 0.10, + 0.04). Phenotypically, IgM NAb titers were moderately correlated with each other (average 0.60), as were IgG NAb titers (average 0.42), but correlations between IgM and IgG NAb titers were weak (average 0.09). Phenotypic correlations of total IgG were moderate with NAb IgG (average 0.46) but weak with NAb IgM (average 0.01). Estimates of genetic correlations among NAb showed similar patterns but with small SE, with estimates averaging 0.76 among IgG NAb, 0.63 among IgM NAb, 0.17 between IgG and IgM NAb, 0.64 between total IgG and IgG NAb, and 0.13 between total IgG and IgM NAb. Phenotypically, pigs that survived had slightly higher levels of NAb and total IgG than pigs that died. Genetically, higher levels of NAb tended to be associated with greater disease resilience based on lower mortality and fewer parenteral antibiotic treatments. Genome-wide association analyses for NAb titers identified several genomic regions, with several candidate genes for immune response. CONCLUSIONS: Levels of NAb in blood of healthy young piglets are heritable and potential genetic indicators of resilience to polymicrobial disease.

The effects of a globin blocker on the resolution of 3'mRNA sequencing data in porcine blood.

BACKGROUND: Gene expression profiling in blood is a potential source of biomarkers to evaluate or predict phenotypic differences between pigs but is expensive and inefficient because of the high abundance of globin mRNA in porcine blood. These limitations can be overcome by the use of QuantSeq 3'mRNA sequencing (QuantSeq) combined with a method to deplete or block the processing of globin mRNA prior to or during library construction. Here, we validated the effectiveness of QuantSeq using a novel specific globin blocker (GB) that is included in the library preparation step of QuantSeq. RESULTS: In data set 1, four concentrations of the GB were applied to RNA samples from two pigs. The GB significantly reduced the proportion of globin reads compared to non-GB (NGB) samples (P = 0.005) and increased the number of detectable non-globin genes. The highest evaluated concentration (C1) of the GB resulted in the largest reduction of globin reads compared to the NGB (from 56.4 to 10.1%). The second highest concentration C2, which showed very similar globin depletion rates (12%) as C1 but a better correlation of the expression of non-globin genes between NGB and GB (r = 0.98), allowed the expression of an additional 1295 non-globin genes to be detected, although 40 genes that were detected in the NGB sample (at a low level) were not present in the GB library. Concentration C2 was applied in the rest of the study. In data set 2, the distribution of the percentage of globin reads for NGB (n = 184) and GB (n = 189) samples clearly showed the effects of the GB on reducing globin reads, in particular for HBB, similar to results from data set 1. Data set 3 (n = 84) revealed that the proportion of globin reads that remained in GB samples was significantly and positively correlated with the reticulocyte count in the original blood sample (P < 0.001). CONCLUSIONS: The effect of the GB on reducing the proportion of globin reads in porcine blood QuantSeq was demonstrated in three data sets. In addition to increasing the efficiency of sequencing non-globin mRNA, the GB for QuantSeq has an advantage that it does not require an additional step prior to or during library creation. Therefore, the GB is a useful tool in the quantification of whole gene expression profiles in porcine blood.

Genome-Wide Analysis of Allele-Specific Expression Patterns in Seventeen Tissues of Korean Cattle (Hanwoo).

The functional hemizygosity could be caused by the MAE of a given gene and it can be one of the sources to affect the phenotypic variation in cattle. We aimed to identify MAE genes across the transcriptome in Korean cattle (Hanwoo). For three Hanwoo family trios, the transcriptome data of 17 tissues were generated in three offspring. Sixty-two MAE genes had a monoallelic expression in at least one tissue. Comparing genotypes among each family trio, the preferred alleles of 18 genes were identified (maternal expression, n = 9; paternal expression, n = 9). The MAE genes are involved in gene regulation, metabolic processes, and immune responses, and in particular, six genes encode transcription factors (FOXD2, FOXM1, HTATSF1, SCRT1, NKX6-2, and UBN1) with tissue-specific expression. In this study, we report genome-wide MAE genes in seventeen tissues of adult cattle. These results could help to elucidate epigenetic effects on phenotypic variation in Hanwoo.

Estimation of pork quality in live pigs using biopsied muscle fibre number composition.

Here, we newly provided the parameters for estimating meat quality in live pigs using the muscle biopsy. The biopsied longissimus thoracis muscle was used to identify the muscle fibre characteristics (MFCs). Of the various MFCs in biopsied muscle, muscle fibre number (MFN) composition showed the greatest correlation with the MFCs in postmortem muscle (P<0.001). Moreover, the pigs cluster groups, based on the biopsied MFN composition, demonstrated statistically significant differences in meat quality traits such as muscle pH, drip loss, and meat colour (P<0.05). Therefore, we conclude that the MFN parameters in live pigs are closely related to the postmortem muscle metabolic rate and ultimately with the quality of meat. We suggest that the higher type I and lower type IIB MFN in biopsied muscle will result in better pork quality.

Identification of the acclimation genes in transcriptomic responses to heat stress of White Pekin duck.

White Pekin duck is an important meat resource in the livestock industries. However, the temperature increase due to global warming has become a serious environmental factor in duck production, because of hyperthermia. Therefore, identifying the gene regulations and understanding the molecular mechanism for adaptation to the warmer environment will provide insightful information on the acclimation system of ducks. This study examined transcriptomic responses to heat stress treatments (3 and 6 h at 35 °C) and control (C, 25 °C) using RNA-sequencing analysis of genes from the breast muscle tissue. Based on three distinct differentially expressed gene (DEG) sets (3H/C, 6H/C, and 6H/3H), the expression patterns of significant DEGs (absolute log2 > 1.0 and false discovery rate < 0.05) were clustered into three responsive gene groups divided into upregulated and downregulated genes. Next, we analyzed the clusters that showed relatively higher expression levels in 3H/C and lower levels in 6H/C with much lower or opposite levels in 6H/3H; we referred to these clusters as the adaptable responsive gene group. These genes were significantly enriched in the ErbB signaling pathway, neuroactive ligand-receptor interaction and type II diabetes mellitus in the KEGG pathways (P < 0.01). From the functional enrichment analysis and significantly regulated genes observed in the enriched pathways, we think that the adaptable responsive genes are responsible for the acclimation mechanism of ducks and suggest that the regulation of phosphoinositide 3-kinase genes including PIK3R6, PIK3R5, and PIK3C2B has an important relationship with the mechanisms of adaptation to heat stress in ducks.

Genomic structure, expression and association study of the porcine FSD2.

The fibronectin type III and SPRY domain containing 2 (FSD2) on porcine chromosome 7 is considered a candidate gene for pork quality, since its two domains, which were present in fibronectin and ryanodine receptor. The fibronectin type III and SPRY domains were first identified in fibronectin and ryanodine receptor, respectively, which are candidate genes for meat quality. The aim of this study was to elucidate the genomic structure of FSD2 and functions of single nucleotide polymorphisms (SNPs) within FSD2 that are related to meat quality in pigs. Using a bacterial artificial chromosome clone sequence, we revealed that porcine FSD2 consisted of 13 exons encoding 750 amino acids. In addition, FSD2 was expressed in heart, longissimus dorsi muscle, psoas muscle, and tendon among 23 kinds of porcine tissues tested. A total of ten SNPs, including four missense mutations, were identified in the exonic region of FSD2, and two major haplotypes were obtained based on the SNP genotypes of 633 Berkshire pigs. Both haplotypes were associated significantly with intramuscular fat content (IMF, P < 0.020) and moisture percentage (MP, P < 0.002). Moreover, haplotype 2 was associated with meat color, affecting yellowness (P = 0.002). These haplotype effects were further supported by the alteration of putative protein structures with amino acid substitutions. Taken together, our results suggest that FSD2 haplotypes are involved in regulating meat quality including IMF, MP, and meat color in pigs, and may be used as meaningful molecular makers to identify pigs with preferable pork quality.

Copy number variations in Hanwoo and Yanbian cattle genomes using the massively parallel sequencing data.

Hanwoo is an indigenous Korean beef cattle breed, and it shared an ancestor with Yanbian cattle that are found in the Northeast provinces in China until the last century. During recent decades, those cattle breeds experienced different selection pressures. Here, we present genome-wide copy number variations (CNVs) by comparing Hanwoo and Yanbian cattle sequencing data. We used ~3.12 and ~3.07 billion sequence reads from Hanwoo and Yanbian cattle, respectively. A total of 901 putative CNV regions (CNVRs) were identified throughout the genome, representing 5,513,340bp. This is a smaller number than has been reported in previous studies, indicating that Hanwoo are genetically close to Yanbian cattle. Of the CNVRs, 53.2% and 46.8% were found to be gains and losses in Hanwoo. Potential functional roles of each CNVR were assessed by annotating all CNVRs and gene ontology (GO) enrichment analysis. We found that 278 CNVRs overlapped with cattle gene-sets (genic-CNVRs) that could be promising candidates to account for economically important traits in cattle. The enrichment analysis indicated that genes were significantly over-represented in GO terms, including developmental process, multicellular organismal process, reproduction, and response to stimulus. These results provide a valuable genomic resource for determining how CNVs are associated with cattle traits.

Effects of intergenic single nucleotide polymorphisms in the fast myosin heavy chain cluster on muscle fiber characteristics and meat quality in Berkshire pigs.

The heterogeneity of myosin heavy chain (MyHC) isoforms is closely related to muscle fiber characteristics, and meat quality in pigs. The fast MyHC gene cluster on porcine chromosome 12 contains MYH2, MYH1, and MYH4, which encode three fast MyHC isoforms expressed in adult skeletal muscle. Here, we identified four intergenic single nucleotide polymorphisms (SNPs) in the gene cluster. Of these, SNP1 and SNP2, which were located near MYH2, were genotyped in 199 Berkshire pigs. SNP1 were significantly associated with the total fiber number (P=0.046) and intramuscular fat contents (P=0.041), and SNP2 had significant effects on type I fiber number/area compositions (P=0.039 and P=0.041, respectively), water holding capacity (drip loss, P=0.045; cooking loss, P=0.001), and meat color (P=0.003). Taken together, our results suggest that the intergenic SNPs near MYH2 can affect muscle fiber formation and meat quality by modifying the regulatory elements of this gene cluster in pigs.