Analyses of Long Non-Coding RNA and mRNA profiling using RNA sequencing in chicken testis with extreme sperm motility

A first complete catalog of highly expressed genes in eight chicken tissues reveals uncharacterized gene families specific for the chicken testis

Based on next-generation sequencing, we established a repertoire of differentially overexpressed genes (DoEGs) in eight adult chicken tissues: the testis, brain, lung, liver, kidney, muscle, heart, and intestine. With 4,499 DoEGs, the testis had the highest number and proportion of DoEGs compared with the seven somatic tissues. The testis DoEG set included the highest proportion of long noncoding RNAs (lncRNAs; 1,851, representing 32% of the lncRNA genes in the whole genome) and the highest proportion of protein-coding genes (2,648, representing 14.7% of the protein-coding genes in the whole genome). The main significantly enriched Gene Ontology terms related to the protein-coding genes were "reproductive process," "tubulin binding," and "microtubule cytoskeleton." Using real-time quantitative reverse transcription-polymerase chain reaction, we confirmed the overexpression of genes that encode proteins already described in chicken sperm [such as calcium binding tyrosine phosphorylation regulated (CABYR), spermatogenesis associated 18 (SPATA18), and CDK5 regulatory subunit associated protein (CDK5RAP2)] but whose testis origin had not been previously confirmed. Moreover, we demonstrated the overexpression of vertebrate orthologs of testis genes not yet described in the adult chicken testis [such as NIMA related kinase 2 (NEK2), adenylate kinase 7 (AK7), and CCNE2]. Using clustering according to primary sequence homology, we found that 1,737 of the 2,648 (67%) testis protein-coding genes were unique genes. This proportion was significantly higher than the somatic tissues except muscle. We clustered the other 911 testis protein-coding genes into 495 families, from which 47 had all paralogs overexpressed in the testis. Among these 47 testis-specific families, eight contained uncharacterized duplicated paralogs without orthologs in other metazoans except birds: these families are thus specific for chickens/birds.NEW & NOTEWORTHY Comparative next-generation sequencing analysis of eight chicken tissues showed that the testis has highest proportion of long noncoding RNA and protein-coding genes of the whole genome. We identified new genes in the chicken testis, including orthologs of known mammalian testicular genes. We also identified 47 gene families in which all the members were overexpressed, if not exclusive, in the testis. Eight families, organized in duplication clusters, were unknown, without orthologs in metazoans except birds, and are thus specific for chickens/birds.

Comprehensive analysis of the differential expression of mRNAs, lncRNAs, and miRNAs in Zi goose testis with high and low sperm mobility

Sperm mobility (SM) is an objective index for measuring sperm motility; however, the mechanisms underlying its regulation in geese remain unclear. The present study sought to elucidate the genetic mechanism underlying SM traits in Zi geese (Anser cygnoides L.). To this end, three successive experiments were performed. In Experiment I, SM was determined in 40 ganders; the 3 ganders with the highest mobility and three with the lowest mobility were assigned to the high and low sperm mobility rank (SMR) groups, respectively. In Experiment II, the differences in fertility between the two SMR groups were assessed within two breeding flocks comprising the selected six ganders from Experiment I and 30 females (each flock had 3 ganders and 15 females). In Experiment III, the testes of the 6 ganders were harvested for histological observation and whole-transcriptome sequencing. Results revealed better fertility, well-developed seminiferous tubules, and abundant mature sperm in the high-SMR-flock compared to those of the low-SMR-flock (89 vs. 81%) (P < 0.05). Differential expression (DE) analysis identified 76 mRNAs, 344 lncRNAs, and 17 miRNAs between the SMR groups, with LOC106049708, XPNPEP3, GNB3, ADCY8, PRKAG3, oha-miR-182-5p, and ocu-miR-10b-5p identified as key mRNAs and miRNAs contributing to SM. Enrichment analysis implicated these DE RNAs in pathways related to ATP binding, cell metabolism, apelin signaling, Wnt signaling, and Adherens junctions. Additionally, competing endogenous RNA (ceRNA) networks comprising 9 DE mRNAs, 17 DE miRNAs, and 169 DE lncRNAs were constructed. Two ceRNA network pathways (LOC106049708-oha-miR-182-5p-MSTRG.2479.6 and PRKAG3-ocu-miR-10b-5p-MSTRG.9047.14) were identified as key regulators of SM in geese. These findings offer crucial insights into the identification of key genes and ceRNA pathways influencing sperm mobility in geese.

Transcriptomic analysis of testis and epididymis tissues from Banna mini-pig inbred line boars with single-molecule long-read sequencing†

In mammals, testis and epididymis are critical components of the male reproductive system for androgen production, spermatogenesis, sperm transportation, as well as sperm maturation. Here, we report single-molecule real-time sequencing data from the testis and epididymis of the Banna mini-pig inbred line (BMI), a promising laboratory animal for medical research. We obtained high-quality full-length transcriptomes and identified 9879 isoforms and 8761 isoforms in the BMI testis and epididymis, respectively. Most of the isoforms we identified have novel exon structures that will greatly improve the annotation of testis- and epididymis-expressed genes in pigs. We also found that 3055 genes (over 50%) were shared between BMI testis and epididymis, indicating widespread expression profiles of genes related to reproduction. We characterized extensive alternative splicing events in BMI testis and epididymis and showed that 96 testis-expressed genes and 79 epididymis-expressed genes have more than six isoforms, revealing the complexity of alternative splicing. We accurately defined the transcribed isoforms in BMI testis and epididymis by combining Pacific Biotechnology Isoform-sequencing (PacBio Iso-Seq) and Illumina RNA Sequencing (RNA-seq) techniques. The refined annotation of some key genes governing male reproduction will facilitate further understanding of the molecular mechanisms underlying BMI male sterility. In addition, the high-confident identification of 548 and 669 long noncoding RNAs (lncRNAs) in these two tissues has established a candidate gene set for future functional investigations. Overall, our study provides new insights into the role of the testis and epididymis during BMI reproduction, paving the path for further studies on BMI male infertility.

Long noncoding RNA profile of the intracranial artery in patients with moyamoya disease

OBJECTIVE

Moyamoya disease (MMD) is a rare cerebrovascular disease characterized by progressive stenosis of the internal carotid artery (ICA) and secondary formation of collateral vessels. Revascularization surgery is performed in patients with MMD to prevent stroke; however, the pathogenesis of MMD remains unknown. Recently, long noncoding RNAs (lncRNAs) have been found to play a key role in gene regulation and are implicated in various vascular diseases. However, the lncRNA expression profile in MMD lesions has not been investigated. In this study the authors aimed to determine the characteristics of lncRNA expression in MMD lesions.

METHODS

