Exploring the gonad transcriptome of two extreme male pigs with RNA-seq

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

Comparative transcriptomic analysis of gonadal development and renewal in the ovoviviparous black rockfish (Sebastes schlegelii)

BACKGROUND

The black rockfish (Sebastes schlegelii) has an ovoviviparous reproductive pattern and long-term sperm storage, resulting in asynchronous gonadal development between the sexes. However, the comprehensive understanding of gonadal development in black rockfish has not yet been achieved. Here, we studied gonadal development and germ cell renewal using histology and RNA-seq.

RESULTS

In this study, RNA-seq was performed on testes and ovaries to characterize key pathways and genes that are active during development and gamete maturation in black rockfish. Differentially expressed genes (DEGs) were identified and annotated in 4 comparisons (F_III vs. F_IV, F_IV vs. F_V, M_III vs. M_IV and M_IV vs. M_V). Based on analysis of DEGs enriched in the testis, 11 and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped to the M_III vs. M_IV group and the M_IV vs. M_V group, respectively. DEGs in ovarian development were also classified into 10 groups according to their biological functions. The expression patterns of the selected genes determined by qPCR were significantly correlated with the RNA-Seq results, supporting the reliability and accuracy of the RNA-Seq analysis. E2 levels showed down regulation from previtellogenesis to mature stage in female and T level showed down regulation from spermatogenesis to regressed stage in the male.

CONCLUSIONS

The categories "intercellular interaction and cytoskeleton", "molecule amplification" and "repair in the cell cycle" were revealed to be crucial in testis development and spermatogenesis, as was the biosynthesis of a series of metabolites. Our results provide comprehensive insight into black rockfish gonadal development and provide a basis for further study of reproductive physiology and molecular biology in ovoviviparity teleosts.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Excavation of Genes Related to the Mining of Growth, Development, and Meat Quality of Two Crossbred Sheep Populations Based on Comparative Transcriptomes

Simple Summary

In this study, we measured the performance parameters of two crossbred sheep breeds, using Masson staining of the muscle tissue, and using the Illumina high-throughput sequencing platform to determine the differentially expressed genes (DEGs) in Dorper (DP) × Small-tailed Han (STH) sheep and Mongolia (MG) × Small-tailed Han sheep (STH). New transcripts of the muscle transcriptome were examined for the first time. DP × STH sheep were superior to MG × STH sheep in terms of meat quality and muscle morphology. In addition, 13 DEGs were found to play important roles in growth, development, and meat quality. The findings of this work may provide valuable resources for future research on muscle development in sheep.

Abstract

Crossbreeding can improve production performance and meat quality in sheep. The objective of this study was to look for genes related to sheep growth, development, and muscle. In this study, Dorper (DP) × Small Tailed Han (STH) sheep and Mongolia (MG) × Small-tailed Han (STH) sheep were used to estimate the productive performance and meat quality in a crossbreed. Subsequently, transcriptome analysis and bioinformatic analysis were performed on the Longissimus dorsi muscles of DP × STH and MG × STH sheep to identify differentially expressed genes (DEGs) related to growth, development, and meat quality. The presence of DEGs was confirmed by real-time PCR (qPCR). Productive performance and meat quality of the DP × STH sheep were better than the MG × STH sheep. Compared to DP × STH, a total of 1445 DEGs were identified in MG × STH sheep (1026 DEG were up-regulated and 419 DEG were down-regulated). Of these, 38 DEGs were related to growth, 161 to development, and 43 to muscle. In addition, 13 co-expressed genes (FGFRL1, SIX1, PLCB1, CRYAB, MYL2, ADIPOQ, GPX1, PPARD, GPC1, CDC42, LOC101106246, IGF1, and LARGE) were identified. The expression of DEGs was consistent with the comparative transcriptome analysis. This work provides genetics resources for future research on muscle development in sheep.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Effects of Chinese wolfberry and Astragalus extract on the antioxidant capacity of Tibetan pig liver

The objective of this study is to determine the effect of Chinese wolfberry (Lycium barbarum) and Astragalus (Astragalus membranaceus) extract (WAE) on the antioxidant capacity of Tibetan pig liver, and discussed the regulatory effect of WAE on the liver antioxidant mechanism. Twelve healthy 120-day-old Tibetan black pigs (35±2 kg) were divided randomly into two groups. The WAE group was fed a basal diet supplemented with 1% WAE for 90 days. The control group was fed the same diet, but without the WAE. We found that liver superoxide dismutase 1 (SOD1) activity (P<0.05), total antioxidative capacity (T-AOC) (P<0.05), and catalase (CAT) activity (P<0.01) significantly increased in the WAE group compared with the control group; malondialdehyde (MDA) content decreased, but this was not significant (P >0.05). Transcriptome sequencing analysis detected 106 differentially expressed genes (DEGs) related to oxidative stress. GO enrichment analysis showed these DEGs were involved in the positive regulation of reactive oxygen metabolism and biosynthesis, process regulation, and regulation of the oxidative stress response. KEGG Pathway enrichment analysis showed they were enriched in the PI3K-Akt, AMPK, Rap1, and peroxisome signaling pathways. The expression levels of key peroxisome biosynthesis genes (e.g., PEX3 and PEX11B) and key antioxidant genes (e.g., CAT and SOD1) were significantly higher in the WAE group than in the control group. The PRDX1 and PRDX5 content also was significantly higher in the WAE group. This study showed that the WAE regulated the antioxidant and anti-stress ability of Tibetan pig liver through a "peroxisome antioxidant-oxidant stress" signaling pathway.

De novo transcriptome assembly from the gonads of a scleractinian coral, Euphyllia ancora: molecular mechanisms underlying scleractinian gametogenesis

BACKGROUND

Sexual reproduction of scleractinians has captured the attention of researchers and the general public for decades. Although extensive ecological data has been acquired, underlying molecular and cellular mechanisms remain largely unknown. In this study, to better understand mechanisms underlying gametogenesis, we isolated ovaries and testes at different developmental phases from a gonochoric coral, Euphyllia ancora, and adopted a transcriptomic approach to reveal sex- and phase-specific gene expression profiles. In particular, we explored genes associated with oocyte development and maturation, spermiogenesis, sperm motility / capacitation, and fertilization.

RESULTS

