Microarray profiling for differential gene expression in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows

Allelic to genome wide perspectives of swine genetic variation to litter size and its component traits

Litter size is a complex and sex limited trait that depends on various biological, managemental and environmental factors. Owing to its low heritability it is inefficaciously selected by traditional methods. However, due to higher heritability of ovulation rate and embryo survival, selection based on component traits of litter size is advocated. QTL analysis and candidate gene approach are among the various supplementary/alternate strategies for selection of litter size. QTL analysis is aimed at identifying genomic regions affecting trait of interest significantly. Candidate gene approach necessitates identification of genes potentially affecting the trait. There are various genes that significantly affect litter size and its component traits viz. ESR, LEP, BF, IGFBP, RBP4, PRLR, CTNNAL1, WNT10B, TCF12, DAZ, and RNF4. These genes affect litter size in a complex interacting manner. Lately, genome wide association study (GWAS) have been utilized to unveil the genetic and biological background of litter traits, and elucidate the genes governing litter size. Favorable SNPs in these genes have been identified and offers a scope for inclusion in selection programs thereby increasing breeding efficiency and profit in pigs. The review provides a comprehensive coverage of investigations carried out globally to unravel the genetic variation in litter size and its component traits in pigs, both at allelic and genome wide level. It offers a current perspective on different strategies including the profiling of candidate genes, QTLs, and genome wide association studies as an aid to efficient selection for litter size and its component traits.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.

Ovarian Circular RNAs Associated with High and Low Fertility in Large White Sows during the Follicular and Luteal Phases of the Estrous Cycle

In this study, the ovarian tissues of Large White pigs were mined for novel circular RNAs (circRNAs), following which, their molecular characteristics and potential mechanisms for fertility regulation were examined. RNA sequencing was used for transcriptome analysis of ovarian follicles and corpora lutea in Large White sows with high (H) and low (L) fertility during the follicular (F) and luteal (L) phases of the estrous cycle. In total, 21,386 circRNA derived from 4535 host genes were identified. Differentially expressed circRNAs were detected in the LH vs. LL (1079) and in the FH vs. FL (1077) comparisons, and their host genes were enriched in steroid biosynthesis and forkhead box O (FOXO), thyroid hormone, cell cycle, and tumor growth factor (TGF)-beta signaling pathways. Protein–protein interaction networks were constructed on the basis of the host genes that were significantly enriched in pathways related to reproductive processes, with AKT3 and PP2CB serving as the hub genes in the networks of the LH vs. LL and FH vs. FL comparisons, respectively. The microRNA (miRNA) binding sites of the differentially expressed circRNAs were predicted, and 128 (LH vs. LL) and 113 (FH vs. FL) circRNA–miRNA pairs were identified. Finally, circRNA–miRNA negative regulatory networks were established on the basis of the gene expression profiles and bioinformatic analyses. In the current study, differentially expressed circRNAs were observed in ovarian tissues between the H and L fertility groups in both F and L phases of the estrous cycle, which suggested roles in pig fertility regulation. These findings provide new clues for elucidating fertility differences in pigs.

NRDR inhibits estradiol synthesis and is associated with changes in reproductive traits in pigs

Cumulus cells secreting steroid hormones have important functions in oocyte development. Several members of the short-chain dehydrogenase/reductase (SDR) family are critical to the biosynthesis of steroid hormones. NADPH-dependent retinol dehydrogenase/reductase ( NRDR), a member of the SDR superfamily, is overexpressed in pig breeds that also show high levels of androstenone. However, the potential functions and regulatory mechanisms of NRDR in pig ovaries have not been reported to date. The present study demonstrated that NRDR is highly expressed in pig ovaries and is specifically located in cumulus granulosa cells. Functional studies showed that NRDR inhibition increased estradiol synthesis. Both pregnant mare serum gonadotropin and human chorionic gonadotropin downregulated the expression of NRDR in pig cumulus granulosa cells. When the relationship between reproductive traits and single-nucleotide polymorphisms (SNPs) of the NRDR gene was examined, we found that two SNPs affected reproductive traits. SNP rs701332503 was significantly associated with a decrease in the total number of piglets born during multiparity, and rs326982309 was significantly associated with an increase in the average birth weight during primiparity. Thus, NRDR has an important role in steroid hormone biosynthesis in cumulus granulosa cells, and NRDR SNPs are associated with changes in porcine reproduction traits.

