Comprehensive analysis of long noncoding RNA and mRNA expression patterns in sheep testicular maturation

Analysis of histomorphology and SERNINA5 gene expression in different regions of epididymis of cattleyak

The molecular mechanism of sterility in cattleyak is still unresolved. The related factors of infertility in cattleyak were studied by tissue section, SERPINA5 gene cloning and bioinformatics technology. Tissue sections of the epididymis showed poorly structured and disorganized epithelial cells in the corpus of the epididymis compared to the caput of the epididymis, while in the cauda part of the epididymis, the extra basal smooth muscle was thinner, the surface of the epithelial lumen was discontinuous and the epithelium was markedly degenerated. The results of gene cloning showed that the coding sequence (CDS) region of the SERPINA5 gene in cattleyak was 1215 bp in length, encoding a total of 404 amino acids, of which the isoleucine content was the highest, accounting for a total of 49 amino acids (12.1%). The results of real-time fluorescence quantitative PCR (qPCR) showed that the expression of the SERPINA5 gene in the epididymis caput in cattleyak was significantly higher than that in the corpus and cauda (P < 0.05), but there were no significant differences between the corpus and cauda. In the current study, histological and bioinformatics analysis, physicochemical properties, and the expression analysis of the SERPINA5 gene in different regions of the epididymis in cattleyak were carried out to explore the biological complications of cattleyak infertility.

A meta-analysis approach for annotation and identification of lncRNAs controlling perirenal fat deposition in suckling lambs

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals.

Genome-wide association study considering genotype-by-environment interaction for productive and reproductive traits using whole-genome sequencing in Nellore cattle

BACKGROUND

The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.

RESULTS

Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.

CONCLUSIONS

GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.

Sperm long non-coding RNAs as markers for ram fertility

It is critical in sheep farming to accurately estimate ram fertility for maintaining reproductive effectiveness and for production profitability. However, there is currently a lack of reliable biomarkers to estimate semen quality and ram fertility, which is hindering advances in animal science and technology. The objective of this study was to uncover long non-coding RNAs (lncRNAs) in sperm from rams with distinct fertility phenotypes. Mature rams were allocated into two groups: high and low fertility (HF; n = 31; 94.5 ± 2.8%, LF; n = 25; 83.1 ± 5.73%; P = 0.028) according to the pregnancy rates sired by the rams (average pregnancy rate; 89.4 ± 7.2%). Total RNAs were isolated from sperm of the highest- and lowest-fertility rams (n = 4, pregnancy rate; 99.2 ± 1.6%, and 73.6 ± 4.4%, respectively) followed by next-generation sequencing of the transcripts. We uncovered 11,209 lncRNAs from the sperm of rams with HF and LF. In comparison to each other, there were 93 differentially expressed (DE) lncRNAs in sperm from the two distinct fertility phenotypes. Of these, 141 mRNAs were upregulated and 134 were downregulated between HF and LF, respectively. Genes commonly enriched for 9 + 2 motile cilium and sperm flagellum were ABHD2, AK1, CABS1, ROPN1, SEPTIN2, SLIRP, and TEKT3. Moreover, CABS1, CCDC39, CFAP97D1, ROPN1, SLIRP, TEKT3, and TTC12 were commonly enriched in flagellated sperm motility and sperm motility. Differentially expressed mRNAs were enriched in the top 16 KEGG pathways. Targets of the differentially expressed lncRNAs elucidate functions in cis and trans manner using the genetic context of the lncRNA locus, and lncRNA sequences revealed 471 mRNAs targets of 10 lncRNAs. This study illustrates the existence of potential lncRNA biomarkers that can be implemented in analyzing the quality of ram sperm and determining the sperm fertility and is used in breeding soundness exams for precision livestock farming to ensure food security on a global scale.

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

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

Comprehensive transcriptomic analysis revealing the regulatory dynamics and networks of the pituitary-testis axis in sheep across developmental stages

Spermatogenesis is a complex process intricately regulated by the hypothalamic-pituitary-testis (HPT) axis. However, research on the regulatory factors governing the HPT axis remains limited. This study addresses this gap by conducting a comprehensive analysis of transcriptomes from the pituitary and testis tissues across various developmental stages, encompassing embryonic day (E120), neonatal period (P0), pre-puberty (P90), and post-puberty day (P270). Utilizing edgeR and WGCNA, we identified stage-specific genes in both the pituitary and testis throughout the four developmental stages. Notably, 380, 242, 34, and 479 stage-specific genes were identified in the pituitary, while 886, 297, 201, and 3,678 genes were identified in the testis. Subsequent analyses unveiled associations between these stage-specific genes and crucial pathways such as the cAMP signaling pathway, GnRH secretion, and male gamete generation. Furthermore, leveraging single-cell data from the pituitary and testis, we identified some signaling pathways involving BMP, HGF, IGF, and TGF-β, highlighting mutual regulation between the pituitary and testis at different developmental stages. This study sheds light on the pivotal role of the pituitary-testis axis in the reproductive process of sheep across four distinct developmental stages. Additionally, it delves into the intricate regulatory networks governing reproduction, offering novel insights into the dynamics of the pituitary-testis axis within the reproductive system.

A novel lncRNA LOC105613571 binding with BDNF in pituitary promotes gonadotropin secretion by AKT/ERK-mTOR pathway in sheep associated with prolificacy

