Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo.

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

Comprehensive variation discovery and recovery of missing sequence in the pig genome using multiple de novo assemblies.

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

Comparison of Y-chromosome-linked TSPY, TSPY2, and PRAMEY genes in Taurus cattle, yaks, and interspecific hybrid bulls.

Domestic yaks (Bos grunniens) and domestic Taurus cattle (Bos taurus) are closely related. An interesting phenomenon in interspecific crossings is male sterility in the F1 hybrid (yattle) and F2 backcross, with no late meiotic cells or spermatids in the seminiferous tubules. The mammalian Y chromosome is crucial for spermatogenesis and male fertility. This study investigated the copy number variations and mRNA of Y-transitional region genes TSPY2 (testis specific protein, Y-linked 2 and testis-specific Y-encoded protein 3-like) and PRAMEY (preferentially expressed antigen in melanoma, Y-linked), and Y-ampliconic region genes TSPY (testis-specific Y-encoded protein 1-like), ZNF280BY (zinc finger protein 280B, Y-linked) and HSFY (heat-shock transcription factor, Y-linked) in mature testes from Taurus cattle, yaks, and yattle. Phylogenetic trees divided 33 copies of TSPY into major 2 types (TSPY-T1 and TSPY-T2), 19 copies of TSPY2 into 2 types (TSPY2-T1 and T2), and 8 copies of PRAMEY into 4 types (PRAMEY-T1 to T4). Searching by the Basic Local Alignment Search Tool of the TSPY2 coding sequences in GenBank revealed that TSPY2 was conserved in Bovidae. The TSPY2-T2 sequences were absent, whereas PRAMEY-T2 and PRAMEY-T4 were amplified on the yak Y chromosome. The average copy numbers of TSPY-T2 and ZNF280BY were significantly different between cattle and yaks. The TSPY-T2, TSPY2, PRAMEY, ZNF280BY, and HSFY genes were uniquely or predominantly expressed in testes. Reverse-transcription quantitative PCR showed that the TSPY-T2, PRAMEY-T2, HSFY, ZNF280BY, protamine 1 (PRM1), and protamine 2 (PRM2) genes were almost not expressed in yattle. The PRM1 and PRM2 genes are used as positive markers for spermatozoa. Thus, our results showed that the genomic structure of the Y-transitional and Y-ampliconic region differed between Taurus cattle and yaks. Dysregulated expression of Y-ampliconic region genes TSPY-T2, HSPY, ZNF280BY, and Y-transitional region gene PRAMEY-T2 may be associated with hybrid male sterility in yattle.

Spontaneous single nucleotide polymorphism in porcine microRNA-378 seed region leads to functional alteration.

Sequence variation in a microRNA (miRNA) seed region can influence its biogenesis and effects on target mRNAs; however, in mammals, few seed region mutations leading to functional alterations have been reported to date. Here, we report the identification of a single nucleotide polymorphism (SNP) with functional consequence located in the seed region of porcine miR-378. In vitro analysis of this rs331295049 A17G SNP showed significantly up-regulated expression of the mature miR-378 (miR-378/G). In silico target prediction indicated that the SNP would modulate secondary structure and result in functional loss affecting >85% of the known target genes of the wild-type miR-378 (miR-378/A), and functional gain affecting >700 new target genes, and dual-luciferase reporter assay verified this result. This report of a SNP in the seed region of miR-378 leads to functional alteration and indicates the potential for substantive functional consequences to the molecular physiology of a mammalian organism.

microRNA profiling in three main stages during porcine spermatogenesis.

BACKGROUND: Spermatogenesis is an intricate biological event wherein an undifferentiated spermatogonium develops into mature sperms. MicroRNAs are a type of single strand small non-coding RNA molecule and are implicated in the regulation of many crucial pathways during cell proliferation, apoptosis, and differentiation. METHOD: Here, we present a comprehensive comparison of miRNA expression profiling in three main stages during porcine spermatogenesis using high-throughput sequencing. RESULTS: We built three small RNA libraries for the testis, the epididymis and the ejaculated sperm from a Landrace boar, and in total obtained 3821 precursor hairpins encoding for 4761 mature miRNAs, of which 23 are miRNA*. Notably, 940 precursor miRNAs produced both the 5'- and 3'- strands as sister pairs, indicating the distinctive expression patterns of germ cell miRNAs. Additionally, 418 out of 710 co-expressed miRNAs were identified as being differentially expressed between libraries (P < 0.001). Apart from the sexual specific X chromosome, many miRNAs were found to be located on chromosome 12, which may play potential roles in spermatogenesis according to the result of synteny analysis with human and mouse. The Gene Ontology and KEGG pathway analysis revealed that the target genes of co-expressed miRNAs were highly involved in the cell cycle process, metal ion binding, modification of plasma membrane, and the p53 signal pathway.