SLA-11 mutations are associated with litter size traits in Large White and Chinese DIV pigs

Litter size is an important economic traits in pigs. SLA-11 gene is a member of SLA (swine leukocyte antigen) complex. In our previous study, the SLA-11 gene was differentially expressed in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows. Here, we identified two mutations (c.754-132 T > C and c.1421 + 38 T > C) in SLA-11 gene and analyzed the associations of two SNPs with litter size traits in Large White (n = 263) and DIV (n = 117) sows. The results showed that in Large White pigs, SLA-11 c.754-132 CC sows produced 0.74 and 0.87 more pigs per litter for TNB and NBA of all parities than did TT sows (p < .05); In DIV pigs, SLA-11 c.754-132 CC sows produced 1.17 more pigs per litter for TNB of all parities than did TC sows (p < .05). In Large White pigs, SLA-11 c.1421 + 38 CC sows produced 0.9 more pigs per litter for TNB of all parities than did TT sows (p < .05), while in DIV pigs SLA-11 c.1421 + 38 CC sows produced 0.84 and 0.7 less pigs per litter for TNB and NBA of all parities than did TT sows (p < .05). Our research indicated that SLA-11 mutations were potential molecular markers for improving the litter size traits in pigs.

Identification and verification of potential piRNAs from domesticated yak testis

PIWI-interacting RNAs (piRNA) are small non-coding RNA molecules expressed in animal germ cells that interact with PIWI family proteins to form RNA–protein complexes involved in epigenetic and post-transcriptional gene silencing of retrotransposons and other genetic elements in germ line cells, including reproductive stem cell self-sustainment, differentiation, meiosis and spermatogenesis. In the present study, we performed high-throughput sequencing of piRNAs in testis samples from yaks in different stages of sexual maturity. Deep sequencing of the small RNAs (18–40 nt in length) yielded 4,900,538 unique reads from a total of 53,035,635 reads. We identified yak small RNAs (18–30 nt) and performed functional characterization. Yak small RNAs showed a bimodal length distribution, with two peaks at 22 nt and >28 nt. More than 80% of the 3,106,033 putative piRNAs were mapped to 4637 piRNA-producing genomic clusters using RPKM. 6388 candidate piRNAs were identified from clean reads and the annotations were compared with the yak reference genome repeat region. Integrated network analysis suggested that some differentially expressed genes were involved in spermatogenesis through ECM–receptor interaction and PI3K-Akt signaling pathways. Our data provide novel insights into the molecular expression and regulation similarities and diversities in spermatogenesis and testicular development in yaks at different stages of sexual maturity.

MicroRNA-144 is regulated by CP2 and decreases COX-2 expression and PGE2 production in mouse ovarian granulosa cells

Mammalian folliculogenesis is a complex process in which primordial follicles develop into pre-ovulatory follicles, followed by ovulation to release mature oocytes. In this study, we explored the role of miR-144 in ovulation. miR-144 was one of the differentially expressed microRNAs, which showed 5.59-fold changes, in pre-ovulatory ovarian follicles between Large White and Chinese Taihu sows detected by Solexa deep sequencing. We demonstrated that overexpression of miR-144 significantly decreased the luciferase reporter activity under the control of the cyclooxygenase-2 (COX-2) or mothers against decapentaplegic homologue 4 (Smad4) 3'-untranslated region (3'-UTR) and suppressed COX-2 and Smad4 expression. In contrast, a miR-144 inhibitor increased COX-2 and Smad4 expression in mouse granulosa cells (mGCs). Meanwhile, Smad4 upregulated COX-2 expression, but this effect was abolished when the mGCs were treated with the transforming growth factor beta signalling pathway inhibitor SB431542. Moreover, luciferase reporter, chromatin immunoprecipitation and electrophoretic mobility shift assay results showed that the transcription factor CP2 upregulated miR-144 expression, which partially contributed to the suppression of COX-2 in mGCs. Both CP2 and miR-144 alter prostaglandin E2 (PGE2) production by regulating COX-2 expression. In addition, miR-144 regulated mGC apoptosis and affected follicular atresia, but these activities did not appear to be through COX-2 and Smad4. Taken together, we revealed an important CP2/miR-144/COX-2/PGE2/ovulation pathway in mGCs.