The pituitary is a vital endocrine organ for synthesis and secretion of gonadotropic hormones (FSH and LH), and the gonadotropin showed fluctuations in animals with different fecundity. Long non-coding RNAs (lncRNAs) have been identified as regulatory factors for the reproductive process. However, the profiles of lncRNAs and their roles involved in sheep fecundity remains unclear. In this study, we performed RNA-sequencing for the sheep pituitary gland associated with different fecundity, and identified a novel candidate lncRNA LOC105613571 targeting BDNF related to gonadotropin secretion. Our results showed that expression of lncRNA LOC105613571 and BDNF could be significantly upregulated by GnRH stimulation in sheep pituitary cells in vitro. Notably, either lncRNA LOC105613571 or BDNF silencing inhibited cell proliferation while promoted cell apoptosis. Moreover, lncRNA LOC105613571 knockdown could also downregulate gonadotropin secretion via inactivation AKT, ERK and mTOR pathway. In addition, co-treatment with GnRH stimulation and lncRNA LOC105613571 or BDNF knockdown showed the opposite effect on sheep pituitary cells in vitro. In summary, BDNF-binding lncRNA LOC105613571 in sheep regulates pituitary cell proliferation and gonadotropin secretion via the AKT/ERK-mTOR pathway, providing new ideas for the molecular mechanisms of pituitary functions.

Molecular characterization, expression patterns and cellular localization of BCAS2 gene in male Hezuo pig

Background

Breast carcinoma amplified sequence 2 (BCAS2) participates in pre-mRNA splicing and DNA damage response, which is implicated in spermatogenesis and meiosis initiation in mouse. Nevertheless, the physiological roles of BCAS2 in the testes of large mammals especially boars remain largely unknown.

Methods

In this study, testes were collected from Hezuo pig at three development stages including 30 days old (30 d), 120 days old (120 d), and 240 days old (240 d). BCAS2 CDS region was firstly cloned using RT-PCR method, and its molecular characteristics were identified using relevant bioinformatics software. Additionally, the expression patterns and cellular localization of BCAS2 were analyzed by quantitative real-time PCR (qRT-PCR), Western blot, immunohistochemistry and immunofluorescence.

Results

The cloning and sequence analysis indicated that the Hezuo pig BCAS2 CDS fragment encompassed 678 bp open reading frame (ORF) capable of encoding 225 amino acid residues, and possessed high identities with some other mammals. The results of qRT-PCR and Western blot displayed that BCAS2 levels both mRNA and protein were age-dependent increased (p < 0.01). Additionally, immunohistochemistry and immunofluorescence results revealed that BCAS2 protein was mainly observed in nucleus of gonocytes at 30 d testes as well as nucleus of spermatogonia and Sertoli cells at 120 and 240 d testes. Accordingly, we conclude that BCAS2 is critical for testicular development and spermatogenesis of Hezuo pig, perhaps by regulating proliferation or differentiation of gonocytes, pre-mRNA splicing of spermatogonia and functional maintenance of Sertoli cells, but specific mechanism still requires be further investigated.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

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 and characterization of unique and common lncRNAs and mRNAs in the pituitary, ovary, and uterus of Hu sheep with different prolificacy

LncRNAs are regarded as regulators in various animal reproductive physiological processes. However, the regulation of lncRNAs in the reproductive organ development of Hu sheep with different prolificacy remains unknown. Herein, numerous tissue-unique and -common differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs), and fecundity-unique DELs and DEGs were identified among different comparison groups at horizontal and vertical levels. Moreover, the tissue-unique and -common, and fecundity-unique female reproduction-associated DEGs and DELs were screened, and the interaction networks were constructed. Furthermore, MSTRG.43442.1 was mainly present in the cytoplasm of tested cells. The key genes ADAMTS1 and DCN were mainly localized in the granulosa cells, pituitary cells and/or endometrial epithelial cells of ovary, pituitary and/or uterus. Overall, this study identified large numbers of unique and common DELs and DEGs in the female reproductive organs of Hu sheep with different prolificacy and provided new insights into understanding the regulation of Hu sheep fecundity.

Evidence for a functional role of Start, a long noncoding RNA, in mouse spermatocytes

A mouse testis-specific long noncoding RNA (lncRNA), Start, is localized in the cytosol of Leydig cells and in the nucleus of pachytene spermatocytes. We previously showed that Start regulates steroidogenesis through controlling the expression of Star and Hsd3b1 genes in Leydig cells, but its function in germ cells was not known. Here we verified that a spermatocyte-specific protease gene, Prss43/Tessp-3, was downregulated in Start-knockout testes. To investigate the transcriptional regulatory activity of Start in spermatocytes, we first performed a series of reporter gene assays using a thymidine kinase promoter in spermatocyte-derived GC-2spd(ts) cells. A 5.4-kb genome sequence encompassing Start exhibited enhancer activity for this promoter, and the activity was decreased by knockdown of Start. Deletion of the Start promoter and replacement of the Start sequence abolished the enhancer activity and, consistently, the activity was detected in further experiments only when Start was actively transcribed. We then examined whether the Prss43/Tessp-3 gene could be a target of Start. A reporter gene assay demonstrated that the 5.4-kb sequence exhibited enhancer activity for a Prss43/Tessp-3 promoter in GC-2spd(ts) cells and that the activity was significantly decreased by knockdown of Start. These results suggest that Start functions in transcriptional activation of the Prss43/Tessp-3 gene in spermatocytes. Given that Start is presumed to regulate steroidogenic genes at the posttranscriptional level in Leydig cells, the function in spermatocytes is a novel role of Start. These findings provide an insight into multifunctionality of lncRNAs in the testis.

