Long noncoding RNAs in spermatogenesis: insights from recent high-throughput transcriptome studies

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.

Oxidative Stress Biomarkers in Male Infertility: Established Methodologies and Future Perspectives

Male fertility can be affected by oxidative stress (OS), which occurs when an imbalance between the production of reactive oxygen species (ROS) and the body's ability to neutralize them arises. OS can damage cells and influence sperm production. High levels of lipid peroxidation have been linked to reduced sperm motility and decreased fertilization ability. This literature review discusses the most commonly used biomarkers to measure sperm damage caused by ROS, such as the high level of OS in seminal plasma as an indicator of imbalance in antioxidant activity. The investigated biomarkers include 8-hydroxy-2-deoxyguanosine acid (8-OHdG), a marker of DNA damage caused by ROS, and F2 isoprostanoids (8-isoprostanes) produced by lipid peroxidation. Furthermore, this review focuses on recent methodologies including the NGS polymorphisms and differentially expressed gene (DEG) analysis, as well as the epigenetic mechanisms linked to ROS during spermatogenesis along with new methodologies developed to evaluate OS biomarkers. Finally, this review addresses a valuable insight into the mechanisms of male infertility provided by these advances and how they have led to new treatment possibilities. Overall, the use of biomarkers to evaluate OS in male infertility has supplied innovative diagnostic and therapeutic approaches, enhancing our understanding of male infertility mechanisms.

Impacts of non-coding RNAs in the pathogenesis of varicocele

Two classes of non-coding RNAs, namely lncRNAs and miRNAs have been reported to be involved in the pathogenesis of varicocele. MIR210HG, MLLT4-AS1, gadd7, and SLC7A11-AS1 are among lncRNAs whose expression has been changed in patients with varicocele in association with the sperm quality. Animal studies have also suggested contribution of NONRATG001060, NONRATG002949, NONRATG013271, NONRATG027523 and NONRATG023747 lncRNAs in this pathology. Meanwhile, expression of some miRNAs, such as miR-210-3p, miR-21, miR-34a, miR-122a, miR-181a, miR-34c and miR-192a has been altered in this condition. Some of these transcripts have the potential to predict the sperm quality. We summarize the impacts of lncRNAs and miRNAs in the pathogenesis of varicocele.

Functions and mechanism of noncoding RNA in regulation and differentiation of male mammalian reproduction

Noncoding RNAs (ncRNAs) are active regulators of a wide range of biological and physiological processes, including the majority of mammalian reproductive events. Knowledge of the biological activities of ncRNAs in the context of mammalian reproduction will allow for a more comprehensive and comparative understanding of male sterility and fertility. In this review, we describe recent advances in ncRNA-mediated control of mammalian reproduction and emphasize the importance of ncRNAs in several aspects of mammalian reproduction, such as germ cell biogenesis and reproductive organ activity. Furthermore, we focus on gene expression regulatory feedback loops including hormones and ncRNA expression to better understand germ cell commitment and reproductive organ function. Finally, this study shows the role of ncRNAs in male reproductive failure and provides suggestions for further research.

LncRNA Functional Screening in Organismal Development

Controversy continues over the functional prevalence of long non-coding RNAs (lncRNAs) despite their being widely investigated in all kinds of cells and organisms. In animals, lncRNAs have aroused general interest from exponentially increasing transcriptomic repertoires reporting their highly tissue-specific and developmentally dynamic expression, and more importantly, from growing experimental evidence supporting their functionality in facilitating organogenesis and individual fitness. In mammalian testes, while a great multitude of lncRNA species are identified, only a minority of them have been shown to be useful, and even fewer have been demonstrated as true requirements for male fertility using knockout models to date. This noticeable gap is attributed to the virtual existence of a large number of junk lncRNAs, the lack of an ideal germline culture system, difficulty in loss-of-function interrogation, and limited screening strategies. Facing these challenges, in this review, we discuss lncRNA functionality in organismal development and especially in mouse testis, with a focus on lncRNAs with functional screening.

Seminal plasma biomarkers for predicting successful sperm retrieval in patients with nonobstructive azoospermia: a narrative review of human studies

BACKGROUND

Non-obstructive azoospermia (NOA) is considered to be the most severe form of male infertility. Before the emergence of surgical testicular sperm extraction and assisted reproductive technology, NOA patients could hardly become biological fathers of their children. However, failure of the surgery could cause physical and psychological harm to patients such as testicular damage, pain, hopeless of fertility and additional cost. Therefore, predicting the successful sperm retrieval (SSR) is so important for NOA patients to make their choice whether to do the surgery or not. Because seminal plasma is secreted by the testes and accessory gonads, it can reflect the spermatogenic environment, making it a preferential choice for SSR valuation. The purpose of this paper is to summarize the available evidence and provide the reader with a broad overview of biomarkers in seminal plasma for SSR prediction.

RESULTS

A total of 15,390 studies were searched from PUBMED, EMBASE, CENTRAL and Web of Science, but only 6615 studies were evaluated after duplications were removed. The abstracts of 6513 articles were excluded because they were irrelevant to the topic. The full texts of 102 articles were obtained, with 21 of them being included in this review. The included studies range in quality from medium to high. In the included articles, surgical sperm extraction methods included conventional testicular sperm extraction (TESE) and microdissection testicular sperm extraction (micro-TESE). Currently, the biomarkers in seminal plasma used to predict SSR are primarily RNAs, metabolites, AMH, inhibin B, leptin, survivin, clusterin, LGALS3BP, ESX1, TEX101, TNP1, DAZ, PRM1 and PRM2.

CONCLUSION

The evidence does not conclusively indicate that AMH and INHB in seminal plasma are valuable to predict the SSR. It is worth noting that RNAs, metabolites and other biomarkers in seminal plasma have shown great potential in predicting SSR. However, existing evidence is insufficient to provide clinicians with adequate decision support, and more prospective, large sample size, and multicenter trials are urgently needed.

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.

Whole-genome sequencing identifies new candidate genes for nonobstructive azoospermia

BACKGROUND

Genetic causes that lead to spermatogenetic failure in patients with nonobstructive azoospermia (NOA) have not been yet completely established.

OBJECTIVE

To identify low-frequency NOA-associated single nucleotide variants (SNVs) using whole-genome sequencing (WGS).

MATERIALS AND METHODS

Men with various types of NOA (n = 39), including samples that had been previously tested with whole-exome sequencing (WES; n = 6) and did not result in diagnostic conclusions. Variants were annotated using the Ensembl Variant Effect Predictor, utilizing frequencies from GnomAD and other databases to provide clinically relevant information (ClinVar), conservation scores (phyloP), and effect predictions (i.e., MutationTaster). Structural protein modeling was also performed.

RESULTS

