Integrated analysis of mRNAs and long noncoding RNAs in the semen from Holstein bulls with high and low sperm motility

Harnessing the value of fertility biomarkers in bull sperm for buck sperm

Efficient and sustainable reproduction and production of cattle and goats are vitally important for ensuring global food security. There is a need for potent biomarkers to accurately evaluate semen quality and predict male fertility. Although there is a reasonable set of biomarkers identified in bull sperm, there is a significant lack of such information in buck sperm along with a lack of transfer of proven technologies in goat reproductive biotechnology. These gaps are important problems because they are preventing advances in fundamental andrology and applied science of goat production. Both cattle and goats are ruminants, and they share significant similarities in their genetics and physiology although subtle differences do exist. This review harnesses the power of utilizing the knowledge developed in bull sperm to generate information on buck sperm fertility markers. These include genomic, functional genomic, epigenomic fertility markers. Revealing molecular underpinnings of such similarity and diversity using systems biology is expected to advance both fundamental and applied andrology of livestock and endangered species.

Transcriptomics of bovine sperm and oocytes

Traditionally, sperm and embryos were studied using microscopy to assess morphology and motility. However, OMICS technologies, especially transcriptomic analysis, are now being used to screen the molecular dynamics of fertility markers at cellular and molecular levels, with high sensitivity. Transcriptomics is the study of the transcriptome - RNA transcripts produced by the genome - using high-throughput methods to understand how the RNAs are expressed. In this review, we have discussed gene contributions to sperm structure and function and their role in fertilization and early embryo development. Further, we identified miRNAs shared by sperm, oocytes, and early embryos and their roles in fertilization and early embryo development.

Development and optimization of an efficient RNA isolation protocol from bovine (Bos indicus) spermatozoa

Achieving the optimum extraction of RNA from spermatozoal cells is crucial for carrying out effective high-throughput analysis regarding its role in fertility and other reproduction processes in Bos indicus. Nevertheless, semen comprises spermatozoa and several other secretions from the male reproductive system, which as well consist of diverse somatic cell types. Therefore, the elimination of somatic cells guarantees the purity of the sperm RNA. In the present study, we tested five different RNA isolation protocols and evaluated them for their yield and purity using spectrophotometer and polymerase chain reaction. Among the five RNA isolation protocols, the Triazol + RNAeasy plus Kit + TCEP method revealed optimum performance. We successfully achieved isolation of spermatozoal RNA without any spermatozoal DNA contamination from Bos indicus spermatozoa that contains approx. 1000 to 10,000 times less RNA as compared to other mammalian somatic cells. RNA quality was assessed using primers Protamine1 (spermatozoal RNA and spermatozoal DNA), CDH1 (epithelial cell), KIT (germ cell) and PTPRC (leukocytes) designed using primer BLAST where there was no product amplified except Prm1 whose product size was specific for spermatozoal RNA. The results of our investigation on RNA isolation procedures indicate that the inclusion of a disulphide reducing agent (TCEP) is crucial for the process of sperm cell lysis.

Single Sperm RNA signatures reveal MicroRNA biomarkers for male subfertility

PURPOSE

To investigate small RNA profiles in sperm, identify stable miRNA patterns unique to sperm, and assess the behavior of consistently expressed miRNAs in sperm from subfertile men compared to fertile controls.

METHODS

The small RNA profiles of single sperm from four proven fertile men were analyzed using Small RNA next-generation sequencing (NGS). Subsequently, a specific set of miRNAs was validated using RT-qPCR on additional sperm samples from 65 subfertile men from an infertility clinic and 30 proven fertile men.

RESULTS

Small RNA sequencing revealed a diverse range of sperm small RNA biotypes, including miRNAs. The mapped read percentage ranged from 22.19% for single sperm to 83.29% for enriched sperm samples used at different RNA concentrations. In single sperm, a smaller proportion of sequences were attributed to piRNAs (2.79%), miRNA (0.94%), tRNA (0.82%), and rRNA (0.47%) compared to enriched sperm samples, where piRNA (41.68%), tRNA (20.31%), miRNA (11.11%), and rRNA (6.54%) were observed. Distinct detection rates and a higher number of detected miRNAs were noted with enriched sperm samples compared to single sperm obtained using either a micromanipulator or microdissection systems. Among the identified miRNAs, 110 were consistently present in all samples. RT-qPCR revealed 15 miRNAs with increased expression and 5 miRNAs with decreased expression in sperm samples from subfertile men compared to proven fertile men. These differentially validated miRNAs were significantly correlated, either positively or negatively, with sperm count, motility, and morphology.

CONCLUSION

The study extensively examines small RNAs in single sperm, identifying sperm-specific miRNAs that could serve as molecular markers to distinguish between subfertile and fertile men in clinical settings.

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.

Transcriptome Analysis Reveals Spermatogenesis-Related CircRNAs and LncRNAs in Goat Spermatozoa