Comprehensive analysis of differentially expressed profiles of mRNA, lncRNA, and miRNA of Yili geese ovary at different egg-laying stages

Background

The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development.

Results

In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145–5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed.

Conclusions

Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08774-4.

Proteome Informatics in Tibetan Sheep (Ovis aries) Testes Suggest the Crucial Proteins Related to Development and Functionality

Testis has an indispensable function in male reproduction of domestic animals. Tibetan sheep (Ovis aries) is a locally adapted breed of sheep raised in the Qinghai-Tibet Plateau, with outsized roles in providing the livelihood for millions of residents. Nevertheless, less is known on how protein expression and their functional roles in developmental testes of such breed limit their use in breeding efforts. In this study, we obtained comprehensive protein profiles from testes of Tibetan sheep at three developmental stages (including pre-puberty, post-puberty, and adulthood) using data-independent acquisition-based proteomic strategy to quantitatively identify the differentially abundant proteins (DAPs) associated with testicular development and function and to unravel the molecular basis of spermatogenesis. A total of 6,221 proteins were differentially expressed in an age-dependent manner. The reliability of the gene expression abundance was corroborated by quantitative PCR and targeted parallel reaction monitoring. These DAPs were significantly enriched to biological processes concerning spermatid development and sperm deformation, mitosis, glycolytic process, cell-cell/extracellular matrix (ECM) junctions, cell proliferation, apoptosis, and migration and to the pathways including, developmental process and sexual reproduction-related (such as VEGF, estrogen, insulin, GnRH, Hippo, PI3K-Akt, mTOR, MAPK, and AMPK), and testicular cell events-related pathways (such as tight/gap/adherens junctions, ECM-receptor interaction, regulation of actin cytoskeleton, glycolysis, cell cycle, and meiosis). Based on these bioinformatics analysis, we constructed four protein–protein interaction network, among which the proteins are involved in mitosis, meiosis, spermiogenesis, and testicular microenvironment, respectively. Altogether, these bioinformatics-based sequencing results suggest that many protein-coding genes were expressed in a development-dependent manner in Tibetan sheep testes to contribute to the testicular cell development and their surrounding microenvironment remodeling at various stages of spermatogenesis. These findings have important implications for further understanding of the mechanisms underlying spermatogenesis in sheep and even other plateau-adapted animals.

Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus)

Background

An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear.

Results

Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG.

Conclusions

These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08679-2.

Comparative analysis of long noncoding RNA and mRNA expression provides insights into adaptation to hypoxia in Tibetan sheep

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA–mRNA and lncRNA–mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS_00306477, TCONS_00306029, TCONS_00139593, TCONS_00293272, and TCONS_00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments.

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

The selection of testis size can improve the reproductive capacity of livestock used for artificial insemination and has been considered as an important strategy for accelerating the breeding process. Although much work has been done to investigate the mechanisms of testis development in various species, there is little information available in regard to the differences in transcriptomic profiling of sheep testes at different developmental stages. In this work, we aimed to identify differentially expressed genes (DEGs) by RNA-Seq in sheep during different growth stages, including 0 month old (infant, M0), 3 months old (puberty, M3), 6 months old (sexual maturity, M6) and 12 months old (body maturity, M12). A total of 4,606 (2,381 up and 2,225 down), 7,500 (4,368 up and 3,132 down), 15 (8 up and seven down) DEGs were identified in M3_vs_M0, M6_vs_M3, and M12_vs_M6 comparison, respectively. Of which, a number of genes were continuously up-regulated and down-regulated with testicular development, including ODF3, ZPBP1, PKDREJ, MYBL1, PDGFA, IGF1, LH, INSL3, VIM, AMH, INHBA, COL1A1, COL1A2, and INHA. GO analysis illustrated that DEGs were mainly involved in testis development and spermatogenesis. KEGG analysis identified several important pathways and verified several reproduction-associated DEGs such as COL1A1, COL1A2, PDGFA, and IGF1. In addition, two gene modules highly associated with testis development and core genes with testis size were identified using weighted gene co-expression network analysis (WGCNA), including hub genes positively associated with testis size such as RANBP9, DNAH17, SPATA4, CIB4 and SPEM1, and those negatively associated with testis size such as CD81, CSK, PDGFA, VIM, and INHBA. This study comprehensively identified key genes related to sheep testicular development, which may provide potential insights for understanding male fertility and better guide in animal breeding.

Rare copy number variants in ASTN2 gene in patients with neurodevelopmental disorders

INTRODUCTION

In humans the normal development of cortical regions depends on the complex interactions between a number of proteins that promote the migrations of neuronal precursors from germinal zones and assembly into neuronal laminae. ASTN2 is one of the proteins implicated in such a complex process. Recently it has been observed that ASTN2 also regulates the surface expression of multiple synaptic proteins resulting in a modulation of synaptic activity. Several rare copy number variants (CNVs) in ASTN2 gene were identified in patients with neurodevelopmental disorders (NDDs) including autism spectrum disorders (ASD), attention deficit-hyperactivity disorders and intellectual disability.

METHODS

By using comparative genomic hybridization array technology, we analyzed the genomic profiles of five patients of three unrelated families with NDDs. Clinical diagnosis of ASD was established according to the Statistical Manual of Mental Disorders, Fifth Edition (APA 2013) criteria.

RESULTS