The authors collected microsamples of the middle cerebral artery (MCA) from patients with MMD (n = 21) and patients with control conditions (n = 11) who underwent neurosurgical treatment. Using microarray experiments, the authors compared the profiles of lncRNA expression in the MCAs of the MMD and control patient groups and identified differentially expressed lncRNAs (fold change > 2, q < 0.05). In addition, the neighboring coding genes, whose transcription can be regulated in cis by the identified differentially expressed lncRNAs, were investigated and Gene Ontology (GO) analysis was applied to predict associated biological functions.

RESULTS

The authors detected 308 differentially expressed lncRNAs (fold change > 2, q < 0.05), including 306 upregulated and 2 downregulated lncRNAs in the MCA from patients with MMD. Regarding the prediction of biological function, GO analyses with possible coding genes whose transcription was regulated in cis by the identified differentially expressed lncRNAs suggested involvement in the antibacterial humoral response, T-cell receptor signaling pathway, positive regulation of cytokine production, and branching involved in blood vessel morphogenesis.

CONCLUSIONS

The profile of lncRNA expression in MMD lesions was different from that in the normal cerebral artery, and differentially expressed lncRNAs were identified. This study provides new insights into the pathophysiology of MMD.

Specific expression of alternatively spliced genes in the turkey (Meleagris gallopavo) reproductive tract revealed their function in spermatogenesis and post-testicular sperm maturation

The tissue-specific profile of alternatively spliced genes (ASGs) and their involvement in reproduction processes characteristic of turkey testis, epididymis, and ductus deferens were investigated for the first time in birds. Deep sequencing of male turkey reproductive tissue RNA samples (n = 6) was performed using Illumina RNA-Seq with 2 independent methods, rMATs and SUPPA2, for differential alternative splicing (DAS) event prediction. The expression of selected ASGs was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. The testis was found to be the site of the highest number of posttranscriptional splicing events within the reproductive tract, and skipping exons were the most frequently occurring class of alternative splicing (AS) among the reproductive tract. Statistical analysis revealed 86, 229, and 6 DAS events in the testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparison, respectively. Alternative splicing was found to be a mechanism of gene expression regulation within the turkey reproduction tract. In testis, modification was observed for spermatogenesis specific genes; the changes in 5' UTR could act as regulator of MEIG1 expression (a player during spermatocytes meiosis), and modification of 3' UTR led to diversification of CREM mRNA (modulator of gene expression related to the structuring of mature spermatozoa). Sperm tail formation can be regulated by changes in the 5' UTR of testicular SLC9A3R1 and gene silencing by producing dysfunctional variants of ODF2 in the testis and ATP1B3 in the epididymis. Predicted differentially ASGs in the turkey reproductive tract seem to be involved in the regulation of spermatogenesis, including acrosome formation and sperm tail formation and binding of sperm to the zona pellucida. Several ASGs were classified as cilia by actin and microtubule cytoskeleton organization. Such genes may play a role in the organization of sperm flagellum and post-testicular motility development. To our knowledge, this is the first functional investigation of alternatively spliced genes associated with tissue-specific processes in the turkey reproductive tract.

Whole transcriptome sequencing reveals core genes related to spermatogenesis in Bactrian camels

Bactrian camels survive and reproduce better in extreme climatic conditions than other domestic animals can. However, the reproductive efficiency of camels under their natural pastoral conditions is low. Several factors affect mammalian reproductive performance, including testicular development, semen quality, libido, and mating ability. Testis is a main reproductive organ of the male and is responsible for producing spermatozoa and hormones. However, our understanding of the expression patterns of the genes in camel testis is minimal. Thus, we performed total RNA-sequencing to investigate the gene expression pattern. As a result, 1,538 differential expressed mRNAs (DEmRNAs), 702 differential expressed long non-coding RNAs (DElncRNAs), and 61 differential expressed microRNAs (DEmiRNAs) were identified between pubertal and adult Bactrian camel testes. Then the genomic features, length distribution, and other characteristics of the lncRNAs and mRNAs in the Bactrian camel testis were investigated. Target genes of the DEmiRNAs and DEmRNAs were further subjected to gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Genes, such as AMHR2, FGF1, ACTL7A, GATA4, WNT4, ID2, LAMA1, IGF1, INHBB, and TLR2, were mainly involved in the TGF-β, PI3K-AKT, Wnt, GnRH, and Hippo signaling pathways which relate to spermatogenesis. Some of the DEmiRNAs were predicted to be associated with numerous DElncRNAs and DEmRNAs through competing endogenous RNA (ceRNA) regulatory network. At last, the candidate genes were validated by RT-qPCR, dual fluorescent reporter gene, and a fluorescence in situ hybridization (FISH) assay. This research provides high-throughput RNA sequencing data of the testes of Bactrian camels across different developmental stages. It lays the foundation for further investigations on lncRNAs, miRNAs, and mRNAs that involved in Bactrian camel spermatogenesis.

Identification of potential candidate genes and regulatory pathways related to reproductive capacity in hypothalamus and pituitarium of male ducks (Anas platyrhynchos) by differential transcriptome analysis

The improvement of reproductive capacity of poultry is important for the poultry industry. The existing studies on reproductive capacity mainly focus on the testis tissue, but few reports on regulationary effect of brain neuroendocrime on reproductive capacity have been available. The hypothalamus-pituitarium-gonad (HPG) axis is an important pathway regulating spermatogenesis and sexual behavior. This study analyzed the gene expression in the hypothalamus and pituitary tissues of male ducks in high-semen-quality group (DH), low-semen-quality group (DL), and non-response group (DN) by RNA-sequencing. A total of 1980 differentially expressed genes (DEGs) were identified, and significantly less DEGs were found in pituitary gland than in hypothalamus. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DEGs were mainly enriched in nerve-related and synapse-related biological processes, mitochondrial inner membrane formation pathway, and ribosome structure pathway. Notably, the neuroactive ligand-receptor interaction pathway significantly enriched in all three comparisons (DH vs. DL, DH vs. DN, and DL vs. DN) was related to different reproductive performance such as semen quality and sexual response. Furthermore, six genes, including POMC, CPLX2, HAPLN2, EGR4, TOX3, and MSH4, were identified as candidate genes regulating reproductive capacity. Our findings provide new insights into the regulation mechanisms underlying the reproductive performance of male poultry, and offer a valuable reference for duck breeding programs aimed at promoting reproductive capacity.

Comparative Transcriptome Analysis Provided a New Insight into the Molecular Mechanisms of Epididymis Regulating Semen Volume in Drakes