Developmental piRNA profiles of the invasive vector mosquito Aedes albopictus

Background

In eukaryotic organisms, Piwi-interacting RNAs (piRNAs) control the activities of mobile genetic elements and ensure genome maintenance. Recent evidence indicates that piRNAs are involved in multiple biological pathways, including transcriptional regulation of protein-coding genes, sex determination and even interactions between host and pathogens. Aedes albopictus is a major invasive species that transmits a number of viral diseases in humans. Ae. albopictus has the largest genome and the highest abundance of repetitive sequences when compared with members that belong to Culicidae with a published genome. Analysis of piRNA profiles will provide a developmental and evolutionary perspective on piRNAs in Ae. albopictus.

Methods

piRNAs were identified and characterized during the development of Ae. albopictus, and piRNA expression patterns in adult males and females as well as sugar-fed females and blood-fed females were compared.

Results

Our results reveal that, despite the large genome size of Ae. albopictus, the piRNA pool of Ae. albopictus (1.2 × 10⁷) is smaller than those of Aedes aegypti (1.7 × 10⁷) and Drosophila melanogaster (1.6 × 10⁷). In Ae. albopictus, piRNAs displayed the highest abundance at the embryo stage and the lowest abundance at the pupal stage. Approximately 50 % of the piRNAs mapped to intergenic regions with no known functions. Approximately 30 % of the piRNAs mapped to repetitive elements, and 77.69 % of these repeat-derived piRNAs mapped to Class I TEs; 45.42 % of the observed piRNA reads originated from piRNA clusters, and most of the top 10 highest expressed piRNA clusters and 100 highest expressed piRNAs from each stage displayed biased expression patterns across the developmental stages. All anti-sense-derived piRNAs displayed a preference for uridine at the 5′ end; however, the sense-derived piRNAs showed adenine bias at the tenth nucleotide position and a typical ping-pong signature, suggesting that the biogenesis of piRNAs was conserved throughout development. Our results also show that 962 piRNAs displayed sex-biased expression, and 522 piRNAs showed higher expression in the blood-fed females than in the sugar-fed females.

Conclusions

Our results suggest that piRNAs, aside from silencing transposable elements in Ae. albopictus, may have a role in other biological pathways.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1815-8) contains supplementary material, which is available to authorized users.

The transcription factor ccaat/enhancer binding protein β (C/EBPβ) and miR-27a regulate the expression of porcine Dickkopf2 (DKK2)

Using Affymetrix porcine Gene-Chip analyses, we found that Dickkopf2 (DKK2), a WNT antagonist, is differentially expressed in pre-ovulatory follicles between Large White and Chinese Taihu sows. This study aims to identify the regulatory factors responsible for DKK2 expression. Deletion fragment and mutation analyses identified DKK2-D3 as the porcine DKK2 core promoter. There were four C/EBPβ binding sites within the DKK2 core promoter. The C allele that results from a spontaneous alteration (DKK2 c.−1130 T > C) in the core promoter was associated with a higher total number born (TNB) and a higher number born alive (NBA) in all parities in a synthetic pig population. This was possibly the result of a change in C/EBPβ binding ability, which was confirmed using chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift assays (EMSA). Moreover, C/EBPβ specifically bound to and activated the DKK2 promoter, as revealed by mutation analysis, overexpression and RNA interference (RNAi) experiments. We also confirmed that miR-27a is a negative regulator of the DKK2 gene using miR-27a overexpression and inhibition experiments and mutation analyses. RTCA xCELLigence experiments showed that miR-27a suppressed Chinese hamster ovary (CHO) cell proliferation by down-regulating DKK2 gene expression. Taken together, our findings suggest that C/EBPβ and miR-27a control DKK2 transcription.