Effects of taurine on metabolomics of bovine mammary epithelial cells under high temperature conditions.

High temperature induces heat stress, adversely affecting the growth and lactation performance of cows. Research has shown the protective effect of taurine against hepatotoxicity both in vivo and in vitro. This study aimed to investigate the effect of taurine on the metabolomics of mammary epithelial cells of dairy cows under high-temperature conditions. Mammary epithelial cells were exposed to 0 mmol/L (HS, control), 8 mmol/L (HT-8), and 32 mmol/L (HT-32) of taurine, then incubated at 42°C for 6 h. Metabolomics analysis was conducted using Liquid Chromatograph Mass Spectrometer (LC-MS). Compared with the HS group, 2,873 and 3,243 metabolites were detected in the HT-8 group in positive and negative ion modes. Among these, 108 and 97 metabolites were significantly upregulated in positive and negative ion modes, while 60 and 166 metabolites were downregulated. Notably, 15 different metabolites such as palmitic acid, adenine and hypoxanthine were screened out in the HT-8 group. Compared with the HS group, 2,873 and 3,243 metabolites were, respectively, detected in the HT-32 group in the positive and negative ion modes. Among those metabolites, 206 metabolites were significantly up-regulated, while 206 metabolites were significantly downregulated in the positive mode. On the other hand, 497 metabolites were significantly upregulated in the negative mode, while 517 metabolites were reported to be downregulated. Noteworthy, 30 distinct metabolites, such as palmitic acid, phytosphingosine, hypoxanthine, nonanoic acid, and octanoic acid, were screened out in the HT-32 group. KEGG enrichment analysis showed that these metabolites were mainly involved in lipid metabolism, purine metabolism and other biological processes. Overall, our study indicates that taurine supplementation alters the metabolites primarily associated with purine metabolism, lipid metabolism and other pathways to alleviate heat stress in bovine mammary epithelial cells.

Promoter hypermethylation of PIWI/piRNA pathway genes associated with diminished pachytene piRNA production in bovine hybrid male sterility.

Hybrid male sterility (HMS) is a postzygotic reproductive isolation mechanism that enforces speciation. A bovine example of HMS is the yattle (also called dzo), an interspecies hybrid of taurine cattle (Bos taurus) and yak (Bos grunniens). The molecular mechanisms underlying HMS of yattle are not well understood. Epigenetic modifications of DNA methylation and P-element induced wimpy testis (PIWI)-interacting RNA (piRNAs) are important regulators in spermatogenesis. In this study, we investigated DNA methylation patterns and piRNA expression in adult testes in hybrid infertile yattle bulls and fertile cattle and yak bulls using whole genome bisulphite-seq and small RNA-seq. Promoter hypermethylation in yattle were associated with DNA methylation involved in gamete generation, piRNA metabolic processes, spermatogenesis, and spermatid development (P < 2.6 × 10-5). Male infertility in yattle was associated with the promoter hypermethylation-associated silencing of PIWI/piRNA pathway genes including PIWIL1, DDX4, PLD6, MAEL, FKBP6, TDRD1 and TDRD5. The downstream effects of silencing these genes were diminished production of 29- to 31- nucleotide pachytene piRNAs in yattle testes. Hypermethylation events at transposable element loci (LINEs, SINEs, and LTRs) were found in yattle. LINE-derived prepachytene piRNAs increased and SINE-derived prepachytene piRNAs were reduced in yattle testes. Our data suggests that DNA methylation affects the PIWI/piRNA pathway and is involved in gene expression and pachytene piRNA production during spermatogenesis in bovine HMS. DNA hypermethylation and disruption of piRNA production contributed to unsuccessful germ cell development that may drive bovine HMS.