We identified new rare CNVs encompassing ASTN2 gene in three unrelated families with different clinical phenotypes of NDDs. In particular, we identified a deletion of about 70 Kb encompassing intron 19, a 186 Kb duplication encompassing the sequence between the 5'-end and the first intron of the gene and a 205 Kb deletion encompassing exons 6-11.

CONCLUSION

The CNVs reported here involve regions not usually disrupted in patients with NDDs with two of them affecting only the expression of the long isoforms. Further studies will be needed to analyze the impact of these CNVs on gene expression regulation and to better understand their impact on the protein function.

Comprehensive Analysis of miRNAs and Target mRNAs between Immature and Mature Testis Tissue in Chinese Red Steppes Cattle

Simple Summary

MicroRNAs are small molecules that can regulate the relative abundance of their target genes by binding to the 3′ untranslated region of the target genes at the post-transcriptional level to affect various biological processes, such as biosynthesis, fat metabolism and proliferation, apoptosis, and cell differentiation. Fertility is one of the most important economic traits in livestock production. Bulls require the continuous production of high-quality spermatozoa in abundance. The quality of semen is an exceptionally important factor affecting the fertilization rate of the dairy cow and is also associated with the increasing conception rate in the process of artificial insemination. Therefore, accurately predicting fertility potential for a semen sample from donor bull for artificial insemination is crucial for consistently high reproductive efficiency. The present study performed a genome-wide sequencing analysis of microRNAs and mRNAs between immature and mature testes of Chinese Red Steppes. These results provide novel candidate microRNAs and functional genes related to bull reproduction traits and the networks between microRNAs and target genes, which will provide a useful genetic mechanism and epigenetic information for marker-assisted selection of bulls with excellent sperm quality in the future.

Abstract

This study aims to screen potential regulators and regulate fecundity networks between microRNAs (miRNAs) and target genes. The bovine testes of immature and mature Chinese Red Steppes were performed by genome-wide analysis of mRNAs and miRNAs. Compared with testicular tissues of newborns, 6051 upregulated genes and 7104 downregulated genes in adult cattle were identified as differentially expressed genes (DEGs). The DEGs were significantly enriched in 808 GO terms (p < 0.05) including male gonad development, male genitalia development, spermatogenesis, and sperm motility. Moreover, DEGs were also significantly enriched in 105 KEGG pathways (p < 0.05), including cGMP-PKG signaling pathway and calcium signaling pathway. To explore the expression of miRNA-regulated gene expression, 896 differentially expressed target genes negatively regulated with the expression levels of 31 differentially expressed miRNAs (DERs) were predicted and analyzed, and a network-integrated analysis was constructed. Furthermore, real-time PCR was performed to verify the expression levels of DEGs and DERs. Our results identified novel candidate DEGs and DERs correlated with male reproduction and intricate regulating networks between miRNAs and genes, which will be valuable for future genetic and epigenetic studies of sperm development and maturity, as well as providing valuable insights into the molecular mechanisms of male fertility and spermatogenesis in cattle.

SMAD2 regulates testicular development and testosterone synthesis in Hu sheep

The SMAD protein family plays crucial roles in reproduction as a downstream target genes of the TGFβ signaling pathway. Many studies have focused on the expression change exploration of SMADs during testicular development and investigation of SMAD2 in hormone synthesis regulation. However, little attention has been given to determining the regulatory mechanism of SMADs in sheep testes. In the present study, we first detected SMAD mRNA expression levels in three-month-old (3 M), six-month-old (6 M), nine-month-old (9 M) and two-year-old (2Y) sheep testes. Different SMADs showed various expression patterns. In addition, the subcellular localization of SMAD2 was also analyzed, and Sertoli cells (SCs), Leydig cells (LCs) and spermatogonia presented mainly positive staining. Protein and nucleic acid sequence alignment showed that the SMAD2 gene was extremely homologous between various species. SMAD2 interference RNA was transfected into sheep LCs to examine the cell proliferation and hormone levels. The testosterone level was significantly decreased, and cell proliferation efficiency presented the same trend (P < 0.05). Moreover, SMAD2 downregulation promoted cell apoptosis (P < 0.05) and changed the cell cycle. In total, our results revealed that downregulating the expression of SMAD2 can effectively inhibit testosterone levels by affecting cell proliferation and apoptosis.

Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep

The tail fat of sheep is the most typical deposited fat, and it can be widely used in human daily life, such as diet, cosmetics, and industrial raw materials. To understand the potential regulatory mechanism of different growth stages of tail fat in Sunite sheep, we performed high-throughput RNA sequencing to characterize the long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles of the sheep tail fat at the age of 6, 18, and 30 months. A total of 223 differentially expressed genes (DEGs) and 148 differentially expressed lncRNAs were found in the tail fat of 6-, 18-, and 30-month-old sheep. Based on functional analysis, we found that fat-related DEGs were mainly expressed at 6 months of age and gradually decreased at 18 and 30 months of age. The target gene prediction analysis shows that most of the lncRNAs target more than 20 mRNAs as their transregulators. Further, we obtained several fat-related differentially expressed target genes; these target genes interact with different differentially expressed lncRNAs at various ages and play an important role in the development of tail fat. Based on the DEGs and differentially expressed lncRNAs, we established three co-expression networks for each comparison group. Finally, we concluded that the development of the sheep tail fat is more active during the early stage of growth and gradually decreases with the increase in age. The mutual regulation of lncRNAs and mRNAs may play a key role in this complex biological process.

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.

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

Simple Summary

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

Abstract

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

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.

Transcriptomic analysis reveals the molecular mechanisms of rumen wall morphological and functional development induced by different solid diet introduction in a lamb model

