Identification of Differentially Expressed Genes Associated with Litter Size in Berkshire Pig Placenta

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Identification of miR-192 target genes in porcine endometrial epithelial cells based on miRNA pull-down

INTRODUCTION

MicroRNAs (miRNAs)-a class of small endogenous non-coding RNAs-are widely involved in post-transcriptional gene regulation of numerous physiological processes. High-throughput sequencing revealed that the miR-192 expression level appeared to be significantly higher in the blood exosomes of sows at early gestation than that in non-pregnant sows. Furthermore, miR-192 was hypothesized to have a regulatory role in embryo implantation; however, the target genes involved in exerting the regulatory function of miR-192 required further elucidation.

METHODS

In the present study, potential target genes of miR-192 in porcine endometrial epithelial cells (PEECs) were identified through biotin-labeled miRNA pull-down; functional and pathway enrichment analysis was performed via gene ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment. Bioinformatic analyses were concurrently used to predict the potential target genes associated with sow embryo implantation. In addition, double luciferase reporter vectors, reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), and Western blot were performed to verify the targeting and regulatory roles of the abovementioned target genes.

RESULTS

A total of 1688 differentially expressed mRNAs were identified via miRNA pull-down. Through RT-qPCR, the accuracy of the sequencing data was verified. In the bioinformatics analysis, potential target genes of miR-192 appeared to form a dense inter-regulatory network and regulated multiple signaling pathways, such as metabolic pathways and the PI3K-Akt, MAPKs, and mTOR signaling pathways, that are relevant to the mammalian embryo implantation process. In addition, CSK (C-terminal Src kinase) and YY1 (Yin-Yang-1) were predicted to be potential candidates, and we validated that miR-192 directly targets and suppresses the expression of the CSK and YY1 genes.

CONCLUSION

We screened 1688 potential target genes of miR-192 were screened, and CSK and YY1 were identified as miR-192 target genes. The outcomes of the present study provide novel insights into the regulatory mechanism of porcine embryo implantation and the identification of miRNA target genes.

Placental structure, function, and mitochondrial phenotype relate to fetal size in each fetal sex in mice†

Fetal growth depends on placental function, which requires energy from mitochondria. Here we investigated whether mitochondrial function in the placenta relates to the growth of the lightest and heaviest fetuses of each sex within the litter of mice. Placentas from the lightest and heaviest fetuses were taken to evaluate placenta morphology (stereology), mitochondrial energetics (high-resolution respirometry), mitochondrial regulators, nutrient transporters, hormone handling, and signaling pathways (qPCR and Western blotting). We found that mitochondrial complex I and II oxygen consumption rate was greater for placentas supporting the lightest female fetuses, although placental complex I abundance of the lightest females and complexes III and V of the lightest males were decreased compared to their heaviest counterparts. Expression of mitochondrial biogenesis (Nrf1) and fission (Drp1 and Fis1) genes was lower in the placenta from the lightest females, whilst biogenesis-related gene Tfam was greater in the placenta of the lightest male fetuses. In addition, placental morphology and steroidogenic gene (Cyp17a1 and Cyp11a1) expression were aberrant for the lightest females, but glucose transporter (Slc2a1) expression was lower in only the lightest males versus their heaviest counterparts. Differences in intra-litter placental phenotype were related to changes in the expression of hormone-responsive (androgen receptor) and metabolic signaling (AMPK, AKT, and PPARγ) pathways. Thus, in normal mouse pregnancy, placental structure, function, and mitochondrial phenotype are differentially responsive to the growth of the female and male fetus. This study may inform the design of sex-specific therapies for placental insufficiency and fetal growth abnormalities with life-long benefits for the offspring.

Neuronatin gene expression levels affect foetal growth and development by regulating glucose transport in porcine placenta

Imprinted genes play important regulatory roles in the growth and development of placentas and foetuses during pregnancy. In a previous study, we found that the imprinted gene Neuronatin (NNAT) is involved in foetal development; NNAT expression was significantly lower in the placentas of piglets that died neonatally compared to the placentas of surviving piglets. However, the function and mechanism of NNAT in regulating porcine placental development is still unknown. In this study, we collected the placentas of high- and low-weight foetuses at gestational day (GD 65, 90), (n = 4-5 litters/GD) to investigate the role of NNAT in regulating foetal growth and development. We found that the mRNA and protein levels of NNAT were significantly higher in the placentas of high-weight than low-weight foetuses. We then overexpressed NNAT in porcine placental trophoblast cell lines (pTr2) and demonstrated that NNAT activated the PI3K-AKT pathway, and further promoted the expression of glucose transporter 1 (GLUT1) and increased cellular calcium ion levels, which improved glucose transport in placental trophoblast cells in vitro. To conclude, our study suggests that NNAT expression impacts porcine foetal development by regulating placental glucose transport.