Using WGS, we revealed potential NOA-associated SNVs, such as: TKTL1, IGSF1, ZFPM2, VCX3A (novel disease causing variants), ESX1, TEX13A, TEX14, DNAH1, FANCM, QRICH2, FSIP2, USP9Y, PMFBP1, MEI1, PIWIL1, WDR66, ZFX, KCND1, KIAA1210, DHRSX, ZMYM3, FAM47C, FANCB, FAM50B (genes previously known to be associated with infertility) and ALG13, BEND2, BRWD3, DDX53, TAF4, FAM47B, FAM9B, FAM9C, MAGEB6, MAP3K15, RBMXL3, SSX3 and FMR1NB genes, which may be involved in spermatogenesis.

DISCUSSION AND CONCLUSION

In this study, we identified novel potential candidate NOA-associated genes in 29 individuals out of 39 azoospermic males. Note that in 5 out of 6 patients subjected previously to WES analysis, which did not disclose potentially causative variants, the WGS analysis was successful with NOA-associated gene findings.

DNAH5 gene and its correlation with linc02220 expression and sperm characteristics

BACKGROUND

Numerous pieces of evidence show that many environmental and genetic factors can cause male infertility. Much research in recent years has investigated the function of long non-coding RNAs (lncRNAs) in fertility. The main objective of the current study was to investigate the expression of Dynein Axonemal Heavy Chain 5 (DNAH5) as a gene that plays an essential role in sperm motility in individuals with asthenozoospermia and terato-asthenozoospermia. Alterations in linc02220 expression (located close to the DNAH5 gene), its action potential in DNAH5 regulating, and the correlation between their expression and normal sperm morphology and motility were also examined.

METHOD AND MATERIAL

This study examined the semen of 31 asthenozoospermia individuals (AZ), 33 terato-asthenozoospermia (TAZ) individuals, and 33 normospermia (NZ) individuals with normal sperm as a control group. The expression levels of DNAH5 and linc02220 in the sperm samples were analyzed by real-time PCR.

RESULTS

Gene expression analysis revealed a significant association between DNAH5 expression and sperm motility and morphology (p < 0.0001). The DNAH5 expression levels in the TAZ and AZ groups were also significantly reduced; however, linc02220 was significantly upregulated in both TAZ and AZ groups compared to the NZ group (p < 0.0001). DNAH5 expression in the TAZ and AZ groups was negatively correlated with linc02220 expression, thus, DNAH5 downregulation was associated with linc02220 overexpression (p < 0.05).

CONCLUSIONS

The gene linc02220 could be a potential regulatory target for DNAH5, and both could affect sperm's normal motility and morphology.

Developmental alterations in DNA methylation during gametogenesis from primordial germ cells to sperm

Because epigenetics is a critical component for gene expression, the hypothesis was tested that DNA methylation alterations are dynamic and continually change throughout gametogenesis to generate the mature sperm. Developmental alterations and stage-specific DNA methylation during gametogenesis from primordial germ cells (PGCs) to mature sperm are investigated. Individual developmental stage germ cells were isolated and analyzed for differential DNA methylation regions (DMRs). The number of DMRs was highest in the first three comparisons with mature PGCs, prospermatogonia, and spermatogonia. The most statistically significant DMRs were present at all stages of development and had variations involving both increases or decreases in DNA methylation. DMR-associated genes were identified and correlated with gene functional categories, pathways, and cellular processes. Observations identified a dynamic cascade of epigenetic changes during development that is dramatic during the early developmental stages. Complex epigenetic alterations are required to regulate genome biology and gene expression during gametogenesis.

•

A dynamic cascade of epigenetic change throughout gametogenesis from PGC to sperm

•

Most dramatic epigenetic alterations in PGC and spermatogenic stem cell stages

•

Different DNA methylation regions between and within stages were identified

•

Complex epigenetic alterations required for gene expression during gametogenesis

Epigenetics in the Uterine Environment: How Maternal Diet and ART May Influence the Epigenome in the Offspring with Long-Term Health Consequences

The societal burden of non-communicable disease is closely linked with environmental exposures and lifestyle behaviours, including the adherence to a poor maternal diet from the earliest preimplantation period of the life course onwards. Epigenetic variations caused by a compromised maternal nutritional status can affect embryonic development. This review summarises the main epigenetic modifications in mammals, especially DNA methylation, histone modifications, and ncRNA. These epigenetic changes can compromise the health of the offspring later in life. We discuss different types of nutritional stressors in human and animal models, such as maternal undernutrition, seasonal diets, low-protein diet, high-fat diet, and synthetic folic acid supplement use, and how these nutritional exposures epigenetically affect target genes and their outcomes. In addition, we review the concept of thrifty genes during the preimplantation period, and some examples that relate to epigenetic change and diet. Finally, we discuss different examples of maternal diets, their effect on outcomes, and their relationship with assisted reproductive technology (ART), including their implications on epigenetic modifications.

A Novel Meiosis-Related lncRNA, Rbakdn, Contributes to Spermatogenesis by Stabilizing Ptbp2

Spermatocyte meiosis is the cornerstone of mammalian production. Thousands of long noncoding RNAs (lncRNAs) have been reported to be functional in various cellular processes, but the function of lncRNAs in meiosis remains largely unknown. Here, we profiled lncRNAs in spermatocytes at stage I of meiosis and identified a testis-specific lncRNA, Rbakdn, as a vital regulator of meiosis. Rbakdn is dynamically expressed during meiosis I, and Rbakdn knockdown inhibits meiosis in vitro. Furthermore, Rbakdn knockdown in testes in mice by intratesticular injection disturbs meiosis, reduces testicular volume, and increases apoptosis of spermatocytes, resulting in vacuolation of the seminiferous tubules. Rbakdn can bind to Ptbp2, an RNA-binding protein that is important in the regulation of the alternative splicing of many genes in spermatogenesis. Rbakdn knockdown leads to a decrease in Ptbp2 through the ubiquitination degradation pathway, indicating that Rbakdn maintains the stability of Ptbp2. In conclusion, our study identified an lncRNA, Rbakdn, with a crucial role in meiosis.

Testis transcriptome profiling identified lncRNAs involved in spermatogenic arrest of cattleyak

Cattleyaks are the crossbred offspring between cattle and yaks, exhibiting the prominent adaptability to the harsh environment as yaks and much higher growth performances than yaks around Qinghai-Tibet plateau. Unfortunately, cattleyak cannot be effectively used in yak breeding due to its male infertility resulted from spermatogenic arrest. In this study, we performed RNA sequencing (RNA-seq) and bioinformatics analysis to determine the expression profiles of long noncoding RNA (lncRNA) from cattleyak and yak testis. A total of 604 differentially expressed (DE) lncRNAs (135 upregulated and 469 downregulated) were identified in cattleyak with respect to yak. Through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, we identified several DE lncRNAs regulating the mitotic cell cycle processes by targeting the genes significantly associated with the mitotic cell cycle checkpoint and DNA damage checkpoint term and also significantly involved in p53 signaling pathway, mismatch repair and homologous recombination pathway (P < 0.05). The reverse transcription PCR (RT-PCR) and quantitative Real-Time PCR (qRT-PCR) analysis of the randomly selected fourteen DE lncRNAs and the seven target genes validated the RNA-seq data and their true expressions during spermatogenesis in vivo. Molecular cloning and sequencing indicated that the testis lncRNAs NONBTAT012170 and NONBTAT010258 presented higher similarity among different cattleyak and yak individuals. The downregulation of these target genes in cattleyak contributed to the abnormal DNA replication and spermatogenic arrest during the S phase of mitotic cell cycle. This study provided a novel insight into lncRNA expression profile changes associated with spermatogenic arrest of cattleyak.