Background

This study aimed to elucidate the molecular mechanisms of solid diet introduction initiating the cellular growth and maturation of rumen tissues and characterize the shared and unique biological processes upon different solid diet regimes.

Methods

Twenty-four Hu lambs were randomly allocated to three groups fed following diets: goat milk powder only (M, n = 8), goat milk powder + alfalfa hay (MH, n = 8), and goat milk powder + concentrate starter (MC, n = 8). At 42 days of age, the lambs were slaughtered. Ruminal fluid sample was collected for analysis of concentration of volatile fatty acid (VFA) and microbial crude protein (MCP). The sample of the rumen wall from the ventral sac was collected for analysis of rumen papilla morphology and transcriptomics.

Results

Compared with the M group, MH and MC group had a higher concentration of VFA, MCP, rumen weight, and rumen papilla area. The transcriptomic results of rumen wall showed that there were 312 shared differentially expressed genes (DEGs) between in “MH vs. M” and “MC vs. M”, and 232 or 796 unique DEGs observed in “MH vs. M” or “MC vs. M”, respectively. The shared DEGs were most enriched in VFA absorption and metabolism, such as peroxisome proliferator-activated receptor (PPAR) signaling pathway, butanoate metabolism, and synthesis and degradation of ketone bodies. Additionally, a weighted gene co-expression network analysis identified M16 (2,052 genes) and M18 (579 genes) modules were positively correlated with VFA and rumen wall morphology. The M16 module was mainly related to metabolism pathway, while the M18 module was mainly associated with signaling transport. Moreover, hay specifically depressed expression of genes involved in cytokine production, immune response, and immunocyte activation, and concentrate starter mainly altered nutrient transport and metabolism, especially ion transport, amino acid, and fatty acid metabolism.

Conclusions

The energy production during VFA metabolism may drive the rumen wall development directly. The hay introduction facilitated establishment of immune function, while the concentrate starter enhanced nutrient transport and metabolism, which are important biological processes required for rumen development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00556-4.

A Testis-Specific Long Noncoding RNA, Start, Is a Regulator of Steroidogenesis in Mouse Leydig Cells

The testis expresses many long noncoding RNAs (lncRNAs), but their functions and overview of lncRNA variety are not well understood. The mouse Prss/Tessp locus contains six serine protease genes and two lncRNAs that have been suggested to play important roles in spermatogenesis. Here, we found a novel testis-specific lncRNA, Start (Steroidogenesis activating lncRNA in testis), in this locus. Start is 1822 nucleotides in length and was found to be localized mostly in the cytosol of germ cells and Leydig cells, although nuclear localization was also observed. Start-knockout (KO) mice generated by the CRISPR/Cas9 system were fertile and showed no morphological abnormality in adults. However, in adult Start-KO testes, RNA-seq and qRT-PCR analyses revealed an increase in the expression of steroidogenic genes such as Star and Hsd3b1, while ELISA analysis revealed that the testosterone levels in serum and testis were significantly low. Interestingly, at 8 days postpartum, both steroidogenic gene expression and testosterone level were decreased in Start-KO mice. Since overexpression of Start in two Leydig-derived cell lines resulted in elevation of the expression of steroidogenic genes including Star and Hsd3b1, Start is likely to be involved in their upregulation. The increase in expression of steroidogenic genes in adult Start-KO testes might be caused by a secondary effect via the androgen receptor autocrine pathway or the hypothalamus-pituitary-gonadal axis. Additionally, we observed a reduced number of Leydig cells at 8 days postpartum. Collectively, our results strongly suggest that Start is a regulator of steroidogenesis in Leydig cells. The current study provides an insight into the overall picture of the function of testis lncRNAs.

Characterization of sheep spermatogenesis through single-cell RNA sequencing

Spermatogenesis is an important biological process in male reproduction. The interaction between male germ cells and somatic cells during spermatogenesis, is necessary for male reproductive activities. This cellular heterogeneity has made it difficult to profile distinct cell types at different stages of development. Here, we present the first comprehensive, unbiased single-cell transcriptomic study of sheep spermatogenesis using 10× genomics single cell sequencing (scRNA-seq). We collected scRNA-seq data from 11 772 cells from the adult sheep testis and identified all known germ cells (including early primary spermatocytes, late primary spermatocytes, round spermatids, elongated spermatids, and sperm), and somatic cells (Sertoli cells and Leydig cells), as well as one somatic cell that unexpectedly contained leukocytes. The functional enrichment analysis indicated that several pathways of cell cycle, gamete generation, protein processing, and mRNA surveillance pathways were significantly enriched in testicular germ cell types, and ribosome pathway was significantly enriched in testicular somatic cell types. Further analysis identified several stage-specific marker genes of sheep germ cells, such as EZH2, SOX18, SCP2, PCNA, and PRKCD. Our research explored for the first time of the changes in the transcription level of various cell types during the process of sheep spermatogenesis, providing new insights for sheep spermatogenesis and spermatogenic cell development.

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.

Long non-coding RNA366.2 controls endometrial epithelial cell proliferation and migration by upregulating WNT6 as a ceRNA of miR-1576 in sheep uterus