Semen volume is an important factor in artificial insemination (AI) of ducks. In drakes, seminal plasma that is produced by the epididymis determines the semen volume. However, the mechanism of epididymis regulating semen volume of drakes remains unclear. Therefore, the aim of the present study was to preliminarily reveal the mechanism regulating the semen volume through comparing the epididymal histomorphology and mRNA expression profiles between drakes with high-volume semen (HVS) and low-volume semen (LVS). Phenotypically, drakes in the HVS group produced more sperm than drakes in the LVS group. In addition, compared with the HVS group, the ductal square of ductuli conjugentes (DC) and dutus epididymidis (DE) in epididymis was significantly smaller in the LVS group, and the lumenal diameter and epithelial thickness of DC/DE were significantly shorter in the LVS group. In transcriptional regulation, 72 different expression genes (DEGs) were identified from the epididymis between HVS and LVS groups. Gene Ontology (GO) analysis indicated that the DEGs were mainly related to hormone secretion, neurotransmitter synthesis/transport, transmembrane signal transduction, transmembrane transporter activity, and nervous system development (p < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis showed that the DEGs were significantly enriched in pathways associated with hormone and neurotransmitter transmission (p < 0.05). In addition, further analysis of the top five pathways enriched by KEGG, nine key candidate genes (including SLC18A2, SNAP25, CACNA1B, GABRG2, DRD3, CAMK2A, NR5A1, and STAR) were identified, which could play a crucial role in the formation of semen. These data provide new insights into the molecular mechanism regulating semen volume of drakes and make feasible the breeding of drakes by semen volume.

Transcriptome sequencing analysis of the role of miR-499-5p and SOX6 in chicken skeletal myofiber specification

MicroRNAs (miRNAs) might play critical roles in skeletal myofiber specification. In a previous study, we found that chicken miR-499-5p is specifically expressed in slow-twitch muscle and that its potential target gene is SOX6. In this study, we performed RNA sequencing to investigate the effects of SOX6 and miR-499-5p on the modulation and regulation of chicken muscle fiber type and its regulatory mechanism. The expression levels of miR-499-5p and SOX6 demonstrated opposing trends in different skeletal muscles and were associated with muscle fiber type composition. Differential expression analysis revealed that miR-499-5p overexpression led to significant changes in the expression of 297 genes in chicken primary myoblasts (CPMs). Myofiber type-related genes, including MYH7B and CSRP3, showed expression patterns similar to those in slow-twitch muscle. According to functional enrichment analysis, differentially expressed genes were mostly associated with muscle development and muscle fiber-related processes. SOX6 was identified as the target gene of miR-499-5p in CPM using target gene mining and luciferase reporter assays. SOX6 knockdown resulted in upregulation of the slow myosin genes and downregulation of fast myosin genes. Furthermore, protein-protein interaction network analysis revealed that MYH7B and RUNX2 may be the direct targets of SOX6. These results indicated that chicken miR-499-5p may promote slow-twitch muscle fiber formation by repressing SOX6 expression. Our study provides a dataset that can be used as a reference for animal meat quality and human muscle disease studies.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Long Noncoding RNAs: Recent Insights into Their Role in Male Infertility and Their Potential as Biomarkers and Therapeutic Targets

Long noncoding RNAs (lncRNAs) are composed of nucleotides located in the nucleus and cytoplasm; these are transcribed by RNA polymerase II and are greater than 200 nt in length. LncRNAs fulfill important functions in a variety of biological processes, including genome imprinting, cell differentiation, apoptosis, stem cell pluripotency, X chromosome inactivation and nuclear transport. As high throughput sequencing technology develops, a substantial number of lncRNAs have been found to be related to a variety of biological processes, such as development of the testes, maintaining the self-renewal and differentiation of spermatogonial stem cells, and regulating spermatocyte meiosis. These indicate that lncRNAs can be used as biomarkers and potential therapeutic targets for male infertility. However, only a few comprehensive reviews have described the role of lncRNAs in male reproduction. In this paper, we summarize recent findings relating to the role of lncRNAs in spermatogenesis, their potential as biomarkers for male infertility and the relationship between reproductive arrest and transgenerational effects. Finally, we suggest specific targets for the treatment of male infertility from the perspective of lncRNAs.

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.

Alterations of mRNA and lncRNA profiles associated with the extracellular matrix and spermatogenesis in goats

Objective

Spermatozoa are produced within the seminiferous tubules after sexual maturity. The expression levels of mRNAs and lncRNAs in testicular tissues are different at each stage of testicular development and are closely related to formation of the extracellular matrix (ECM) and spermatogenesis. Therefore, we set out to study the expression of lncRNAs and mRNAs during the different developmental stages of the goat testis.

Methods

We constructed 12 RNA libraries using testicular tissues from goats aged 3, 6, and 12 months, and studied the functions of mRNAs and lncRNAs using the gene ontogeny (GO) and Kyoto encyclopedia of genes and genomes (KEGG) databases. Relationships between differentially expressed genes (DEGs) were analyzed by lncRNA-mRNA co-expression network and protein-protein interaction network (PPI). Finally, the protein expression levels of matrix metalloproteinase 2 (MMP2), insulin-like growth factor 2 (IGF2), and insulin-like growth factor-binding protein 6 (IGFBP6) were detected by western blotting.

Results

We found 23, 8, and 135 differentially expressed lncRNAs and 161, 12, and 665 differentially expressed mRNAs that were identified between 3 vs 6, 6 vs 12, and 3 vs 12 months, respectively. GO, KEGG, and PPI analyses showed that the differential genes were mainly related to the ECM. Moreover, MMP2 was a hub gene and co-expressed with the lncRNA TCONS-0002139 and TCONS-00093342. The results of quantitative reverse-transcription polymerase chain reaction verification were consistent with those of RNA-seq sequencing. The expression trends of MMP2, IGF2, and IGFBP6 protein were the same as that of mRNA, which all decreased with age. IGF2 and MMP2 were significantly different in the 3 vs 6-month-old group (p<0.05).

Conclusion

These results improve our understanding of the molecular mechanisms involved in sexual maturation of the goat testis.

Integrated Analysis of Long Non-Coding RNA and mRNA Expression Profiles in Testes of Calves and Sexually Mature Wandong Bulls (Bos taurus)

The mRNAs and long non-coding RNAs axes are playing a vital role in the regulating of post-transcriptional gene expression. Thereby, elucidating the expression pattern of mRNAs and long non-coding RNAs underlying testis development is crucial. In this study, mRNA and long non-coding RNAs expression profiles were investigated in 3-month-old calves and 3-year-old mature bulls' testes by total RNA sequencing. Additionally, during the gene level analysis, 21,250 mRNAs and 20,533 long non-coding RNAs were identified. As a result, 7908 long non-coding RNAs (p-adjust < 0.05) and 5122 mRNAs (p-adjust < 0.05) were significantly differentially expressed between the distinct age groups. In addition, gene ontology and biological pathway analyses revealed that the predicted target genes are enriched in the lysine degradation, cell cycle, propanoate metabolism, adherens junction and cell adhesion molecules pathways. Correspondingly, the RT-qPCR validation results showed a strong consistency with the sequencing data. The source genes for the mRNAs (CCDC83, DMRTC2, HSPA2, IQCG, PACRG, SPO11, EHHADH, SPP1, NSD2 and ACTN4) and the long non-coding RNAs (COX7A2, COX6B2, TRIM37, PRM2, INHBA, ERBB4, SDHA, ATP6VOA2, FGF9 and TCF21) were found to be actively associated with bull sexual maturity and spermatogenesis. This study provided a comprehensive catalog of long non-coding RNAs in the bovine testes and also offered useful resources for understanding the differences in sexual development caused by the changes in the mRNA and long non-coding RNA interaction expressions between the immature and mature stages.

Long Non-coding RNAs: Mechanisms, Experimental, and Computational Approaches in Identification, Characterization, and Their Biomarker Potential in Cancer

Long non-coding RNAs are diverse class of non-coding RNA molecules >200 base pairs of length having various functions like gene regulation, dosage compensation, epigenetic regulation. Dysregulation and genomic variations of several lncRNAs have been implicated in several diseases. Their tissue and developmental specific expression are contributing factors for them to be viable indicators of physiological states of the cells. Here we present an comprehensive review the molecular mechanisms and functions, state of the art experimental and computational pipelines and challenges involved in the identification and functional annotation of lncRNAs and their prospects as biomarkers. We also illustrate the application of co-expression networks on the TCGA-LIHC dataset for putative functional predictions of lncRNAs having a therapeutic potential in Hepatocellular carcinoma (HCC).

Transcriptomes of an Array of Chicken Ovary, Intestinal, and Immune Cells and Tissues

While the chicken (Gallus gallus) is the most consumed agricultural animal worldwide, the chicken transcriptome remains understudied. We have characterized the transcriptome of 10 cell and tissue types from the chicken using RNA-seq, spanning intestinal tissues (ileum, jejunum, proximal cecum), immune cells (B cells, bursa, macrophages, monocytes, spleen T cells, thymus), and reproductive tissue (ovary). We detected 17,872 genes and 24,812 transcripts across all cell and tissue types, representing 73% and 63% of the current gene annotation, respectively. Further quantification of RNA transcript biotypes revealed protein-coding and lncRNAs specific to an individual cell/tissue type. Each cell/tissue type also has an average of around 1.2 isoforms per gene, however, they all have at least one gene with at least 11 isoforms. Differential expression analysis revealed a large number of differentially expressed genes between tissues of the same category (immune and intestinal). Many of these differentially expressed genes in immune cells were involved in cellular processes relating to differentiation and cell metabolism as well as basic functions of immune cells such as cell adhesion and signal transduction. The differential expressed genes of the different segments of the chicken intestine (jejunum, ileum, proximal cecum) correlated to the metabolic processes in nutrient digestion and absorption. These data should provide a valuable resource in understanding the chicken genome.

Analysis of long non-coding RNAs in epididymis of cattleyak associated with male infertility

Cattleyak (CY), is a cross breed between cattle and yak (YK), which display equal adaptability to the harsh environment as YK and much higher performances than YK. However, the CY is female fertile and male sterile. Previous studies were conducted on testes tissues to investigate the mechanism of male infertility in CY. There is no systematic research on genes, especially lncRNAs between CY and YK epididymis. In this study, Illumina Hiseq was performed to profile the epididymis transcriptome (lncRNA and mRNA) of CY and YK. In total 18859 lncRNAs were identified, from which lincRNAs 12458, antisense lncRNAs 2345, intronic lncRNAs 3101, and sense lncRNAs 955 respectively. We have identified 345 DE lncRNAs and 3008 DE mRNAs between YK and CY epididymis. Thirteen DEGs were validated by quantitative real-time PCR. Combing with DEG, 14 couples of lncRNAs and their target genes were both DE, and 6 of them including CCDC39, KCNJ16, NECTIN2, MRPL20, PSMC4, and DEFB112 show their potential infertility-related terms such as cellular motility, sperm maturation, sperm storage, cellular junction, folate metabolism, and capacitation. On the other hand, several down-regulated genes such as DEFB124, DEFB126, DEFB125, DEFB127, DEFB129, CES5A, TKDP1, CST3, RNASE9 and CD52 in CY compared to YK were involved in the immune response and sperm maturation. Therefore, comprehensive analysis for lncRNAs and their target genes may enhance our understanding of the molecular mechanisms underlying the process of sperm maturation in CY and may provide important resources for further research.

Transcriptome analysis of turkey (Meleagris gallopavo) reproductive tract revealed key pathways regulating spermatogenesis and post-testicular sperm maturation

The application of transcriptomics to the study of the reproductive tract in male turkeys can significantly increase our current knowledge regarding the specifics of bird reproduction. To characterize the complex transcriptomic changes that occur in the testis, epididymis, and ductus deferens, deep sequencing of male turkey RNA samples (n = 6) was performed, using Illumina RNA-Seq. The obtained sequence reads were mapped to the turkey genome, and relative expression values were calculated to analyze differentially expressed genes (DEGs). Statistical analysis revealed 1,682; 2,150; and 340 DEGs in testis/epididymis, testis/ductus deferens, and epididymis/ductus deferens comparisons, respectively. The expression of selected genes was validated using quantitative real-time reverse transcriptase-polymerase chain reaction. Bioinformatics analysis revealed several potential candidate genes involved in spermatogenesis, spermiogenesis and flagellum formation in the testis, and in post-testicular sperm maturation in the epididymis and ductus deferens. In the testis, genes were linked with the mitotic proliferation of spermatogonia and the meiotic division of spermatocytes. Histone ubiquitination and protamine phosphorylation were shown to be regulatory mechanisms for nuclear condensation during spermiogenesis. The characterization of testicular transcripts allowed a better understanding of acrosome formation and development and flagellum formation, including axoneme structures and functions. Spermatozoa motility during post-testicular maturation was linked to the development of flagellar actin filaments and biochemical processes, including Ca2+ influx and protein phosphorylation/dephosphorylation. Spermatozoa quality appeared to be controlled by apoptosis and antioxidant systems in the epididymis and ductus deferens. Finally, genes associated with reproductive system development and morphogenesis were identified. To the best of our knowledge, this is the first genome-wide functional investigation of genes associated with tissue-specific processes in turkey reproductive tract. A catalog of genes worthy of further studies to understand the avian reproductive physiology and regulation was provided.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

AFB1 Induced Transcriptional Regulation Related to Apoptosis and Lipid Metabolism in Liver of Chicken

Aflatoxin B1 (AFB1) leads to a major risk to poultry and its residues in meat products can also pose serious threat to human health. In this study, after feeding 165-day-old Roman laying hens for 35 days, the toxic effects of aflatoxin B1 at different concentrations were evaluated. The purpose of this study was to explore the mechanism of liver toxicosis responses to AFB1. We found that highly toxic group exposure resulted in liver fat deposition, increased interstitial space, and hepatocyte apoptosis in laying hens. Furthermore, a total of 164 differentially expressed lnRNAs and 186 differentially expressed genes were found to be highly correlated (Pearson Correlation Coefficient > 0.80, p-value < 0.05) by sequencing the transcriptome of control (CB) and highly toxic group (TB3) chickens. We also identify 29 differentially expressed genes and 19 miRNAs that have targeted regulatory relationships. Based on the liver cell apoptosis and fatty liver syndrome that this research focused on, we found that the highly toxic AFB1 led to dysregulation of the expression of PPARG and BCL6. They are cis-regulated by TU10057 and TU45776, respectively. PPARG was the target gene of gga-miR-301a-3p, gga-miR-301b-3p, and BCL6 was the target gene of gga-miR-190a-3p. In summary, highly toxic AFB1 affects the expression levels of protein-coding genes and miRNAs in the liver of Roman layer hens, as well as the expression level of long non-coding RNA in the liver, which upregulates the expression of PPARG and downregulates the expression of Bcl-6. Our study provides information on possible genetic regulatory networks in AFB1-induced hepatic fat deposition and hepatocyte apoptosis.

Genome-wide analysis of differentially expressed profiles of mRNAs, lncRNAs and circRNAs in chickens during Eimeria necatrix infection

BACKGROUND

Eimeria necatrix, the most highly pathogenic coccidian in chicken small intestines, can cause high morbidity and mortality in susceptible birds and devastating economic losses in poultry production, but the underlying molecular mechanisms in interaction between chicken and E. necatrix are not entirely revealed. Accumulating evidence shows that the long-non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are key regulators in various infectious diseases. However, the expression profiles and roles of these two non-coding RNAs (ncRNAs) during E. necatrix infection are still unclear.

METHODS

The expression profiles of mRNAs, lncRNAs and circRNAs in mid-segments of chicken small intestines at 108 h post-infection (pi) with E. necatrix were analyzed by using the RNA-seq technique.

RESULTS

After strict filtering of raw data, we putatively identified 49,183 mRNAs, 818 lncRNAs and 4153 circRNAs. The obtained lncRNAs were classified into four types, including 228 (27.87%) intergenic, 67 (8.19%) intronic, 166 (20.29%) anti-sense and 357 (43.64%) sense-overlapping lncRNAs; of these, 571 were found to be novel. Five types were also predicted for putative circRNAs, including 180 exonic, 54 intronic, 113 antisense, 109 intergenic and 3697 sense-overlapping circRNAs. Eimeria necatrix infection significantly altered the expression of 1543 mRNAs (707 upregulated and 836 downregulated), 95 lncRNAs (49 upregulated and 46 downregulated) and 13 circRNAs (9 upregulated and 4 downregulated). Target predictions revealed that 38 aberrantly expressed lncRNAs would cis-regulate 73 mRNAs, and 1453 mRNAs could be trans-regulated by 87 differentially regulated lncRNAs. Additionally, 109 potential sponging miRNAs were also identified for 9 circRNAs. GO and KEGG enrichment analysis of target mRNAs for lncRNAs, and sponging miRNA targets and source genes for circRNAs identified associations of both lncRNAs and circRNAs with host immune defense and pathogenesis during E. necatrix infection.

CONCLUSIONS

To the best of our knowledge, the present study provides the first genome-wide analysis of mRNAs, lncRNAs and circRNAs in chicken small intestines infected with E. necatrix. The obtained data will offer novel clues for exploring the interaction mechanisms between chickens and Eimeria spp.

Analysis of Long Noncoding RNA and mRNA Expression Profiles of Testes with High and Low Sperm Motility in Domestic Pigeons (Columba livia)

Sperm motility is one of the most important indicators in assessing semen quality, and it is used to evaluate poultry fertility. Many long noncoding RNAs (lncRNAs) and mRNAs are involved in regulating testis development and spermatogenesis. In this study, we employed RNA sequencing to analyse the testis transcriptome (lncRNA and mRNA) of ten pigeons with high and low sperm motility. In total, 46,117 mRNAs and 17,463 lncRNAs were identified, of which 2673 mRNAs and 229 lncRNAs (P < 0.05) were significantly differentially expressed (DE) between the high and low sperm motility groups. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analysis showed that target genes of DE lncRNAs and DE mRNAs were related to calcium ion binding, ATP binding, and spermatogenesis. Moreover, we found that UBB, a target gene of lncRNA MSTRG.7787.5, was involved in germ cell development. Our study provided a catalogue of lncRNAs and mRNAs associated with sperm motility, and they deserve further study to deepen the understanding of biological processes in the pigeon testis.

Revealing stage-specific expression patterns of long noncoding RNAs along mouse spermatogenesis

The discovery of a large number of long noncoding RNAs (lncRNAs), and the finding that they may play key roles in different biological processes, have started to provide a new perspective in the understanding of gene regulation. It has been shown that the testes express the highest amount of lncRNAs among different vertebrate tissues. However, although some studies have addressed the characterization of lncRNAs along spermatogenesis, an exhaustive analysis of the differential expression of lncRNAs at its different stages is still lacking. Here, we present the results for lncRNA transcriptome profiling along mouse spermatogenesis, employing highly pure flow sorted spermatogenic stage-specific cell populations, strand-specific RNAseq, and a combination of up-to-date bioinformatic pipelines for analysis. We found that the vast majority of testicular lncRNA genes are expressed at post-meiotic stages (i.e. spermiogenesis), which are characterized by extensive post-transcriptional regulation. LncRNAs at different spermatogenic stages shared common traits in terms of transcript length, exon number, and biotypes. Most lncRNAs were lincRNAs, followed by a high representation of antisense (AS) lncRNAs. Co-expression analyses showed a high correlation along the different spermatogenic stage transitions between the expression patterns of AS lncRNAs and their overlapping protein-coding genes, raising possible clues about lncRNA-related regulatory mechanisms. Interestingly, we observed the co-localization of an AS lncRNA and its host sense mRNA in the chromatoid body, a round spermatids-specific organelle that has been proposed as a reservoir of RNA-related regulatory machinery. An additional, intriguing observation is the almost complete lack of detectable expression for Y-linked testicular lncRNAs, despite that a high number of lncRNA genes are annotated for this chromosome.

Photoperiodic effect on the testicular transcriptome in broiler roosters

Information about the effects of photoperiod on the testicular transcriptome of broiler roosters is limited. The aim of the present study was to explore the effect of different photoperiodic regimes on gene expression in the testes of broiler breeder roosters. One hundred and twenty Arbor Acres broiler breeder roosters aged 20 weeks were assigned to one of three groups (n = 40) and subjected to different photoperiodic regimes: control (CTR; 12.5 L:11.5 D), short day (SD; 8 L:16 D) and long day (LD; 16 L:8 D). After 4 weeks, the testes of 10 randomly selected birds from each group were dissected, sliced and haematoxylin-eosin stained. The testicular transcriptome of roosters from the SD and LD groups was determined by RNA sequencing (RNA-Seq), and the results were confirmed using quantitative real-time PCR. The seminiferous tubule area and sperm count increased significantly with the prolongation of photoperiod (p < .01). Additionally, the RNA-Seq results indicated that 387 genes were upregulated and 1,052 genes were downregulated in the LD group compared with those in the SD group. Several crucial genes involved in rooster testicular development and reproduction were also screened, including heat shock proteins 90, extracellular regulated protein kinases 1, phosphatidylinositol 3-kinase, adenosine 5'-monophosphate -activated protein kinase, BCL-6 and Smad3. The differentially expressed genes were enriched in the mammalian targets of rapamycin (mTOR), forkhead box (FoxO), transforming growth factor beta (TGF-β) and insulin signalling pathway. In conclusion, a 16 hr photoperiod for 4 weeks increased the seminiferous tubule duct area and promoted spermatogenesis in the rooster's testicles, and the mTOR, FoxO, TGF-β and insulin signalling pathways may be involved.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.

RESULTS

A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.

CONCLUSION

The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Genome-wide analysis of mRNAs and lncRNAs in Mycoplasma bovis infected and non-infected bovine mammary gland tissues

Mycoplasma bovis (M. bovis) causes diseases such as arthritis, pneumonia, abortion, and mastitis, leading to great losses in the bovine dairy industries. RNA types such as messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) play significant roles in regulating the immune responses triggered by bacteria. The expression profiles of mRNA and lncRNA as they occur in bovine mammary gland tissues infected with M. bovis are still not well understood. To illuminate this issue, transcription analysis of mRNA and LncRNAs were conducted on the mammary gland tissues belonging to Holstein cattle infected and not infected with M. bovis. The analysis revealed 1310 differentially expressed mRNAs and 57 differentially expressed lncRNAs in the bovine mammary gland tissues infected and not infected with M. bovis. In addition, 392 novel lncRNAs were detected, 19 of which were differentially expressed. Gene ontology analysis reveals that differentially expressed mRNAs and lncRNAs play significant roles in such vital biological pathways as metabolic pathways, T-cell receptor signaling, TGF-beta signaling, pathways in cancer, PI3K-Akt signaling, NF-kappa B signaling, mTOR signaling, and apoptosis, including in the immune response to cancer. Based on our literature review, this study is the first genome-wide lncRNA research conducted on bovine mammary gland tissues infected with M. bovis. Our results provide bovine mammary gland lncRNA and mRNA resources to understand their roles in the regulation of the immune response against the agent M. bovis in bovine mammary gland tissues.