Long non-coding RNA profile in banana shrimp, Fenneropenaeus merguiensis and the potential role of lncPV13 in vitellogenesis

The long non-coding RNAs (lncRNAs) have been known to play important roles in several biological processes as well as in reproduction. This study aimed to identify lncRNA in ovary female banana shrimp, Fenneropenaeus merguiensis, and investigate the potential role of lncPV13 in the vitellogenesis. After the in silico identification of the ovarian transcriptome, a total of 24,733 putative lncRNAs were obtained, and only 147 putative lncRNAs were significantly differentially expressed among the ovarian development stages. To validate the in silico identification of lncRNAs, the 16 lncRNAs with the highest differential expression in the transcriptome analysis were evaluated by RT-qPCR. The 6 lncRNAs showed higher expression levels in the mature stage than in the previtellogenic stage and were found in several tissues such as in eyestalks, brains, thoracic ganglia, gills, and muscle. Furthermore, most candidate lncRNAs were amplifiable in Litopenaeus vannamei's and Penaeus monodon's DNA but not in Macrobrachium rosenbergii's DNA, suggesting some lncRNAs are expressed in a species-specific manner among penaeid shrimp. In this study, the lncPV13 was investigated for its vitellogenin regulating function by RNA interference. The result indicates that the lncPV13 expression was suppressed in the ovary on day 7 after the injection of double-stranded RNA specific to lncPV13 (dslncPV13), while vitellogenin (Vg) expression was significantly decreased. In contrast, the gonad inhibiting hormone (GIH) expression was significantly increased in the lncPV13 knockdown shrimp. However, the oocyte proliferation was not significantly different between control and lncPV13 knockdown shrimp. This suggests that lncPV13 regulate Vg synthesis through GIH inhibition. Finally, our findings provide lncRNA information and potential lncRNAs involved in the vitellogenesis of female banana shrimp.

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.

LncRNAs induce oxidative stress and spermatogenesis by regulating endoplasmic reticulum genes and pathways

Oligozoospermia or low sperm count is a leading cause of male infertility worldwide. Despite decades of work on non-coding RNAs (ncRNAs) as regulators of spermatogenesis, fertilization, and male fertility, the literature on the function of long non-coding RNAs (lncRNAs) in human oligozoospermia is scarce. We integrated lncRNA and mRNA sequencing data from 12 human normozoospermic and oligozoospermic samples and comprehensively analyzed the function of differentially expressed lncRNAs (DE lncRNAs) and mRNAs (DE mRNAs) in male infertility. The target genes of DE lncRNAs were identified using a Gaussian graphical model. Gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were primarily enriched in protein transport and localization to the endoplasmic reticulum (ER). The lncRNA–mRNA co-expression network revealed cis- and trans-regulated target genes of lncRNAs. The transcriptome data implicated DE lncRNAs and DE mRNAs and their target genes in the accumulation of unfolded proteins in sperm ER, PERK-EIF2 pathway-induced ER stress, oxidative stress, and sperm cell apoptosis in individuals with oligozoospermia. These findings suggest that the identified lncRNAs and pathways could serve as effective therapeutic targets for male infertility.

A panel of extracellular vesicle long noncoding RNAs in seminal plasma for predicting testicular spermatozoa in nonobstructive azoospermia patients

STUDY QUESTION

Whether the testis-specific extracellular vesicle (EV) long noncoding RNAs (lncRNAs) in seminal plasma could be utilized to predict the presence of testicular spermatozoa in nonobstructive azoospermia (NOA) patients?

SUMMARY ANSWER

Our findings indicate that the panel based on seminal plasma EV lncRNAs was a sensitive and specific method in predicting the presence of testicular spermatozoa and may improve clinical decision-making of NOA.

WHAT IS KNOWN ALREADY

The adoption of sperm retrieval techniques, especially microdissection testicular sperm extraction (mTESE), in combination with ICSI has revolutionized treatment for NOA. However, there are no precise and noninvasive methods for predicting whether there are testicular spermatozoa in NOA patients before mTESE.

STUDY DESIGN, SIZE, DURATION

RNA sequencing was performed on seminal plasma EVs from 6 normozoospermic men who underwent IVF due to female factor and 5 idiopathic NOA patients who failed to obtain testicular spermatozoa by mTESE and were diagnosed as having Sertoli cell-only syndrome by postoperative pathology. A biomarker panel of lncRNAs was constructed and verified in 96 NOA patients who underwent mTESE. Decision-making process was established based on the panel in seminal plasma EVs from 45 normozoospermia samples, 43 oligozoospermia samples, 62 cryptozoospermia samples, 96 NOA samples.

PARTICIPANTS/MATERIALS, SETTING, METHODS

RNA sequencing was done to examine altered profiles of EV lncRNAs in seminal plasma. Furthermore, a panel consisting of EV lncRNAs was established and evaluated in training set and validation sets.

MAIN RESULTS AND THE ROLE OF CHANCE

A panel consisting of nine differentially expressed testis-specific lncRNAs, including LOC100505685, SPATA42, CCDC37-DT, GABRG3-AS1, LOC440934, LOC101929088 (XR_927561.2), LOC101929088 (XR_001745218.1), LINC00343 and LINC00301, was established in the training set and the AUC was 0.986. Furthermore, the AUC in the validation set was 0.960. Importantly, the panel had a unique advantage when compared with models based on serum hormones from the same group of NOA cases (AUC, 0.970 vs 0.723; 0.959 vs 0.687, respectively). According to the panel of lncRNAs, a decision-making process was established, that is when the score of an NOA case exceeds 0.532, sperm retrieval surgery may be recommended.

LIMITATIONS, REASONS FOR CAUTION

In the future, the sample size needs to be further expanded. Meanwhile, the regulatory functions and mechanism of lncRNAs in spermatogenesis also need to be elucidated.

WIDER IMPLICATIONS OF THE FINDINGS