Long non-coding RNAs (lncRNAs) play an important regulatory role in mammalian fecundity. Currently, most studies are primarily concentrated on ovarian lncRNAs, ignoring the influence of uterine lncRNAs on the fecundity of female sheep. In this study, we found a higher density of uterine glands and endometrial microvessel density (MVD) in high prolificacy group of Hu sheep compared to low prolificacy groups (p < 0.05) as well as an increased level of serum placental growth factor (PLGF). Hundreds of differentially expressed (DE) lncRNAs were identified in Hu sheep with different fecundity by RNA sequencing (RNA-seq), and their targets were enriched in some signaling pathways involved in endometrial functions, such as the estrogen signaling pathway, nuclear factor kappa B (NF-κB) signaling pathway, oxytocin signaling pathway, and Wnt signaling pathway. Furthermore, the underlying mechanisms of competitive endogenous RNA (ceRNA) of lncRNA366.2-miR-1576- WNT6 were determined by bioinformatics analysis. Functionally, our results indicated that lncRNA366.2 promoted endometrial epithelial cell (EEC) proliferation, migration, and growth factor expression by sponging miR-1576 to upregulate WNT6 expression and activate the Wnt/β-catenin pathway. Taken together, our research indicated the regulatory mechanism of the lncRNA366.2-miR-1576-WNT6 in EEC proliferation and migration. Furthermore, this study provides a new theoretical reference for the identification of candidate genes related to fecundity.

Genome-wide identification and comparison of differentially expressed profiles of miRNAs and lncRNAs with associated ceRNA networks in the gonads of Chinese soft-shelled turtle, Pelodiscus sinensis

Background

The gonad is the major factor affecting animal reproduction. The regulatory mechanism of the expression of protein-coding genes involved in reproduction still remains to be elucidated. Increasing evidence has shown that ncRNAs play key regulatory roles in gene expression in many life processes. The roles of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in reproduction have been investigated in some species. However, the regulatory patterns of miRNA and lncRNA in the sex biased expression of protein coding genes remains to be elucidated. In this study, we performed an integrated analysis of miRNA, messenger RNA (mRNA), and lncRNA expression profiles to explore their regulatory patterns in the female ovary and male testis of Pelodiscus sinensis.

Results

We identified 10,446 mature miRNAs, 20,414 mRNAs and 28,500 lncRNAs in the ovaries and testes, and 633 miRNAs, 11,319 mRNAs, and 10,495 lncRNAs showed differential expression. A total of 2814 target genes were identified for miRNAs. The predicted target genes of these differentially expressed (DE) miRNAs and lncRNAs included abundant genes related to reproductive regulation. Furthermore, we found that 189 DEmiRNAs and 5408 DElncRNAs showed sex-specific expression. Of these, 3 DEmiRNAs and 917 DElncRNAs were testis-specific, and 186 DEmiRNAs and 4491 DElncRNAs were ovary-specific. We further constructed complete endogenous lncRNA-miRNA-mRNA networks using bioinformatics, including 103 DEmiRNAs, 636 DEmRNAs, and 1622 DElncRNAs. The target genes for the differentially expressed miRNAs and lncRNAs included abundant genes involved in gonadal development, including Wt1, Creb3l2, Gata4, Wnt2, Nr5a1, Hsd17, Igf2r, H2afz, Lin52, Trim71, Zar1, and Jazf1.

Conclusions

In animals, miRNA and lncRNA as master regulators regulate reproductive processes by controlling the expression of mRNAs. Considering their importance, the identified miRNAs, lncRNAs, and their targets in P. sinensis might be useful for studying the molecular processes involved in sexual reproduction and genome editing to produce higher quality aquaculture animals. A thorough understanding of ncRNA-based cellular regulatory networks will aid in the improvement of P. sinensis reproductive traits for aquaculture.

A novel testis-specific long noncoding RNA, Tesra, activates the Prss42/Tessp-2 gene during mouse spermatogenesis†

The progression of spermatogenesis is precisely controlled by meiotic stage-specific genes, but the molecular mechanism for activation of such genes is still elusive. Here we found a novel testis-specific long noncoding RNA (lncRNA), Tesra, that was specifically expressed in the mouse testis at the Prss/Tessp gene cluster on chromosome 9. Tesra was transcribed downstream of Prss44/Tessp-4, starting within the gene, as a 4435-nucleotide transcript and developmentally activated at a stage similar to that for Prss/Tessp genes. By in situ hybridization, Tesra was found to be localized in and around germ cells and Leydig cells, being consistent with biochemical data showing its existence in cytoplasmic, nuclear, and extracellular fractions. Based on the finding of more signals in nuclei of pachytene spermatocytes, we explored the possibility that Tesra plays a role in transcriptional activation of Prss/Tessp genes. By a ChIRP assay, the Tesra transcript was found to bind to the Prss42/Tessp-2 promoter region in testicular germ cells, and transient overexpression of Tesra significantly activated endogenous Prss42/Tessp-2 expression and increased Prss42/Tessp-2 promoter activity in a reporter construct. These findings suggest that Tesra activates the Prss42/Tessp-2 gene by binding to the promoter. Finally, we investigated whether Tesra co-functioned with enhancers adjacent to another lncRNA, lncRNA-HSVIII. In the Tet-on system, Tesra transcription significantly increased activity of one enhancer, but Tesra and the enhancer were not interdependent. Collectively, our results proposed a potential function of an lncRNA, Tesra, in transcriptional activation and suggest a novel relationship between an lncRNA and an enhancer.

Exploring the roles of fecundity-related long non-coding RNAs and mRNAs in the adrenal glands of small-tailed Han Sheep

BACKGROUND