Mammalian spermatozoa comprises both coding and non-coding RNAs, which are traditionally believed to be a residual of spermatogenesis. The differential expression level of spermatozoal RNAs is also observed between fertile and infertile human, thereby anticipated as potential molecular marker of male fertility. This study investigated the transcriptome profile of goat (Capra hircus) spermatozoa. The sperm transcriptome was analyzed by three different methods viz. RLM-RACE, long-read RNA sequencing (RNAseq) in Nanopore™ platform, and short-read RNAseq in Illumina™ platform. The Illumina™ sequencing discovered 16,604 transcripts with 357 mRNAs having FPKM (fragments per kilobase per million mapped reads) of more than five. The spermatozoal RNA suite included mRNA (94%), rRNA (3%), miscRNA (1%), circRNA (1%), miRNA (1%), etc. This study also predicted circRNAs (127), lncRNAs (655), and imprinted genes (160) that have potential role in male reproduction. The gene ontology analysis revealed the involvement of spermatozoal RNA in regulating male meiosis (TET3, STAT5B), capacitation (ACRBP, CATSPER4), sperm motility (GAS8, TEKT2), spermatogenesis (ADAMTS2, CREB3L4), etc. The spermatozoal RNA were also associated with different biological pathways viz. Wnt signaling pathway, cAMP signaling pathway, AMPK signaling pathway, and MAPK signaling pathways having potential role in spermatogenesis. Overall, this study enlightened the suite of spRNA transcripts in goat and their relevance in male fertility for diagnostic approach.

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.

Coordinated alternation of DNA methylation and alternative splicing of PBRM1 affect bovine sperm structure and motility

DNA methylation and gene alternative splicing drive spermatogenesis. In screening DNA methylation markers and transcripts related to sperm motility, semen from three pairs of full-sibling Holstein bulls with high and low motility was subjected to reduced representation bisulphite sequencing. A total of 948 DMRs were found in 874 genes (gDMRs). Approximately 89% of gDMR-related genes harboured alternative splicing events, including SMAD2, KIF17, and PBRM1. One DMR in exon 29 of PBRM1 with the highest 5mC ratio was found, and hypermethylation in this region was related to bull sperm motility. Furthermore, alternative splicing events at exon 29 of PBRM1 were found in bull testis, including PBRM1-complete, PBRM1-SV1 (exon 28 deletion), and PBRM1-SV2 (exons 28-29 deletion). PBRM1-SV2 exhibited significantly higher expression in adult bull testes than in newborn bull testes. In addition, PBRM1 was localized to the redundant nuclear membrane of bull sperm, which might be related to sperm motility caused by sperm tail breakage. Therefore, the hypermethylation of exon 29 may be associated with the production of PBRM1-SV2 in spermatogenesis. These findings indicated that DNA methylation alteration at specific loci could regulate gene splicing and expression and synergistically alter sperm structure and motility.

Effect of arginine-induced motility and capacitation on RNA population in goat spermatozoa

INTRODUCTION

In vitro capacitation is essential in assisted reproductive technologies (ART) for embryo production. Recently, arginine has been proven to enhance capacitation in mammalian spermatozoa. However, the detailed mechanism of action of arginine remains elusive.

AIM

This study investigated the effect of arginine-induced capacitation and motility enhancement on the spermatozoal RNA (spRNA) population in goats.

MATERIAL AND METHODS

Goat spermatozoa were treated with arginine for up to six hours and compared with non-treated or PHE (penicillamine, hypotaurine, and epinephrine)-treated spermatozoa at different intervals (0, 1, 2, 4, and 6 hours). Sperm parameters, including viability, individual motility, capacitation, acrosome reaction, and ROS production, were evaluated. The spRNA population was analyzed by short-read RNA sequencing (RNA-seq).

RESULTS

The percentage of capacitated (73.21 ± 4.22%) and acrosome reacted (18.35 ± 0.56%) spermatozoa was highest in arginine treatment, while PHE treatment showed the highest percentage (79.82 ± 4.31%) of motile spermatozoa from 0 to 4 hours of incubation. RNA-seq analysis identified 1,321 differentially expressed genes (DEGs) in arginine-treated spermatozoa compared to the control. The PGK2, RNASE10, ODF1, and ROPN1L genes involved in sperm motility and ACR, DKKL1, KCNJ11, and PRND genes involved in the capacitation process were upregulated in arginine-treated spermatozoa. The DEGs regulate sperm capacitation-related cAMP-PKA, PI3-Akt, calcium, and MAPK signaling pathways.

CONCLUSION

The arginine-induced capacitation and enhanced sperm motility were associated with the upregulation of several genes involved in sperm motility and capacitation pathways. The comparative study also suggests that arginine may be used in lieu of PHE for motility enhancement and in vitro capacitation of goat spermatozoa.