When the score of our panel is below 0.532, subjecting the NOA patients to ineffective surgical interventions may not be recommended due to poor sperm retrieval rate.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the National Natural Science Foundation of China (81871110, 81971314 and 81971759); the Guangdong Special Support Plan-Science and Technology Innovation Youth Top Talents Project (2016TQ03R444); the Science and Technology Planning Project of Guangdong Province (2016B030230001 and 201707010394); the Key Scientific and Technological Program of Guangzhou City (201604020189); the Pearl River S&T Nova Program of Guangzhou (201806010089); the Transformation of Scientific and Technological Achievements Project of Sun Yat-sen University (80000-18843235) and the Youth Teacher Training Project of Sun Yat-sen University (17ykpy68 and 18ykpy09). There are no competing interests related to this study.

TRIAL REGISTRATION NUMBER

N/A.

RNA N6-methyladenosine modification, spermatogenesis, and human male infertility

RNA N6-methyladenosine (m6A) modification is one of the main forms of posttranscriptional modification, and its dysregulation is involved in a series of pathological processes. RNA m6A regulators, which mediate dynamic RNA m6A modification, are expressed in almost all types of testicular cells, including spermatogenetic cells and somatic cells. Cumulative studies have found that knockout of RNA m6A regulators in the testis leads to abnormal metabolism of the target mRNAs, which eventually causes spermatogenetic disorders and infertility. To date, a role for dysregulated RNA m6A modification in human male infertility remains elusive; however, dysregulated expression of RNA m6A regulators in abnormal human semen samples, including oligospermia, asthenozoospermia and azoospermia, has been found. Therefore, we speculate that abnormal RNA m6A methylation may be an important mechanism of male infertility. In this review, we summarize the recent findings regarding the spatiotemporal expression of RNA m6A regulators in the testes, mechanisms of RNA m6A modification in spermatogenesis and the relation between dysregulated RNA m6A regulators and human male infertility. In addition, we also discuss future directions in studying the molecular mechanism of male infertility and exploring their clinical applications from the viewpoint of RNA m6A modification.

Transcriptomics of Meiosis in the Male Mouse

Molecular studies of meiosis in mammals have been long relegated due to some intrinsic obstacles, namely the impossibility to reproduce the process in vitro, and the difficulty to obtain highly pure isolated cells of the different meiotic stages. In the recent years, some technical advances, from the improvement of flow cytometry sorting protocols to single-cell RNAseq, are enabling to profile the transcriptome and its fluctuations along the meiotic process. In this mini-review we will outline the diverse methodological approaches that have been employed, and some of the main findings that have started to arise from these studies. As for practical reasons most studies have been carried out in males, and mostly using mouse as a model, our focus will be on murine male meiosis, although also including specific comments about humans. Particularly, we will center on the controversy about gene expression during early meiotic prophase; the widespread existing gap between transcription and translation in meiotic cells; the expression patterns and potential roles of meiotic long non-coding RNAs; and the visualization of meiotic sex chromosome inactivation from the RNAseq perspective.

Insights and clinical potential of proteomics in understanding spermatogenesis

Introduction: With the worldwide decline on male fertility potential, the importance of the insight of the spermatogenic process has been increasing. In recent years, proteomic methodologies have been applied to seminal fluid of infertile men to search for infertility potential biomarkers. However, to understand the spermatogenic event and to search for treatment to spermatogenic impairment, comparative analysis of testicular proteomics is considered a powerful methodology.Areas covered: Herein, we present a critical overview of the studies addressing proteomic alterations in the development of spermatogenesis during puberty, as well as during the different phases of the spermatogenic event. The comparative studies of the proteomic testicular profile of men with and without spermatogenic impairment are also discussed and key proteins and pathways involved highlighted.Expert opinion: The usage of whole human testicular tissue with its heterogeneous cellular composition makes proteome data interpretation particularly challenging. This may be minimized by controlled experiments involving the collection of testicular tissue and sperm from the same individuals, integrated in a clinically characterized cohort of healthy and infertile men. The analysis of specific subcellular proteomes can add more information to the proteomic puzzle, opening new treatment possibilities for infertile/subfertile men.

Deletion of lncRNA5512 has no effect on spermatogenesis and reproduction in mice

Long non-coding (lnc) RNAs are a series of RNAs longer than 200 nucleotides that do not code for protein products. Whole-genome expression profiles of lncRNAs suggest that they play important roles in spermatogenesis because they are particularly abundant in testes. However, most of their characteristics and functions remain unclear. The aim of this study was to define the function of lncRNA5512, which is abundant in spermatocytes and round spermatids, in mouse fertility invivo. To investigate this we generated lncRNA5512-knockout mice by clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) 9 technology. Knockout mice showed normal spermatogenesis and fertility, and had no detectable abnormalities. This indicates that lncRNA5512 does not affect mouse fertility despite its high expression in the testes. Its specific localisation in spermatocytes and round spermatids suggests that it could be a useful marker for the identification of spermatocytes and round spermatids in mouse testes.

Panoramic transcriptome analysis and functional screening of long noncoding RNAs in mouse spermatogenesis

Long noncoding RNAs (lncRNAs) have emerged as diverse functional regulators involved in mammalian development; however, large-scale functional investigation of lncRNAs in mammalian spermatogenesis in vivo is lacking. Here, we delineated the global lncRNA expression landscape in mouse spermatogenesis and identified 968 germ cell signature lncRNAs. By combining bioinformatics and functional screening, we identified three functional lncRNAs (Gm4665, 1700027A15Rik, and 1700052I22Rik) that directly influence spermatogenesis in vivo. Knocking down Gm4665 hampered the development of round spermatids into elongating spermatids and disrupted key spermatogenic gene expression. Mechanistically, lncRNA Gm4665 localized in the nucleus of round spermatids and occupied the genomic regulatory region of important spermatogenic genes including Ip6k1 and Akap3. These findings provide a valuable resource and framework for future functional analysis of lncRNAs in spermatogenesis and their potential roles in other biological processes.

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.

RNA-Seq analysis reveals functionally relevant coding and non-coding RNAs in crossbred bull spermatozoa

•

RNA-Seq analysis was done to characterize the transcriptome of crossbred bull spermatozoa.

•

Among the 13,814 transcripts detected, 431 had FPKM > 1 and 13,673 had FPKM > 0 or < 1.

•

Coding and non-coding RNAs account for 13,145 (95.15%) and 152 (1.1%), respectively.

•

Sperm transcripts were mainly related to ribosome, oxidative phosphorylation and spliceosome pathways.

•

qPCR analysis showed individual variations in transcriptional abundance of selected genes.