Meat and nutritional quality comparison of purebred and crossbred pigs.

Crossbreeding is an effective method of improving the efficiency and profit of production in commercial pig operations. To understand the effect of crossbreeding on meat and nutrient quality, a combination including three purebred (Duroc, D; Landrace, L; Yorkshire, Y) and two crossbred pig lines (Landrace × Yorkshire, LY; Duroc × (Landrace × Yorkshire), DLY) frequently used internationally were studied. The results showed that meat from the LY and DLY crosses had lower values for lightness L24h∗, shear force and epinephrine and higher values for drip loss, C18:1, insulin, glucagon and monounsaturated fatty acids than D, L and Y pigs. Moreover, LY had higher values for post mortem pH and lower values for a* and b* than the purebreds. In contrast, DLY had lower values for pH and higher values for a* and b* than the purebreds. Meat quality-related gene analysis showed that the CAST, IGF2 and MC4R gene expression levels in the LY and DLY pigs were significantly higher than those in the D, L and Y pigs. These results indicate that crossbreeding can alter the meat quality, nutritive value, energy metabolism and gene expression of pigs. Future research should focus on microRNA expression and DNA methylation that regulate gene expression and thus affect the meat quality.

A mutation in porcine pre-miR-15b alters the biogenesis of MiR-15b\16-1 cluster and strand selection of MiR-15b.

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in translational regulation of the messenger RNA molecules. Sequence variations in the genes encoding miRNAs could influence their biogenesis and function. MiR-15b plays an important role in cellular proliferation, apoptosis and the cell cycle. Here, we report the identification of a C58T mutation in porcine pre-miR-15b. Through in vitro and in vivo experiments, we determined that this mutation blocks the transition from pri-miRNA to pre-miRNA, alters the strand selection between miR-15b-5p and miR-15b-3p, and obstructs biogenesis of the downstream miR-16-1. These results serve to highlight the importance of miRNA mutations and their impacts on miRNA biogenesis.

Dynamic gene expression profiles during postnatal development of porcine subcutaneous adipose.

A better understanding of the control of lipogenesis is of critical importance for both human and animal physiology. This requires a better knowledge of the changes of gene expression during the process of adipose tissue development. Thus, the objective of the current study was to determine the effects of development on subcutaneous adipose tissue gene expression in growing and adult pigs. Here, we present a comprehensive investigation of mRNA transcriptomes in porcine subcutaneous adipose tissue across four developmental stages using digital gene expression profiling. We identified 3,274 differential expressed genes associated with oxidative stress, immune processes, apoptosis, energy metabolism, insulin stimulus, cell cycle, angiogenesis and translation. A set of universally abundant genes (ATP8, COX2, COX3, ND1, ND2, SCD and TUBA1B) was found across all four developmental stages. This set of genes may play important roles in lipogenesis and development. We also identified development-related gene expression patterns that are linked to the different adipose phenotypes. We showed that genes enriched in significantly up-regulated profiles were associated with phosphorylation and angiogenesis. In contrast, genes enriched in significantly down-regulated profiles were related to cell cycle and cytoskeleton organization, suggesting an important role for these biological processes in adipose growth and development. These results provide a resource for studying adipose development and promote the pig as a model organism for researching the development of human obesity, as well as being used in the pig industry.

Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars.

We report the sequencing at 131× coverage, de novo assembly and analyses of the genome of a female Tibetan wild boar. We also resequenced the whole genomes of 30 Tibetan wild boars from six major distributed locations and 18 geographically related pigs in China. We characterized genetic diversity, population structure and patterns of evolution. We searched for genomic regions under selection, which includes genes that are involved in hypoxia, olfaction, energy metabolism and drug response. Comparing the genome of Tibetan wild boar with those of neighboring Chinese domestic pigs further showed the impact of thousands of years of artificial selection and different signatures of selection in wild boar and domestic pig. We also report genetic adaptations in Tibetan wild boar that are associated with high altitudes and characterize the genetic basis of increased salivation in domestic pig.

An atlas of DNA methylomes in porcine adipose and muscle tissues.