Genetic Variations and Differential DNA Methylation to Face Contrasted Climates in Small Ruminants: An Analysis on Traditionally-Managed Sheep and Goats

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (F_ST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

Liver transcriptome profiling and functional analysis of intrauterine growth restriction (IUGR) piglets reveals a genetic correction and sexual-dimorphic gene expression during postnatal development

BACKGROUND

Intrauterine growth restriction (IUGR) remains a major problem associated with swine production. Thus, understanding the physiological changes of postnatal IUGR piglets would aid in improving growth performance. Moreover, liver metabolism plays an important role in the growth and survival of neonatal piglets.

RESULTS

By profiling the transcriptome of liver samples on postnatal Days 1, 7, and 28, our study focused on characterizing the growth, function, and metabolism in the liver of IUGR neonatal piglets. Our study demonstrates that the livers of IUGR piglets were associated with a series of complications, including inflammatory stress and immune dysregulation; cytoskeleton and membrane structure disorganization; dysregulated transcription events; and abnormal glucocorticoid metabolism. In addition, the abnormal liver function index in the serum [alanine aminotransferase (ALT), aspartate aminotransferase (AST), and total protein (TP)], coupled with hepatic pathological and ultrastructural morphological changes are indicative of liver damage and dysfunction in IUGR piglets. Moreover, these results reveal the sex-biased developmental dynamics between male and female IUGR piglets, and that male IUGR piglets may be more sensitive to disrupted metabolic homeostasis.

CONCLUSIONS

These observations provide a detailed reference for understanding the mechanisms and characterizations of IUGR liver functions, and suggest that the potential strategies for improving the survival and growth performance of IUGR offspring should consider the balance between postnatal catch-up growth and adverse metabolic consequences. In particular, sex-specific intervention strategies should be considered for both female and male IUGR piglets.

RNA-seq reveals downregulated osteochondral genes potentially related to tibia bacterial chondronecrosis with osteomyelitis in broilers

Background

Bacterial chondronecrosis with osteomyelitis (BCO) develops in the growth plate (GP) of the proximal femur and tibia and is initiated by damage to the less mineralized chondrocytes followed by colonization of opportunistic bacteria. This condition affects approximately 1% of all birds housed, being considered one of the major causes of lameness in fast growing broilers. Although several studies have been previously performed aiming to understand its pathogenesis, the molecular mechanisms involved with BCO remains to be elucidated. Therefore, this study aimed to generate a profile of global differential gene expression involved with BCO in the tibia of commercial broilers, through RNA sequencing analysis to identity genes and molecular pathways involved with BCO in chickens.

Results

Our data showed 192 differentially expressed (DE) genes: 63 upregulated and 129 downregulated in the GP of the tibia proximal epiphysis of BCO-affected broilers. Using all DE genes, six Biological Processes (BP) were associated with bone development (connective tissue development, cartilage development, skeletal system development, organ morphogenesis, system development and skeletal system morphogenesis). The analyses of the upregulated genes did not indicate any significant BP (FDR < 0.05). However, with the downregulated genes, the same BP were identified when using all DE genes in the analysis, with a total of 26 coding genes explaining BCO in the tibia: ACAN, ALDH1A2, CDH7, CHAD, CHADL, COL11A1, COMP, CSGALNACT1, CYR61, FRZB, GAL3ST1, HAPLN1, IHH, KIF26B, LECT1, LPPR1, PDE6B, RBP4A, SERINC5, SFRP1, SOX8, SOX9, TENM2, THBS1, UCHL1 and WFIKKN2. In addition, seven transcription factors were also associated to BCO: NFATC2, MAFB, HIF1A-ARNT, EWSR1-FLI1, NFIC, TCF3 and NF-KAPPAB.

Conclusions

Our data show that osteochondral downregulated genes are potential molecular causes of BCO in broilers, and the bacterial process seems to be, in fact, a secondary condition. Sixteen genes responsible for bone and cartilage formation were downregulated in BCO-affected broilers being strong candidate genes to trigger this disorder.

N6-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis

Circular RNAs (circRNAs) have been implicated in cancer progression through largely unknown mechanisms. Herein, we identify an N⁶-methyladenosine (m⁶A) modified circRNA, circNSUN2, frequently upregulated in tumor tissues and serum samples from colorectal carcinoma (CRC) patients with liver metastasis (LM) and predicts poorer patient survival. The upregulated expression of circNSUN2 promotes LM in PDX metastasis models in vivo and accelerates cancer cells invasion in vitro. Importantly, N⁶-methyladenosine modification of circNSUN2 increases export to the cytoplasm. By forming a circNSUN2/IGF2BP2/HMGA2 RNA-protein ternary complex in the cytoplasm, circNSUN2 enhances the stability of HMGA2 mRNA to promote CRC metastasis progression. Clinically, the upregulated expressions of circNSUN2 and HMGA2 are more prevalent in LM tissues than in primary CRC tissues. These findings elucidate that N⁶-methyladenosine modification of circNSUN2 modulates cytoplasmic export and stabilizes HMGA2 to promote CRC LM, and suggest that circNSUN2 could represent a critical prognostic marker and/or therapeutic target for the disease.

Serial gene co-expression network approach to mine biological meanings from integrated transcriptomes of the porcine endometrium during estrous cycle

The estrous cycle is a complex process regulated by several hormones. To understand the dynamic changes in gene expression that takes place in the swine endometrium during the estrous cycle relative to the day of estrus onset, we performed RNA-sequencing analysis on days 0, 3, 6, 9, 12, 15, and 18, resulting in the identification of 4495 differentially expressed genes (DEGs; Q ≤ 0.05 and |log₂FC| ≥ 1) at various phases in the estrous cycle. These DEGs were integrated into multiple gene co-expression networks based on different fold changes and correlation coefficient (R²) thresholds and a suitable network, which included 899 genes (|log2FC| ≥ 2 and R² ≥ 0.99), was identified for downstream analyses based on the biological relevance of the Gene Ontology (GO) terms enriched. The genes in this network were partitioned into 6 clusters based on the expression pattern. Several GO terms including cell cycle, apoptosis, hormone signaling, and lipid biosynthetic process were found to be enriched. Furthermore, we found 15 significant KEGG pathways, including cell adhesion molecules, cytokine-cytokine receptor signaling, steroid biosynthesis, and estrogen signaling pathways. We identified several genes and GO terms to be stage-specific. Moreover, the identified genes and pathways extend our understanding of porcine endometrial regulation during estrous cycle and will serve as a good resource for future studies.

Genetic Association between Swine Leukocyte Antigen Class II Haplotypes and Reproduction Traits in Microminipigs

The effects of swine leukocyte antigen (SLA) molecules on numerous production and reproduction performance traits have been mainly reported as associations with specific SLA haplotypes that were assigned using serological typing methods. In this study, we intended to clarify the association between SLA class II genes and reproductive traits in a highly inbred population of 187 Microminipigs (MMP), that have eight different types of SLA class II haplotypes. In doing so, we compared the reproductive performances, such as fertility index, gestation period, litter size, and number of stillbirth among SLA class II low resolution haplotypes (Lrs) that were assigned by a polymerase chain reaction-sequence specific primers (PCR-SSP) typing method. Only low resolution haplotypes were used in this study because the eight SLA class II high-resolution haplotypes had been assigned to the 14 parents or the progenitors of the highly inbred MMP herd in a previous publication. The fertility index of dams with Lr-0.13 was significantly lower than that of dams with Lr-0.16, Lr-0.17, Lr-0.18, or Lr-0.37. Dams with Lr-0.23 had significantly smaller litter size at birth than those with Lr-0.17, Lr-0.18, or Lr-0.37. Furthermore, litter size at weaning of dams with Lr-0.23 was also significantly smaller than those dams with Lr-0.16, Lr-0.17, Lr-0.18, or Lr-0.37. The small litter size of dams with Lr-0.23 correlated with the smaller body sizes of these MMPs. These results suggest that SLA class II haplotypes are useful differential genetic markers for further haplotypic and epistatic studies of reproductive traits, selective breeding programs, and improvements in the production and reproduction performances of MMPs.

Intra-tumor heterogeneity for endometrial cancer and its clinical significance

BACKGROUND

Management of tumors has become more complex owing to tumor heterogeneity. Fewer studies have been performed on intra-tumor heterogeneity of endometrial cancer (EC) until now. Therefore, it is of great clinical value to explore the intra-tumor heterogeneity of EC based on clinical features and gene expression profiles.

METHODS

A total of 1688 patients with EC were screened and 114 patients were finally selected, including specimens from 84 patients with primary EC without relapse (PE) and the paired metastases (P-M) specimens, as well as specimens from 30 patients with primary EC with relapse (RPE) and the paired relapsed EC (P-RE) specimens. Microarray and RNA-seq were used to detect gene expression of EC samples. Clinicopathological characteristics and molecular data were compared between PE and P-M groups and between RPE and P-RE groups to explore the intra-tumor heterogeneity of EC.

RESULTS

The clinical intra-tumor spatial heterogeneity of pathological type, grade, ER status, and PR status between PE and P-M were 17.9%, 13.1%, 28.6%, and 28.6%, respectively. The clinical intra-tumor spatiotemporal heterogeneity of pathological type, grade, ER status, and PR status between RPE and P-RE were 16.7%, 33.3%, 25.0%, and 37.5%, respectively. Cluster analysis sorts EC samples based on progression type of lesion and their pathological type. There were differentially expressed genes between PE and P-M and between RPE and P-RE, of which gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis were mainly enriched in cell proliferation, the p53 signaling pathway, etc. CONCLUSIONS:: Clinical and molecular data showed that there was spatiotemporal heterogeneity in intra-tumor of EC, which may add to the complexity of diagnosis and therapeutics for EC. Considering the intra-tumor heterogeneity, sequential chemotherapy and precision medicine may be a more suitable treatment plan for EC.

Predicting prognosis of endometrioid endometrial adenocarcinoma on the basis of gene expression and clinical features using Random Forest

Traditional clinical features are not sufficient to accurately judge the prognosis of endometrioid endometrial adenocarcinoma (EEA). Molecular biological characteristics and traditional clinical features are particularly important in the prognosis of EEA. The aim of the present study was to establish a predictive model that considers genes and clinical features for the prognosis of EEA. The clinical and RNA sequencing expression data of EEA were derived from samples from The Cancer Genome Atlas (TCGA) and Peking University People's Hospital (PKUPH; Beijing, China). Samples from TCGA were used as the training set, and samples from the PKUPH were used as the testing set. Variable selection using Random Forests (VSURF) was used to select the genes and clinical features on the basis of TCGA samples. The RF classification method was used to establish the prediction model. Kaplan-Meier curves were tested with the log-rank test. The results from this study demonstrated that on the basis of TCGA samples, 11 genes and the grade were selected as the input features. In the training set, the out-of-bag (OOB) error of RF model-1, which was established using the '11 genes', was 0.15; the OOB error of RF model-2, which was established using the 'grade', was 0.39; and the OOB error of RF model-3, established using the '11 genes and grade', was 0.15. In the testing set, the classification accuracy of RF model-1, model-2 and model-3 was 71.43, 66.67 and 80.95%, respectively. In conclusion, to the best of our knowledge, the VSURF was used to select features relevant to EEA prognosis, and an EEA predictive model combining genes and traditional features was established for the first time in the present study. The prediction accuracy of the RF model on the basis of the 11 genes and grade was markedly higher than that of the RF models established by either the 11 genes or grade alone.

Characteristics of Pig Carcass and Primal Cuts Measured by the Autofom Ⅲ Depend on Seasonal Classification

The objective of this study was to investigate slaughtering performance, carcass grade, and quantitative traits of cuts according to seasonal influence by each month in pigs slaughtered in livestock processing complex (LPC) slaughterhouse in Korea, 2017. A total of 267,990 LYD (Landrace×Yorkshire×Duroc) pig data were used in this study. Results of slaughter heads, sex distribution, carcass weight, backfat thickness, grading class, total weight, and fat and lean meat percentages of each cut predicted by Autofom Ⅲ were obtained each month. The number of slaughtered pigs was the highest in early and late fall but the lowest in midsummer. Only in midsummer that the number of females was higher than that of castrates. During 2017, carcass weight was the lowest in late summer. Backfat thickness was in the range of 21–22 mm. In mid and late spring, pigs showed high 1+ grade ratio (37.05% and 36.15%, respectively). For traits of 11 cuts predicted by Autofom Ⅲ, porkbelly showed lower total weight, lean weight, and fat weight in midsummer to early fall but higher lean meat percentage compared to other seasons. Weights of deboned neck, loin, and lean meat were the highest in midfall compared to other seasons (p<0.05). In conclusion, characteristics of slaughtering, grading, and economic traits of pigs seemed to be highly seasonal. They were influenced by seasons. Results of this study could be used as basic data to develop seasonal specified management ways to improve pork production.

An Epigenome-Wide DNA Methylation Map of Testis in Pigs for Study of Complex Traits

Epigenetic changes are important for understanding complex trait variation and inheritance in pigs that are also a valuable biomedical model for human health research. Testis is the main organ for reproduction and boar taint in pigs; however, there have been no studies to-date on adult pig testis epigenome. The main objective of this study was to establish a genome-wide DNA methylation map of pig testis that would help identify candidate epigenetic biomarkers and methylated genes for complex traits such as male reproduction, fertility or boar taint. Reduced Representation Bisulfite Sequencing (RRBS) was used to study methylation levels of cytosine in nine pig testis samples. The results showed that genome-wide methylation status of nine samples overlapped greatly and their variation among pigs were low. The methylation levels of promoter, exon, intron, cytosine and guanine dinucleotide (CpG) islands and CpG island shores regions were 0.15, 0.47, 0.55, 0.39, and 0.53, respectively. Cytosines binding to CpG islands showed different methylation levels between exon and intron regions. All methylation levels of CpG islands were lower than CpG island shores in different genic features. The distribution of 12,738 differentially methylated cytosines (DMCs) within CpG islands, CpG island shores and other regions was 36.86, 21.65, and 41.49%, respectively, and was 0.33, 1.71, 5.95, and 92.01% in promoter, exon, intron and intergenic regions, respectively. Methylation levels of DMCs in promoter, exon and intron regions were significantly different between CpG islands and CpG island shores (P < 0.05). A total of 898 genes with 2089 DMCs were enriched in 112 Gene Ontology (GO) terms. Fifteen methylated genes from our study were associated with fertility or boar taint traits. Our analysis revealed the methylation patterns in different genic features and CpG island regions of testis in pigs, and summarized several candidate genes associated with DMCs and the involved GO terms. These findings are helpful to understand the relationship between DNA methylation and genic CpG islands, to provide candidate epigenetic regions or biomarkers for pig production and welfare and for translational epigenomic studies that use pigs as an animal model for human research.

Analysis of differential gene expression of the transgenic pig with overexpression of PGC1α in muscle

In order to better understand the key regulatory mechanisms of PGC1α in muscle fiber type transition, the RNA-seq was used to compare the change of gene expression in gastrocnemius muscles between wild type pigs and transgenic pigs with overexpression of PGC1α gene in muscle. 371 differentially expressed genes (P ≤ 0.05 and Ratio ≥ 2), including 184 up-regulated genes and 187 down-regulated genes, were identified. Five main signaling pathways including metabolic pathways, ECM-receptor interaction, PPAR signaling pathway, adipocytokine signaling pathway and insulin signaling pathway, were authenticated using KEGG pathway analysis. Our results indicate that the fat metabolism pathway plays an important role in the transformation of muscle fiber types regulated by PGC1α.

The identification of differentially expressed genes between extremes of placental efficiency in maternal line gilts on day 95 of gestation

Background

Placental efficiency (PE) describes the relationship between placental and fetal weights (fetal wt/placental wt). Within litters, PE can vary drastically, resulting in similarly sized pigs associated with differently sized placentas, up to a 25% weight difference. However, the mechanisms enabling the smaller placenta to grow a comparable littermate are unknown. To elucidate potential mechanisms, morphological measurements and gene expression profiles in placental and associated endometrial tissues of high PE and low PE feto-placental units were compared. Tissue samples were obtained from eight maternal line gilts during gestational day 95 ovario-hysterectomies. RNA was extracted from tissues of feto-placental units with the highest and lowest PE in each litter and sequenced.

Results

Morphological measurements, except placental weight, were not different (P > 0.05) between high and low PE. No DEG were identified in the endometrium and 214 DEG were identified in the placenta (FDR < 0.1), of which 48% were upregulated and 52% were downregulated. Gene ontology (GO) analysis revealed that a large percentage of DEG were involved in catalytic activity, binding, transporter activity, metabolism, biological regulation, and localization. Four GO terms were enriched in the upregulated genes and no terms were enriched in the downregulated genes (FDR < 0.05). Eight statistically significant correlations (P < 0.05) were identified between the morphological measurements and DEG.

Conclusion

Morphological measures between high and low PE verified comparisons were of similarly sized pigs grown on different sized placentas, and indicated that any negative effects of a reduced placental size on fetal growth were not evident by day 95. The identification of DEG in the placenta, but absence of DEG in the endometrium confirmed that the placenta responds to the fetus. The GO analyses provided evidence that extremes of PE are differentially regulated, affecting components of placental transport capacity like nutrient transport and blood flow. However, alternative GO terms were identified, indicating the complexity of the relationship between placental and fetal weights. These findings support the use of PE as a marker of placental function and provide novel insight into the genetic control of PE, but further research is required to make PE production applicable.

Factors contributing to the variation in placental efficiency on days 70, 90, and 110 of gestation in gilts

Variations in placental efficiency (PE), a measure of grams of fetus produced per gram of placenta, were initially researched between swine breeds, where increased PE was associated with larger litters. Placental efficiency was also found to vary greatly within production herds and individual litters; however, the use of PE as a selection tool has been debated. Nonetheless, PE is an index of feto-placental adaptation and may help identify compensatory mechanisms that maintain fetal growth when placental size is reduced, potentially providing an opportunity to address production concerns like low birth weights and preweaning survival. Since the nutrient transport capacity of the placenta largely depends on vasculature and nutrient transporter abundance, the objectives of this experiment were to 1) determine the mRNA expression of genes encoding nutrient transporters in the placenta and adjacent endometrium, and 2) evaluate if a relationship existed between PE and vascular density and/or nutrient transporters. Gilts (n = 19) were ovario-hysterectomized on day 70, 90, or 110 of gestation to collect placental and adjacent endometrial samples. The mean litter size was 11.1. Placental efficiency increased (P < 0.0001) throughout the end of gestation, while the range of PE increased from day 70 to 90 and was reduced on day 110 (P < 0.0001). Placental efficiency and placental weight were negatively correlated throughout gestation (70 d, r = -0.83, P < 0.0001; 90 d, r = -0.81, P < 0.0001; 110 d, r = -0.44, P < 0.0007), but the negative correlation between PE and fetal weight was not maintained as gestation progressed (70 d, r = -0.58, P < 0.0001; 90 d, r = -0.36, P < 0.0005; 110 d, r = 0.09, P = 0.51). Based on conditional effects plots, variations in PE were associated with alterations in amino acid transporter expression in the placenta (SLC7A7, SLC3A1) and endometrium (SLC7A1) on day 70. On day 90, PE had a positive relationship with placental expression of a glucose transporter (SLC2A3), and on day 110 PE was positively related to placental vascular density. The results suggest utero-placental adaptations occur as a compensation for reduced placental size to meet the increasing nutrient demands of the growing fetus during late gestation in swine. Furthermore, nutrient requirements differ for individual feto-placental units on a given day; therefore, optimizing nutrient availability during late gestation may improve fetal growth and survival.