Sperm, which are believed to be transcriptionally and translationally inactive, synthesize RNA and proteins before there is gradual disappearance of the ribosome during chromatin compaction. Sperm transfer several functionally relevant transcripts to the oocyte, controlling maternal-zygotic transition and embryonic development. The present study was undertaken to profile and analyze sperm transcripts comprehensively using Next Generation Ribonucleic acid sequencing technology in Holstein Friesian x Tharparkar crossbred bulls. The results from global transcriptomic profiling revealed transcripts for 13,814 genes; of which 431 transcripts were expressed with >1 FPKM and 13,383 transcripts were expressed with >0 or <1 FPKM. The abundant mRNA transcripts of crossbred bull sperm were PRM1 and HMGB4. Gene ontology of transcripts with>1 FPKM revealed there was a major involvement in the structural constituent of ribosomes and translation. Results from pathway enrichment indicated the connection between ribosome, oxidative phosphorylation and spliceosome pathways and the transcripts of crossbred bull spermatozoa. The transcriptional abundance of selected genes, validated using RT-qPCR, indicated significant variations between bulls. Collectively, it may be inferred that the transcripts in crossbred bull sperm were heavily implicated in functions such as the structural constituent of ribosomes and translation, and pathways such as ribosome, oxidative phosphorylation and spliceosome. Further studies using larger sample sizes are required to understand the possible implications of transcriptomic variations on semen quality and fertility.

Transcript Isoforms of SLC7A11-AS1 Are Associated With Varicocele-Related Male Infertility

Oxidative stress is one of the crucial mediators of varicocele-related male infertility. Recently, roles of long noncoding RNAs (lncRNAs) in oxidative stress have begun to emerge, however, little is known about their role in male infertility. The aim of this study was to determine the role of lncRNA SLC7A11-AS1 in varicocele-related male infertility. Through a high-throughput bioinformatics investigation, we predicted that lncRNA SLC7A11-AS1 might be involved in this type of infertility. The reactive oxygen species (ROS) levels and expression levels of SLC7A11-AS1 isoforms were evaluated in ejaculated spermatozoa of 25 infertile men with varicocele and 17 fertile individuals as control. Isoform 6 of SLC7A11-AS1 that showed a significant elevation in infertile men with varicocele relative to the fertile group was overexpressed in testicular-derived carcinoma cell lines (NT2 and NCCIT) followed by assessment of ROS, glutathione (GSH), lipid peroxidation, and cell viability. Overexpression of SLC7A11-AS1 isoform 6 in NT2 and NCCIT cell lines resulted in a significant downregulation of SLC7A11 gene expression, which consequently decreased GSH levels and concomitantly increased ROS levels and enhanced lipid peroxidation, which jeopardized cell survival and promoted cell death. Our finding revealed a potential role of oxidative-related lncRNAs in the pathophysiology of male infertility associated with varicocele.

Investigating the regulatory function of the ANO1-AS2 on the ANO1 gene in infertile men with asthenozoospermia and terato-asthenozoospermia

Long non-coding RNAs (lncRNAs) have a particular expression in the testicular tissue and exhibit a regulatory function on the reproduction system. ANO1-AS2 (linc02584), as an lncRNA is located near the anoctamin1 (ANO1) gene. ANO1 is an important component of the transmembrane system exhibiting expression modifications in the idiopathic infertile men. Therefore, the present study was conducted to investigate the relationship between ANO1-AS2 and ANO1 gene expression with sperm motility and morphology in the patients with asthenozoospermia (AZ) and terato- asthenozoospermia (TAZ). The study population included 32 patients with AZ, 35 patients with TAZ, and 34 people with normozoospermia (NZ, control). The expression levels of ANO1 gene and ANO1-AS2 in the spermatozoa were measured by the quantitative real-time polymerase chain reaction (PCR). Docking analysis was performed to investigate the interactions of the ANO1 gene promoter and intermediate elements with ANO1-AS2. ANO1 gene expression was significantly (P < 0.05) downregulated in the patients however; ANO1-AS2 expression was significantly upregulated (P < 0.05). The subsequent analysis confirmed the inverse correlation between ANO1 and ANO1-AS2. ANO1 gene expression level was significantly positively correlated with sperm motility and morphology (P < 0.05). Moreover, ANO1-AS2 expression showed an inverse correlation with sperm motility and morphology (P < 0.05). Docking analysis confirmed that ANO1-AS2 could stably interact with ANO1 gene promoter. In conclusion, ANO1-AS2 is likely to downregulate the ANO1 gene by interacting with ANO1 gene promoter, which can influence the sperm motility and morphology.

Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides that are transcribed from non-coding loci yet undergo biosynthesis similar to coding mRNAs. The disproportional number of lincRNAs expressed in testes suggests that lincRNAs are important during gametogenesis, but experimental evidence has implicated very few lincRNAs in this process. We took advantage of the relatively limited number of lincRNAs in the genome of the nematode Caenorhabditis elegans to systematically analyse the functions of lincRNAs during meiosis. We deleted six lincRNA genes that are highly and dynamically expressed in the C. elegans gonad and tested the effects on central meiotic processes. Surprisingly, whereas the lincRNA deletions did not strongly impact fertility, germline apoptosis, crossovers, or synapsis, linc-4 was required for somatic growth. Slower growth was observed in linc-4-deletion mutants and in worms depleted of linc-4 using RNAi, indicating that linc-4 transcripts are required for this post-embryonic process. Unexpectedly, analysis of worms depleted of linc-4 in soma versus germline showed that the somatic role stems from linc-4 expression in germline cells. This unique feature suggests that some lincRNAs, like some small non-coding RNAs, are required for germ-soma interactions.

Hypoxia-related long noncoding RNAs are associated with varicocele-related male infertility

One of the main molecular causes that contributes to varicocele-related male infertility is excess production of reactive oxygen species (ROS). It is believed that hypoxia is an important stimulator of ROS in this condition. Recently, the significant roles of long non-coding RNAs (lncRNAs) in hypoxia response have emerged. Despite the investigation of hypoxia, there is scant information about the role of hypoxia-responding lncRNAs in varicocele-related male infertility. In the present study, we deduced eight hypoxia-responding lncRNAs based on high-throughput RNA sequencing data from two Gene Expression Omnibus (GEO) datasets. We used qRT-PCR to assess the expression levels of some of these lncRNAs in 42 ejaculated spermatozoa samples from 25 infertile men with varicocele and 17 fertile men as controls. We identified significant increases in expression levels of hypoxia-related lncRNAs, MIR210HG and MLLT4-AS1 in ejaculated spermatozoa of infertile men with varicocele. These lncRNAs also showed significant positive correlations with ROS levels and meaningful negative correlations with sperm parameters (count and motility). Besides, in silico studies identified several hypoxia response elements (HREs) within selected lncRNAs promoters. Delineation of hypoxia-related lncRNAs in varicocele-related infertility provides a valuable insight into male infertility.

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.

Harnessing the full potential of reproductive genetics and epigenetics for male infertility in the era of "big data"

The complexity of male reproductive impairment has hampered characterization of the underlying genetic causes of male infertility. However, in the last 20 years, more powerful and affordable tools to interrogate the genetic and epigenetic determinants of male infertility have accelerated the number of new discoveries in the characterization of male infertility. With this explosion of new data, integration in a systems-based approach-including complete phenotypic information-to male infertility is imperative. We briefly review the current understanding of genetic and epigenetic causes of male infertility and how findings may be translated into a practical component for the diagnosis and treatment of male infertility.