Long non-coding RNAs (lncRNAs) can play important roles in uterine and ovarian functions. However, little researches have been done on the role of lncRNAs in the adrenal gland of sheep. Herein, RNA sequencing was used to compare and analyze gene expressions in adrenal tissues between follicular phases and luteal phases in FecBBB (MM) and FecB++ (WW) sheep, respectively, and differentially expressed lncRNAs and genes associated with reproduction were identified.

RESULTS

In MM sheep, 38 lncRNAs and 545 mRNAs were differentially expressed in the adrenal gland between the luteal and follicular phases; In WW sheep, 513 differentially expressed lncRNAs and 2481 mRNAs were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that differentially expressed lncRNAs and their target genes are mainly involved in the circadian rhythm, the mitogen activated protein kinase, thyroid, ovarian steroidogenesis and transforming growth factor beta signaling pathways. Differentially expressed lncRNAs can regulate reproduction by modulating genes involved in these signaling pathways and biological processes. Specifically, XLOC_254761, XLOC_357966, 105,614,839 and XLOC_212877 targeting CREB1, PER3, SMAD1 and TGFBR2, respectively, appear to play key regulatory roles.

CONCLUSION

These results broaden our understanding of lncRNAs in adrenal gland of sheep and provide new insights into the molecular mechanisms underlying sheep reproduction.

Genome-Wide Analysis and Function Prediction of Long Noncoding RNAs in Sheep Pituitary Gland Associated with Sexual Maturation

Long noncoding RNA (lncRNA) plays a crucial role in the hypothalamic-pituitary-testis (HPT) axis associated with sheep reproduction. The pituitary plays a connecting role in the HPT axis. However, little is known of their expression pattern and potential roles in the pituitary gland. To explore the potential lncRNAs that regulate the male sheep pituitary development and sexual maturation, we constructed immature and mature sheep pituitary cDNA libraries (three-month-old, TM, and nine-month-old, NM, respectively, n = 3) for lncRNA and mRNA high-throughput sequencing. Firstly, the expression of lncRNA and mRNA were comparatively analyzed. 2417 known lncRNAs and 1256 new lncRNAs were identified. Then, 193 differentially expressed (DE) lncRNAs and 1407 DE mRNAs were found in the pituitary between the two groups. Moreover, mRNA-lncRNA interaction network was constructed according to the target gene prediction of lncRNA and functional enrichment analysis. Five candidate lncRNAs and their targeted genes HSD17B12, DCBLD2, PDPK1, GPX3 and DLL1 that enriched in growth and reproduction related pathways were further filtered. Lastly, the interaction of candidate lncRNA TCONS_00066406 and its targeted gene HSD17B12 were validated in in vitro of sheep pituitary cells. Our study provided a systematic presentation of lncRNAs and mRNAs in male sheep pituitary, which revealed the potential role of lncRNA in male reproduction.

Comprehensive Analysis of Differentially Expressed Profiles of mRNA, lncRNA, and circRNA in the Uterus of Seasonal Reproduction Sheep

Photoperiod is one of the important factors leading to seasonal reproduction of sheep. However, the molecular mechanisms underlying the photoperiod regulation of seasonal reproduction remain poorly understood. In this study, we compared the expression profiles of mRNAs, lncRNAs, and circRNAs in uterine tissues from Sunite sheep during three different photoperiods, namely, the short photoperiod (SP), short transfer to long photoperiod (SLP), and long photoperiod (LP). The results showed that 298, 403, and 378 differentially expressed (DE) mRNAs, 171, 491, and 499 DE lncRNAs, and 124, 270, and 400 DE circRNAs were identified between SP and LP, between SP and SLP, and between LP and SLP, respectively. Furthermore, functional enrichment analysis showed that the differentially expressed RNAs were mainly involved in the GnRH signaling pathway, thyroid hormone synthesis, and thyroid hormone signaling pathway. In addition, co-expression networks of lncRNA–mRNA were constructed based on the correlation analysis between the differentially expressed RNAs. Our study provides new insights into the expression changes of RNAs in different photoperiods, which might contribute to understanding the molecular mechanisms of seasonal reproduction in sheep.

Gene expression patterns and protein cellular localization suggest a novel role for DAZL in developing Tibetan sheep testes

Deleted in azoospermia-like (DAZL) is essential for mammalian spermatogenesis as it regulates proliferation, development, maturation and functional maintenance of male germ cells. Its expression and regulation vary with different species or in the same animal at different developmental stages, and despite its importance, very little is known about its roles in sheep, especially Tibetan sheep. To investigate the expression patterns and regulatory roles of DZAL in Tibetan sheep testis, testicular tissue was isolated from sheep at three crucial development stages: 3 months old, 1 year old and 3 years old. Using quantitative real-time PCR and Western blot, we found that DAZL mRNA first decreased and then increased with advancing age, while DAZL protein exhibited an opposite expression pattern, with first increased and subsequently decreased levels. Immunohistochemistry and immunofluorescence revealed that DAZL protein was located predominantly in the cytoplasm of Leydig cells and in both the cytoplasm and nucleus of spermatids. ELISA indicated that testosterone content within developing testes was first enhanced and then declined. Our results, taken together, demonstrate, for the first time, that DAZL gene is involved in Tibetan sheep spermatogenesis by regulating the development of spermatids in post-pubertal rams, along with a novel role in functional maintenance of Leydig cells in postnatal rams.

Functions and mechanism of noncoding RNA in the somatic cells of the testis