Development and validation of most efficient RNA isolation method from buffalo bull spermatozoa

BACKGROUND

Being highly fragmented and low in concentration, isolation of good quality RNA from sperm cells is a big challenge. Attempts have been made to evaluate various sperm RNA isolation methods from purified buffalo bull sperm cells.

METHODS

Both, non-membrane and membrane-based methods have been evaluated for isolating RNA from Murrah buffalo sperms and compared for their respective efficacies. The traditional TRIzol, TRIzol-heat lysed (H-TRIzol) and cocktail of TCEP-RLT lysis buffer (Qiagen RNeasy mini kit)-TRIzol (C-TRIzol) based isopropanol isolation methods have been evaluated.

RESULTS

H-TRIzol yielded best results among conventional methods. The combined T-RLT RNA isolation protocol yielded best quality and quantity compared to other membrane-based methods, due to high lytic property of cocktail of lysis reagents, necessary for complete breakdown of sperm membrane and RNA binding membrane for RNA isolation. Combined lysis performed by treatment with RLT-T and T-RLT differing in order of reagents used were also evaluated. T-RLT combination giving better results compared to RLT-T due to high gDNA contamination and membrane clogging in later protocol steps.

CONCLUSION

Overall, in terms of total RNA quantity and quality per million spermatozoa, the heat-lysed TRIzol method (H-TRIzol) performs best among RNA separation techniques employed and is also quite easy to perform. This comparative evaluation of sperm RNA isolation protocols can be useful in deciding the best protocol for isolation of good quality and high concentration sperm RNA from buffalo semen, for transcriptome and other downstream studies.

QTLs and Candidate Genes Associated with Semen Traits in Merino Sheep

Ram semen traits play a significant role in conception outcomes, which in turn may influence reproductive efficiency and the overall productivity and profitability of sheep enterprises. Since hundreds of ewes may be inseminated from a single ejaculate, it is important to evaluate semen quality prior to use in sheep breeding programs. Given that semen traits have been found to be heritable, genetic variation likely contributes to the variability observed in these traits. Identifying such genetic variants could provide novel insights into the molecular mechanisms underlying variability in semen traits. Therefore, this study aimed to identify quantitative trait loci (QTLs) associated with semen traits in Merino sheep. A genome-wide association study (GWAS) was undertaken using 4506 semen collection records from 246 Merino rams collected between January 2002 and May 2021. The R package RepeatABEL was used to perform a GWAS for semen volume, gross motility, concentration, and percent post-thaw motility. A total of 35 QTLs, located on 16 Ovis aries autosomes (OARs), were significantly associated with either of the four semen traits in this study. A total of 89, 95, 33, and 73 candidate genes were identified, via modified Bonferroni, within the QTLs significantly associated with volume, gross motility, concentration, and percent post-thaw motility, respectively. Among the candidate genes identified, SORD, SH2B1, and NT5E have been previously described to significantly influence spermatogenesis, spermatozoal motility, and high percent post-thaw motility, respectively. Several candidate genes identified could potentially influence ram semen traits based on existing evidence in the literature. As such, validation of these putative candidates may offer the potential to develop future strategies to improve sheep reproductive efficiency. Furthermore, Merino ram semen traits are lowly heritable (0.071-0.139), and thus may be improved by selective breeding.

Harnessing male germline epigenomics for the genetic improvement in cattle

Sperm is essential for successful artificial insemination in dairy cattle, and its quality can be influenced by both epigenetic modification and epigenetic inheritance. The bovine germline differentiation is characterized by epigenetic reprogramming, while intergenerational and transgenerational epigenetic inheritance can influence the offspring's development through the transmission of epigenetic features to the offspring via the germline. Therefore, the selection of bulls with superior sperm quality for the production and fertility traits requires a better understanding of the epigenetic mechanism and more accurate identifications of epigenetic biomarkers. We have comprehensively reviewed the current progress in the studies of bovine sperm epigenome in terms of both resources and biological discovery in order to provide perspectives on how to harness this valuable information for genetic improvement in the cattle breeding industry.

Integrating sperm cell transcriptome and seminal plasma metabolome to analyze the molecular regulatory mechanism of sperm motility in Holstein stud bulls