Testicular biopsies microarray analysis reveals circRNAs are involved in the pathogenesis of non-obstructive azoospermia

Circular RNAs (circRNAs) have been reported to be involved in many diseases. But there is no report on circRNAs in non-obstructive azoospermia (NOA). The purpose of this paper is to explore the circular RNA expression profile and potential functions of circRNAs in NOA patients. We first preformed circRNA expression profiling analysis using a circRNA microarray in testicular samples from NOA and obstructive azoospermia (OA) patients. CircRNAs were validated by qRT-PCR. Bioinformatics analysis were used to construct the ceRNA network. GO and KEGG enrichment analysis were performed by using DAVID. Microarray analysis identified 82 differentially expressed circRNAs in NOA specimens. The differential expression of hsa_circRNA_402130, hsa_circRNA_072697, hsa_circRNA_030050, hsa_circRNA_100812 and hsa_circRNA_406168 was confirmed by qRT-PCR. Enrichment analysis revealed the association of hsa_circRNA_402130 and hsa_circRNA_072697 with multiple signaling pathways. The data indicated that circRNAs were significantly dysregulated in NOA specimens and might involve in the pathogenesis of NOA.

Two novel testis-specific long noncoding RNAs produced by <i>1700121C10Rik</i> are dispensable for male fertility in mice

Testis-specific genes are prone to affect spermatogenesis or sperm fertility, and thus may play pivotal roles in male reproduction. However, whether a gene really affects male reproduction in vivo needs to be confirmed using a gene knock-out (KO) model, a ‘gold standard’ method. Increasing studies have found that some of the evolutionarily conserved testis-enriched genes are not essential for male fertility. In this study, we report that 1700121C10Rik, a previously uncharacterized gene, is specifically expressed in the testis and produces two long noncoding RNAs (lncRNAs) in mouse: Transcript 1 and Transcript 2. qRT-PCR, northern blotting, and in situ hybridization revealed that expression of both the lncRNAs commenced at the onset of sexual maturity and was predominant in round and elongating spermatids during spermiogenesis. Moreover, we found different subcellular localization of Transcript 1 and Transcript 2 that was predominant in the cytoplasm and nucleus, respectively. 1700121C10Rik-KO mouse model disrupting Transcript 1 and Transcript 2 expression was generated by CRISPR/Cas9 to determine their role in male reproduction. Results showed that 1700121C10Rik-KO male mice were fully fertile with approximately standard testis size, testicular histology, sperm production, sperm morphology, sperm motility, and induction of acrosome reaction. Thus, we conclude that both the testis-specific 1700121C10Rik-produced lncRNAs are dispensable for male fertility in mice under standard laboratory conditions.

Mechanisms regulating mammalian spermatogenesis and fertility recovery following germ cell depletion

Mammalian spermatogenesis is a highly complex multi-step process sustained by a population of mitotic germ cells with self-renewal potential known as spermatogonial stem cells (SSCs). The maintenance and regulation of SSC function are strictly dependent on a supportive niche that is composed of multiple cell types. A detailed appreciation of the molecular mechanisms underpinning SSC activity and fate is of fundamental importance for spermatogenesis and male fertility. However, different models of SSC identity and spermatogonial hierarchy have been proposed and recent studies indicate that cell populations supporting steady-state germline maintenance and regeneration following damage are distinct. Importantly, dynamic changes in niche properties may underlie the fate plasticity of spermatogonia evident during testis regeneration. While formation of spermatogenic colonies in germ-cell-depleted testis upon transplantation is a standard assay for SSCs, differentiation-primed spermatogonial fractions have transplantation potential and this assay provides readout of regenerative rather than steady-state stem cell capacity. The characterisation of spermatogonial populations with regenerative capacity is essential for the development of clinical applications aimed at restoring fertility in individuals following germline depletion by genotoxic treatments. This review will discuss regulatory mechanisms of SSCs in homeostatic and regenerative testis and the conservation of these mechanisms between rodent models and man.

A long noncoding RNA binding to QKI-5 regulates germ cell apoptosis via p38 MAPK signaling pathway

Spermatogenesis is the complex process of male germline development and requires coordinated interactions by multiple gene products that undergo strict developmental regulations. Increasing evidence has suggested that a number of long noncoding RNAs (lncRNAs) may function as important regulatory molecules in various physiological and pathological processes by binding to specific proteins. Here, we identified a subset of QKI-5-binding lncRNAs in the mouse testis through the integrated analyses of RNA immunoprecipitation (RIP)-microarray and biological verification. Among the lncRNAs, we revealed that NONMMUT074098.2 (Lnc10), which was highly expressed in the spermatogonia and spermatocytes of the testis, interacted with QKI-5. Furthermore, Lnc10 depletion promoted germ cell apoptosis via the activation of p38 MAPK, whereas the simultaneous knockdown of QKI-5 could rescue the apoptotic phenotype and the activation of p38 MAPK, which were induced by the loss of Lnc10. These data indicated that the Lnc10-QKI-5 interaction was associated with the regulatory roles of QKI-5 and that the Lnc10-QKI-5 interaction inhibited the regulation of QKI-5 on the downstream p38 MAPK signaling pathway. Additionally, we functionally characterized the biological roles of Lnc10 and found that the knockdown of Lnc10 promoted the apoptosis of spermatogenic cells in vivo; this suggested that Lnc10 had an important biological role in mouse spermatogenesis. Thus, our study provides a potential strategy to investigate the biological significance of lncRNA-RBP interactions during male germline development.

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

Long noncoding RNAs (LncRNAs) have been identified as important regulators of testis development; however, their expression patterns and roles in sheep are not yet clear. Thus, we used stranded specific RNA-seq to profile the testis transcriptome (lncRNAs and mRNAs) in premature and mature sheep. Hormone levels and the testis index were examined, and histological analyses were performed at five stages of testis development, 5-day-old (D5), 3-month-old (3M), 6-month-old (6M), 9-month-old (9M), and 2-year-old (2Y), the results of which indicate a significant difference in hormone levels and testis morphometries between the 3M and 9M stages (P < 0.05). Based on the comparison between 3M and 9M samples, we found 1,118 differentially expressed (DE) lncRNAs and 7,253 DE mRNAs in the testes, and qRT-PCR analysis showed that the results correlated well with the transcriptome data. Furthermore, we constructed lncRNA-protein-coding gene interaction networks. Forty-seven DE lncRNA-targeted genes enriched for male reproduction were obtained by cis- and trans-acting; 51 DE lncRNAs and 45 cis-targets, 2 DE lncRNAs and 2 trans-targets were involved in this network. Of these, 5 lncRNAs and their targets, PRKCD, NANOS3, SERPINA5, and CYP19A1, were enriched for spermatogenesis and male gonad development signaling pathways. We further examined the expression levels of 5 candidate lncRNAs and their target genes during testis development. Lastly, the interaction of lncRNA TCONS__00863147 and its target gene PRKCD was validated in vitro in sheep Leydig cells. This study provides a valuable resource for further study of lncRNA function in sheep testis development and spermatogenesis.

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.