1.6 billion raw reads were obtained from 24 gonadal samples. De novo assembly of trimmed reads, and elimination of contigs derived from symbiotic dinoflagellates (Symbiodiniaceae) and other organisms yielded a reference E. ancora gonadal transcriptome of 35,802 contigs. Analysis of 4 developmental phases identified 2023 genes that were differentially expressed during oogenesis and 678 during spermatogenesis. In premature/mature ovaries, 631 genes were specifically upregulated, with 538 in mature testes. Upregulated genes included those involved in gametogenesis, gamete maturation, sperm motility / capacitation, and fertilization in other metazoans, including humans. Meanwhile, a large number of genes without homology to sequences in the SWISS-PROT database were also observed among upregulated genes in premature / mature ovaries and mature testes.

CONCLUSIONS

Our findings show that scleractinian gametogenesis shares many molecular characteristics with that of other metazoans, but it also possesses unique characteristics developed during cnidarian and/or scleractinian evolution. To the best of our knowledge, this study is the first to create a gonadal transcriptome assembly from any scleractinian. This study and associated datasets provide a foundation for future studies regarding gametogenesis and differences between male and female colonies from molecular and cellular perspectives. Furthermore, our transcriptome assembly will be a useful reference for future development of sex-specific and/or stage-specific germ cell markers that can be used in coral aquaculture and ecological studies.

Transcriptome analysis reveals the genetic basis of skeletal muscle glycolytic potential based on a pig model

Glycolytic potential (GP) calculated based on glucose, glycogen, glucose-6-phosphate, and lactate contents is a critical factor for multiple meat quality characteristics. However, the genetic basis of glycolytic metabolism is still unclear. In this study, we constructed six RNA-Seq libraries using longissimus dorsi (LD) muscles from pigs divergent for GP phenotypic values and generated the whole genome-wide gene expression profiles. Furthermore, we identified 25,880 known and 220 novel genes from these skeletal muscle libraries, and 222 differentially expressed genes (DEGs) between the higher and lower GP groups. Notably, we found that the Lactate dehydrogenase B (LDHB) and Fructose-2, 6-biphosphatase 3 (PFKFB3) expression levels were higher in the higher GP group than the lower GP group, and positively correlated with GP and lactic acid (LA), and reversely correlated with pH value at 45 min postmortem (pH_45min). Besides, LDHB and PFKFB3 expression were positively correlated with drip loss measured at 48 h postmortem (DL_48h) and drip loss measured at 24 h postmortem (DL_24h). Collectively, we identified a serial of DEGs as the potential key candidate genes affecting GP and found that LDHB and PFKFB3 are closely related to GP and GP-related traits. Our results lay a solid basis for in-depth studies of the regulatory mechanisms on GP and GP-related traits in pigs.

Genomic analyses provide insights into breed-of-origin effects from purebreds on three-way crossbred pigs

Crossbreeding is widely used aimed at improving crossbred performance for poultry and livestock. Alleles that are specific to different purebreds will yield a large number of heterozygous single-nucleotide polymorphisms (SNPs) in crossbred individuals, which are supposed to have the power to alter gene function or regulate gene expression. For pork production, a classic three-way crossbreeding system of Duroc × (Landrace × Yorkshire) (DLY) is generally used to produce terminal crossbred pigs with stable and prominent performance. Nonetheless, little is known about the breed-of-origin effects from purebreds on DLY pigs. In this study, we first estimated the distribution of heterozygous SNPs in three kinds of three-way crossbred pigs via whole genome sequencing data originated from three purebreds. The result suggested that DLY is a more effective strategy for three-way crossbreeding as it could yield more stably inherited heterozygous SNPs. We then sequenced a DLY pig family and identified 95, 79, 132 and 42 allele-specific expression (ASE) genes in adipose, heart, liver and skeletal muscle, respectively. Principal component analysis and unrestricted clustering analyses revealed the tissue-specific pattern of ASE genes, indicating the potential roles of ASE genes for development of DLY pigs. In summary, our findings provided a lot of candidate SNP markers and ASE genes for DLY three-way crossbreeding system, which may be valuable for pig breeding and production in the future.

Sex and the TEs: transposable elements in sexual development and function in animals

Transposable elements are endogenous DNA sequences able to integrate into and multiply within genomes. They constitute a major source of genetic innovations, as they can not only rearrange genomes but also spread ready-to-use regulatory sequences able to modify host gene expression, and even can give birth to new host genes. As their evolutionary success depends on their vertical transmission, transposable elements are intrinsically linked to reproduction. In organisms with sexual reproduction, this implies that transposable elements have to manifest their transpositional activity in germ cells or their progenitors. The control of sexual development and function can be very versatile, and several studies have demonstrated the implication of transposable elements in the evolution of sex. In this review, we report the functional and evolutionary relationships between transposable elements and sexual reproduction in animals. In particular, we highlight how transposable elements can influence expression of sexual development genes, and how, reciprocally, they are tightly controlled in gonads. We also review how transposable elements contribute to the organization, expression and evolution of sexual development genes and sex chromosomes. This underscores the intricate co-evolution between host functions and transposable elements, which regularly shift from a parasitic to a domesticated status useful to the host.

Genome-Wide Analysis of Allele-Specific Expression Patterns in Seventeen Tissues of Korean Cattle (Hanwoo)

The functional hemizygosity could be caused by the MAE of a given gene and it can be one of the sources to affect the phenotypic variation in cattle. We aimed to identify MAE genes across the transcriptome in Korean cattle (Hanwoo). For three Hanwoo family trios, the transcriptome data of 17 tissues were generated in three offspring. Sixty-two MAE genes had a monoallelic expression in at least one tissue. Comparing genotypes among each family trio, the preferred alleles of 18 genes were identified (maternal expression, n = 9; paternal expression, n = 9). The MAE genes are involved in gene regulation, metabolic processes, and immune responses, and in particular, six genes encode transcription factors (FOXD2, FOXM1, HTATSF1, SCRT1, NKX6-2, and UBN1) with tissue-specific expression. In this study, we report genome-wide MAE genes in seventeen tissues of adult cattle. These results could help to elucidate epigenetic effects on phenotypic variation in Hanwoo.

Gene expression dynamics during the gonocyte to spermatogonia transition and spermatogenesis in the domestic yak

Background

Spermatogenesis is a cellular differentiation process that includes three major events: mitosis of spermatogonia, meiosis of spermatocytes and spermiogenesis. Steady-state spermatogenesis relies on functions of spermatogonial stem cells (SSCs). Establishing and maintaining a foundational SSC pool is essential for continued spermatogenesis in mammals. Currently, our knowledge about SSC and spermatogenesis is severely limited in domestic animals.