Identification and Analysis of Regulatory Elements in Porcine Bone Morphogenetic Protein 15 Gene Promoter

Bone morphogenetic protein 15 (BMP15) is secreted by the mammalian oocytes and is indispensable for ovarian follicular development, ovulation, and fertility. To determine the regulation mechanism of BMP15 gene, the regulatory sequence of porcine BMP15 was investigated in this study. The cloned BMP15 promoter retains the cell-type specificity, and is activated in cells derived from ovarian tissue. The luciferase assays in combination with a series of deletion of BMP15 promoter sequence show that the -427 to -376 bp region of BMP15 promoter is the primary regulatory element, in which there are a number of transcription factor binding sites, including LIM homeobox 8 (LHX8), newborn ovary homeobox gene (NOBOX), and paired-like homeodomain transcription factor 1 (PITX1). Determination of tissue-specific expression reveals that LHX8, but not PITX1 and NOBOX, is exclusively expressed in pig ovary tissue and is translocated into the cell nuclei. Overexpression of LHX8 in Chinese hamster ovary (CHO) cells could significantly promote BMP15 promoter activation. This study confirms a key regulatory element that is located in the proximal region of BMP15 promoter and is regulated by the LHX8 factor.

Proteomic Analysis of Fetal Ovary Reveals That Ovarian Developmental Potential Is Greater in Meishan Pigs than in Yorkshire Pigs

Time-dependent expression of functional proteins in fetal ovaries is important to understand the developmental process of the ovary. This study was carried out to enhance our understanding of the developmental process of porcine fetal ovaries and to better address the differences in fetal ovary development of local and foreign pigs. The objective of the present study is to test the expression of key proteins that regulate the growth and development of fetal ovaries in Meishan and Yorkshire porcine breeds by using proteomics technology. Six Meishan and 6 Yorkshire pregnant gilts were used in this experiment. Fetal ovaries were obtained from Yorkshire and Meishan gilts on days 55 and 90 of the gestation period. Using 2D-DIGE (two dimensional-difference in gel electrophoresis) analysis, the results showed that there are about 1551 and 1400 proteins in gilt fetal ovaries on days 55 and 90, respectively of the gestation. Using MALDI TOF-TOF MS analysis, 27 differentially expressed proteins were identified in the fetal ovaries of the 2 breeds on day 55 of gestation, and a total of 18 proteins were identified on day 90 of gestation. These differentially expressed proteins were involved in the regulation of biological processes (cell death, stress response, cytoskeletal proteins) and molecular functions (enzyme regulator activity). We also found that alpha-1-antitrypsin, actin, vimentin, and PP2A proteins promote the formation of primordial follicles in the ovaries of Yorkshire pigs on day 55 of gestation while low expression heat shock proteins and high expression alpha-fetoproteins (AFP) may promote Meishan fetal ovarian follicular development on day 90 of gestation. These findings provide a deeper understanding of how reduced expression of heat shock proteins and increased expression of AFP can significantly reduce the risk of reproductive disease in obese Meishan sows. Our study also shows how these proteins can increase the ovulation rate and may be responsible for the low reproductive efficiency reported in other obese breeds. The ovarian developmental potential was found to be greater in Meishan pigs than in Yorkshire pigs.

The novel porcine gene early growth response 4 (Egr4) is differentially expressed in the ovaries of Erhualian and Pietrain pigs