LncRNA Gm2044 promotes 17β-estradiol synthesis in mpGCs by acting as miR-138-5p sponge

Long noncoding RNAs (lncRNAs) have been demonstrated to play vital roles in mammalian reproduction. Our previous research revealed that lncRNA Gm2044 is highly expressed in mouse spermatocytes and regulates male germ cell function. The gene annotation database BioGPS shows that Gm2044 is not only highly expressed in testicular tissue but also in ovarian tissue, which suggests that Gm2044 may be involved in female reproductive development. In this study, we confirmed that lncRNA Gm2044 promotes 17β-estradiol synthesis in mouse pre-antral follicular granulosa cells (mpGCs). Furthermore, bioinformatics methods, western blot, and the luciferase assay proved that Gm2044 functions as a miR-138-5p sponge to inhibit the direct target of miR-138-5p, Nr5a1, which enhances 17β-estradiol synthesis through cyp19a1 activation. Taken together, our results provide an insight into the mechanistic roles of lncRNA Gm2044 for 17β-estradiol synthesis by acting as competing-endogenous RNAs to modulate the function of mpGCs. Studying the potential lncRNAs, which regulate estradiol release, will be beneficial for the diagnosis and treatment of steroid hormone-related disease.

Non-coding RNA regulation in reproduction: Their potential use as biomarkers

Non-coding RNAs (ncRNAs) are crucial regulatory elements in most biological processes and reproduction is also controlled by them. The different types of ncRNAs, as well as the high complexity of these regulatory pathways, present a complex scenario; however, recent studies have shed some light on these questions, discovering the regulatory function of specific ncRNAs on concrete reproductive biology processes. This mini review will focus on the role of ncRNAs in spermatogenesis and oogenesis, and their potential use as biomarkers for reproductive diseases or for reproduction success.

Insight into bovine (Bos indicus) spermatozoal whole transcriptome profile

Mature spermatozoa harbor both coding and non-coding type of RNAs which regulates spermatogenesis, fertilization and early development. Characterization of bovine sperm transcriptome can provide more insight into the molecular mechanisms involved in these processes. Here, we have analyzed whole transcriptome profile of Bos indicus spermatozoa to access the global RNA expression. RNA-Seq analysis identified 14,306 genes expressed with FPKM >0, while 405 genes expressed when threshold increased to FPKM >5. Functional annotations showed that sperm transcripts were associated with molecular processes (translation, ribosomal small and large subunit assembly) and cellular components (cytosolic small and large ribosomal subunit and membranes) related to known sperm functions at fertilization and spermatogenesis. The RNA-Seq data was validated using droplet digital PCR where both highly abundant gene viz. RN7SL1 and less abundant gene viz. ZFP280B were validated. This study may provide future directions in reproductive biology of Bos indicus.

Transcriptome profiling of fetal Klinefelter testis tissue reveals a possible involvement of long non-coding RNAs in gonocyte maturation

In humans, the most common sex chromosomal disorder is Klinefelter syndrome (KS), caused by the presence of one or more extra X-chromosomes. KS patients display a varying adult phenotype but usually present with azoospermia due to testicular degeneration, which accelerates at puberty. The timing of the germ cell loss and whether it is caused by dysgenetic fetal development of the testes is not known. We investigated eight fetal KS testes and found a marked reduction in MAGE-A4-positive pre-spermatogonia compared with testes from 15 age-matched controls, indicating a failure of the gonocytes to differentiate into pre-spermatogonia. Transcriptome analysis by RNA-sequencing of formalin-fixed, paraffin-embedded testes originating from four fetal KS and five age-matched controls revealed 211 differentially expressed transcripts in the fetal KS testis. We found a significant enrichment of upregulated X-chromosomal transcripts and validated the expression of the pseudoautosomal region 1 (PAR1) gene, AKAP17A. Moreover, we found enrichment of long non-coding RNAs in the KS testes (e.g. LINC01569 and RP11-485F13.1). In conclusion, our data indicate that the testicular phenotype observed among adult men with KS is initiated already in fetal life by failure of the gonocyte differentiation into pre-spermatogonia, which could be due to aberrant expression of long non-coding RNAs.

Identification of programmed cell death 1 and its ligand in the testicular tissue of mice

PROBLEM

This study aims to determine the expression and localization of programmed cell death 1 (PD-1) and programmed cell death 1 ligand 1 (PD-L1) in the testes of mice at different developmental stages.

METHOD OF STUDY

By means of RT-qPCR, Western blot and immunofluorescence, the expression and localization of PD-1 and PD-L1 were detected in the testicular tissues of mice at different postnatal times: P7, P14, P21, P28, P35, and adulthood. Meanwhile, the level of soluble PD-L1 (sPD-L1) was evaluated by ELISA in the testicular interstitial fluid (IF) of the adult mice, culture supernatants of TM4 cell lines (Sertoli cells lines), and primary Sertoli cells at P14.

RESULTS

Pd-1 mRNA levels were unexpectedly low. From P7 to P21, there was limited PD-1 protein detected while PD-1 was evident at P28 and afterward at significantly higher levels than at P14 and P21 (P < 0.05). Despite being found in the interstitial area at P7, P14, and P21, PD-1 was also detected in the germ cells of the seminiferous tubules after P28. Pd-l1 mRNA exhibited age-related changes, peaking at P21, while PD-L1 protein was constitutively expressed at any stage, specifically localized in the nucleus of Sertoli cells. Moreover, the level of sPD-L1 in IF was significantly higher than that in the culture supernatants of both TM4 and primary Sertoli cells at P14.

CONCLUSIONS

PD-1 and PD-L1 were present in the testicular tissue of adult mice. The expression and localization of PD-1 fluctuated with age, and PD-1 was mainly localized to advanced germ cells, suggesting that it may play a role in spermiogenesis. PD-L1 was constitutively expressed in the nucleus of Sertoli cells, which could secrete sPD-L1 into the testicular interstitial space and thus may be involved in testicular immune privilege.

Long non-coding RNA: its evolutionary relics and biological implications in mammals: a review