Identification of Differentially Expressed Genes in Porcine Ovaries at Proestrus and Estrus Stages Using RNA-Seq Technique

Estrus is an important factor for the fecundity of sows, and it is involved in ovulation and hormone secretion in ovaries. To better understand the molecular mechanisms of porcine estrus, the expression patterns of ovarian mRNA at proestrus and estrus stages were analyzed using RNA sequencing technology. A total of 2,167 differentially expressed genes (DEGs) were identified (P ≤ 0.05, |log2  Ratio| ≥ 1), of which 784 were upregulated and 1,383 were downregulated in the estrus compared with the proestrus group. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the cellular process, single-organism process, cell and cell part, and binding and metabolic process. In addition, a pathway analysis showed that these DEGs were significantly enriched in 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including cell adhesion molecules, ECM-receptor interaction, and cytokine-cytokine receptor interaction. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) confirmed the differential expression of 10 selected DEGs. Many of the novel candidate genes identified in this study will be valuable for understanding the molecular mechanisms of the sow estrous cycle.

Effect of dietary chitosan oligosaccharide supplementation on the pig ovary transcriptome

Fecundity improvement is one of the most important economic traits for the swine industry. In this study, we identified 486 differentially expressed genes associated with sow prolificacy from COS administrated sow ovaries by RNA-seq.