Results

In the present study, we examined transcriptomes of testes from domestic yaks at four different stages (3, 5, 8 and 24 months of age) and attempted to identify genes that are associated with key developmental events of spermatogenesis. Histological analyses showed that the most advanced germ cells within seminiferous tubules of testes from 3, 5, 8 and 24 months old yaks were gonocytes, spermatogonia, spermatocytes and elongated spermatids, respectively. RNA-sequencing (RNA-seq) analyses revealed that 11904, 4381 and 2459 genes were differentially expressed during the gonocyte to spermatogonia transition, the mitosis to meiosis transition and the meiosis to post-meiosis transition. Further analyses identified a list of candidate genes than may regulate these important cellular processes. CXCR4, a previously identified SSC niche factor in mouse, was one of the up-regulated genes in the 5 months old yak testis. Results of immunohistochemical staining confirmed that CXCR4 was exclusively expressed in gonocytes and a subpopulation of spermatogonia in the yak testis.

Conclusions

Together, these findings demonstrated histological changes of postnatal testis development in the domestic yak. During development of spermatogonial lineage, meiotic and haploid germ cells are supported by dynamic transcriptional regulation of gene expression. Our transcriptomic analyses provided a list of candidate genes that potentially play crucial roles in directing the establishment of SSC and spermatogenesis in yak.

Electronic supplementary material

The online version of this article (10.1186/s40104-019-0360-7) contains supplementary material, which is available to authorized users.

GROWTH AND DEVELOPMENT SYMPOSIUM: STEM AND PROGENITOR CELLS IN ANIMAL GROWTH: Long noncoding RNAs in adipogenesis and adipose development of meat animals12

Sequencing technology, especially next-generation RNA sequencing, has greatly facilitated the identification and annotation of long noncoding RNAs (lncRNAs). In mammals, a large number of lncRNAs have been identified, which regulate various biological processes. An increasing number of lncRNAs have been identified which could function as key regulators of adipogenesis (adipocyte formation), a key step of the development of adipose tissue. Because proper adipose tissue development is a key factor affecting animal growth efficiency, lean/fat ratio, and meat quality, summarizing the roles and recent advances of lncRNAs in adipogenesis is needed in order to develop strategies to effectively manage fat deposition. In this review, we updated lncRNAs contributed to the regulation of adipogenesis, focusing on their roles in fat development of farm animals.

Comparative transcriptomic analysis of hypothalamus-pituitary-liver axis in bighead carp (Hypophthalmichthys nobilis) with differential growth rate

Background

Growth rate is one of the most important features for aquaculture species and deciphering its regulation mechanism has great significance both in genetics and in economics. Hypothalamus-pituitary growth axis (HP growth axis) or neuro-endocrine axis plays a vital role in growth regulation in different aquaculture animals.

Results

In this study, the HP and liver transcriptomes of two female groups (H and L) with phenotypically extreme growth rate were sequenced using RNA-Seq. A total of 30,524 and 22,341 genes were found expressed in the two tissues, respectively. The average expression levels for the two tissues were almost the same, but the median differed significantly. A differential expression analysis between H and L groups identified 173 and 204 differentially expressed genes (DEGs) in HP and liver tissue, respectively. Pathway analysis revealed that DEGs in HP tissue were enriched in regulation of cell proliferation and angiogenesis while in liver tissue these genes were overrepresented in sterol biosynthesis and transportation. Genomic overlapping analyses found that 4 and 5 DEGs were within growth-related QTL in HP and liver tissue respectively. A deeper analysis of these 9 genes revealed 3 genes were functionally linked to the trait of interest. The expression of 2075 lncRNAs in HP tissue and 1490 in liver tissue were also detected, and some of lncRNAs were highly expressed in the two tissues.

Conclusions

Above all, the results of the present study greatly contributed to the knowledge of the regulation of growth and then assisted the design of new selection strategies for bighead carp with improved growth-related traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5691-4) contains supplementary material, which is available to authorized users.

Integrative Analyses of Long Non-coding RNA and mRNA Involved in Piglet Ileum Immune Response to Clostridium perfringens Type C Infection

Long non-coding RNAs (lncRNAs) have been shown to play important roles in regulating host immune and inflammatory responses to bacterial infection. Infection with Clostridium perfringens (C. perfringens), a food-borne zoonotic pathogen, can lead to a series of inflammatory diseases in human and piglet, greatly challenging the healthy development of global pig industry. However, the roles of lncRNAs involved in piglet immune response against C. perfringens type C infection remain unknown. In this study, the regulatory functions of ileum lncRNAs and mRNAs were investigated in piglet immune response to C. perfringens type C infection among resistance (IR), susceptibility (IS) and sham-inoculation (control, IC) groups. A total of 480 lncRNAs and 3,669 mRNAs were significantly differentially expressed, the differentially expressed lncRNAs and mRNAs in the IR and IS groups were enriched in various pathways of ABC transporters, olfactory transduction, PPAR signaling pathway, chemokine signaling pathway and Toll-like receptor signaling pathway, involving in regulating piglet immune responses and resistance during infection. There were 212 lncRNAs and 505 target mRNAs found to have important association with C. perfringens infectious diseases, furthermore, 25 dysregulated lncRNAs corresponding to 13 immune-related target mRNAs were identified to play potential roles in defense against bacterial infection. In conclusion, the results improve our understanding on the characteristics of lncRNAs and mRNAs on regulating host immune response against C. perfringens type C infection, which will provide a reference for future research into exploring C. perfringens-related diseases in human.

A RNA-Seq Analysis to Describe the Boar Sperm Transcriptome and Its Seasonal Changes

Understanding the molecular basis of cell function and ultimate phenotypes is crucial for the development of biological markers. With this aim, several RNA-seq studies have been devoted to the characterization of the transcriptome of ejaculated spermatozoa in relation to sperm quality and fertility. Semen quality follows a seasonal pattern and decays in the summer months in several animal species. The aim of this study was to deeply profile the transcriptome of the boar sperm and to evaluate its seasonal changes. We sequenced the total and the short fractions of the sperm RNA from 10 Pietrain boars, 5 collected in summer and 5 five sampled in winter, and identified a complex and rich transcriptome with 4,436 coding genes of moderate to high abundance. Transcript fragmentation was high but less obvious in genes related to spermatogenesis, chromatin compaction and fertility. Short non-coding RNAs mostly included piwi-interacting RNAs, transfer RNAs and microRNAs. We also compared the transcriptome of the summer and the winter ejaculates and identified 34 coding genes and 7 microRNAs with a significantly distinct distribution. These genes were mostly related to oxidative stress, DNA damage and autophagy. This is the deepest characterization of the boar sperm transcriptome and the first study linking the transcriptome and the seasonal variability of semen quality in animals. The annotation described here can be used as a reference for the identification of markers of sperm quality in pigs.

