Differences in X-chromosome transcriptional activity and cholesterol metabolism between placentae from swine breeds from Asian and Western origins

Integrated analysis of lncRNA, miRNA and mRNA expression profiles reveals regulatory pathways associated with pig testis function

Long noncoding RNA (lncRNA) and microRNA (miRNA) are known to play pivotal roles in mammalian testicular function and spermatogenesis. However, their impact on porcine male reproduction has yet to be well unraveled. Here, we sequenced and identified lncRNA and miRNA expressed in the testes of Chinese indigenous Banna mini-pig inbred line (BMI) and introduced Western Duroc (DU) and Large White (LW) pigs. By pairwise comparison (BMI vs DU, BMI vs LW, and DU vs LW), we found the gene expression differences in the testes between Chinese local pigs and introduced Western commercial breeds were more striking than those between introduced commercial breeds. Furthermore, we found 1622 co-differentially expressed genes (co-DEGs), 122 co-differentially expressed lncRNAs (co-DELs), 39 co-differentially expressed miRNAs (co-DEMs) in BMI vs introduced commercial breeds (DU and LW). Functional analysis revealed that these co-DEGs and co-DELs/co-DEMs target genes were enriched in male sexual function pathways, including MAPK, AMPK, TGF-β/Smad, Hippo, NF-kappa B, and PI3K/Akt signaling pathways. Additionally, we established 10,536 lncRNA-mRNA, 11,248 miRNA-mRNA pairs, and 62 ceRNA (lncRNA-miRNA-mRNA) networks. The ssc-miR-1343 had the most interactive factors in the ceRNA network, including 20 mRNAs and 3 lncRNAs, consisting of 56 ceRNA pairs. These factors played extremely important roles in the regulation of testis function as key nodes in the interactive regulatory network. Our results provide insight into the functional roles of lncRNAs and miRNAs in porcine testis and offer a valuable resource for understanding the differences between Chinese indigenous and introduced Western pigs.

Transcriptome and methylome sequencing reveals altered long non-coding RNA genes expression and their aberrant DNA methylation in equine sarcoids

Recent publications confirmed that long non-coding RNAs (lncRNAs) perform an essential function in gene-specific transcription regulation. Nevertheless, despite its important role, lncRNA has not yet been described in equine sarcoids, the skin neoplasia of horses. Therefore, the aim of this study is to deepen the knowledge about lncRNA expression in the pathogenesis of equine sarcoids and provide new insight into the regulatory function of lncRNA in the bovine papillomavirus-dependent neoplasia of horse dermal tissues. RNA sequencing (RNA-seq) data from 12 equine sarcoid samples and the corresponding controls were reanalyzed in this study. A total of 3396 differentially expressed (DE) lncRNAs and 128 DElncRNA-DE genes (DEGs) pairs were identified. Differentially expressed lncRNAs predicted target genes were enriched in pathways associated with inter alia the extracellular matrix disassembly and cancer pathways. Furthermore, methylation data from the same samples were integrated into the analysis, and 12 DElncRNAs were described as potentially disturbed by aberrant methylation. In conclusion, this study presents novel data about lncRNA's role in the pathogenesis of equine sarcoids.

Identification and Characterization of lncRNA and mRNA in Testes of Landrace and Hezuo Boars

Simple Summary

Precocious puberty is an excellent reproductive trait in domestic animals, which can generate higher breeding benefits in livestock production. However, regulators associated with this sexual maturation process remain largely unknown. Chinese Hezuo (HZ) boars are known for their early sexual maturity. In this work, the characteristics of precocious puberty in HZ pigs were confirmed by histological analysis, and some important long noncoding RNA (lncRNA) and mRNA were identified in the testes of immature (30-day-old) and mature (120-day-old) HZ boars, which could play a key role in precocious puberty. These results will provide a theoretical basis for further research on the regulatory mechanism of precocious puberty, which is important for accelerating the breeding process of highly fertile animals.

Abstract