The central dogma of gene expression propounds that DNA is transcribed to mRNA and finally gets translated into protein. Only 2–3% of the genomic DNA is transcribed to protein-coding mRNA. Interestingly, only a further minuscule part of genomic DNA encodes for long non-coding RNAs (lncRNAs) which are characteristically more than 200 nucleotides long and can be transcribed from both protein-coding (e.g. H19 and TUG1) as well as non-coding DNA by RNA polymerase II. The lncRNAs do not have open reading frames (with some exceptions), 3`-untranslated regions (3’-UTRs) and necessarily these RNAs lack any translation-termination regions, however, these can be spliced, capped and polyadenylated as mRNA molecules. The flexibility of lncRNAs confers them specific 3D-conformations that eventually enable the lncRNAs to interact with proteins, DNA or other RNA molecules via base pairing or by forming networks. The lncRNAs play a major role in gene regulation, cell differentiation, cancer cell invasion and metastasis and chromatin remodeling. Deregulation of lncRNA is also responsible for numerous diseases in mammals. Various studies have revealed their significance as biomarkers for prognosis and diagnosis of cancer. The aim of this review is to overview the salient features, evolution, biogenesis and biological importance of these molecules in the mammalian system.

Noncoding RNAs: Potential players in the self-renewal of mammalian spermatogonial stem cells

Spermatogonial stem cells (SSCs), a unique population of male germ cells with self-renewal ability, are the foundation for maintenance of spermatogenesis throughout the life of the male. Although many regulatory molecules essential for SSC self-renewal have been identified, the fundamental mechanism underlying how SSCs acquire and maintain their self-renewal activity remains largely to be elucidated. In recent years, many types of noncoding RNAs (ncRNAs) have been suggested to regulate the SSC self-renewal through multiple ways, indicating ncRNAs play crucial roles in SSC self-renewal. In this paper, we mainly focus on four types of ncRNAs including microRNA, long ncRNA, piwi-interacting RNA, as well as circular RNAs, and reviewed their potential roles in SSC self-renewal that discovered recently to help us gain a better understanding of molecular mechanisms by which ncRNAs perform their function in regulating SSC self-renewal.

New transcriptomic tools to understand testis development and functions

The testis plays a central role in the male reproductive system - secreting several hormones including male steroids and producing male gametes. A complex and coordinated molecular program is required for the proper differentiation of testicular cell types and maintenance of their functions in adulthood. The testicular transcriptome displays the highest levels of complexity and specificity across all tissues in a wide range of species. Many studies have used high-throughput sequencing technologies to define the molecular dynamics and regulatory networks in the testis as well as to identify novel genes or gene isoforms expressed in this organ. This review intends to highlight the complementarity of these transcriptomic studies and to show how the use of different sequencing protocols contribute to improve our global understanding of testicular biology.

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.

LncRNA Gm2044 highly expresses in spermatocyte and inhibits Utf1 translation by interacting with Utf1 mRNA

Spermatogenesis is a process which includes the following phases: spermatogonial stem cell proliferation and differentiation, spermatogonia, spermatocyte, spermatid and mature sperm. Spermatogenic failure is the important factor resulting in male infertility. Recent studies showed that long noncoding RNA (lncRNA) have been found to be involved in the regulation of male reproduction. However, lncRNA associated with spermatogenesis and their mechanisms of action are unclear. The aim of this study is to explore the role and molecular mechanism of lncRNA in spermatogenesis. LncRNA microarray of germ cells and bioinformatic analysis showed lncRNA Gm2044 may play potential roles in spermatogenesis. The expression level of RNA and protein were analyzed by RT-qPCR and western blotting, respectively. The interaction of lncRNA with mRNA was detected by RNA pull down and cellular proliferation was measured using CCK-8 reagent. Testis-enriched lncRNA Gm2044 is abundant in mouse spermatocytes. Gm2044 can suppress the translation of adjacent spermatogenesis-related gene Utf1 by interacting with Utf1 mRNA. Furthermore, the proliferation of mouse spermatogonia GC-1 cell line and spermatocyte GC-2 cell line was inhibited by Gm2044.

CONCLUSION

LncRNA Gm2044 was identified to inhibit Utf1 mRNA translation and play important roles in spermatogenesis.

Different expression patterns of sperm motility-related genes in testis of diploid and tetraploid cyprinid fish†

Sperm motility is an important standard to measure the fertility of male. In our previous study, we found that the diploid spermatozoa from allotetraploid hybrid (4nAT) had longer durations of rapid and slow progressive motility than haploid spermatozoa from common carp (COC). In this study, to explore sperm motility-related molecular mechanisms, we compared the testis tissues transcriptomes from 2-year-old male COC and 4nAT. The RNA-seq data revealed that 2985 genes were differentially expressed between COC and 4nAT, including 2216 upregulated and 769 downregulated genes in 4nAT. Some differentially expressed genes, such as tubulin genes, dynein, axonemal, heavy chain(dnah) genes, mitogen-activated protein kinase(mapk) genes, tektin 4, FOX transcription factors, proteasome genes, and ubiquitin carboxyl-terminal hydrolase(uchl) genes, are involved in the regulation of cell division, flagellar and ciliary motility, gene transcription, cytoskeleton, energy metabolism, and the ubiquitin–proteasome system, suggesting that these genes were related to sperm motility of the 4nAT. We confirmed the differential expression of 12 such genes in 4nAT by quantitative PCR. By western blotting, we also confirmed increased expression of Uchl3 in 4nAT testis. In addition, we identified 1915 and 2551 predicted long noncoding RNA (lncRNA) transcripts from testis tissue transcriptomes of COC and 4nAT, respectively. Of these, 1575 lncRNAs were specifically expressed in 4nAT and 939 were specifically expressed in COC. This study provides insights into the transcriptome profile of testis tissues from diploid and tetraploid, which are useful for research on regulatory mechanisms behind sperm motility in male polyploidy.

1700108J01Rik and 1700101O22Rik are mouse testis-specific long non-coding RNAs

Long non-coding RNAs (lncRNAs; > 200 nucleotides in length) have attracted attention as fine-tuners of gene expression. However, little is known about the cell- and stage-specific expression pattern and function of lncRNAs in spermatogenesis. The purpose of this study was to identify mouse testis-associated lncRNAs using a combination of computational and experimental approaches. We first used the FANTOM5 database to survey lncRNA expression in the mouse testis and performed reverse transcription quantitative polymerase chain reaction (real-time PCR) and in situ hybridization (ISH) analyses. In silico analysis showed that most of the highly expressed lncRNAs in the adult mouse testis were testis-specific lncRNAs and were expressed at and following the initiation of spermatogenesis. We selected the antisense lncRNA 1700108J01Rik and long intergenic non-coding RNA 1700101O22Rik from the most highly expressed lncRNAs in the adult testis for further analysis. Real-time PCR analysis confirmed that 1700108J01Rik and 1700101O22Rik were specifically expressed in the testis. ISH analysis revealed that the two mouse-testis-specific lncRNAs were expressed exclusively in testicular germ cells in meiotic prophase and the round spermatid stage, which coincide with the period of transcriptional reactivation during spermatogenesis. The cytoplasmic distribution of these lncRNAs revealed by ISH suggests their involvement in post-transcriptional gene regulation rather than in epigenetic or transcriptional regulation. Our data provide new insight into testis-associated lncRNAs that will be useful in expression and functional studies of spermatogenesis.