Transcriptome analysis to identify long non coding RNA (lncRNA) and characterize their functional role in back fat tissue of pig

Long non coding RNAs (lncRNA) have been previously found to be involved in important cellular activities like epigenetics, implantation, cell growth etc. in pigs. However, comprehensive analysis of lncRNA in back fat tissues at different developmental stages in pigs is still lacking. In this study we conducted transcriptome analysis in the back fat tissue of a F1 crossbred Korean Native Pig (KNP) × Yorkshire Pig to identify lncRNA. We investigated their role in 16 pigs at two different growth stages; stage 1 (10 weeks, n = 8) and stage 2 (26 weeks, n = 8). After quality assessment of sequencing reads, we got a total of 1,641,165 assembled transcripts out of eight paired end read from each stage. Among them, 6808 lncRNA transcripts were identified by filtering on the basis of multiple parameters like read length ≥ 200 nucleotides, exon numbers ≥2, FPKM ≥0.5, coding potential score < 0 etc. PFAM and RFAM were used to filter out all possible protein coding genes and housekeeping RNAs respectively. A total of 103 lncRNAs and 1057 mRNAs were found to be differentially expressed (DE) between the two stages (|log2FC| > 2, q < 0.05). We also identified 306 genes located around 100 kb upstream and 234 genes downstream around these DE lncRNA transcripts. The expression of top eleven DE lncRNAs (COL4A6, LY7S, MYH2, OXCT1, SMPDL3A, TMEM182, TTC36, RFOOOO4, RFOOO15, RFOOO45, CADM2) had been validating by qRT-PCR. Pathway and GO terms analysis showed that, positive regulation of biosynthetic process, Wnt signaling pathway, cellular protein modification process, and positive regulation of nitrogen compound were differentially enriched. Our results suggested that, KEGG pathways such as protein digestion and absorption, Arrhythmogenic right ventricular cardiomyopathy (ARVC) to be significantly enriched in both DE lncRNAs as well as DE mRNAs and involved in back fat tissues development. It also suggests that, identified lncRNAs are involved in regulation of important adipose tissues development pathways.

Survey of allele specific expression in bovine muscle

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3’UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

Comparative study on seasonal hair follicle cycling by analysis of the transcriptomes from cashmere and milk goats

Guard hair and cashmere undercoat are developed from primary and secondary hair follicle, respectively. Little is known about the gene expression differences between primary and secondary hair follicle cycling. In this study, we obtained RNA-seq data from cashmere and milk goats grown at four different seasons. We studied the differentially expressed genes (DEGs) during the yearly hair follicle cycling, and between cashmere and milk goats. WNT, NOTCH, MAPK, BMP, TGFβ and Hedgehog signaling pathways were involved in hair follicle cycling in both cashmere and milk goat. However, Milk goat DEGs between different months were significantly more than cashmere goat DEGs, with the largest difference being identified in December. Some expression dynamics were confirmed by quantitative PCR and western blot, and immunohistochemistry. This study offers new information sources related to hair follicle cycling in milk and cashmere goats, which could be applicable to improve the wool production and quality.

Investigation of allele-specific expression of genes involved in adipogenesis and lipid metabolism suggests complex regulatory mechanisms of PPARGC1A expression in porcine fat tissues

Background

The expression of genes involved in regulating adipogenesis and lipid metabolism may affect economically important fatness traits in pigs. Allele-specific expression (ASE) reflects imbalance between allelic transcript levels and can be used to identify underlying cis-regulatory elements. ASE has not yet been intensively studied in pigs. The aim of this investigation was to analyze the differential allelic expression of four genes, PPARA, PPARG, SREBF1, and PPARGC1A, which are involved in the regulation of fat deposition in porcine subcutaneous and visceral fat and longissimus dorsi muscle.

Results

Quantification of allelic proportions by pyrosequencing revealed that both alleles of PPARG and SREBF1 are expressed at similar levels. PPARGC1A showed the greatest ASE imbalance in fat deposits in Polish Large White (PLW), Polish Landrace and Pietrain pigs; and PPARA in PLW pigs. Significant deviations of mean PPARGC1A allelic transcript ratio between cDNA and genomic DNA were detected in all tissues, with the most pronounced difference (p < 0.001) in visceral fat of PLW pigs. To search for potential cis-regulatory elements affecting ASE in the PPARGC1A gene we analyzed the effects of four SNPs (rs337351686, rs340650517, rs336405906 and rs345224049) in the promoter region, but none were associated with ASE in the breeds studied. DNA methylation analysis revealed significant CpG methylation differences between samples showing balanced (allelic transcript ratio ≈1) and imbalanced allelic expression for CpG site at the genomic position in chromosome 8 (SSC8): 18527678 in visceral fat (p = 0.017) and two CpG sites (SSC8:18525215, p = 0.030; SSC8:18525237, p = 0.031) in subcutaneous fat.

Conclusions

Our analysis of differential allelic expression suggests that PPARGC1A is subjected to cis-regulation in porcine fat tissues. Further studies are necessary to identify other regulatory elements localized outside the PPARGC1A proximal promoter region.

Electronic supplementary material

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

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species

Background

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.

Results

Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.

Conclusions

LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5037-7) contains supplementary material, which is available to authorized users.

Genome Wide Identification of Novel Long Non-coding RNAs and Their Potential Associations With Milk Proteins in Chinese Holstein Cows

Long non-coding RNAs (lncRNAs) have emerged as a novel class of regulatory molecules involved in various biological processes. However, their role in milk performance is unknown. Here, whole transcriptome RNA sequencing was used to generate the lncRNA transcriptome profiles in mammary tissue samples from 6 Chinese Holstein cows with 3 extremely high and 3 low milk protein percentage phenotypes. In this study, 6,450 lncRNA transcripts were identified through 5 stringent steps and filtration by coding potential. In total, 31 lncRNAs and 18 novel genes were identified to be differentially expressed in high milk protein samples (HP) relative to low milk protein samples (LP), respectively. Differentially expressed lncRNAs were selected to predict target genes through bioinformatics analysis, followed by the integration of differentially expressed mRNA data, gene function, gene ontology (GO) and pathway, genome wide association study (GWAS) and quantitative trait locus (QTL) information, as well as network analysis to further characterize potential interactions. Several lncRNAs were found (such as XLOC_059976) that could be used as candidate markers for milk protein content prediction. This is the first study to perform global expression profiling of lncRNAs and mRNAs related to milk protein traits in dairy cows. These results provide important information and insights into the synthesis of milk proteins, and potential targets for the future improvement of milk quality.