Considering that artificial insemination is the most widely used assisted reproductive technique in the dairy industry, the semen quality of bulls is very important for selecting excellent stud bulls. Sperm motility is one of the important traits of semen quality, and related genes may be regulated by environmental factors. Seminal plasma can affect sperm cell transcriptome and further affect sperm motility through exosome or other processes. However, the molecular regulation mechanism of bull sperm motility has not been studied by combining the sperm cell transcriptome with seminal plasma metabolome. The number of motile sperm per ejaculate (NMSPE) is an integrated indicator for assessing sperm motility in stud bulls. In the present study, we selected 7 bulls with higher NMSPE (5,698.55 million +/- 945.40 million) as group H and 7 bulls with lower NMSPE (2,279.76 million +/- 1,305.69 million) as group L from 53 Holstein stud bulls. The differentially expressed genes (DEGs) in sperm cells were evaluated between the two groups (H vs. L). We conducted gene co-expression network analysis (WGCNA) on H and L groups of bulls, as well as two monozygotic twin Holstein bulls with different NMSPE values, to screen candidate genes for NMSPE. The regulatory effect of seminal plasma metabolome on the candidate genes of NMSPE was also investigated. A total of 1,099 DEGs were identified in the sperm cells of H and L groups. These DEGs were primarily concentrated in energy metabolism and sperm cell transcription. The significantly enriched Kyoto encyclopedia of genes and genomes (KEGG) pathways of the 57 differential metabolites were the aminoacyl-tRNA biosynthesis pathway and vitamin B6 metabolism pathway. Our study discovered 14 genes as the potential candidate markers for sperm motility, including FBXO39. We observed a broad correlation between transcriptome of sperm cells and seminal plasma metabolome, such as three metabolites, namely, mesaconic acid, 2-coumaric acid, and 4-formylaminoantipyrine, might regulate FBXO39 expression through potential pathways. The genes related to seminal plasma metabolites expressed in sperm cells are not only located near the quantitative trait loci of reproductive traits, but also enriched in the genome-wide association study signal of sire conception rate. Collectively, this study was the first to investigate the interplays among transcriptome of sperm cells and seminal plasma metabolome from Holstein stud bulls with different sperm motility.

Unraveling the harmful effect of oxidative stress on male fertility: A mechanistic insight

Male infertility is a widely debated issue that affects males globally. There are several mechanisms involved. Oxidative stress is accepted to be the main contributing factor, with sperm quality and quantity affected by the overproduction of free radicals. Excess reactive oxygen species (ROS) cannot be controlled by the antioxidant system and, thus, potentially impact male fertility and hamper sperm quality parameters. Mitochondria are the driving force of sperm motility; irregularities in their function may lead to apoptosis, alterations to signaling pathway function, and, ultimately, compromised fertility. Moreover, it has been observed that the prevalence of inflammation may arrest sperm function and the production of cytokines triggered by the overproduction of ROS. Further, oxidative stress interacts with seminal plasma proteomes that influence male fertility. Enhanced ROS production disturbs the cellular constituents, particularly DNA, and sperms are unable to impregnate the ovum. Here, we review the latest information to better understand the relationship between oxidative stress and male infertility, the role of mitochondria, the cellular response, inflammation and fertility, and the interaction of seminal plasma proteomes with oxidative stress, as well as highlight the influence of oxidative stress on hormones; collectively, all of these factors are assumed to be important for the regulation of male infertility. This article may help improve our understanding of male infertility and the strategies to prevent it.

Sperm Transcriptome Analysis Accurately Reveals Male Fertility Potential in Livestock

Nowadays, selection of superior male candidates in livestock as a source of frozen semen based on sperm quality at the cellular level is not considered accurate enough for predicting the potential of male fertility. Sperm transcriptome analysis approaches, such as messenger RNA levels, have been shown to correlate with fertility rates. Using this technology in livestock growth has become the principal method, which can be widely applied to predict male fertility potential in the livestock industry through the analysis of the sperm transcriptome. It provides the gene expression to validate the function of sperm in spermatogenesis, fertilization, and embryo development, as the parameters of male fertility. This review proposes a transcriptomic analysis approach as a high-throughput method to predict the fertility potential of livestock more accurately in the future.

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

The Novel Competing Endogenous Long Noncoding RNA SM2 Regulates Gonadotropin Secretion in the Hu Sheep Anterior Pituitary by Targeting the Oar-miR-16b/TGF-β/SMAD2 Signaling Pathway

Pituitary gonadotropins play a pivotal role in reproduction. Long noncoding RNAs (lncRNAs) have been identified as important regulators in the hypothalamic−pituitary−ovarian (HPO) axis associated with reproduction. However, the contributions of lncRNAs to pituitary gonadotropin secretion remain largely unknown. Therefore, this work was performed to uncover the functional mechanisms of the novel lncRNA TCONS_00083279 (lncRNA SM2) and its potential targeting pathway oar-miR-16b/TGF-beta/SMAD2, which is associated with gonadotropin secretion in sheep pituitary cells. In the present study, the lncRNA SM2 showed high expression levels in the sheep pituitary gland, and it was located in both the nucleus and the cytoplasm of pituitary cells. lncRNA SM2 knockdown inhibited pituitary cell proliferation and FSH and LH secretion. The function of the lncRNA SM2 was sponged by oar-miR-16b, and this regulated the growth and gonadotropin secretion of pituitary cells by modulating SMAD2, as shown by the dual-luciferase reporter assay. FSH and LH levels were both upregulated by SMAD2 overexpression. Moreover, the levels of the lncRNA SM2, SMAD2 and TGFR1, as well as FSH and LH, in sheep pituitary cells increased significantly under gonadotropin-releasing hormone (GnRH) stimulation (p < 0.05). This work illustrates that the lncRNA SM2 regulates gonadotropin secretion in the Hu sheep anterior pituitary by targeting the oar-miR-16b/TGF-β/SMAD2 signaling pathway, providing a valuable resource for understanding the molecular mechanisms underlying sheep reproduction.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