Digital gene expression profiling and validation study highlight Cyclin F as an important regulator for sperm motility of chickens

In poultry industry, around 5 to 12% roosters were eliminated from the breeding program because of low sperm motility. Relatively few studies have been directed toward understanding and explaining the genetics mechanisms involved in sperm motility regulation in chickens. In the present study, digital gene expression (DGE) profiling and bioinformation analysis were used to explore the globally differentially expressed genes (DEG) in the testis of low sperm motility and high sperm motility roosters. Further validation study of key candidate genes was also performed. The DGE identified 652 DEGs, including 473 up-regulated and 179 down-regulated genes in the low sperm motility testis. Those DEGs were enriched on 21 terms of biological process category, 10 terms of cellular component category, including motile cilium, and 13 terms of molecular function category including microtubule motor activity and ATP binding. The kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis indicated that these DEGs were involved in the FoxO signaling pathway and insulin resistance pathway. Quantitative real time PCR (qRT-PCR) studies of 8 DEGs were used to validate the DGE results. A key candidate gene Cyclin F (CCNF) was extremely low expressed in the low sperm motility testis (log2 ratio (low sperm motility/high sperm motility) = -5.23). The CCNF gene silencing in the chicken DF-1 cell line induced the reduced cell activity and proliferation. In summary, the present study provides insight into the potential genetic regulation of sperm motility and highlighted the underlying pathways (Insulin resistance and FoxO signaling pathways) and important candidate genes such as CCNF.