The early growth response 4 (Egr4) gene plays a critical role in human and mouse fertility. In the present study, Affymetrix microarray gene chips were used to evaluate differential gene expression in the ovaries between Erhualian and Pietrain pigs. In all, 487 and 573 transcripts were identified with significantly higher and lower expression, respectively, in Erhualian compared with Pietrain sows. The Egr4 gene, one of the differentially expressed genes, was cloned and its genomic structure was analysed. Egr4 expression is increased 120-fold in ovaries from Erhualian sows. The full-length cDNA of porcine Egr4 was obtained by in silico cloning and 5' rapid amplification of cDNA ends. The gene consists of two exons and its predicted protein contains a Cys2His2 zinc finger structure. The porcine transcript is alternatively spliced by exon sequence deletion, producing two different mRNAs differing at the 5' end of Exon 2. Egr4 transcripts were detected in the central nervous system, including the cerebrum, cerebellum, hypothalamus and pituitary gland, and were highly expressed in the ovary. The Egr4 gene was evaluated as a candidate gene for porcine reproductivity. To investigate the role of Egr4 in the ovary, Egr4 was knocked down using short interference (si) RNA in porcine granulosa cells. Knockdown of Egr4 using siRNA effectively inhibited Egr4 mRNA and protein expression and knockdown significantly affected the expression of Bax, P450arom, P450scc, Egr1, Egr2, and Egr3. In conclusion, these observations establish an important role for Egr4 in the porcine ovary.

Transcription factor C/EBPβ and 17β-estradiol promote transcription of the porcine p53 gene

The tumor protein 53 (p53) gene played a crucial role in maternal reproduction except its classic roles in maintaining genomic stability and preventing tumorigenesis. However, little is known concerning the regulatory elements which control the expression of p53 gene. In this study, we predicted two binding sites (-490/-477 and -405/-392) of transcription factor CCAAT/enhancer binding protein beta (C/EBPβ) within the core promoter (-985/-273) determined by promoter deletion analysis, and discovered that the second site (-405/-392) was important for p53 promoter activity by site-directed mutagenesis. Then the binding of C/EBPβ to the p53 promoter was identified by electrophoretic mobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP). Moreover, evidence from C/EBPβ overexpression and RNAi studies showed C/EBPβ regulated p53 promoter activity and endogenous p53 expression. Meanwhile, we observed p53 mRNA at the peak in 10(-6)mol/L 17β-estradiol treated cells for 24h via enhancing its core promoter activity. Taken together, our study indicates that C/EBPβ and 17β-estradiol are the essential regulatory factors for p53 transcription.

Associations of TCF12, CTNNAL1 and WNT10B gene polymorphisms with litter size in pigs

In previous research, several WNT signaling pathway genes including transcription factor 12 (TCF12), catenin alpha-like protein 1 (CTNNAL1) and wingless-type MMTV integration site family, member 10B (WNT10B) were differentially expressed in PMSG-hCG stimulated preovulatory ovarian follicles of Large White and Chinese Taihu sows. In the present research, these three genes were selected as the candidate genes for litter size traits in pigs. Four mutations (TCF12 c.-201+65 G>A, TCF12 c.-200-300 G>A, CTNNAL1 c.1878 G>C and WNT10B c.*12 C>T) were detected in eleven pig populations, and results indicated CTNNAL1 c.1878 G and WNT10B c.*12 C were the major alleles in all tested pig populations, while TCF12 c.-201+65 A and TCF12 c.-200-300 A were the major alleles in several Chinese native pig breeds. Association analysis of four mutations with litter size in Large White and DIV pigs showed that both the signficant differences of total number born (TNB) and number born alive (NBA) among three genotypes and the significance of additive effects appeared at TCF12 c.-200-300 G>A and CTNNAL1 c.1878 G>C loci, suggesting these two mutations might be reliable markers for pig selection and breeding.

Identification of the promoter region and genetic mutations of the porcine GALP gene

Galanin-like peptide (GALP) gene, encoding a member of the galanin family of neuropeptides involved in reproduction, was differentially expressed in PMSG-hCG stimulated pre-ovulatory ovarian follicles of Chinese Taihu and Large White sows in our previous study. In the present study, promoter region and genetic mutations of the porcine GALP gene were determined. A 1,322 bp contig in 5'-flanking region was predicted to contain 5 potential transcription promoters by Neural Network Promoter Prediction version 2.2. 5'-deletion expression in both CHO and hela cells showed that there were a negative regulatory element at -852 to -803 bp and a positive regulatory element at -1,318 to -1,269 bp. Comparative sequence analyses of Chinese Taihu and Large White GALP gene sequence revealed the c.*27C>G mutation in the 3'-UTR and the c.88-1225C>G mutation in intron 1, which can be detected by HhaI and AluI PCR-RFLP, respectively. The association analysis with litter size traits showed that at both loci CC and GG genotypes were different for NBA for all parities in DIV pigs (P < 0.05). However, two SNPs were not in significant linkage disequilibrium analyzed using SHEsis online software, and could be used in pig breeding individually.