ncRNAs are involved in numerous biological processes by regulating gene expression and cell stability. Studies have shown that ncRNAs also contribute to spermatogenesis. Leydig cells (LCs) and Sertoli cells (SCs) are somatic cells of the testis that support spermatogenesis and are vital to male fertility. In this review, we summarized the findings from studies on ncRNAs in SCs and LCs. In SCs, ncRNAs play key roles in phagocytosis, immunoprotection and development of SCs. In LCs, ncRNAs are involved in steroidogenesis, in particular production of testosterone as well as development of LCs. Here, we discuss the possible target genes and functions of ncRNAs in both types of cells. These ncRNAs regulate the expression of target genes or mRNA coding sequence regions, resulting in a chain reaction that influences cell function. In addition, microRNAs, lncRNAs, piRNA-like RNAs (pilRNAs) and natural antisense transcripts (NATs) are discussed in this review. In summary, we suggest that these ncRNAs might act in coordination to control spermatogenesis and maintain the environmental homeostasis of the testis.

LncNONO-AS regulates AR expression by mediating NONO

Spermatogenesis and healthy testicular development are prerequisites for male reproductive function. Androgen receptor (AR), an important receptor in testicular sertoli cells, is involved in androgen specific response and its dysfunction will lead to abnormal sperm development, resulting in male infertility. NONO (non-POU-domain-containing octamer binding protein) can act as a coactivator to enhance the transcription of AR, while AR may be regulated by NONO in testicular sertoli cells. LncRNAs are involved in almost every step of spermatogenesis. However, there are few studies focus on the relationship between lncRNAs and spermatogenesis in goat testis. Therefore, in this research, high throughput sequencing and bioinformatics analysis were performed on testicular tissues of Dazu black goats at different stages of development to obtain the target NONO lncRNA. It's called lncNONO-AS. This study further explored the biological functions of lncRNA through RNA pull down, overexpression, interference, fluorescence quantification, Western blot and other techniques on the basis of in vitro culture of testis sertoli cells, and we got the following results: The gene expression levels of NONO and AR in lncNONO-AS overexpression group were significantly higher than that in the empty vector group (P < 0.01). Compared with the untreated negative control group, the expression of NONO decreased from 1.00 to 0.68 (P < 0.01), and the expression of AR decreased from 1.01 to 0.34 (P < 0.01). The results showed that lncNONO-AS could regulate the expression of AR by mediating the expression of NONO.

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.

Digital gene expression analyses of mammary glands from meat ewes naturally infected with clinical mastitis

Clinical mastitis in sheep has gravely restrained production performance for a long time. Knowledge of mechanisms of its pathogenesis and resistance in meat sheep mammary gland with clinical mastitis are not yet understood, especially for clinical mastitis caused by natural infection. In this work, RNA-sequencing was firstly used to screen the differentially expressed genes (DEGs) in clinical mastitic mammary tissues (CMMTs) when compared with healthy mammary tissues (HMTs) from meat sheep flocks. We identified 420 DEGs including 316 upregulated and 104 downregulated genes in CMMTs. Gene ontology annotation revealed these DEGs were mainly engaged in immune response and inflammation response. Pathway enrichment showed they were primarily enriched in pathways relevant to inflammation, immune response and metabolism. Alternative splicing analysis showed most common differential splicing genes in CMMTs and HMTs were implicated in immune response. Immunostaining for three immune response-related proteins encoded by DEGs were mainly observed in mammary epithelium from both CMMTs and HMTs, and their positive signals were more intensive in CMMTs than those in HMTs. These findings provide experimental basis and reference for further researching the molecular genetic mechanisms, particularly immune defence mechanisms, of sheep mammary gland during clinical mastitis.

Identification and characterization of mRNAs and lncRNAs in the uterus of polytocous and monotocous Small Tail Han sheep (Ovis aries)

Background

Long non-coding RNAs (lncRNAs) regulate endometrial secretion and uterine volume. However, there is little research on the role of lncRNAs in the uterus of Small Tail Han sheep (FecB++). Herein, RNA-seq was used to comparatively analyze gene expression profiles of uterine tissue between polytocous and monotocous sheep (FecB++) in follicular and luteal phases.

Methods

To identify lncRNA and mRNA expressed in the uterus, the expression of lncRNA and mRNA in the uterus of Small Tail Han sheep (FecB++) from the polytocous group (n = 6) and the monotocous group (n = 6) using RNA-sequencing and real-time polymerase chain reaction (RT-PCR). Identification of differentially expressed lncRNAs and mRNAs were performed between the two groups and two phases . Gene ontology (GO) and pathway enrichment analyses were performed to analyze the biological functions and pathways for the differentially expressed mRNAs. LncRNA-mRNA co-expression network was constructed to further analyses the function of related genes.

Results

In the follicular phase, 473 lncRNAs and 166 mRNAs were differentially expressed in polytocous and monotocous sheep; in the luteal phase, 967 lncRNAs and 505 mRNAs were differentially expressed in polytocous and monotocous sheep. GO and KEGG enrichment analysis showed that the differentially expressed lncRNAs and their target genes are mainly involved in ovarian steroidogenesis, retinol metabolism, the oxytocin signaling pathway, steroid hormone biosynthesis, and the Foxo signaling pathway. Key lncRNAs may regulate reproduction by regulating genes involved in these signaling pathways and biological processes. Specifically, UGT1A1, LHB, TGFB1, TAB1, and RHOA, which are targeted by MSTRG.134747, MSTRG.82376, MSTRG.134749, MSTRG.134751, and MSTRG.134746, may play key regulatory roles. These results offer insight into molecular mechanisms underlying sheep prolificacy.