Sperm Transcripts Associated With Odorant Binding and Olfactory Transduction Pathways Are Altered in Breeding Bulls Producing Poor-Quality Semen

Spermatozoa carries a reservoir of mRNAs regulating sperm functions and fertilizing potential. Although it is well recognized that a considerable proportion of high genetic merit breeding bulls produce poor-quality semen, the transcriptomic alterations in spermatozoa from such bulls are not understood. In the present study, comparative high-throughput transcriptomic profiling of spermatozoa from good and poor-quality semen-producing bulls was carried out to identify the transcripts associated with semen quality. Using next-generation sequencing (NGS), we identified 11,632 transcripts in Holstein Friesian bull spermatozoa; after total hit normalization, a total of 544 transcripts were detected, of which 185 transcripts were common to both good and poor-quality semen, while 181 sperm transcripts were unique to good quality semen, and 178 transcripts were unique to poor-quality semen. Among the co-expressed transcripts, 31 were upregulated, while 108 were downregulated, and 46 were neutrally expressed in poor-quality semen. Bioinformatics analysis revealed that the dysregulated transcripts were predominantly involved in molecular function, such as olfactory receptor activity and odor binding, and in biological process, such as detection of chemical stimulus involved in sensory perception, sensory perception of smell, signal transduction, and signal synaptic transmission. Since a majority of the dysregulated transcripts were involved in the olfactory pathway (85% of enriched dysregulated genes were involved in this pathway), the expression of selected five transcripts associated with this pathway (OR2T11, OR10S1, ORIL3, OR5M11, and PRRX1) were validated using real-time qPCR, and it was found that their transcriptional abundance followed the same trend as observed in NGS; the sperm transcriptional abundance of OR2T11 and OR10S1 differed significantly (p < 0.05) between good and poor-quality semen. It is concluded that poor-quality semen showed altered expression of transcripts associated with olfactory receptors and pathways indicating the relationship between olfactory pathway and semen quality in bulls.

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

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

LncRNAs in domesticated animals: from dog to livestock species

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

Cellular and Molecular Insights Into the Etiology of Subfertility/Infertility in Crossbred Bulls (Bos taurus × Bos indicus): A Review

Crossbreeding of indigenous cattle (Bos indicus) with improved (Bos taurus) breeds gained momentum and economic relevance in several countries to increase milk production. While production performance of the crossbred offspring is high due to hybrid vigor, they suffer from a high incidence of reproductive problems. Specifically, the crossbred males suffer from serious forms of subfertility/infertility, which can have a significant effect because semen from a single male is used to breed several thousand females. During the last two decades, attempts have been made to understand the probable reasons for infertility in crossbred bulls. Published evidence indicates that testicular cytology indices, hormonal concentrations, sperm phenotypic characteristics and seminal plasma composition were altered in crossbred compared to purebred males. A few recent studies compared crossbred bull semen with purebred bull semen using genomics, transcriptomics, proteomics and metabolomics; molecules potentially associated with subfertility/infertility in crossbred bulls were identified. Nevertheless, the precise reason behind the poor quality of semen and high incidence of sub-fertility/infertility in crossbred bulls are not yet well defined. To identify the underlying etiology for infertility in crossbred bulls, a thorough understanding of the magnitude of the problem and an overview of the prior art is needed; however, such systematically reviewed information is not available. Therefore, the primary focus of this review is to compile and analyze earlier findings on crossbred bull fertility/infertility. In addition, the differences between purebred and crossbred males in terms of testicular composition, sperm phenotypic characteristics, molecular composition, environmental influence and other details are described; future prospects for research on crossbred males are also outlined.

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.

Sperm Functional Genome Associated With Bull Fertility

Bull fertility is an important economic trait in sustainable cattle production, as infertile or subfertile bulls give rise to large economic losses. Current methods to assess bull fertility are tedious and not totally accurate. The massive collection of functional data analyses, including genomics, proteomics, metabolomics, transcriptomics, and epigenomics, helps researchers generate extensive knowledge to better understand the unraveling physiological mechanisms underlying subpar male fertility. This review focuses on the sperm phenomes of the functional genome and epigenome that are associated with bull fertility. Findings from multiple sources were integrated to generate new knowledge that is transferable to applied andrology. Diverse methods encompassing analyses of molecular and cellular dynamics in the fertility-associated molecules and conventional sperm parameters can be considered an effective approach to determine bull fertility for efficient and sustainable cattle production. In addition to gene expression information, we also provide methodological information, which is important for the rigor and reliability of the studies. Fertility is a complex trait influenced by several factors and has low heritability, although heritability of scrotal circumference is high and that it is a known fertility maker. There is a need for new knowledge on the expression levels and functions of sperm RNA, proteins, and metabolites. The new knowledge can shed light on additional fertility markers that can be used in combination with scrotal circumference to predict the fertility of breeding bulls. This review provides a comprehensive review of sperm functional characteristics or phenotypes associated with bull fertility.

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.

