Jianfeng Guo

Transcriptomic Analysis of the Porcine Endometrium during Embryo Implantation

In pigs, successful embryo implantation is an important guarantee for producing litter size, and early embryonic loss occurring on day 12–30 of gestation critically affects the potential litter size. The implantation process is regulated by the expression of numerous genes, so comprehensive analysis of the endometrium is necessary. In this study, RNA sequencing (RNA-Seq) technology is used to analyze endometrial tissues during early pregnancy. We investigated the changes of gene expression between three stages (day 12, 18, and 25) by multiple comparisons. There were 1557, 8951, and 2345 differentially expressed genes (DEGs) revealed between the different periods of implantation. We selected several genes for validation by the use of quantitative real-time RT-PCR. Bioinformatic analysis of differentially expressed genes in the endometrium revealed a number of biological processes and pathways potentially involved in embryo implantation in the pig, most noticeably cell proliferation, regulation of immune response, interaction of cytokine-cytokine receptors, and cell adhesion. These results showed that specific gene expression patterns reflect the different functions of the endometrium in three stages (maternal recognition, conceptus attachment, and embryo implantation). This study identified comprehensive transcriptomic profile in the porcine endometrium and thus could be a foundation for targeted studies of genes and pathways potentially involved in abnormal endometrial receptivity and embryo loss in early pregnancy.

Selection of Reference Genes for Gene Expression Studies in Porcine Whole Blood and Peripheral Blood Mononuclear Cells under Polyinosinic:Polycytidylic Acid Stimulation

Investigating gene expression of immune cells of whole blood or peripheral blood mononuclear cells (PBMC) under polyinosinic:polycytidylic acid (poly I:C) stimulation is valuable for understanding the immune response of organism to RNA viruses. Quantitative real-time PCR (qRT-PCR) is a standard method for quantification of gene expression studies. However, the reliability of qRT-PCR data critically depends on proper selection of reference genes. In the study, using two different analysis programs, geNorm and NormFinder, we systematically evaluated the gene expression stability of six candidate reference genes (GAPDH, ACTB, B2M, RPL4, TBP, and PPIA) in samples of whole blood and PBMC with or without poly I:C stimulation. Generally, the six candidate genes performed a similar trend of expression stability in the samples of whole blood and PBMC, but more stably expressed in whole blood than in PBMC. geNorm ranked B2M and PPIA as the best combination for gene expression normalization, while according to NormFinder, TBP was ranked as the most stable reference gene, followed by B2M and PPIA. Comprehensively considering the results from the two programs, we recommended using the geometric mean of the three genes, TBP, PPIA and B2M, to normalize the gene expression of whole blood and PBMC with poly I:C stimulation. Our study is the first detailed survey of the gene expression stability in whole blood and PBMC with or without poly I:C stimulation and should be helpful for investigating the molecular mechanism involved in porcine whole blood and PBMC in response to poly I:C stimulation.

Phylogenetic relationships of pig breeds from Shandong province of China and their influence by modern commercial breeds by analysis of mitochondrial DNA sequences

Mitochondrial cytochrome b (Cytb) and partial control region (CR) DNA sequence of 165 individuals out of eight pig breeds (Laiwu Black, Dapulian Black, Licha Black, Yimeng Black, Yantai Black, Wulian Black, Luyan White and Changwei White) from Shandong province of China were sequenced to determine phylogenetic relationships and to what extent they were affected by modern commercial breeds (Landrace, Yorkshire and Duroc) in maternal lineage. The eight breeds from Shandong had lower maternal evolutionary divergence than the three commercial breeds. Additionally, they shared the same haplotype (T2) of Cytb. All these results suggest they may have originated from closely related maternal ancestors, if not from a single ancestor. The fact that haplotypes of Dapulian Black and Laiwu Black formed their own subclades, while those of other breeds of Shandong dispersed in the Asian clade of the tree based on CR, indicates that Dapulian Black and Laiwu Black are indigenous breeds with relatively independent maternal lineage, while the others have extensively gene flow with other breeds. Except in Luyan White, no European haplotypes were found in pig breeds from Shandong, demonstrating the maternal introgression from the modern commercial pigs has had no or very little impact on the breeds from Shandong. Comparatively, CR is more suitable to analyze the phylogenetic relationships of closely related groups, whereas the four SNPs of Cytb are useful to determine whether modern commercial breeds (Duroc, Landrace and sometimes Yorkshire) were used as dam in the breeding process of new breeds or synthetic lines in China.