Genome-wide re-sequencing and transcriptome analysis reveal candidate genes associated with the pendulous comb phenotype in domestic chickens

To determine the causative variations associated with two chicken comb phenotypes, pendulous comb (PC) or upright comb (UC), two pooled genomic DNA samples from PC and UC chickens were re-sequenced by Next-Generation Sequencer, and genome-wide Single nucleotide polymorphisms (SNPs) were detected. Using three selective sweep approaches, F_ST , θπ, and Tajima's D, with top 5% window values serving as the threshold, a total of 84 positively selective genes (PSGs) were identified. There were no SNPs in exons of the PSGs with significant differences in allele frequencies between the two comb phenotype groups. Then, 515 differentially expressed genes (DEGs) between the PC and UC were identified by RNA-seq. Three genes including CD36 (CD36 molecule), ADAMTSL3 (ADAMTS-like 3), and AOX1 (aldehyde oxidases 1) are overlapped between PSGs and DEGs. After genotyping seven candidate SNPs in the regulatory regions of the three overlapping genes in 120 chickens from two other breeds, two variants (rs14607046 and rs731818051) in the regulatory regions of AOX1 and ADAMTSL3 were found to have significant differences in allele frequency between the PC and UC, suggesting that the two variants may be causative mutations for PC. Overall, our study shed light on the genetic basis underlying the PC phenotype in chickens.