It is evident that epigenetic factors, especially DNA methylation, have essential roles in obesity development. Here, using pig as a model, we investigate the systematic association between DNA methylation and obesity. We sample eight variant adipose and two distinct skeletal muscle tissues from three pig breeds living within comparable environments but displaying distinct fat level. We generate 1,381 Gb of sequence data from 180 methylated DNA immunoprecipitation libraries, and provide a genome-wide DNA methylation map as well as a gene expression map for adipose and muscle studies. The analysis shows global similarity and difference among breeds, sexes and anatomic locations, and identifies the differentially methylated regions. The differentially methylated regions in promoters are highly associated with obesity development via expression repression of both known obesity-related genes and novel genes. This comprehensive map provides a solid basis for exploring epigenetic mechanisms of adipose deposition and muscle growth.

Repertoire of porcine microRNAs in adult ovary and testis by deep sequencing.

BACKGROUND: MicroRNAs (miRNAs), a large family of short endogenous RNAs known to post-transcriptionally repress gene expression, participate in the regulation of almost every cellular process. Changes in miRNA expression are associated with many pathologies. Ovarian folliculogenesis and testicular spermatogenesis are complex and coordinated biological processes, in which tightly regulated expression and interaction of a multitude of genes could be regulated by these miRNAs. Identification and preliminary characterization of gonad-specific miRNAs would be a prerequisite for a thorough understanding of the role that miRNA-mediated posttranscriptional gene regulation plays in mammalian reproduction. METHOD: Here, we present the identification of a repertoire of porcine miRNAs in adult ovary and testis using deep sequencing technology. A bioinformatics pipeline was developed to distinguish authentic mature miRNA sequences from other classes of small RNAs represented in the sequencing data. RESULTS: Using this approach, we detected 582 precursor hairpins (pre-miRNAs) encoding for 732 mature miRNAs, of which 673 are unique. Statistically, 224 unique miRNAs (out of 673, 33.28%) were identified which had significant differential expression (DE) between ovary and testis libraries (P < 0.001). Most of DE miRNAs located on the X chromosome (X-linked miRNAs) (24 out of 34, 70.59%) significantly up-regulated in ovary versus testis (P < 0.001). Predictably, X-linked miRNAs are expressed in a testis-preferential or testis-specific pattern. To explore the potential for co-expression among genomic location clusters of X-linked miRNAs, we surveyed the relationship between the distance separating miRNA loci and the coordinate expression patterns of 32 high confidence X-linked miRNAs in seven normal pig tissues using the real-time quantitative PCR (q-PCR) approach. Our results show that proximal pairs of miRNAs are generally co-expressed implying that miRNAs within 50 kb of genomic bases are typically derived from a common transcript. CONCLUSIONS: The present study characterizes the miRNA transcriptome of adult porcine gonads, with an emphasis on the co-expression patterns of X-linked miRNAs. Our report should facilitate studies of the organ-specific reproductive roles of miRNAs.

MicroRNAome of porcine pre- and postnatal development.

The domestic pig is of enormous agricultural significance and valuable models for many human diseases. Information concerning the pig microRNAome (miRNAome) has been long overdue and elucidation of this information will permit an atlas of microRNA (miRNA) regulation functions and networks to be constructed. Here we performed a comprehensive search for porcine miRNAs on ten small RNA sequencing libraries prepared from a mixture of tissues obtained during the entire pig lifetime, from the fetal period through adulthood. The sequencing results were analyzed using mammalian miRNAs, the precursor hairpins (pre-miRNAs) and the first release of the high-coverage porcine genome assembly (Sscrofa9, April 2009) and the available expressed sequence tag (EST) sequences. Our results extend the repertoire of pig miRNAome to 867 pre-miRNAs (623 with genomic coordinates) encoding for 1,004 miRNAs, of which 777 are unique. We preformed real-time quantitative PCR (q-PCR) experiments for selected 30 miRNAs in 47 tissue-specific samples and found agreement between the sequencing and q-PCR data. This broad survey provides detailed information about multiple variants of mature sequences, precursors, chromosomal organization, development-specific expression, and conservation patterns. Our data mining produced a broad view of the pig miRNAome, consisting of miRNAs and isomiRs and a wealth of information of pig miRNA characteristics. These results are prelude to the advancement in pig biology as well the use of pigs as model organism for human biological and biomedical studies.