MicroRNA Transcriptome of Poly I:C-Stimulated Peripheral Blood Mononuclear Cells Reveals Evidence for MicroRNAs in Regulating Host Response to RNA Viruses in Pigs

MicroRNAs (miRNAs) are one family of small noncoding RNAs that function to modulate the activity of specific mRNA targets in animals. To understand the role of miRNAs in regulating genes involved in the host immune response to RNA viruses, we profiled and characterized the miRNAs of swine peripheral blood mononuclear cells (PBMC) stimulated with poly I:C, a synthetic dsRNA analog, by miRNA-sequencing (miRNA-seq). We identified a total of 905 miRNAs, of which 503 miRNAs were firstly exploited herein with no annotation in the latest miRBase 21.0. Expression analysis demonstrated that poly I:C stimulation can elicit significantly differentially expressed (DE) miRNAs in Dapulian (n = 20), one Chinese indigenous breed, as well as Landrace (n = 23). By integrating the mRNA expression profiles of the same sample with miRNA profiles, we carried out function analyses of the target genes of these DE miRNAs, with the results indicating that target genes were most enriched in some immune-related pathways and gene ontology (GO) terms, suggesting that DE miRNAs play an important role in the regulation of host to poly I:C stimulation. Furthermore, we also detected 43 and 61 significantly DE miRNAs between the two breeds in the control sample groups and poly I:C stimulation groups, respectively, which may be involved in regulation of the different characteristics of the two breeds. This study describes for the first time the PBMC miRNA transcriptomic response to poly I:C stimulation in pigs, which not only contributes to a broad view of the pig miRNAome but improves our understanding of miRNA function in regulating host immune response to RNA viruses.

Transcriptomic Analysis Identifies Candidate Genes and Gene Sets Controlling the Response of Porcine Peripheral Blood Mononuclear Cells to Poly I:C Stimulation

Polyinosinic-polycytidylic acid (poly I:C), a synthetic dsRNA analog, has been demonstrated to have stimulatory effects similar to viral dsRNA. To gain deep knowledge of the host transcriptional response of pigs to poly I:C stimulation, in the present study, we cultured and stimulated peripheral blood mononuclear cells (PBMC) of piglets of one Chinese indigenous breed (Dapulian) and one modern commercial breed (Landrace) with poly I:C, and compared their transcriptional profiling using RNA-sequencing (RNA-seq). Our results indicated that poly I:C stimulation can elicit significantly differentially expressed (DE) genes in Dapulian (g = 290) as well as Landrace (g = 85). We also performed gene set analysis using the Gene Set Enrichment Analysis (GSEA) package, and identified some significantly enriched gene sets in Dapulian (g = 18) and Landrace (g = 21). Most of the shared DE genes and gene sets were immune-related, and may play crucial rules in the immune response of poly I:C stimulation. In addition, we detected large sets of significantly DE genes and enriched gene sets when comparing the gene expression profile between the two breeds, including control and poly I:C stimulation groups. Besides immune-related functions, some of the DE genes and gene sets between the two breeds were involved in development and growth of various tissues, which may be correlated with the different characteristics of the two breeds. The DE genes and gene sets detected herein provide crucial information towards understanding the immune regulation of antiviral responses, and the molecular mechanisms of different genetic resistance to viral infection, in modern and indigenous pigs.

Effects of L-tryptophan, Fructan, and Casein on Reducing Ammonia, Hydrogen Sulfide, and Skatole in Fermented Swine Manure