Development of a porcine (Sus scofa) embryo-specific microarray: array annotation and validation

BACKGROUND

The domestic pig is an important livestock species and there is strong interest in the factors that affect the development of viable embryos and offspring in this species. A limited understanding of the molecular mechanisms involved in early embryonic development has inhibited our ability to fully elucidate these factors. Next generation deep sequencing and microarray technologies are powerful tools for delineation of molecular pathways involved in the developing embryo.

RESULTS

Here we present the development of a porcine-embryo-specific microarray platform created from a large expressed sequence tag (EST) analysis generated by Roche/454 next-generation sequencing of cDNAs constructed from critical stages of in vivo or in vitro porcine preimplantation embryos. Two cDNA libraries constructed from in vitro and in vivo produced preimplantation porcine embryos were normalized and sequenced using 454 Titanium pyrosequencing technology. Over one million high-quality EST sequences were obtained and used to develop the EMbryogene Porcine Version 1 (EMPV1) microarray composed of 43,795 probes. Based on an initial probe sequence annotation, the EMPV1 features 17,409 protein-coding, 473 pseudogenes, 46 retrotransposed, 2,359 non-coding RNA, 4,121 splice variants in 2,862 genes and a total of 12,324 Novel Transcript Regions (NTR). After re-annotation, the total unique genes increased from 11,961 to 16,281 and 1.9% of them belonged to a large olfactory receptor (OR) gene family. Quality control on the EMPV1 was performed and revealed an even distribution of ten clusters of spiked-in control spots and array to array (dye-swap) correlation was 0.97.

CONCLUSIONS

Using next-generation deep sequencing we have produced a large EST dataset to allow for the selection of probe sequences for the development of the EMPV1 microarray platform. The quality of this embryo-specific array was confirmed with a high-level of reproducibility using current Agilent microarray technology. With more than an estimated 20,000 unique genes represented on the EMPV1, this platform will provide the foundation for future research into the in vivo and in vitro factors that affect the viability of porcine embryos, as well as the effects of these factors on the live offspring that result from these embryos.

Coordinated miRNA/mRNA expression profiles for understanding breed-specific metabolic characters of liver between Erhualian and large white pigs

MicroRNAs (miRNAs) are involved in the regulation of various metabolic processes in the liver, yet little is known on the breed-specific expression profiles of miRNAs in coordination with those of mRNAs. Here we used two breeds of male newborn piglets with distinct metabolic characteristics, Large White (LW) and Erhualian (EHL), to delineate the hepatic expression profiles of mRNA with microarray and miRNAs with both deep sequencing and microarray, and to analyze the functional relevance of integrated miRNA and mRNA expression in relation to the physiological and biochemical parameters. EHL had significantly lower body weight and liver weight at birth, but showed elevated serum levels of total cholesterol (TCH), high-density lipoprotein cholesterol (HDLC) and low-density lipoprotein cholesterol (LDLC), as well as higher liver content of cholesterol. Higher serum cortisol and lower serum insulin and leptin were also observed in EHL piglets. Compared to LW, 30 up-regulated and 18 down-regulated miRNAs were identified in the liver of EHL, together with 298 up-regulated and 510 down-regulated mRNAs (FDR<10%). RT-PCR validation of some differentially expressed miRNAs (DEMs) further confirmed the high-throughput data analysis. Using a target prediction algorithm, we found significant correlation between the up-regulated miRNAs and down-regulated mRNAs. Moreover, differentially expressed genes (DEGs), which are involved in proteolysis, were predicted to be mediated by DEMs. These findings provide new information on the miRNA and mRNA profiles in porcine liver, which would shed light on the molecular mechanisms underlying the breed-specific traits in the pig, and may serve as a basis for further investigation into the biological functions of miRNAs in porcine liver.