Hippocampus-specific regulation of long non-coding RNA and mRNA expression in germ-free mice

Gut microbiota can affect multiple brain functions and cause behavioral alterations through the microbiota-gut-brain axis. In our previous study, we found that the absence of gut microbiota can influence the expression of microRNAs and mRNAs in the hippocampal region of the germ-free (GF) mice. Long non-coding RNAs (lncRNAs) are increasingly being recognized as an important functional transcriptional regulator in the brain. In the present study, we aim to identify possible biological pathways and functional networks for lncRNA-associated transcript of the gut microbiota in relation to the brain function. The profiles of lncRNA and mRNA from specific pathogen-free (SPF), colonized GF (CGF), and GF mice were generated using the Agilent Mouse LncRNA Array v2.0. Differentially expressed (DE) lncRNAs and mRNAs were identified, and lncRNA target genes were also predicted. Ingenuity pathway analysis (IPA) was performed to analyze related signaling pathways and biological functions associated with these dysregulated mRNAs and target genes. Validation with quantitative real-time PCR was performed on several key genes. Compared with SPF mice a total of 2230 DE lncRNAs were found in GF mice. Among these, 1355 were upregulated and 875 were downregulated. After comparing the target genes of DE lncRNAs with mRNA datasets, 669 overlapping genes were identified. IPA core analyses revealed that most of these genes were highly associated with cardiac hypertrophy, nuclear factors of activated T cells (NFAT) gonadotropin-releasing hormone (GnRH), calcium, and cAMP-response element-binding protein (CREB) signaling pathways. Additionally, mRNA expression levels of APP, CASP9, IGFBP2, PTGDS, and TGFBR2 genes that are involved in central nervous system functions were significantly changed in the GF mouse hippocampus. Through this study, for the first time, we describe the effect of gut microbiota on the hippocampal lncRNA regulation. This will help in enhancing the overall knowledge about microbiota-gut-brain axis.