The effects of daily dietary Bacillus subtilis (Bs), and adding L-tryptophan, fructan, or casein to fecal fermentation broths were investigated as means to reduce the production of noxious gas during manure fermentation caused by ammonia, hydrogen sulfide (H2S), and 3-methylindole (skatole). Eighty swine (50.0±0.5 kg) were equally apportioned to an experimental group given Bs in daily feed, or a control group without Bs. After 6 weeks, fresh manure was collected from both groups for fermentation studies using a 3×3 orthogonal array, in which tryptophan, casein, and fructan were added at various concentrations. After fermentation, the ammonia, H2S, L-tryptophan, skatole, and microflora were measured. In both groups, L-tryptophan was the principle additive increasing skatole production, with significant correlation (r = 0.9992). L-tryptophan had no effect on the production of ammonia, H2S, or skatole in animals fed Bs. In both groups, fructan was the principle additive that reduced H2S production (r = 0.9981). Fructan and Bs significantly interacted in H2S production (p = 0.014). Casein was the principle additive affecting the concentration of ammonia, only in the control group. Casein and Bs significantly interacted in ammonia production (p = 0.039). The predominant bacteria were Bacillus spp. CWBI B1434 (26%) in the control group, and Streptococcus alactolyticus AF201899 (36%) in the experimental group. In summary, daily dietary Bs reduced ammonia production during fecal fermentation. Lessening L-tryptophan and increasing fructan in the fermentation broth reduced skatole and H2S.

Genome-wide association study for intramuscular fat content in Chinese Lulai black pigs

Objective Intramuscular fat (IMF) content plays an important role in meat quality. Identification of single nucleotide polymorphisms (SNPs) and genes related to pig IMF, especially using pig populations with high IMF content variation, can help to establish novel molecular breeding tools for optimizing IMF in pork and unveil the mechanisms that underlie fat metabolism. Methods We collected muscle samples of 453 Chinese Lulai black pigs, measured IMF content by Soxhlet petroleum-ether extraction method, and genotyped genome-wide SNPs using GeneSeek Genomic Profiler Porcine HD BeadChip. Then a genome-wide association study was performed using a linear mixed model implemented in the GEMMA software. Results A total of 43 SNPs were identified to be significantly associated with IMF content by the cutoff p<0.001. Among these significant SNPs, the greatest number of SNPs (n = 19) were detected on Chr.9, and two linkage disequilibrium blocks were formed among them. Additionally, 17 significant SNPs are mapped to previously reported quantitative trait loci (QTLs) of IMF and confirmed previous QTLs studies. Forty-two annotated genes centering these significant SNPs were obtained from Ensembl database. Overrepresentation test of pathways and gene ontology (GO) terms revealed some enriched reactome pathways and GO terms, which mainly involved regulation of basic material transport, energy metabolic process and signaling pathway. Conclusion These findings improve our understanding of the genetic architecture of IMF content in pork and facilitate the follow-up study of fine-mapping genes that influence fat deposition in muscle.

Correlation analyses of CpG island methylation of cluster of differentiation 4 protein with gene expression and T lymphocyte subpopulation traits

Objective Cluster of differentiation 4 protein (CD4) gene is an important immune related gene which plays a significant role in T cell development and host resistance during viral infection. Methods In order to unravel the relationship of CpG island methylation level of CD4 gene with its gene expression and T lymphocyte subpopulation traits, we used one typical Chinese indigenous breed (Dapulian, DP) and one commercial breed (Landrace), then predicted the CpG island of CD4 gene, determined the methylation status of CpG sites by bisulfite sequencing polymerase chain reaction (BSP), and carried out the correlation analyses of methylation frequencies of CpG sites with mRNA expression and T lymphocyte subpopulation traits. Results There was one CpG island predicted in the upstream −2 kb region and exon one of porcine CD4 gene, which located 333 bp upstream from the start site of gene and contained nine CpG sites. The correlation analysis results indicated that the methylation frequency of CpG_2 significantly correlated with CD4 mRNA expression in the DP and Landrace combined population, though it did not reach significance level in DP and Landrace separately. Additionally, 15 potential binding transcription factors (TFs) were predicted within the CpG island, and one of them (Jumonji) contained CpG_2 site, suggesting that it may influence the CD4 gene expression through the potential binding TFs. We also found methylation frequency of CpG_2 negatively correlated with T lymphocyte subpopulation traits CD4+CD8−CD3−, CD4−CD8+CD3− and CD4+/CD8+, and positively correlated with CD4−CD8+CD3+ and CD4+CD8+CD3+ (for all correlation, p<0.01) in DP and Landrace combined population. Thus, the CpG_2 was a critical methylation site for porcine CD4 gene expression and T lymphocyte subpopulation traits. Conclusion We speculated that increased methylation frequency of CpG_2 may lead to the decreased expression of CD4, which may have some kind of influence on T lymphocyte subpopulation traits and the immunity of DP population.