DNA methylation patterns and gene expression associated with litter size in Berkshire pig placenta

Increasing litter size is of great interest to the pig industry. DNA methylation is an important epigenetic modification that regulates gene expression, resulting in livestock phenotypes such as disease resistance, milk production, and reproduction. We classified Berkshire pigs into two groups according to litter size and estimated breeding value: smaller (SLG) and larger (LLG) litter size groups. Genome-wide DNA methylation and gene expression were analyzed using placenta genomic DNA and RNA to identify differentially methylated regions (DMRs) and differentially expressed genes (DEGs) associated with litter size. The methylation levels of CpG dinucleotides in different genomic regions were noticeably different between the groups, while global methylation pattern was similar, and excluding intergenic regions they were found the most frequently in gene body regions. Next, we analyzed RNA-Seq data to identify DEGs between the SLG and LLG groups. A total of 1591 DEGs were identified: 567 were downregulated and 1024 were upregulated in LLG compared to SLG. To identify genes that simultaneously exhibited changes in DNA methylation and mRNA expression, we integrated and analyzed the data from bisulfite-Seq and RNA-Seq. Nine DEGs positioned in DMRs were found. The expression of only three of these genes (PRKG2, CLCA4, and PCK1) was verified by RT-qPCR. Furthermore, we observed the same methylation patterns in blood samples as in the placental tissues by PCR-based methylation analysis. Together, these results provide useful data regarding potential epigenetic markers for selecting hyperprolific sows.