Association of single nucleotide polymorphism at long non-coding RNA 8138.1 with duration of fertility in egg-laying hens

A previous genome-wide transcriptional analysis identified long non-coding RNA 8138.1 (lncRNA8138.1) as a candidate gene related to hen duration of the fertility (DF) trait. LncRNA8138.1 gene response to growth factor and reproductive system development suggests it has a vital role in reproduction. In this study, we investigated the lncRNA8138.1 gene sequence in a population of egg-laying hens. The sequence analysis of the lncRNA8138.1 gene containing about 1.6 k nucleotides (nt) was observed with four single nucleotide polymorphisms (SNPs) and 7 nt indel including r.4937159A > G; r.4937219T > C; r.4937258G > C; r.4937318C > G and g.4937319_4937325delinsTGTGTGG. Next, the genomic DNAs from laying hen populations were subjected to polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) to detect a region of 457 bp carrying lncRNA8138.1 r.4937159A > G substitution. Further inspection of the region containing r.4937159A > G mutation revealed three genotypes viz., AA, AG, and GG were observed with respective frequencies of 0.106, 0.607, and 0.287 in laying hen population 1 (P1) (n = 1, 042) and respective frequencies of 0.176, 0.708, and 0.116 in laying hen population 2 (P2) (n = 826). Moreover, to further examining the frequencies of r.4937159A > G genotypes in P1 and P2, and their additive and dominance effects; r.4937159A > G locus was significantly associated with DF-trait in both P1 and P2 (EN: the number of eggs, FN: the number of fertile eggs after a single AI), and DN (the number of days post-insemination until last fertile egg). In testing for additive and dominance effects, additive effect was significant (P < 0.05) in both P1 and P2 for DF-trait, and the dominance effect was significant (P < 0.05) for EN and FN traits, suggesting that r.4937159A > G polymorphism is a potential biomarker for DF-trait. However, the identified novel r.4937159A > G mutation and others require further investigation to confirm phenotypic causality and potential genetic relationships with reproductive traits. Overall, our findings suggest the significance of genetic variation in long non-coding RNAs may assist in future breeding programs to improve selection for prolonged DF-trait.

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

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

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

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

Genome-wide identification and characterization of long non-coding RNAs during postnatal development of rabbit adipose tissue

Background

The rabbit is widely used as an important experimental model for biomedical research, and shows low adipose tissue deposition during growth. Long non-coding RNAs (lncRNAs) are associated with adipose growth, but little is known about the function of lncRNAs in the rabbit adipose tissue.

Methods

Deep RNA-sequencing and comprehensive bioinformatics analyses were used to characterize the lncRNAs of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed (DE) lncRNAs were identified at the three growth stages by DESeq. The cis and trans prediction ways predicted the target genes of the DE lncRNAs. To explore the function of lncRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the candidate genes.

Results

A total of 991,157,544 clean reads were generated after RNA-Seq of the three growth stages, of which, 30,353 and 107 differentially expressed (DE) lncRNAs were identified. Compared to the protein-coding transcripts, the rabbit lncRNAs shared some characteristics such as shorter length and fewer exons. Cis and trans target gene prediction revealed, 43 and 64 DE lncRNAs respectively, corresponding to 72 and 20 protein-coding genes. GO enrichment and KEGG pathway analyses revealed that the candidate DE lncRNA target genes were involved in oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, and other adipose growth-related pathways. Six DE lncRNAs were randomly selected and validated by q-PCR.

Conclusions

This study is the first to profile the potentially functional lncRNAs in the adipose tissue growth in rabbits, and contributes to our understanding of mammalian adipogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12944-018-0915-1) contains supplementary material, which is available to authorized users.

Transcriptome sequencing reveals key potential long non-coding RNAs related to duration of fertility trait in the uterovaginal junction of egg-laying hens

Duration of fertility, (DF) is an important functional trait in poultry production and lncRNAs have emerged as important regulators of various process including fertility. In this study we applied a genome-guided strategy to reconstruct the uterovaginal junction (UVJ) transcriptome of 14 egg-laying birds with long- and short-DF (n = 7); and sought to uncover key lncRNAs related to duration of fertility traits by RNA-sequencing technology. Examination of RNA-seq data revealed a total of 9977 lncRNAs including 2576 novel lncRNAs. Differential expression (DE) analysis of lncRNA identified 223 lncRNAs differentially expressed between the two groups. DE-lncRNA target genes prediction uncovered over 200 lncRNA target genes and functional enrichment tests predict a potential function of DE-lncRNAs. Gene ontology classification and pathway analysis revealed 8 DE-lncRNAs, with the majority of their target genes enriched in biological functions such as reproductive structure development, developmental process involved in reproduction, response to cytokine, carbohydrate binding, chromatin organization, and immune pathways. Differential expression of lncRNAs and target genes were confirmed by qPCR. Together, these results significantly expand the utility of the UVJ transcriptome and our analysis identification of key lncRNAs and their target genes regulating DF will form the baseline for understanding the molecular functions of lncRNAs regulating DF.

Epigenetic reprogramming during spermatogenesis and male factor infertility