Transcriptome signature of liver tissue with divergent mutton odour and flavour using RNA deep sequencing

Mutton consumption is less popular in many Asian countries including Indonesia, whose consumers often complain about the unpleasant flavour and odour of the meat. The main causes of mutton odour are the two compounds of branched chain fatty acid (BCFA): methylnonanoic (MNA), phenol, 3-methyl (MP), 4-methylnonanoic (MNA) and 4-ethyloctanoic (EOA) present in all the adipose tissue; and the 3-methylindole (MI) or skatole and indole, which are originated from pastoral diets. It is crucial to understand the genetic mechanism of mutton odour and flavour (MOF) to select sheep for lower BCFA and indole thus reduce the unpleasant flavour of meat. The aim of the present study was to investigate transcriptome profiling in liver tissue with divergent MOF using RNA deep sequencing. Liver tissues from higher (n = 3) and lower (n = 3) MOF sheep were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample ranged from 21.37 to 25.37 million. Approximately 103 genes were differentially expressed (DEGs) with significance level of p-adjusted value <0.05. Among them, 60 genes were up-regulated, and 43 were down-regulated (p < 0.01, FC > 1.5) in higher MOF group. Differentially regulated genes in high MOF liver samples were enriched in biological processes such as cellular response to chemical stimulus and endogenous stimulus; cellular components such as such as basement membrane and extracellular matrix; and molecular functions such as haeme binding and oxidoreductase activity. Among the DEGs, metabolic phase I related genes belonging to the cytochrome P450 CYP2A6 were dominantly expressed. Additionally, phase II conjugation genes including UDP glucuronosyltransferases UGT2B18, sulfotransferase SULT1C1, and glutathione S-transferase GSTM1 were identified. The dominant candidate genes for SOF could be cytochrome P450, sodium-channel protein, transmembrane protein, glutathione transferase, UDP glucuronosyltransferases and sulfotransferase. Pathway analysis identified steroid hormone biosynthesis and chemical carcinogenesis by cytochrome P450 pathways which may play important roles in MOF-related molecules metabolism. This work highlighted potential genes and gene-networks that may affect meat off flavour and odour in sheep.

Long non-coding RNAs involved in the regulatory network during porcine pre-implantation embryonic development and iPSC induction

Long non-coding RNAs (lncRNA) play a key role in the orchestration of transcriptional regulation during development and many other cellular processes. The importance of the regulatory co-expression network was highlighted in the identification of the mechanism of these processes in humans and mice. However, elucidation of the properties of porcine lncRNAs involved in the regulatory network during pre-implantation embryonic development and fibroblast reprogramming to induced pluripotent stem cell (iPSC) has been limited to date. Using a weighted gene co-expression network analysis, we constructed the regulatory network and determined that the novel lncRNAs were functionally involved in key events of embryonic development during the pre-implantation period; moreover, reprogramming could be delineated by a small number of potentially functional modules of co-expressed genes. These findings indicate that lncRNAs may be involved in the transcriptional regulation of zygotic genome activation, first lineage segregation and somatic reprogramming to pluripotency. Furthermore, we performed a conservation and synteny analysis with the significant lncRNAs involved in these vital events and validated the results via experimental assays. In summary, the current findings provide a valuable resource to dissect the protein coding gene and lncRNA regulatory networks that underlie the progressive development of embryos and somatic reprogramming.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

A comprehensive transcriptome analysis of skeletal muscles in two Polish pig breeds differing in fat and meat quality traits

Pork is the most popular meat in the world. Unfortunately, the selection pressure focused on high meat content led to a reduction in pork quality. The present study used RNA-seq technology to identify metabolic process genes related to pork quality traits and fat deposition. Differentially expressed genes (DEGs) were identified between pigs of Pulawska and Polish Landrace breeds for two the most important muscles (semimembranosus and longissimus dorsi). A total of 71 significant DEGs were reported: 15 for longissimus dorsi and 56 for semimembranosus muscles. The genes overexpressed in Pulawska pigs were involved in lipid metabolism (APOD, LXRA, LIPE, AP2B1, ENSSSCG00000028753 and OAS2) and proteolysis (CST6, CTSD, ISG15 and UCHL1). In Polish Landrace pigs, genes playing a role in biological adhesion (KIT, VCAN, HES1, SFRP2, CDH11, SSX2IP and PCDH17), actin cytoskeletal organisation (FRMD6, LIMK1, KIF23 and CNN1) and calcium ion binding (PVALB, CIB2, PCDH17, VCAN and CDH11) were transcriptionally more active. The present study allows for better understanding of the physiological processes associated with lipid metabolism and muscle fiber organization. This information could be helpful in further research aiming to estimate the genetic markers.

Whole-genome analysis of structural variations between Xiang pigs with larger litter sizes and those with smaller litter sizes

To gain a better knowledge of structural variations (SVs) in Xiang pig, we used next-generation sequencing to analyze the Xiang pigs with larger (XL) or smaller litter sizes (XS). Our analysis yielded 28,040 putative SVs in the Xiang pig. These SVs distributed throughout all of chromosomes. Some functional regions including exons and untranslated regions were less varied than introns and intergenic regions. We detected 4637 and 4119 specific SVs, which contained 1697 and 1582 genes in XL and XS group, respectively. These genes were mainly enriched in the well-known pathways involved in development and reproduction processes. Population validation was carried out on 50 SVs candidates using PCR method in 144 Xiang pig crowds. All of 50 SVs were confirmed by PCR method and 14 SVs were associated with the litter size of Xiang pigs. These results may be helpful for the elucidation of growth and reproduction regulation in Xiang pig.

Systems genomics study reveals expression quantitative trait loci, regulator genes and pathways associated with boar taint in pigs