Data mining in networks of differentially expressed genes during sow pregnancy

Small to moderate gains in Pig fertility can mean large returns in overall efficiency, and developing methods to improve it is highly desirable. High fertility rates depend on completion of successful pregnancies. To understand the molecular signals associated with pregnancy in sows, expression profiling experiments were conducted to identify differentially expressed genes in ovary and myometrium at different pregnancy periods using the Affymetrix Porcine GeneChip^TM. A total of 974, 1800, 335 and 710 differentially expressed transcripts were identified in the myometrium during early pregnancy (EP) and late pregnancy (LP), and in the ovary during EP and LP, respectively. Self-Organizing Map (SOM) clusters indicated the differentially expressed genes belonged to 7 different functional groups. Based on BLASTX searches and Gene Ontology (GO) classifications, 129 unique genes closely related to pregnancy showed differential expression patterns. GO analysis also indicated that there were 21 different molecular function categories, 20 different biological process categories, and 8 different cellular component categories of genes differentially expressed during sow pregnancy. Gene regulatory network reconstruction provided us with an interaction model of known genes such as insulin-like growth factor 2 (IGF2) gene, estrogen receptor (ESR) gene, retinol-binding protein-4 (RBP4) gene, and several unknown candidate genes related to reproduction. Several pitch point genes were selected for association study with reproduction traits. For instance, DPPA5 g.363 T>C was found to associate with litter born weight at later parities in Beijing Black pigs significantly (p < 0.05). Overall, this study contributes to elucidating the mechanism underlying pregnancy processes, which maybe provide valuable information for pig reproduction improvement.

Discovery of potential piRNAs from next generation sequences of the sexually mature porcine testes

Piwi-interacting RNAs (piRNAs), a new class of small RNAs discovered from mammalian testes, are involved in transcriptional silencing of retrotransposons and other genetic elements in germ line cells. In order to identify a full transcriptome set of piRNAs expressed in the sexually mature porcine testes, small RNA fractions were extracted and were subjected to a Solexa deep sequencing. We cloned 6,913,561 clean reads of Sus Scrofa small RNAs (18-30 nt) and performed functional characterization. Sus Scrofa small RNAs showed a bimodal length distribution with two peaks at 21 nt and 29 nt. Then from 938,328 deep-sequenced small RNAs (26-30 nt), 375,195 piRNAs were identified by a k-mer scheme and 326 piRNAs were identified by homology searches. All piRNAs predicted by the k-mer scheme were then mapped to swine genome by Short Oligonucleotide Analysis Package (SOAP), and 81.61% of all uniquely mapping piRNAs (197,673) were located to 1124 defined genomic regions (5.85 Mb). Within these regions, 536 and 501 piRNA clusters generally distributed across only minus or plus genomic strand, 48 piRNA clusters distributed on two strands but in a divergent manner, and 39 piRNA clusters distributed on two strands in an overlapping manner. Furthermore, expression pattern of 7 piRNAs identified by homology searches showed 5 piRNAs displayed a ubiquitous expression pattern, although 2 piRNAs were specifically expressed in the testes. Overall, our results provide new information of porcine piRNAs and their specific expression pattern in porcine testes suggests that piRNAs have a role in regulating spermatogenesis.

Exploring the application of high-throughput genomics technologies in the field of maternal-embryo communication

Deciphering the complex molecular dialogue between the maternal tract and embryo is crucial to increasing our understanding of pregnancy failure, infertility problems and in the modulation of embryo development, which has consequences through adulthood. High-throughput genomic technologies have been applied to look for a holistic view of the molecular interactions occurring during this dialogue. Among these technologies, microarrays have been widely used, being one of the most popular tools in maternal-embryo communication. Today, next generation sequencing technologies are dwarfing the capabilities of microarrays. The application of these new technologies has broadened to almost all areas of genomics research, because of their massive sequencing capacity. We review the current status of high-throughput genomic technologies and their application to maternal-embryo communication research. We also survey next generation technologies and their huge potential in many research areas. Given the diversity of unanswered questions in the field of maternal-embryo communication and the wide range of possibilities that these technologies offer, here we discuss future perspectives on the use of these technologies to enhance maternal-embryo research.