Chinese HZ boars are typical plateau miniature boars characterized by precocious puberty, which is closely related to testicular development and spermatogenesis. Accumulating evidence indicates that lncRNA is involved in the testicular development and regulation of spermatogenesis. However, little is known about the lncRNA precocious regulation in testicular development and spermatogenesis on early sexual maturity of HZ boars. Thus, we investigated the expression and characterization of lncRNA and mRNA in 30-day-old and 120-day-old HZ boar testes using transcriptome to explore precocious puberty. Landrace (LC) boar was treated as the control. Histological analyses indicated that HZ boar underwent puberty development at an earlier stage than LC boar and had achieved sexual maturity at 120 days old. RNA-Seq yielded a total of 187 lncRNAs and 984 mRNAs; these molecules were identified as possible candidates for precocious puberty. GO terms and KEGG pathways enrichment analyses revealed that the differentially expressed lncRNA and their targeted genes were involved in metabolic pathways regulating testis development and spermatogenesis, such as the PI3K-Akt, TGF-beta and Wnt pathways. Further screening, some lncRNA (such as LOC102166140, LOC110259451, and MSTRG.15011.2), and mRNA (such as PDCL2, HSD17B4, SHCBP1L, CYP21A2, and SPATA3) were found to be possibly associated with precocious puberty, which would add to our understanding of the molecular regulatory mechanisms of precocious puberty. This study provided valuable information for further study of the role of lncRNA and mRNA in the process of precocious puberty.

Transcriptomic and ChIP-seq Integrative Analysis Reveals Important Roles of Epigenetically Regulated lncRNAs in Placental Development in Meishan Pigs

The development of the placental fold, which increases the maternal–fetal interacting surface area, is of primary importance for the growth of the fetus throughout the whole pregnancy. However, the mechanisms involved remain to be fully elucidated. Increasing evidence has revealed that long non-coding RNAs (lncRNAs) are a new class of RNAs with regulatory functions and could be epigenetically regulated by histone modifications. In this study, 141 lncRNAs (including 73 up-regulated and 68 down-regulated lncRNAs) were identified to be differentially expressed in the placentas of pigs during the establishment and expanding stages of placental fold development. The differentially expressed lncRNAs and genes (DElncRNA-DEgene) co-expression network analysis revealed that these differentially expressed lncRNAs (DElncRNAs) were mainly enriched in pathways of cell adhesion, cytoskeleton organization, epithelial cell differentiation and angiogenesis, indicating that the DElncRNAs are related to the major events that occur during placental fold development. In addition, we integrated the RNA-seq (RNA sequencing) data with the ChIP-seq (chromatin immunoprecipitation sequencing) data of H3K4me3/H3K27ac produced from the placental samples of pigs from the two stages (gestational days 50 and 95). The analysis revealed that the changes in H3K4me3 and/or H3K27ac levels were significantly associated with the changes in the expression levels of 37 DElncRNAs. Furthermore, several H3K4me3/H3K27ac-lncRNAs were characterized to be significantly correlated with genes functionally related to placental development. Thus, this study provides new insights into understanding the mechanisms for the placental development of pigs.

No imprinted XIST expression in pigs: biallelic XIST expression in early embryos and random X inactivation in placentas

Dosage compensation, which is achieved by X-chromosome inactivation (XCI) in female mammals, ensures balanced X-linked gene expression levels between the sexes. Although eutherian mammals commonly display random XCI in embryonic and adult tissues, imprinted XCI has also been identified in extraembryonic tissues of mouse, rat, and cow. Little is known about XCI in pigs. Here, we sequenced the porcine XIST gene and identified an insertion/deletion mutation between Asian- and Western-origin pig breeds. Allele-specific analysis revealed biallelic XIST expression in porcine ICSI blastocysts. To investigate the XCI pattern in porcine placentas, we performed allele-specific RNA sequencing analysis on individuals from reciprocal crosses between Duroc and Rongchang pigs. Our results were the first to reveal that random XCI occurs in the placentas of pigs. Next, we investigated the H3K27me3 histone pattern in porcine blastocysts, showing that only 17-31.8% cells have attained XCI. The hypomethylation status of an important XIST DMR (differentially methylated region) in gametes and early embryos demonstrated that no methylation is pre-deposited on XIST in pigs. Our findings reveal that the XCI regulation mechanism in pigs is different from that in mice and highlight the importance of further study of the mechanisms regulating XCI during early porcine embryo development.

Analysis of Long Non-Coding RNA and mRNA Expression Profiling in Immature and Mature Bovine (Bos taurus) Testes

Testis development and spermatogenesis are strictly regulated by numbers of genes and non-coding genes. However, long non-coding RNAs (lncRNAs) as key regulators in multitudinous biological processes have not been systematically identified in bovine testes during sexual maturation. In this study, we comprehensively analyzed lncRNA and mRNA expression profiling of six bovine testes at 3 days after birth and 13 months by RNA sequencing. 23,735 lncRNAs and 22,118 mRNAs were identified, in which 540 lncRNAs (P-value < 0.05) and 3,525 mRNAs (P-adjust < 0.05) were significantly differentially expressed (DE) between two stages. Correspondingly, the results of RT-qPCR analysis showed well correlation with the transcriptome data. Moreover, GO and KEGG enrichment analyses showed that DE genes and target genes of DE lncRNAs were enriched in spermatogenesis. Furthermore, we constructed lncRNA–gene interaction networks; consequently, 15 DE lncRNAs and 12 cis-target genes were involved. The target genes (SPATA16, TCF21, ZPBP, PACRG, ATP8B3, COMP, ACE, and OSBP2) were found associated with bovine sexual maturation. In addition, the expression of lncRNAs and cis-target genes was detected in bovine Leydig cells, Sertoli cells, and spermatogonia. Our study identified and analyzed lncRNAs and mRNAs in testis tissues, suggesting that lncRNAs may regulate testis development and spermatogenesis. Our findings provided new insights for further investigation of biological function in bovine lncRNA.

Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.