Boar taint is an offensive odour and/or taste from a proportion of non-castrated male pigs caused by skatole and androstenone accumulation during sexual maturity. Castration is widely used to avoid boar taint but is currently under debate because of animal welfare concerns. This study aimed to identify expression quantitative trait loci (eQTLs) with potential effects on boar taint compounds to improve breeding possibilities for reduced boar taint. Danish Landrace male boars with low, medium and high genetic merit for skatole and human nose score (HNS) were slaughtered at ~100 kg. Gene expression profiles were obtained by RNA-Seq, and genotype data were obtained by an Illumina 60K Porcine SNP chip. Following quality control and filtering, 10,545 and 12,731 genes from liver and testis were included in the eQTL analysis, together with 20,827 SNP variants. A total of 205 and 109 single-tissue eQTLs associated with 102 and 58 unique genes were identified in liver and testis, respectively. By employing a multivariate Bayesian hierarchical model, 26 eQTLs were identified as significant multi-tissue eQTLs. The highest densities of eQTLs were found on pig chromosomes SSC12, SSC1, SSC13, SSC9 and SSC14. Functional characterisation of eQTLs revealed functions within regulation of androgen and the intracellular steroid hormone receptor signalling pathway and of xenobiotic metabolism by cytochrome P450 system and cellular response to oestradiol. A QTL enrichment test revealed 89 QTL traits curated by the Animal Genome PigQTL database to be significantly overlapped by the genomic coordinates of cis-acting eQTLs. Finally, a subset of 35 cis-acting eQTLs overlapped with known boar taint QTL traits. These eQTLs could be useful in the development of a DNA test for boar taint but careful monitoring of other overlapping QTL traits should be performed to avoid any negative consequences of selection.

Gene expression profile changes in the jejunum of weaned piglets after oral administration of Lactobacillus or an antibiotic

The small intestine plays an essential role in the health and well-being of animals. Previous studies have shown that Lactobacillus has a protective effect on intestinal morphology, intestinal epithelium integrity and appropriate maturation of gut-associated tissues. Here, gene expression in jejunum tissue of weaned piglets was investigated by RNA-seq analysis after administration of sterile saline, Lactobacillus reuteri, or an antibiotic (chlortetracycline). In total, 401 and 293 genes were significantly regulated by chlortetracycline and L. reuteri, respectively, compared with control treatment. Notably, the HP, NOX1 and GPX2 genes were significantly up-regulated in the L. reuteri group compared with control, which is related to the antioxidant ability of this strain. In addition, the expression of genes related to arachidonic acid metabolism and linoleic acid metabolism enriched after treatment with L. reuteri. The fatty acid composition in the jejunum and colon was examined by GC-MS analysis and suggested that the MUFA C18:1n9c, and PUFAs C18:2n6c and C20:4n6 were increased in the L. reuteri group, verifying the GO enrichment and KEGG pathway analyses of the RNA-seq results. The results contribute to our understanding of the probiotic activity of this strain and its application in pig production.

RNA-Seq Analyses Identify Frequent Allele Specific Expression and No Evidence of Genomic Imprinting in Specific Embryonic Tissues of Chicken

Epigenetic and genetic cis-regulatory elements in diploid organisms may cause allele specific expression (ASE) – unequal expression of the two chromosomal gene copies. Genomic imprinting is an intriguing type of ASE in which some genes are expressed monoallelically from either the paternal allele or maternal allele as a result of epigenetic modifications. Imprinted genes have been identified in several animal species and are frequently associated with embryonic development and growth. Whether genomic imprinting exists in chickens remains debatable, as previous studies have reported conflicting evidence. Albeit no genomic imprinting has been reported in the chicken embryo as a whole, we interrogated the existence or absence of genomic imprinting in the 12-day-old chicken embryonic brain and liver by examining ASE in F1 reciprocal crosses of two highly inbred chicken lines (Fayoumi and Leghorn). We identified 5197 and 4638 ASE SNPs, corresponding to 18.3% and 17.3% of the genes with a detectable expression in the embryonic brain and liver, respectively. There was no evidence detected of genomic imprinting in 12-day-old embryonic brain and liver. While ruling out the possibility of imprinted Z-chromosome inactivation, our results indicated that Z-linked gene expression is partially compensated between sexes in chickens.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

Genome-wide analysis of DNA methylation in obese, lean, and miniature pig breeds

DNA methylation is a crucial epigenetic modification involved in diverse biological processes. There is significant phenotypic variance between Chinese indigenous and western pig breeds. Here, we surveyed the genome-wide DNA methylation profiles of blood leukocytes from three pig breeds (Tongcheng, Landrace, and Wuzhishan) by methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in gene body regions and repetitive sequences. LINE/L1 and SINE/tRNA-Glu were the predominant methylated repeats in pigs. The methylation level in the gene body regions was higher than in the 5' and 3' flanking regions of genes. About 15% of CpG islands were methylated in the pig genomes. Additionally, 2,807, 2,969, and 5,547 differentially methylated genes (DMGs) were identified in the Tongcheng vs. Landrace, Tongcheng vs. Wuzhishan, and Landrace vs. Wuzhishan comparisons, respectively. A total of 868 DMGs were shared by the three contrasts. The DMGs were significantly enriched in development- and metabolism-related biological processes and pathways. Finally, we identified 32 candidate DMGs associated with phenotype variance in pigs. Our research provides a DNA methylome resource for pigs and furthers understanding of epigenetically regulated phenotype variance in mammals.

Histological and transcriptome analyses of testes from Duroc and Meishan boars

Meishan boars are known for their early sexual maturity. However, they exhibit a significantly smaller testicular size and a reduced proportion of Sertoli cells and daily sperm production compared with Duroc boars. The testes of Duroc and Meishan boars at 20, 75 and 270 days of age were used for histological and transcriptome analyses. Haematoxylin-eosin staining was conducted to observe histological structure of the testes in Duroc and Meishan boars at different ages. Although spermatogenesis occurred prior to 75 days in Meishan boars, the number of spermatogonia and Sertoli cells in Meishan boars were less than in Duroc boars at adulthood. The diameters of the seminiferous tubules of the testes differed significantly during the initiation of development of the seminiferous tubules between the two breeds. We obtained differentially expressed functional genes and analysed seven pathways involved in male sexual maturity and spermatogenesis using RNA-seq. We also detected four main alternative splicing events and many single nucleotide polymorphisms from testes. Eight functionally important genes were validated by qPCR, and Neurotrophin 3 was subjected to quantification and cellular localization analysis. Our study provides the first transcriptome evidence for the differences in sexual function development between Meishan and Duroc boars.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Improvement of genome assembly completeness and identification of novel full-length protein-coding genes by RNA-seq in the giant panda genome