Mitochondrial Reactive Oxygen Species (ROS) Production Alters Sperm Quality

Besides ATP production, mitochondria are key organelles in several cellular functions, such as steroid hormone biosynthesis, calcium homoeostasis, intrinsic apoptotic pathway, and the generation of reactive oxygen species (ROS). Despite the loss of the majority of the cytoplasm occurring during spermiogenesis, mammalian sperm preserves a number of mitochondria that rearrange in a tubular structure at the level of the sperm flagellum midpiece. Although sperm mitochondria are destroyed inside the zygote, the integrity and the functionality of these organelles seem to be critical for fertilization and embryo development. The aim of this review was to discuss the impact of mitochondria-produced ROS at multiple levels in sperm: the genome, proteome, lipidome, epigenome. How diet, aging and environmental pollution may affect sperm quality and offspring health—by exacerbating oxidative stress—will be also described.

A systems biology framework integrating GWAS and RNA-seq to shed light on the molecular basis of sperm quality in swine

BACKGROUND

Genetic pressure in animal breeding is sparking the interest of breeders for selecting elite boars with higher sperm quality to optimize ejaculate doses and fertility rates. However, the molecular basis of sperm quality is not yet fully understood. Our aim was to identify candidate genes, pathways and DNA variants associated to sperm quality in swine by analysing 25 sperm-related phenotypes and integrating genome-wide association studies (GWAS) and RNA-seq under a systems biology framework.

RESULTS

By GWAS, we identified 12 quantitative trait loci (QTL) associated to the percentage of head and neck abnormalities, abnormal acrosomes and motile spermatozoa. Candidate genes included CHD2, KATNAL2, SLC14A2 and ABCA1. By RNA-seq, we identified a wide repertoire of mRNAs (e.g. PRM1, OAZ3, DNAJB8, TPPP2 and TNP1) and miRNAs (e.g. ssc-miR-30d, ssc-miR-34c, ssc-miR-30c-5p, ssc-miR-191, members of the let-7 family and ssc-miR-425-5p) with functions related to sperm biology. We detected 6128 significant correlations (P-value ≤ 0.05) between sperm traits and mRNA abundances. By expression (e)GWAS, we identified three trans-expression QTL involving the genes IQCJ, ACTR2 and HARS. Using the GWAS and RNA-seq data, we built a gene interaction network. We considered that the genes and interactions that were present in both the GWAS and RNA-seq networks had a higher probability of being actually involved in sperm quality and used them to build a robust gene interaction network. In addition, in the final network we included genes with RNA abundances correlated with more than four semen traits and miRNAs interacting with the genes on the network. The final network was enriched for genes involved in gamete generation and development, meiotic cell cycle, DNA repair or embryo implantation. Finally, we designed a panel of 73 SNPs based on the GWAS, eGWAS and final network data, that explains between 5% (for sperm cell concentration) and 36% (for percentage of neck abnormalities) of the phenotypic variance of the sperm traits.

CONCLUSIONS

By applying a systems biology approach, we identified genes that potentially affect sperm quality and constructed a SNP panel that explains a substantial part of the phenotypic variance for semen quality in our study and that should be tested in other swine populations to evaluate its relevance for the pig breeding sector.

Knockdown of Long Non-Coding RNA AFAP1-AS1 Promoted Viability and Suppressed Death of Cardiomyocytes in Response to I/R In Vitro and In Vivo