Transcriptomic Analysis of Oenococcus oeni SD-2a Response to Acid Shock by RNA-Seq

Oenococcus oeni can be applied to conduct malolactic fermentation (MLF), but also is the main species growing naturally in wine. Due to the high stress tolerance, it is an interesting model for investigating acid response mechanisms. In this study, the changes in the transcriptome of O.oeni SD-2a during the adaptation period have been studied. RNA-seq was introduced for the transcriptomic analysis of O. oeni samples treated with pH 4.8 and pH 3.0 at 0 and 1 h, respectively. Gene ontology (GO) and Kyoto encyclopedia of genes and genome (KEGG) were performed to compare the transcriptome data between different treatments. From GO analysis, the majority of differentially expressed genes (DEGs) (pH 3.0_1 h-VS-pH 4.8_1 h, pH 3.0_1 h-VS-pH 4.8_0 h, and pH 4.8_1 h-VS-pH 4.8_0 h) were found to be involved in the metabolic process, catalytic activity, cellular process, and binding. KEGG analysis reveals that the most functional gene categories affected by acid are membrane transport, amino acid metabolism and carbohydrate metabolism. Some genes, like the heat shock protein Hsp20, malate transporter and malate permease, were also over-expressed in response to acid stress. In addition, a considerable proportion of gene indicate a significantly different expression in this study, are novel, which needs to be investigated further. These results provide a new viewpoint and crucial resource on the acid stress response in O. oeni.