High-quality and complete gene models are the basis of whole genome analyses. The giant panda (Ailuropoda melanoleuca) genome was the first genome sequenced on the basis of solely short reads, but the genome annotation had lacked the support of transcriptomic evidence. In this study, we applied RNA-seq to globally improve the genome assembly completeness and to detect novel expressed transcripts in 12 tissues from giant pandas, by using a transcriptome reconstruction strategy that combined reference-based and de novo methods. Several aspects of genome assembly completeness in the transcribed regions were effectively improved by the de novo assembled transcripts, including genome scaffolding, the detection of small-size assembly errors, the extension of scaffold/contig boundaries, and gap closure. Through expression and homology validation, we detected three groups of novel full-length protein-coding genes. A total of 12.62% of the novel protein-coding genes were validated by proteomic data. GO annotation analysis showed that some of the novel protein-coding genes were involved in pigmentation, anatomical structure formation and reproduction, which might be related to the development and evolution of the black-white pelage, pseudo-thumb and delayed embryonic implantation of giant pandas. The updated genome annotation will help further giant panda studies from both structural and functional perspectives.

Identification of Genes Related to Growth and Lipid Deposition from Transcriptome Profiles of Pig Muscle Tissue

Transcriptome profiles established using high-throughput sequencing can be effectively used for screening genome-wide differentially expressed genes (DEGs). RNA sequences (from RNA-seq) and microRNA sequences (from miRNA-seq) from the tissues of longissimus dorsi muscle of two indigenous Chinese pig breeds (Diannan Small-ear pig [DSP] and Tibetan pig [TP]) and two introduced pig breeds (Landrace [LL] and Yorkshire [YY]) were examined using HiSeq 2000 to identify and compare the differential expression of functional genes related to muscle growth and lipid deposition. We obtained 27.18 G clean data through the RNA-seq and detected that 18,208 genes were positively expressed and 14,633 of them were co-expressed in the muscle tissues of the four samples. In all, 315 DEGs were found between the Chinese pig group and the introduced pig group, 240 of which were enriched with functional annotations from the David database and significantly enriched in 27 Gene Ontology (GO) terms that were mainly associated with muscle fiber contraction, cadmium ion binding, response to organic substance and contractile fiber part. Based on functional annotation, we identified 85 DEGs related to growth traits that were mainly involved in muscle tissue development, muscle system process, regulation of cell development, and growth factor binding, and 27 DEGs related to lipid deposition that were mainly involved in lipid metabolic process and fatty acid biosynthetic process. With miRNA-seq, we obtained 23.78 M reads and 320 positively expressed miRNAs from muscle tissues, including 271 known pig miRNAs and 49 novel miRNAs. In those 271 known miRNAs, 20 were higher and 10 lower expressed in DSP-TP than in LL-YY. The target genes of the 30 miRNAs were mainly participated in MAPK, GnRH, insulin and Calcium signaling pathway and others involved cell development, growth and proliferation, etc. Combining the DEGs and the differentially expressed (DE) miRNAs, we drafted a network of 46 genes and 18 miRNAs for regulating muscle growth and a network of 15 genes and 16 miRNAs for regulating lipid deposition. We identified that CAV2, MYOZ2, FRZB, miR-29b, miR-122, miR-145-5p and miR-let-7c, etc, were key genes or miRNAs regulating muscle growth, and FASN, SCD, ADORA1, miR-4332, miR-182, miR-92b-3p, miR-let-7a and miR-let-7e, etc, were key genes or miRNAs regulating lipid deposition. The quantitative expressions of eight DEGs and seven DE miRNAs measured with real-time PCR certified that the results of differential expression genes or miRNAs were reliable. Thus, 18,208 genes and 320 miRNAs were positively expressed in porcine longissimus dorsi muscle. We obtained 85 genes and 18 miRNAs related to muscle growth and 27 genes and 16 miRNAs related to lipid deposition, which provided new insights into molecular mechanism of the economical traits in pig.

Optimization of next-generation sequencing transcriptome annotation for species lacking sequenced genomes

Next‐generation sequencing methods, such as RNA‐seq, have permitted the exploration of gene expression in a range of organisms which have been studied in ecological contexts but lack a sequenced genome. However, the efficacy and accuracy of RNA‐seq annotation methods using reference genomes from related species have yet to be robustly characterized. Here we conduct a comprehensive power analysis employing RNA‐seq data from Drosophila melanogaster in conjunction with 11 additional genomes from related Drosophila species to compare annotation methods and quantify the impact of evolutionary divergence between transcriptome and the reference genome. Our analyses demonstrate that, regardless of the level of sequence divergence, direct genome mapping (DGM), where transcript short reads are aligned directly to the reference genome, significantly outperforms the widely used de novo and guided assembly‐based methods in both the quantity and accuracy of gene detection. Our analysis also reveals that DGM recovers a more representative profile of Gene Ontology functional categories, which are often used to interpret emergent patterns in genomewide expression analyses. Lastly, analysis of available primate RNA‐seq data demonstrates the applicability of our observations across diverse taxa. Our quantification of annotation accuracy and reduced gene detection associated with sequence divergence thus provides empirically derived guidelines for the design of future gene expression studies in species without sequenced genomes.

MicroRNA-23b Promotes Avian Leukosis Virus Subgroup J (ALV-J) Replication by Targeting IRF1

Avian leukosis virus subgroup J (ALV-J) can cause several different leukemia-like proliferative diseases in the hemopoietic system of chickens. Here, we investigated the transcriptome profiles and miRNA expression profiles of ALV-J-infected and uninfected chicken spleens to identify the genes and miRNAs related to ALV-J invasion. In total, 252 genes and 167 miRNAs were differentially expressed in ALV-J-infected spleens compared to control uninfected spleens. miR-23b expression was up-regulated in ALV-J-infected spleens compared with the control spleens, and transcriptome analysis revealed that the expression of interferon regulatory factor 1 (IRF1) was down-regulated in ALV-J-infected spleens compared to uninfected spleens. A dual-luciferase reporter assay showed that IRF1 was a direct target of miR-23b. miR-23b overexpression significantly (P = 0.0022) decreased IRF1 mRNA levels and repressed IRF1-3′-UTR reporter activity. In vitro experiments revealed that miR-23b overexpression strengthened ALV-J replication, whereas miR-23b loss of function inhibited ALV-J replication. IRF1 overexpression inhibited ALV-J replication, and IRF1 knockdown enhanced ALV-J replication. Moreover, IRF1 overexpression significantly (P = 0.0014) increased IFN-β expression. In conclusion, these results suggested that miR-23b may play an important role in ALV-J replication by targeting IRF1.