Long non-coding RNA (lncRNA) plays a pivotal role in the development of myocardial infarction (MI). The aim of this study was to investigate the effects of lncRNA actin filament-associated protein 1 antisense RNA 1 (AFAP1-AS1) on cell cycle, proliferation, and apoptosis. RT-qPCR was used to detect the expression levels of AFAP1-AS1, miR-512-3p, and reticulon 3 (RTN3) in rat model of I/R. The simulated MI environment was constructed. MTT assay and flow cytometry were used to detect changes in cardiomyocyte viability and cell cycle/apoptosis after MI by AFAP1-AS1 silencing or RTN3 silencing. The targeting relationship of miR-512-3p and AFAP1-AS1 and RTN3 in cardiomyocytes was verified by dual luciferase reporter assay. The expression levels of AFAP1-AS1 and RTN3 were significantly upregulated in a rat model of LAD ligation (or MI) ligation, while the expression level of miR-512-3p was significantly reduced. Overexpressed AFAP1-AS1 and RTN3 promoted cardiomyocyte apoptosis and inhibited cardiomyocyte proliferation. MiR-512-3p was a direct target of AFAP1-AS1, and RTN3 was a direct target of miR-512-3p. AFAP1-AS1 promoted the progression of MI by targeting miR-512-3p. AFAP1-AS1 promoted the progression of MI by modulating the miR-512-3p/RTN3 axis. AFAP1-AS1 may be a potential therapy target for MI. Graphical Abstract The role of AFAP1-AS1 in regulating MI injury in vivo. (A) Effect of AFAP1-AS1 in MI injury in vivo. (B) The mRNA level of RTN3 in MI injury in vivo. (C) The protein level of RTN3 in MI injury in vivo. (D) Effect of miR-512-3p in MI model group. (E) TUNEL assay. *P < 0.05, **P < 0.01 vs the sham group; #P < 0.05, ##P < 0.01 vs the MI group.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm's response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Long non-coding RNAs (lncRNAs) in spermatogenesis and male infertility

BACKGROUND

Long non-coding RNAs (lncRNAs) have a size of more than 200 bp and are known to regulate a host of crucial cellular processes like proliferation, differentiation and apoptosis by regulating gene expression. While small noncoding RNAs (ncRNAs) such as miRNAs, siRNAs, Piwi-interacting RNAs have been extensively studied in male germ cell development, the role of lncRNAs in spermatogenesis remains largely unknown.

OBJECTIVE

In this article, we have reviewed the biology and role of lncRNAs in spermatogenesis along with the tools available for data analysis.

RESULTS AND CONCLUSIONS

Till date, three microarray and four RNA-seq studies have been undertaken to identify lncRNAs in mouse testes or germ cells. These studies were done on pre-natal, post-natal, adult testis, and different germ cells to identify lncRNAs regulating spermatogenesis. In case of humans, five RNA-seq studies on different germ cell populations, including two on sperm, were undertaken. We compared three studies on human germ cells to identify common lncRNAs and found 15 lncRNAs (LINC00635, LINC00521, LINC00174, LINC00654, LINC00710, LINC00226, LINC00326, LINC00494, LINC00535, LINC00616, LINC00662, LINC00668, LINC00467, LINC00608, and LINC00658) to show consistent differential expression across these studies. Some of the targets of these lncRNAs included CENPB, FAM98B, GOLGA6 family, RPGR, TPM2, GNB5, KCNQ10T1, TAZ, LIN28A, CDKN2B, CDKN2A, CDKN1A, CDKN1B, CDKN1C, EZH2, SUZ12, VEGFA genes. A lone study on human male infertility identified 9879 differentially expressed lncRNAs with three (lnc32058, lnc09522, and lnc98497) of them showing specific and high expression in immotile sperm in comparison to normal motile sperm. A few lncRNAs (Mrhl, Drm, Spga-lncRNAs, NLC1-C, HongrES2, Tsx, LncRNA-tcam1, Tug1, Tesra, AK015322, Gm2044, and LncRNA033862) have been functionally validated for their roles in spermatogenesis. Apart from rodents and humans, studies on sheep and bull have also identified lncRNAs potentially important for spermatogenesis. A number of these non-coding RNAs are strong candidates for further research on their roles in spermatogenesis.

Marker-assisted selection vis-à-vis bull fertility: coming full circle-a review

Bull fertility is considered an indispensable trait, as far as farm economics is concerned since it is the successful conception in a cow that provides calf crop, along with the ensuing lactation. This ensures sustainability of a dairy farm. Traditionally, bull fertility did not receive much attention by the farm managers and breeding animals were solely evaluated based on phenotypic predictors, namely, sire conception rate and seminal parameters in bull. With the advent of the molecular era in animal breeding, attempts were made to unravel the genetic complexity of bull fertility by the identification of genetic markers related to the trait. Marker-Assisted Selection (MAS) is a methodology that aims at utilizing the genetic information at markers and selecting improved populations for important traits. Traditionally, MAS was pursued using a candidate gene approach for identifying markers related to genes that are already known to have a physiological function related to the trait but this approach had certain shortcomings like stringent criteria for significance testing. Now, with the availability of genome-wide data, the number of markers identified and variance explained in relation to bull fertility has gone up. So, this presents a unique opportunity to revisit MAS by selection based on the information of a large number of genome-wide markers and thus, improving the accuracy of selection.

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.

Review: Long non-coding RNA in livestock

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

Identification and quantification of coding and long non-coding RNAs in stallion spermatozoa separated by density

BACKGROUND

It is not unusual for stallions to have fertility problems. For many, artificial insemination with more dense spermatozoa (isolated by density gradient centrifugation) results in greater pregnancy rates compared with the rates when using unfractionated spermatozoa. RNAs in spermatozoa delivered to the oocyte at conception are required for embryo development. Novel molecular assays of spermatozoa that reflect function are needed to predict the fertility of stallions.