Src homology domain-containing phosphatase 2 suppresses cellular senescence in glioblastoma

BACKGROUND

Epidermal growth factor receptor (EGFR) signalling is frequently altered during glioblastoma de novo pathogenesis. An important downstream modulator of this signal cascade is SHP2 (Src homology domain-containing phosphatase 2).

METHODS

We examined the The Cancer Genome Atlas (TCGA) database for SHP2 mutations. We also examined the expression of a further 191 phosphatases in the TCGA database and used principal component and comparative marker analysis available from the Broad Institute to recapitulate the TCGA-defined subgroups and identify the specific phosphatases defining each subgroup. We identified five siRNAs from two independent commercial sources that were reported by the vendor to be pre-optimised in their specificity of SHP2 silencing. The specificity and physiological effects of these siRNAs were tested using an in vitro glioma model.

RESULTS

TCGA data demonstrate SHP2 to be mutated in 2% of the glioblastoma multiforme's studied. Both mutations identified in this study are likely to be activating mutations. We found that the four subgroups of GBM as defined by TCGA differ significantly with regard to the expression level of specific phosphatases as revealed by comparative marker analysis. Surprisingly, the four subgroups can be defined solely on the basis of phosphatase expression level by principal component analysis. This result suggests that critical phosphatases are responsible for the modulation of specific molecular pathways within each subgroup. Src homology domain-containing phosphatase 2 constitutes one of the 12 phosphatases that define the classical subgroup. We confirmed the biological significance by siRNA knockdown of SHP2. All five siRNAs tested reduced SHP2 expression by 70-100% and reduced glioblastoma cell line growth by up to 80%. Profiling the established molecular targets of SHP2 (ERK1/2 and STAT3) confirmed specificity of these siRNAs. The loss of cell viability induced by SHP2 silencing could not be explained by a significant increase in apoptosis alone as demonstrated by terminal deoxyribonucleotidyl transferase-mediated nick-end labelling and propidium iodide staining. Src homology domain-containing phosphatase 2 silencing, however, did induce an increase in β-galactosidase staining. Propidium iodide staining also showed that SHP2 silencing increases the population of glioblastoma cells in the G1 phase of the cell cycle and reduces the population of such cells in the G2/M- and S-phase.

CONCLUSION

Src homology domain-containing phosphatase 2 promotes the growth of glioblastoma cells by suppression of cellular senescence, a phenomenon not described previously. Selective inhibitors of SHP2 are commercially available and may be considered as a strategy for glioblastoma therapy.

Analysis of muscle and ovary transcriptome of Sus scrofa: assembly, annotation and marker discovery

Pig (Sus scrofa) is an important organism for both agricultural and medical purpose. This study aims to investigate the S. scrofa transcriptome by the use of Roche 454 pyrosequencing. We obtained a total of 558 743 and 528 260 reads for the back-leg muscle and ovary tissue each. The overall 1 087 003 reads give rise to 421 767 341 bp total residues averaging 388 bp per read. The de novo assemblies yielded 11 057 contigs and 60 270 singletons for the back-leg muscle, 12 204 contigs and 70 192 singletons for the ovary and 18 938 contigs and 102 361 singletons for combined tissues. The overall GC content of S. scrofa transcriptome is 42.3% for assembled contigs. Alternative splicing was found within 4394 contigs, giving rise to 1267 isogroups or genes. A total of 56 589 transcripts are involved in molecular function (40 916), biological process (38 563), cellular component (35 787) by further gene ontology analyses. Comparison analyses showed that 336 and 553 genes had significant higher expression in the back-leg muscle and ovary each. In addition, we obtained a total of 24 214 single-nucleotide polymorphisms and 11 928 simple sequence repeats. These results contribute to the understanding of the genetic makeup of S. scrofa transcriptome and provide useful information for functional genomic research in future.