Genome-wide differential mRNA expression profiles in follicles of two breeds and at two stages of estrus cycle of gilts

Estrus expression by gilts and sows is hereditable and important for heat detection. To better understand the molecular biological mechanisms of estrus expression in gilts, the mRNA expression profiles of follicular tissue from Large White gilts in diestrus (LD, n = 3) and estrus (LE, n = 3), and Chinese indigenous Mi gilts in diestrus (MD, n = 2) and estrus (ME, n = 3) were investigated using RNA sequencing. We detected 122,804-335,295 SNPs, 6,140-14,947 InDel and 12 types of AS events (39.57% TSS, 34.90% TTS) in 11 samples. A total of 2,838 differentially expressed genes (DEGs) were found in LD vs MD, LE vs ME, LE vs LD, or ME vs MD comparisons. Two DEGs (ACP5 and PIGS) were observed in all comparisons. Two new genes (ENSSSCG00000028235 and ENSSSCG00000021903) were exclusively expressed in Mi and Large White gilts, respectively. Bioinformatics analyses indicate that these DEGs are involved in single-organism process, catalytic activity, cell adhesion and enriched in ECM-receptor interaction, olfactory transduction, ovarian steroidogenesis, steroid biosynthesis and CAMs signaling pathways. These results of RNA-Seq have provided important information for screening the key functional genes or molecular markers of estrus expression in gilts.