OBJECTIVES

To describe and compare the RNA populations in more dense and less dense spermatozoa from stallions.

MATERIALS AND METHODS

Spermatozoa from five stallions were separated into more dense and less dense populations by density gradient centrifugation. Complementary DNA libraries were made from each of the ten total RNA samples after ribosomal RNA removal. Next-generation sequencing characterized the RNA populations in more and less dense spermatozoa. Quantitative reverse transcription-PCR was used to confirm differential expression of selected RNAs.

RESULTS

Stallion spermatozoa contain 11 215 RNAs, with the most prevalent RNA being a 1492 base long non-coding RNA. The levels of 159 RNAs were greater in more dense spermatozoa, while levels of seven other RNAs were greater in less dense spermatozoa. Quantitative reverse transcription-PCR confirmed the threefold greater levels of solute carrier family 26 member 8 (SLC26A8) mRNA in less dense spermatozoa, and sixfold and threefold greater expression levels of the SCP2 sterol binding domain containing 1 (SCP2D1) and spermatogenesis-associated protein 31D1 (SPATA31D1) mRNAs in more dense spermatozoa, respectively.

DISCUSSION AND CONCLUSION

We identified 11 215 RNAs in stallion spermatozoa and 166 with differential expression between more dense and less dense fractions. Many prevalent RNAs were also found in bull, boar, and human spermatozoa. Many differentially expressed RNAs are known to be testis- or spermatozoa-specific. Our results may lead to identification of an RNA population in spermatozoa that is optimal for establishing successful pregnancies.

Genome-wide methylation and transcriptome of blood neutrophils reveal the roles of DNA methylation in affecting transcription of protein-coding genes and miRNAs in E. coli-infected mastitis cows

Background

Neutrophils are the first effectors of inflammatory response triggered by mastitis infection, and are important defense cells against pathogenic Escherichia coli (E. coli). DNA methylation, as a critical epigenetic mechanism for regulating gene function, is involved in bovine mastitis.

Results

In this study, we sequenced the blood neutrophils of healthy and E. coli-infected mastitic half-sib cows for the overall DNA methylation levels using transcriptome sequencing and reduced representation bisulfite sequencing. The methylation levels in the mastitis cows (MCs) were decreased compared with healthy cows (HCs). A total of 494 differentially methylated regions were identified, among which 61 were up-methylated and 433 were down-methylated (MCs vs. HCs). The expression levels of 1094 differentially expressed genes were up-regulated, and 245 genes were down-regulated. Twenty-nine genes were found in methylation and transcription data, among which seven genes’ promoter methylation levels were negatively correlated with expression levels, and 11 genes were differentially methylated in the exon regions. The bisulfite sequencing PCR and quantitative real-time PCR validation results demonstrated that the promoter methylation of CITED2 and SLC40A1 genes affected differential expression. The methylation of LGR4 exon 5 regulated its own alternative splicing. The promoter methylation of bta-miR-15a has an indirect effect on the expression of its target gene CD163. The CITED2, SLC40A1, and LGR4 genes can be used as candidates for E. coli-induced mastitis resistance.

Conclusions

This study explored the roles of DNA methylation in affecting transcription of protein-coding genes and miRNAs in E. coli-induced mastitis, thereby helping explain the function of DNA methylation in the pathogenesis of mastitis and provided new target genes and epigenetic markers for mastitis resistance breeding in dairy cattle.

CircRNA Role and circRNA-Dependent Network (ceRNET) in Asthenozoospermia

The role of circRNA in reproduction is under investigation. CircRNAs are expressed in human testis, spermatozoa (SPZ), and seminal plasma. Their involvement in embryo development has also been suggested. Asthenozoospermia, a common cause of male infertility, is characterized by reduced or absent sperm motility in fresh ejaculate. While abnormal mitochondrial function, altered sperm tail, and genomic causes have been deeply investigated, the epigenetic signature of asthenozoospermic derived SPZ still remains unexplored. CircRNAs may take part in the repertoire of differentially expressed molecules in infertile men. Considering this background, we carried out a circRNA microarray, identifying a total of 9,138 transcripts, 22% of them novel based and 83.5% with an exonic structure. Using KEGG analysis, we evaluated the circRNA contribution in pathways related to mitochondrial function and sperm motility. In order to discriminate circRNAs with a differential expression in SPZ with differential morphological parameters, we separated sperm cells by Percoll gradient and analyzed their differential circRNA payload. A bioinformatic approach was then utilized to build a circRNA/miRNA/mRNA network. With the aim to demonstrate a dynamic contribution of circRNAs to the sperm epigenetic signature, we verified their modulation as a consequence of an oral amino acid supplementation, efficacious in improving SPZ motility.