Identification of Reproduction-Related Gene Polymorphisms Using Whole Transcriptome Sequencing in the Large White Pig Population

Recent developments in high-throughput sequencing techniques have enabled large-scale analysis of genetic variations and gene expression in different tissues and species, but gene expression patterns and genetic variations in livestock are not well-characterized. In this study, we have used high-throughput transcriptomic sequencing of the Finnish Large White to identify gene expression patterns and coding polymorphisms within the breed in the testis and oviduct. The main objective of this study was to identify polymorphisms within genes that are highly and specifically expressed in male and/or female reproductive organs. The differential expression (DE) analysis underlined 1234 genes highly expressed in the testis and 1501 in the oviduct. Furthermore, we used a novel in-house R-package hoardeR for the identification of novel genes and their orthologs, which underlined 55 additional DE genes based on orthologs in the human, cow, and sheep. Identification of polymorphisms in the dataset resulted in a total of 29,973 variants, of which 10,704 were known coding variants. Fifty-seven nonsynonymous SNPs were present among genes with high expression in the testis and 67 were present in the oviduct, underlining possible influential genes for reproduction traits. Seven genes (PGR, FRAS1, TCF4, ADAT1, SPAG6, PIWIL2, and DNAH8) with polymorphisms were highlighted as reproduction-related based on their biological function. The expression and SNPs of these genes were confirmed using RT-PCR and Sanger sequencing. The identified nonsynonymous mutations within genes highly expressed in the testis or oviduct provide a list of candidate genes for reproduction traits within the pig population and enable identification of biomarkers for sow and boar fertility.

ALDB: a domestic-animal long noncoding RNA database

Background

Long noncoding RNAs (lncRNAs) have attracted significant attention in recent years due to their important roles in many biological processes. Domestic animals constitute a unique resource for understanding the genetic basis of phenotypic variation and are ideal models relevant to diverse areas of biomedical research. With improving sequencing technologies, numerous domestic-animal lncRNAs are now available. Thus, there is an immediate need for a database resource that can assist researchers to store, organize, analyze and visualize domestic-animal lncRNAs.

Results

The domestic-animal lncRNA database, named ALDB, is the first comprehensive database with a focus on the domestic-animal lncRNAs. It currently archives 12,103 pig intergenic lncRNAs (lincRNAs), 8,923 chicken lincRNAs and 8,250 cow lincRNAs. In addition to the annotations of lincRNAs, it offers related data that is not available yet in existing lncRNA databases (lncRNAdb and NONCODE), such as genome-wide expression profiles and animal quantitative trait loci (QTLs) of domestic animals. Moreover, a collection of interfaces and applications, such as the Basic Local Alignment Search Tool (BLAST), the Generic Genome Browser (GBrowse) and flexible search functionalities, are available to help users effectively explore, analyze and download data related to domestic-animal lncRNAs.

Conclusions

ALDB enables the exploration and comparative analysis of lncRNAs in domestic animals. A user-friendly web interface, integrated information and tools make it valuable to researchers in their studies. ALDB is freely available from http://res.xaut.edu.cn/aldb/index.jsp.

Differential gene expression in ovaries of Qira black sheep and Hetian sheep using RNA-Seq technique

The Qira black sheep and the Hetian sheep are two local breeds in the Northwest of China, which are characterized by high-fecundity and low-fecundity breed respectively. The elucidation of mRNA expression profiles in the ovaries among different sheep breeds representing fecundity extremes will helpful for identification and utilization of major prolificacy genes in sheep. In the present study, we performed RNA-seq technology to compare the difference in ovarian mRNA expression profiles between Qira black sheep and Hetian sheep. From the Qira black sheep and the Hetian sheep libraries, we obtained a total of 11,747,582 and 11,879,968 sequencing reads, respectively. After aligning to the reference sequences, the two libraries included 16,763 and 16,814 genes respectively. A total of 1,252 genes were significantly differentially expressed at Hetian sheep compared with Qira black sheep. Eight differentially expressed genes were randomly selected for validation by real-time RT-PCR. This study provides a basic data for future research of the sheep reproduction.

Systematic identification and characterization of long intergenic non-coding RNAs in fetal porcine skeletal muscle development

Long intergenic non-coding RNAs (lincRNAs) play important roles in many cellular processes. Here, we present the first systematic identification and characterization of lincRNAs in fetal porcine skeletal muscle. We obtained a total of 55.02 million 90-bp paired-end reads and assembled 54,550 transcripts using cufflinks. We developed a pipeline to identify 570 multi-exon lincRNAs by integrating a set of previous approaches. These putative porcine lincRNAs share many characteristics with mammalian lincRNAs, such as a relatively short length, small number of exons and low level of sequence conservation. We found that the porcine lincRNAs were preferentially located near genes mediating transcriptional regulation rather than those with developmental functions. We further experimentally analyzed the features of a conserved mouse lincRNA gene and found that isoforms 1 and 4 of this lincRNA were enriched in the cell nucleus and were associated with polycomb repressive complex 2 (PRC2). Our results provide a catalog of fetal porcine lincRNAs for further experimental investigation of the functions of these genes in the skeletal muscle developmental process.

Transcriptomic analysis of Rongchang pig brains and livers

Recent developments in high-throughput RNA sequencing (RNA-seq) technology have led to a dramatic impact on our understanding of the structure and expression profiles of the mammalian transcriptome. To gain insights into the usefulness of swine production and biomedical model, the transcriptome profiling of Rongchang pig brains and livers was characterized using RNA-seq technology to uncover functional candidate molecules. In the study, total RNAs from brains and livers of Rongchang pig were sequenced and 8.6Gb sequencing data was obtained. This analysis revealed tissue specificity through the identification of 5575 and 4600 differentially expressed genes (DEGs) in brains and livers, respectively and the functional analysis of DEGs. Furthermore, 83 neuropeptide gene transcripts, 69 neuropeptide receptor gene transcripts, 10 pro-neuropeptide convertase gene transcripts and many other neuropeptide related protein gene transcripts were identified. Totally, the major characteristics of the transcriptional profiles of Rongchang pig brains and livers were present.

Integrated analysis of miRNA and mRNA expression profiles in development of porcine testes

Gene expression profile in the development of porcine testes investigates the intricate physiological process in pig testis development and spermatogenesis.

Current transcriptomics in pig immunity research

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.