Genome-Wide Identification of Histone Modifications Involved in Placental Development in Pigs

Development of placental folds is a critical event affecting placental function in pigs because it can increase surface area for improvement in capillary density as gestation advances. However, the molecular mechanisms of the event are not well defined. Histone modifications have important roles in gene regulation. To investigate their effects on regulation of genes controlling porcine placental development, RNA-seq and ChIP-seq of porcine placental tissues from gestational days 50 (establishment stage of placental folds) and 95 (expanding stage of placental folds) were carried out in this study. The differentially expressed genes were identified and of which the down- and up-regulated genes are related to endoplasmic reticulum (ER) stress and angiogenesis, respectively. In addition, we mapped the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), which are associated with transcriptional activation. A number of differential modification regions between the 2 gestational stages were identified and majority of them are those with increased signals of H3K4me3 (14,576 out of 16,931). Furthermore, we observed that the increase of H3K4me3 is significantly correlated with the elevated expression levels of the neighboring genes, and notably, these genes were enriched in pathways related to blood vessel formation and microvascular permeability. Taken together, the findings suggest important roles of histone modifications on placental remolding in response to developmental changes.

Improved analysis of anthocyanins and vitamin C in blue-purple potato cultivars

Methods were optimized for extraction and quantification of anthocyanins (ACY) and vitamin C in potatoes. Acidified aqueous methanol (70%) was the optimal extraction solvent and freeze-drying significantly improved the extraction yield of ACY. The content of ACY varied widely in five potato cultivars from 0.42 to 3.18mg/g dry weight, with the latter being the highest value found in the Finnish cultivar 'Synkeä Sakari'. Compared with dithiothreitol (DTT), tris(2-carboxyethyl) phosphine hydrochloride (TCEP) was more efficient in reducing dehydroascorbic acid (DHA) to ascorbic acid (AA) and for quantifying the content of total ascorbic acid (TAA). For raw potatoes, quantification of TAA after treatment with TCEP was more reliable than a direct analysis of AA, whereas AA can be analyzed directly in steam-cooked samples. The TAA contents in the three potato cultivars were around 0.30-0.35mg/g dry weight. The loss of AA in steam cooking was 24%.

Polymorphism in the intron 20 of porcine O-linked N-acetylglucosamine transferase

Objective

O-linked N-acetylglucosamine (O-GlcNAc) transferase (OGT) catalyzes the addition of O-GlcNAc and GlcNAcylation has extensive crosstalk with phosphorylation to regulate signaling and transcription. Pig OGT is located near the region of chromosome X that affects follicle stimulating hormone level and testes size. The objective of this study was to find the variations of OGT between European and Chinese pigs.

Methods

Pigs were tested initially for polymorphism in OGT among European and Chinese pigs by polymerase chain reaction and sequencing at the U.S. Meat Animal Research Center (USMARC). The polymorphism was also determined in an independent population of pigs including European and Chinese Meishan (ME) breeds at the National Institute of Animal Science (NIAS, RDA, Korea).

Results

The intron 20 of OGT from European and Chinese pigs was 514 and 233 bp, respectively, in the pigs tested initially. They included 1 White composite (WC) boar and 7 sows (2 Minzu×WC, 2 Duroc [DU]×WC, 2 ME×WC, 1 Fengzing×WC) at USMARC. The 281-bp difference was due to an inserted 276-bp element and GACTT in European pigs. When additional WC and ME boars, the grandparents that were used to generate the 1/2ME×1/2WC parents, and the 84 boars of 16 litters from mating of 1/2ME×1/2WC parents were analyzed, the breeds of origin of X chromosome quantitative trait locus (QTL) were confirmed. The polymorphism was determined in an independent population of pigs including DU, Landrace, Yorkshire, and ME breeds at NIAS. OGT was placed at position 67 cM on the chromosome X of the USMARC swine linkage map.

Conclusion