Infertility is an often devastating diagnosis encountered by around one in six couples who are trying to conceive. Moving away from the long-held belief that infertility is primarily a female issue, it is now recognised that half, if not more, of these cases may be due to male factors. Recent evidence has suggested that epigenetic abnormalities in chromatin dynamics, DNA methylation or sperm-borne RNAs may contribute to male infertility. In light of advances in deep sequencing technologies, researchers have been able to increase the coverage and depth of sequencing results, which in turn has allowed more comprehensive analyses of spermatozoa chromatin dynamics and methylomes and enabled the discovery of new subsets of sperm RNAs. This review examines the most current literature related to epigenetic processes in the male germline and the associations of aberrant modifications with fertility and development.

Identification of molecular pathways and candidate genes associated with cocks' comb size trait by genome-wide transcriptome analysis

The comb of the male is an important secondary sexual characteristic. Although quantitative trait loci (QTLs) related to comb size have been identified, molecular mechanisms underlying this trait remain mostly unknown. In this study, RNA sequencing (RNA-seq) was employed to compare whole transcriptomic differences between two groups of Partridge Shank chickens that are divergent in comb sizes. A total of 563 differentially expressed genes (DEGs) were identified, including 277 up-regulated and 286 down-regulated DEGs. According to the animal QTL database, eight DEGs including BMP2 and CHADL matching the reported QTLs were associated with the comb size. Functional annotation analysis revealed that DEGs were involved in cell communication and calcium signaling. Protein-protein interaction network analysis showed that STK32A, PIK3R1, EDN1, HSPA5, and HSPA8 have an impact on comb growth. Moreover, potential alternative splicing events and single nucleotide polymorphisms were also identified. Our data provide a source for identifying genes and pathways with functions critical to comb size and accelerate studies involving molecular mechanisms of this sexual ornament.

Advances of Long Noncoding RNAs-mediated Regulation in Reproduction

OBJECTIVE

Advances in genomics and molecular biology have led to the discovery of a large group of uncharacterized long noncoding RNAs (lncRNAs). Emerging evidence indicated that many lncRNAs function in multiple biological processes and its dysregulation often causes diseases. Recent studies suggested that almost all regulatory lncRNAs interact with biological macromolecules such as DNA, RNA, and protein. LncRNAs regulate gene expression mainly on three levels, including epigenetic modification, transcription, and posttranscription, through DNA methylation, histone modification, and chromatin remodeling. LncRNAs can also affect the development of diseases and therefore be used to diagnose and treat diseases. With new sequencing and microarray techniques, hundreds of lncRNAs involved in reproductive disorders have been identified, but their functions in these disorders are undefined.

DATA SOURCES

This review was based on articles published in PubMed databases up to July 10, 2017, with the following keywords: "long noncoding RNAs", "LncRNA", "placentation", and "reproductive diseases".

STUDY SELECTION

Original articles and reviews on the topics were selected.

RESULTS

LncRNAs widely participate in various physiological and pathological processes as a new class of important regulatory factors. In spermatogenesis, spermatocytes divide and differentiate into mature spermatozoa. The whole process is elaborately regulated by the expression of phase-specific genes that involve many strains of lncRNAs. Literature showed that lncRNA in reproductive cumulus cells may contribute to the regulation of oocyte maturation, fertilization, and embryo development.

CONCLUSIONS

LncRNA has been found to play a role in the development of reproduction. Meanwhile, we reviewed the studies on how lncRNAs participate in reproductive disorders, which provides a basis for the study of lncRNA in reproduction regulation.

Regulatory role of long non-coding RNAs during reproductive disease

Long non-coding RNA (lncRNA) is a group of RNAs with broad biogenesis, which are longer than 200 nt and highly conserved in their secondary and tertiary structures. lncRNA that broadly participates in varied physiological processes in organisms has abundant biological function and can regulate expression of target genes at transcriptional, post-transcriptional and epigenetic levels. LncRNAs can also affect the development of diseases, and therefore be used to diagnose and treat diseases. With new sequencing and microarray techniques, hundreds of lncRNAs involved in reproductive disorders have been identified, but their functions in these disorders are undefined. In this paper, we reviewed the studies on how lncRNAs participate in the development of reproductive disorders, hoping our outcome can instruct the future study and provide new biomarkers and therapies for reproductive disorders.

Dissection of Myogenic Differentiation Signatures in Chickens by RNA-Seq Analysis

A series of elaborately regulated and orchestrated changes in gene expression profiles leads to muscle growth and development. In this study, RNA sequencing was used to profile embryonic chicken myoblasts and fused myotube transcriptomes, long non-coding RNAs (lncRNAs), and messenger RNAs (mRNAs) at four stages of myoblast differentiation. Of a total of 2484 lncRNA transcripts, 2288 were long intergenic non-coding RNAs (lincRNAs) and 198 were antisense lncRNAs. Additionally, 1530 lncRNAs were neighboring 2041 protein-coding genes (<10 kb upstream and downstream) and functionally enriched in several pathways related to skeletal muscle development that have been extensively studied, indicating that these genes may be in cis-regulatory relationships. In addition, Pearson’s correlation coefficients demonstrated that 990 lncRNAs and 7436 mRNAs were possibly in trans-regulatory relationships. These co-expressed mRNAs were enriched in various developmentally-related biological processes, such as myocyte proliferation and differentiation, myoblast differentiation, and myoblast fusion. The number of transcripts (906 lncRNAs and 4422 mRNAs) differentially expressed across various stages declined with the progression of differentiation. Then, 4422 differentially expressed genes were assigned to four clusters according to K-means analysis. Genes in the K1 cluster likely play important roles in myoblast proliferation and those in the K4 cluster were likely associated with the initiation of myoblast differentiation, while genes in the K2 and K3 clusters were likely related to myoblast fusion. This study provides a catalog of chicken lncRNAs and mRNAs for further experimental investigations and facilitates a better understanding of skeletal muscle development.

Long noncoding RNAs regulate Wnt signaling during feather regeneration

Long noncoding RNAs (lncRNAs) are non-protein coding transcripts that are involved in a broad range of biological processes. Here, we examined the functional roles of lncRNAs in feather regeneration. RNA-seq profiling of the regenerating feather blastema revealed that the Wnt signaling is among the most active pathways during feather regeneration, with the Wnt ligands and their inhibitors showing distinct expression patterns. Co-expression analysis identified hundreds of lncRNAs with similar expression patterns to either the Wnt ligands (the Lwnt group) or their downstream target genes (the Twnt group). Among these, we randomly picked two lncRNAs in the Lwnt group, and three lncRNAs in the Twnt group to validate their expression and function. Members in the Twnt group regulated feather regeneration and axis formation, whereas members in the Lwnt group showed no obvious phenotype. Further analysis confirmed that the three Twnt group members inhibit Wnt signal transduction and at the same time are down-stream target genes of this pathway. Our results suggested that the feather regeneration model can be utilized to systematically annotate the functions of lncRNAs in the chicken genome.