There was complete concordance with the insertion in European pigs at USMARC and NIAS. This polymorphism could be a useful marker to identify the breed of origin of X chromosome QTL in pigs produced by crossbreeding Chinese and European pigs.

The effect of maternal obesity on fatty acid transporter expression and lipid metabolism in the full-term placenta of lean breed swine

This study was conducted to evaluate the influence of back-fat thickness (BF), at mating of sows, on the maternal and newborn circulating lipids, expression of placental fatty acids (FA) transporters and lipid accumulation in placenta. Full-term placentas were obtained by vaginal delivery from BFI (9-14 mm; n = 37), BFII (15-19 mm; n = 43) and BFIII (20-27 mm; n = 38) sows according to BF at mating, and frozen placental sections were analysed for fat accumulation. Blood samples were collected from the sows of day 105 pregnancy and from cord blood at delivery. mRNA and protein expression levels were evaluated with real-time RT-PCR and Western blotting. Our results demonstrated that BFII females had significantly increased litter weight and placental efficiency, decreased maternal triglyceride (TG) and non-esterified fatty acids (NEFA) levels, decreased maternal IL-6, TNFα and leptin levels compared to BFIII females (p < .05). BFIII sows were associated with significantly decreased newborn TG levels, increased newborn glucose, IL-6 and TNFα levels compared to BFI or BFII sows (p < .05). BFI and BFII females had significantly decreased placental TG, NEFA and cholesterol (CHOL) contents compared to BFIII females (p < .05). Moreover, decreased CD36, FATP1, FABP4, and FABP1 mRNA and protein and FATP4 protein expression, and increased LPL activity were also observed in BFIII group compared with BFII group (p < .05). PPARγ mRNA and protein and lipogenic genes such as SREBP-1c, ACSL1, ACCα, FAS and SCD mRNA expression were downregulated or upregulated, respectively, in the placentas of BFIII sows compared to BFI or BFII sows (p < .05). Overall, this study demonstrated that there is no advantage, in terms of litter live size, litter weight and placental FA transport and metabolism, in performing the mating of sows with BF>19 mm.

Systematic Analysis of Long Noncoding RNAs in the Senescence-accelerated Mouse Prone 8 Brain Using RNA Sequencing

Long noncoding RNAs (lncRNAs) may play an important role in Alzheimer's disease (AD) pathogenesis. However, despite considerable research in this area, the comprehensive and systematic understanding of lncRNAs in AD is still limited. The emergence of RNA sequencing provides a predictor and has incomparable advantage compared with other methods, including microarray. In this study, we identified lncRNAs in a 7-month-old mouse brain through deep RNA sequencing using the senescence-accelerated mouse prone 8 (SAMP8) and senescence-accelerated mouse resistant 1 (SAMR1) models. A total of 599,985,802 clean reads and 23,334 lncRNA transcripts were obtained. Then, we identified 97 significantly upregulated and 114 significantly downregulated lncRNA transcripts from all cases in SAMP8 mice relative to SAMR1 mice. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses revealed that these significantly dysregulated lncRNAs were involved in regulating the development of AD from various angles, such as nerve growth factor term (GO: 1990089), mitogen-activated protein kinase signaling pathway, and AD pathway. Furthermore, the most probable AD-associated lncRNAs were predicted and listed in detail. Our study provided the systematic dissection of lncRNA profiling in SAMP8 mouse brain and accelerated the development of lncRNA biomarkers in AD. These attracting biomarkers could provide significant insights into AD therapy in the future.

Difference in expression patterns of placental cholesterol transporters, ABCA1 and SR-BI, in Meishan and Yorkshire pigs with different placental efficiency

Cholesterol is a key cell membrane component and precursor of steroid hormones. The maternal cholesterol is an important exogenous cholesterol source for the developing embryos and its transportation is mediated by ABCA1 and SR-BI. Here we reported that during the peri-implantation period in pigs, ABCA1 was expressed by uterine luminal epithelium (LE) and interestingly, its expression was more abundantly in LE on mesometrial side of uterus. However, SR-BI was expressed primarily by LE, glandular epithelial cells (GE) and trophoblast cells (Tr). During the placentation period, the expression levels of ABCA1 and SR-BI proteins at epithelial bilayer and placental areolae were significantly higher in Chinese Meishan pigs compared to Yorkshire pigs. Consisitently, mRNA levels of HMGCR, the rate-limiting enzyme for cholesterol synthesis, were significantly higher in Meishan placentas than in Yorkshire placentas. Our findings revealed the routes of transplacental cholesterol transport mediated by ABCA1 and SR-BI in pigs and indicated that ABCA1 related pathway may participate in anchoring the conceptus to the mesometrial side of uterus. Additionally, an ABCA1 dependent compensatory mechanism related to the placental efficiency in response to the smaller placenta size in Meishan pigs was suggested.

Genetic determinants of pig birth weight variability

Background

Piglet birth weight variability, a trait also known as the within-litter homogeneity of birth weight, reflects the sow’s prolificacy, because it is positively genetically correlated with preweaning mortality but negatively correlated with the mean growth of piglets during sucking. In addition, the maternal additive genetic variance and heritability has been found exist for this trait, thus, reduction in the variability of piglet birth weight to improve the sow prolificacy is possible by selective breeding.

Results

We performed a genome wide association study (GWAS) in 82 sows with extreme standard deviation of birth weights within the first parity to identify significant SNPs, and finally 266 genome-wide significant SNPs (p < 0.01) were identified. These SNPs were mainly enriched on chromosome 7, 1, 13, 14, 15 and 18. We further scanned genes of the top 50 SNPs with the lowest p values and found some genes involved in plasma glucose homeostasis (GLP1R) and lipid metabolism as well as maternal-fetal lipid transport (AACS, APOB, OSBPL10 and LRP1B) which may contribute to the birth weight variability trait.

Conclusions

Birth weight variability trait has a low heritability. It is not easy to get significant signal by GWAS using small sample size. Herein, we identified some candidate chromosome regions especially chromosome 7 and suggested five genes which may provide some information for the further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0309-6) contains supplementary material, which is available to authorized users.

The Expression Pattern of MicroRNAs and the Associated Pathways Involved in the Development of Porcine Placental Folds That Contribute to the Expansion of the Exchange Surface Area

The development of the microscopically folded structure of the diffuse epitheliochorial placenta in pigs is important because it expands the surface area for maternal-fetal exchange, resulting in an increase in placental efficiency. To better understand the regulatory mechanisms involved in this process, we characterized miRNA expression profiles in porcine placentas during the initiation and establishment of placental fold development. A total of 42 miRNAs were found to be differentially expressed, and their putative target genes were predicted using four target prediction programs. Following a comparative analysis with published gene expression pattern data obtained from porcine placentas in the corresponding stages of placental fold development, only those genes that were negatively correlated with miRNA expression were retained for further function and pathway enrichment analysis. The results showed that the up-regulated miRNAs were associated mainly with extracellular matrix remodeling and tissue morphogenesis, while the down-regulated miRNAs were related to cell proliferation and signal transduction. Furthermore, we provide evidence that miR-130b may facilitate the expression of HPSE, which has been reported to be a regulator of the folding of the pig placenta, by suppressing the expression of PPARG. In addition, we also reveal that the miRNA-target pairs expressed in the pig placenta may trigger the degradation of the stromal matrix and basement membrane (miR-29a-COL1A2, COL3A1, and LAMC1) and regulate trophoblast epithelial cell adherens junctions (the miR-200 family and miR-205-ZEB2-CDH1) and proliferation (miR-17-92 cluster-HBP1 and ULK1). Taken together, these results indicate that miRNAs and related pathways may have potential roles in porcine placental fold development.

Identification of the porcine XIST gene and its differential CpG methylation status in male and female pig cells

XIST, a long non-coding RNA, plays an important role in triggering X chromosome inactivation in eutherians, and is used extensively for qualifying stem cells and cloned embryos. However, a porcine XIST has not yet been thoroughly identified despite its biological importance in a wide variety of research fields. Here, we present a full-length porcine XIST sequence assembled using known sequences (GenBank), RNA-Seq data (NCBI SRA), and PCR/sequencing. The proposed porcine XIST gene model encodes a 25,215-bp transcript consisting of 7 exons, including two conserved and two porcine-specific repeat regions. Transcription covering the entire XIST region was observed specifically in female cells, but not in male cells. We also identified eight transcription starting sites (TSSs) and evaluated CpG methylation patterns in the upstream (+2.0 kb) and downstream (−2.0 kb) regions. Sixty-seven CG di-nucleotides identified in the target region were considered to be candidate CpG sites, and were enriched in the following two regions: −284 to +53 bp (13 sites) and +285 to +1,727 bp (54 sites) from the selected TSS. Male 5` region of XIST (64.5 sites, 96.26%) had a higher level of CpG methylation than female DNA (33.4 sites, 49.85%). Taken together, our results revealed that the porcine XIST gene is expressed exclusively in female cells, which is influenced by the lower level of CpG methylation in the putative promoter region compared with male cells. The porcine XIST presented in this study represents a useful tool for related research areas such as porcine embryology and stem cell biology.