Go to Top Go to Bottom
Anim Biosci > Volume 35(4); 2022 > Article
Chen, Miao, Xu, Pu, Li, Han, Mao, and Ma: Alterations of mRNA and lncRNA profiles associated with the extracellular matrix and spermatogenesis in goats



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


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


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


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


Mammalian testicular development and spermatogenesis are complex processes regulated by the transcriptome; they are dynamic and staged [1]. Testis-specific genes play a crucial role in male reproduction by influencing spermatogenesis and fertility [2]. Optimum fertility levels are critical to individual enterprises and to the livestock industry as a whole [3]. Past breeding efforts demonstrate that genetic selection is a highly successful tool for improving livestock populations [4]. Therefore, a comprehensive understanding of the molecular mechanisms at play within testicular tissues during goat sexual maturation will significantly promote the future success of goat breeding.
As an important part of the seminiferous tubule wall, the extracellular matrix (ECM) directly affects spermatogenesis [5]. Matrix metalloproteinases (MMPs) are a group of endopeptidases that function in the degradation of collagen in the ECM and during the regulation of homeostasis in the ECM [6]. They play an essential role in spermatogenesis [7] and remodeling by interacting with insulin-like growth factor-binding proteins (IGFBPs) and insulin-like growth factor (IGF) [8]. The matrix metalloproteinase 2 (MMP2) is an extracellular zinc protease that plays a vital role in testicular maturation and spermatogenesis [9]. In addition, IGF2 and IGFBP6 can bind with the ECM to affect a variety of cell functions [10].
Long non-coding RNAs (lncRNAs) are highly expressed non-coding RNAs that play an essential role in regulating gene expression in organisms [11], including during spermatogenesis [12]. For instance, in the process of spermatocyte division and differentiation into mature sperm, many lncRNAs are involved in synchronising the expression of specific genes to achieve regulatory purposes [13]. lncRNAs have been studied in the testicular tissues of livestock and poultry, including chicken [14], sheep [15], and pigs [16]. However, the comparative analysis of lncRNAs and mRNAs in prepubertal, pubertal, and post-pubertal testicular tissues of goats is unclear.
We used RNA-seq to evaluate the differential expression profiles of lncRNAs and mRNAs in goat testicular tissues to explore the molecular mechanisms of ECM formation and spermatogenesis. Our purpose was to both enrich our understanding of the goat genome and also to provide a basis for further research on the function of candidate genes. We used RNA-seq to research lncRNA and mRNA expression levels in the testicular tissue of goats aged 3, 6, and 12 months; these ages represent the three growth stages of prepubertal, pubertal, and post-pubertal, respectively. In addition, we further studied the protein expressions of MMP2, IGF2, and IGFBP6 by western blotting. Our research may provide new insights into the molecular mechanisms of testicular development in the goat and provide a basis for further exploration of goat breeding, spermatogenesis, and marker screening, among other roles.


Ethics statement

All experiments involving goats were conducted in strict compliance with the relevant guidelines set by the Ethics Committee of Guizhou Animal Husbandry and Veterinary Research Institute (JXYJS-20190312). All tests were conducted in accordance with the relevant guidelines and regulations formulated by the Ministry of Agriculture, People’s Republic of China.

Sample collection and RNA isolation

The 3, 6, and 12-month-old male Guizhou Black Goats were obtained from a goat farm in Weining County (Guizhou, China). All experimental animals were euthanized; the anesthesia program gave 2.5% thiopental sodium general anesthesia at a dose of 10 mg/kg, followed by electrocution. Four healthy goats were selected for each age stage. Tissue samples for the extraction of RNA and protein were removed from testicular tissue and immediately frozen in liquid nitrogen, transported to the laboratory, and stored in the refrigerator at −80°C.
Each frozen tissue sample was ground into a fine powder in liquid nitrogen, and the total RNA was extracted with Trizol reagent according to the manufacturer’s instructions. Then the concentration, purity, and integrity of total RNA were detected. RNA concentration was measured using a Qubit RNA Assay Kit in a Qubit 2.0 Flurometer (Life Technologies, Carlsbad, CA, USA). RNA purity was checked using a NanoPhotometer spectrophotometer (IMPLEN, Palo Alto, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the bioanalyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA).

Library preparation and RNA sequencing

After qualifying, rRNA removal, library preparation, and examination of library quality was performed. After passing the quality inspection, computer sequencing was undertaken. The raw reads obtained by sequencing were sequentially removed from the reads with adapter; thereafter, clean reads were obtained after those reads with a proportion of undefinable base information more significant than 10% and low-quality reads.

Identification and prediction of differentially expressed genes

The reference genome index (https://www.ncbi.nlm.nih.gov/genome/?term=Capra+aegagrus+hircus) was built using Bowtie (v2.0.6) [17], and paired-end clean reads were aligned to the reference genome using TopHat (v2.0.9) [18]. Basic screening mainly included the selection of transcripts with length ≥200 bp and exon number ≥2; the reads of each transcript with minimum coverage ≥3 were calculated by StringTie; similar or identical transcripts were screened by comparing them with known non-lncRNAs and non-mRNAs using Cuffcompare in goats; by comparing with known mRNAs and using class_code information in Cuffcompare results, candidate transcripts of different types were screened (lincRNA, intronic lncRNA, anti-sense lncRNA, and sense-overlapping lncRNA). We applied four approaches to the analysis of Coding potential: the Coding-Non-Coding-Index (CNCI) [19]; the Coding Potential Calculator (CPC) [20]; Pfam [21], and the Coding Potential Assessment Tool (CPAT) [22].

Target gene prediction and enrichment analysis

The target genes of lncRNAs were predicted by cis- or trans-action. Cis-action screened out the protein coding genes adjacent to lncRNA (upstream and downstream 100 K) as its target genes. Trans-action target genes were predicted by selecting the Pearson correlation coefficient (PCC) of lncRNAs and protein-coding genes between samples >0.6, and comparing the base complementary coordination relationship between lncRNAs and mRNAs.
The FPKMs (Reads Per Kilobase of exon model per Million mapped reads) of lncRNA and coding genes in each sample were calculated using Cuffdiff (v2.1.1) [23]. The FPKM value was used to estimate the expression levels of genes. Gene ontogeny (GO) enrichment analysis of differentially expressed genes (DEG) was conducted using the goseq R package [24]. Kyoto encyclopedia of genes and genomes (KEGG) pathways were detected using KOBAS (KEGG Orthology Based Annotation System) [25]. Corrected p-values <0.05 were considered significantly enriched by DEGs.

PPI network and lncRNA–mRNA co-expression network analysis

The protein–protein interaction (PPI) network based on DEGs was significantly enriched in GO terms between 3 vs 6 months of age. The PPI network was obtained through an online analysis tool (https://www.string-db.org/). By calculating the expression correlation between lncRNAs and genes encoded by custom scripts, and using WGCNA to cluster genes from different samples, standard expression modules were found [26]. Subsequently, their functions were analyzed using functional enrichment analysis. The co-expression of lncRNAs-mRNAs was analyzed by calculating PCC between the encoding gene and the specific expression level of lncRNAs. The differentially expressed lncRNAs were correlated with the predicted cis- or trans-acting target mRNAs. The visualization of gene interaction was realized using Cytoscape software (v3.3.0) (The Cytoscape Consortium, USA).

Quantitative reverse-transcription polymerase chain reaction validation

We identified 11 differentially expressed mRNAs and five differentially expressed lncRNAs. Primers were designed and detected using the NCBI Pick Primers and Primer-BLAST tools, respectively. Primer information is listed in Table 1. The 2−ΔΔCt quantitative method was used to calculate the relative expression level of mRNAs and lncRNAs [27], and normalization took place using the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH). The total RNA of goat testicular tissue was isolated by TRIzol (Invitrogen Life Technologies, Carlsbad, CA, USA). cDNAs were synthesized from RNA using a PrimeScript RT Kit with gDNA Eraser (Takara, Dalian, China). Quantitative polymerase chain reaction (qPCR) was performed with SYBR Green Master mix (Roche Applied Science, Mannheim, Germany). The 20 μL reaction solution contained 9 μL SYBR, 2 μL cDNA, 1 μL each of the forward and reverse primers, and 7 μL ddH2O. The reaction system consisted of a holding stage of 5 min at 50°C, and 10 min at 95°C, 40 cycles of 95°C for 15 s, 60°C for 30 s, and 72°C for 30 s; and a final stage at 95°C for 15 s, 60°C for 30 s, and 95°C for 15 s.

Western blotting analyses

The testicular tissue samples were washed three times with phosphate buffered saline (Servicebio, China), then cut into small pieces and place in an equalizer with 12,000 g of homogenate for 10 min. The supernatant was collected, and the protein concentration was measured using a BCA protein concentration assay kit according to the manufacturer’s instructions (Servicebio, China) and normalized to GAPDH levels. Transfer condition: 200 mA transfer for 1 h. Primary antibodies against IGF2, IGFBP6, and MMP2 (diluted 1:500, 1:1,000 and 1:1,000; cat. no: ab170304, ab109765, and ab97779, respectively [Abcam, Cambridge, UK]), and GAPDH (diluted 1:30,000; cat. no: GB12002, Proteintech Group, Inc., Wuhan, China) were incubated at 4°C for 3 h. Subsequently, the goat anti-rabbit horseradish peroxidase immunoglobulin G secondary antibody (diluted 1:5,000; Servicebio, cat. no: GB23303, China) was incubated at room temperature for 30 min. Images were acquired using a Perform Electrochemiluminescence kit (Servicebio, China) following the manufacturer’s instructions. For image analysis, films were scanned, organized, and desaturated using Adobe Photoshop (Adobe, Santa Clara, CA, USA). The optical density of the target band was analyzed using the Alpha processing system (Alpha Innotech, Shanghai, China).

Statistical analysis

The original data was organized in Microsoft Excel 2010. The one-way analysis of variance (age factor) analytical method was performed to compare means using SPSS 21.0 software (SPSS, Chicago, USA). All results were expressed as mean±standard deviation. Similar letters indicate no significant difference (p>0.05), while different letters indicate a significant difference (p<0.05).


Testis transcriptome characterization and lncRNA identification

We constructed 12 cDNA libraries using total testicular RNA from four individuals from each of the three specific stages of growth and development in goats. Each sample contained approximately 7.52 Gb of high-quality sequence data. The data has been deposited into the NCBI database (NCBI Accession No: PRJNA613301). Pairwise comparison analysis of different groups showed that 161, 12, and 665 differentially expressed mRNAs were present in the testicular tissues of goats between 3 vs 6, 6 vs 12, and 3 vs 12 months, respectively. In addition, there was a common differential gene, MMP2, among the three groups of differential genes. In total, we intersected the results of CPC/CNCI/PFAM/CPAT and obtained lncRNAs for subsequent analysis. Results including 26,380 lincRNAs (76.64%), 3,996 anti-sense lncRNAs (11.61%), 3,235 intronic lncRNAs (9.40%), and 808 sense lncRNAs (2.35%) (Figure 1).

Differentially expressed gene analysis

A total of 161 (105 up, 56 down), 12 (9 up, 3 down), and 665 (512 up, 153 down) differentially expressed mRNAs were identified. Furthermore, 23 (5 up, 18 down), 8 (5 up, 3 down), and 135 (33 up, 102 down) differentially expressed lncRNAs were present in 3 vs 6, 6 vs 12, and 3 vs 12-month-old goat testes, respectively. The differentially expressed mRNAs and lncRNAs were clustered together owing to their similar expression profiles, which are shown in Figure 2A and Figure 2B as heat maps, respectively.

Functional enrichment analysis

According to GO analysis results of DEG, a total of 21 GO terms were significantly enriched between 3 vs 6 months of age (Figure 3A). In comparison, there was no significant enrichment in GO terms between 6 vs 12 months of age (Figure 3B). Based on GO analysis of biological processes, we found that these DEGs were significantly enriched in three GO terms (reactive oxygen species metabolic process, superoxide metabolic process, and transcription initiation DNA-dependent). In addition, according to the GO analysis result of molecular function, the transcription initiation factor activity and ECM structural constituent were also significantly enriched. Moreover, the level of significant enrichment in cell composition was the largest, with 16 GO terms, including transcription factor complex, DNA-directed RNA polymerase II, membrane-enclosed lumen, organelle lumen, intracellular organelle lumen, ECM, proteinaceous ECM, and extracellular region (Figure 3A).
The significantly enriched KEGG pathway analysis of mRNA revealed that many pathways were related to spermatogenesis; as an example, the ECM-receptor interaction signaling pathway was present. KEGG pathway analysis of lncRNA cis- and trans-target genes revealed that many significant enrichment pathways were consistent with the results of mRNA enrichment (Figure 4); including ECM-receptor interaction signaling pathway, focal adhesion, and ribosome.

PPI network and Co-expression network construction

PPI network analysis of DEGs from significantly enriched GO terms. There were 38 edges and 16 nodes in the network (Figure 5). According to the results of the PPI network analysis, MMP2 was the hub gene (degree = 29). In addition, among these DEGs, most were related to the ECM (Figure 5). PPI network analysis of DEGs was significantly enriched in GO terms. Line color represented the different types of interaction evidence.
According to the co-expression network, TCONS-0002139 and TCONS-00093342 were co-expressed with MMP2 and collagen type IV alpha 4 chain (COL4A4). In addition, lncRNA TCONS_00093342 was co-expressed with IGFBP6, and lncRNA TCONS-0002139 was also co-expressed with COL4A1 (Figure 6).

Quantitative polymerase chain reaction validation

To test the reliability of our RNA-seq sequencing results, 11 differentially expressed mRNAs and five differentially expressed lncRNAs were used in validation through quantitative reverse-transcription PCR (qRT-PCR). The relative fold changes in the qRT-PCR assay showed that the qRT-PCR verification results were consistent with the RNA-seq sequencing results (Tables 2, 3).

Western blotting validation

Western blotting results showed that the protein expression trend of IGF2, IGFBP6, and MMP2 was consistent with the mRNA expression trend, with the highest expression levels at 3 months of age followed by a gradual decrease with increasing age. In addition, IGF2 and MMP2 were significantly different between 3 vs 6 months of age (p<0.05), with no significant difference between 6 vs 12 months of age (p>0.05) (Figure 7).


lncRNAs, as new regulatory molecules in cell development, have attracted wide attention [28]. One reason is that lncRNAs may help regulate gene expression at transcriptional and post-transcriptional levels through genetic and epigenetic mechanisms [29]. In addition, they play a vital role in male fertility and spermatogenesis [30]. In recent years, reports employing lncRNA-seq technology have investigated the testicular tissue of ruminants at different developmental stages. For example, Yang et al [31] report a comparative analysis and study on lncRNAs and mRNAs during the maturation of sheep testis, which provides a valuable resource for further research regarding the functions of lncRNAs in ovine testicular development. Gao et al [32] conducted a comparative analysis of mRNAs and lncRNAs in prepubertal and post-pubertal bovine testicular tissues, which provides new information for further study of the biological functions of bovine lncRNAs.
There are also related studies regarding the expression of lncRNAs in different developmental stages of goat testis. Bo et al [33] employed RNA-seq in their study of the lncRNA of testicular tissues in prepubertal and pubertal goats. Their findings indicate that lncRNAs regulate different modes of spermatogenesis and testicular growth, and provide a new perspective for the analysis of lncRNA expression and age-related changes in goat testis [33]. However, the current study was limited to prepubertal and pubertal testicular tissue in goats, and lacked post-pubertal research. Therefore, we analyzed mRNA and lncRNA expression in the testicular tissues of goats aged 3, 6, and 12 months. We found the number of differentially expressed mRNAs was greater than that of differentially expressed lncRNAs. In addition, the number of DEGs between prepubertal and pubertal was higher than those between post-pubertal and pubertal. These findings were consistent with the findings in sheep testicles [34]. This phenomenon suggests that goat testicular tissue has a more complex gene regulation mechanism prior to puberty.
To further study the biological processes involved in testicular development and spermatogenesis, we functionally classified the differentially expressed mRNAs in the testicular tissues of goats at different developmental stages. The primary function of DEGs was elucidated by GO analysis. According to the GO term results, there was an evident and significant enrichment phenomenon between 3 vs 6 months of age. In addition, many of these GO terms were associated with the ECM, especially in terms of cellular components (Figure 3). Subsequently, PPI analysis was performed based on the significantly enriched GO terms at the age of 3 vs 6 months. Analysis results showed that these genes were mainly related to the ECM, including ECM structural constituent, ECM organization, and ECM-receptor interaction. In addition, according to PPI network, MMP2 was the hub gene (Figure 4). MMP2 was not only associated with the ECM, but also with the extracellular region (GO: 0005576). It is known that the ECM plays an essential role in forming the seminiferous tubule wall [35]. This phenomenon indicated that the age of 3 to 6 months was an essential period during formation of the extracellular region of goat testicular tissue, and that it is a critical stage in seminiferous tubule formation.
In the 3 vs 6 months GO terms, the number of DEGs in the extracellular region was the largest (Figure 3A). These genes included MMP2, IGF2, IGFBP6, and COL4A1. It is known that the testicular ECM is essential for the movement of germ cells through the blood-testis-barrier (BTB) during spermatogenesis; it is also known that proteins in the ECM modulate BTB dynamics through cytokines [36]. As a cytokine, IGF2 plays an essential regulatory role in the ECM [37]. In addition, many GO terms related to spermatogenesis were also enriched, including regulation of cell growth, IGF binding, cell growth, growth factor binding, and regulation of growth (Figure 3). These results further suggest that these genes regulate testicular development and spermatogenesis in goats through the ECM.
KEGG analysis of lncRNAs targets genes and DEGs revealed many pathways related to the ECM and spermatogenesis (Figure 5). These included focal adhesion, cell adhesion molecules, and ECM-receptor interaction. For example, the adhesion complex acts as an anchor junction between the germ cell epithelium and tissue ECM, regulating tight junctions between the Sertoli cells and germ cells as sperm cells pass through the germinal epithelium [36]. The association of lncRNAs with functionally annotated mRNAs can be achieved by co-expression networks [38]. We constructed a co-expression network of lncRNAs and their target genes to explore the potential regulatory mechanism of lncRNAs-mRNAs. The results of co-expression analysis showed that there were 11 mRNAs both co-expressing with lncRNA Tcons_00021399 and Tcons_0009334, including MMP2 (Figure 6). In addition, COL4A1 was co-expressed with Tcons_0021399, and IGFBP6 was co-expressed with Tcons_0009334. According to previous studies, these genes are particularly associated with the ECM in spermatogenesis; for example, MMP2 is indispensable in the degradation and remodeling of the ECM [39] and COL4A1 is closely associated with ECM recombination [40].
The remodeling process depends on the interaction between cells and the ECM [41]. IGFBP6 is involved in the IGF regulatory pathway, which regulates the proliferation and differentiation of undifferentiated spermatogonia [42]. The IGF system is involved in various cellular biological functions, inducing the production of the ECM [43]. In addition, the expression of these genes is regulated by lncRNAs. For example, LINC01128 plays a regulatory role in cell invasion, migration, and proliferation through the LINC01128/miR-299-3p/MMP2 axis [44]. Based on the above analysis, we speculated that the lncRNAs Tcons_00021399 and Tcons_0009334 may play an essential role in regulating molecular mechanisms related to ECM recombination or spermatogenesis in goat testis.
According to our results, the protein expression levels of IGF2, IGFBP6, and MMP2 were highest at 3 months of age, and then decreased with age (Table 3; Figure 8). IGF2 and MMP2 were significantly different between 3 vs 6 months of age (p<0.05). These results are similar to those found in the development of testicular tissue in mice [45]. In addition, IGF2 expression levels in the liver, kidney, and heart of male mice are significantly down-regulated with age [46], which was consistent with our findings in goat testicular tissues. These results further demonstrate the reliability of our transcriptome analysis from a protein perspective. However, the molecular mechanism between these lncRNAs and mRNAs needs further study.


We studied mRNA and lncRNA expression levels using RNA-seq in the testicular tissues of goats at 3, 6, and 12 months of age. Several differentially expressed lncRNAs and mRNAs associated with the ECM and spermatogenesis were obtained. This study provides a high-quality resource for future goat transcriptome studies, especially for Guizhou Black goats. In addition, this study improves the comparative understanding of the molecular mechanisms of testicular development during goat sexual maturation.



We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.


This study was supported by the Science and Technology Program of Guizhou Province ([2020]4009(002)); China Agriculture Research System of MOF and MARA; Team Cultivation Project of Guizhou Animal Husbandry and Veterinary Research Institute (No. 03[2018]).

Figure 1
Venn diagram and pie chart of differentially expressed genes. (A) Pie chart of differentially expressed mRNAs. (B) Venn diagram according to four tools used for coding protein analysis. (C) Statistical analysis of different lncRNAs. 3M, 3 months old; 6M, 6 months old; 12M, 12 months old.
Figure 2
Hierarchical heat map of a differential expression cluster graph. (A) Clustering heat map of differentially expressed mRNAs. (B) Clustering heat map of differentially expressed lncRNAs. Red indicates highly expressed genes and green indicates lowly expressed genes. 3M, 3 months old; 6M, 6 months old; 12M, 12 months old.
Figure 3
The top 30 GO terms listed in the GO analysis of differentially expressed mRNAs. (A) 3 vs 6 months of age. (B) 6 vs 12 months of age. * Indicates significant enrichment (p<0.05). GO, gene ontology.
Figure 4
KEGG analysis of 3 vs 12 month old differentially expressed genes and lncRNA target genes. (A) Differentially expressed genes. (B) Differentially expressed lncRNA cis target genes. (C) Differentially expressed lncRNA trans target genes; KEGG, Kyoto encyclopedia of genes and genomes.
Figure 5
PPI network analysis of differentially expressed genes significantly enriched in GO terms. Line color represents different types of interaction evidence. PPI, protein–protein interaction; GO, gene ontology.
Figure 6
The lncRNA-mRNA co-expression network. The size of a node in the interaction network is proportional to the degree of the node. The thickness of inter-node lines is directly proportional to the value of Pearson’s correlation coefficient. Thicker lines indicate higher values.
Figure 7
Relative abundance of IGF2, IGFBP6, and MMP2 protein in goat testes. Each sample was assessed three times. The results are represented by means±standard deviation. GAPDH was used as a loading control. 3M, 3 months old; 6M, 6 months old; 12M, 12 months old; IGF2, insulin-like growth factor 2; IGFBP6, insulin-like growth factor-binding protein 6; MMP2, matrix metalloproteinase 2; GAPDH, glyceraldehyde-3-phosphate dehydrogenase. A,B Different letters indicate significant differences (p<0.05); similar letters indicate that the difference is not significant (p>0.05).
Table 1
List of the primers used in quantitative reverse-transcription polymerase chain reaction
NCBI Accession No Gene Primer sequence (5′–3′) Tm Product length
60 270
60 184
60 97
60 180
60 179
58 127
60 102
59 244
58 194
58 170
59 155
59 209
60 113
59 268
59 122
61 231
61 131
Table 2
List of lncRNA-seq and lncRNA-qPCR results
Item Gene 3M 6M 12M
lncRNA-seq TCONS_00020249 16.403±0.903A 10.934±0.959B 9.692±1.077B
TCONS_00021399 13.982±1.176A 7.393±1.923B 5.979±1.298B
TCONS_00093342 12.928±1.289A 6.098±1.183B 4.329±0.168B
TCONS_00124262 5.290±0.886A 1.301±0.379B 2.401±0.893AB
TCONS_00070285 11.145±1.282A 5.203±0.700B 5.595±1.067B
lncRNA-qPCR TCONS_00020249 1.003±0.015A 0.596±0.002B 0.599±0.028B
TCONS_00021399 1.001±0.003A 0.844±0.002B 0.492±0.002C
TCONS_00093342 1.002±0.011A 0.300±0.048B 0.327±0.002B
TCONS_00124262 0.998±0.007A 0.308±0.003C 0.501±0.004B
TCONS_00070285 1.001±0.057A 0.667±0.002B 0.403±0.007C

qPCR, quantitative polymerase chain reaction; 3M, 3 months old; 6M, 6 months old; 12M, 12 months old.

A–C Different letters in the same row indicate significant differences (p<0.05).

Table 3
List of mRNA-seq and mRNA-qPCR results
Item Gene 3M 6M 12M
mRNA-seq SPARCL1 90.890±4.409A 52.901±3.092B 44.064±3.533B
CRISP2 442.548±40.613B 646.077±30.079A 674.702±33.187A
HSD3B7 7.596±0.656A 4.008±0.364B 3.592±0.390B
PMP22 7.644±0.660A 3.664±0.602B 2.624±0.371B
COL4A1 22.911±2.870A 10.757±1.118B 8.111±1.363B
COL4A4 69.471±7.089A 27.079±1.636B 20.927±2.475B
HHIPL1 6.534±0.361A 3.731±0.410B 2.822±0.382B
AMH 8.897±0.975A 5.167±0.124B 4.944±0.488B
MMP2 9.052±0.875A 4.427±0.422B 2.602±0.206B
IGFBP6 130.894±9.796A 77.288±7.249B 69.930±9.074B
IGF2 7.918±0.921A 3.719±0.571B 4.114±0.481B
mRNA-qPCR SPARCL1 0.997±0.008A 0.376±0.011B 0.362±0.015B
CRISP2 0.997±0.009B 1.270±0.009A 0.485±0.006C
HSD3B7 1.002±0.009A 0.118±0.005C 0.148±0.004B
PMP22 1.000±0.012A 0.491±0.026B 0.350±0.002C
COL4A1 1.003±0.015A 0.031±0.001B 0.041±0.001B
COL4A4 1.003±0.055A 0.596±0.024B 0.576±0.010B
HHIPL1 1.003±0.010A 0.031±0.001B 0.046±0.002B
AMH 1.002±0.008A 0.020±0.000B 0.011±0.000B
MMP2 1.030±0.081A 0.520±0.317A 0.086±0.002B
IGFBP6 1.001±0.047A 0.950±0.143A 0.486±0.017A
IGF2 1.021 ±0.069A 0.519±0.024B 0.073±0.002C

qPCR, quantitative polymerase chain reaction; 3M, 3 months old; 6M, 6 months old; 12M, 12 months old.

A–C Different letters in the same row indicate significant differences (p<0.05).


1. Zhu Z, Li C, Yang S, et al. Dynamics of the transcriptome during human spermatogenesis: predicting the potential key genes regulating male gametes generation. Sci Rep 2016; 6:19069 https://doi.org/10.1038/srep19069
crossref pmid pmc
2. Li C, Shen C, Shang X, et al. Two novel testis-specific long noncoding RNAs produced by 1700121C10Rik are dispensable for male fertility in mice. J Reprod Dev 2020; 66:57–65. https://doi.org/10.1262/jrd.2019-104
crossref pmid
3. Zeng F, Chen Y, Guo C, et al. Analysis of differentially abundant proteins related to boar fertility in seminal plasma using iTRAQ-based quantitative proteomics. J Proteomics 2021; 236:104120 https://doi.org/10.1016/j.jprot.2021.104120
crossref pmid
4. Cole JB, Dürr JW, Nicolazzi EL. Invited review: The future of selection decisions and breeding programs: What are we breeding for, and who decides? J Dairy Sci 2021; 104:5111–24. https://doi.org/10.3168/jds.2020-19777
crossref pmid
5. Mayerhofer A, Walenta L, Mayer C, Eubler K, Welter H. Human testicular peritubular cells, mast cells and testicular inflammation. Andrologia 2018; 50:e13055 https://doi.org/10.1111/and.13055
crossref pmid
6. Nosrati R, Kheirouri S, Ghodsi R, Ojaghi H. The effects of zinc treatment on matrix metalloproteinases: a systematic review. J Trace Elem Med Biol 2019; 56:107–15. https://doi.org/10.1016/j.jtemb.2019.08.001
crossref pmid
7. Matuszczak E, Komarowska M, Sankiewicz A, et al. Plasma concentration of MMP-1 and MMP-2 in boys with cryptorchidism and its lack of correlation with INSL3 and inhibin B. Scand J Clin Lab Invest 2019; 79:412–8. https://doi.org/10.1080/00365513.2019.1637534
crossref pmid
8. Wathes D, Cheng Z, Fenwick MA, Fitzpatrick R, Patton J. Influence of energy balance on the somatotrophic axis and matrix metalloproteinase expression in the endometrium of the postpartum dairy cow. Reproduction (Cambridge, England) 2011; 141:269–81. https://doi.org/10.1530/rep-10-0177
crossref pmid pmc
9. Voit-Ostricki L, Lovas S, Watts CR. Conformation and domain movement analysis of human matrix metalloproteinase-2: role of associated Zn2+ and Ca2+ ions. Int J Mol Sci 2019; 20:4194 https://doi.org/10.3390/ijms20174194
crossref pmc
10. Clemmons DR. Role of IGF-binding proteins in regulating IGF responses to changes in metabolism. J Mol Endocrinol 2018; 61:T139–T69. https://doi.org/10.1530/jme-18-0016
crossref pmid
11. Batista P, Chang H. Long noncoding RNAs: cellular address codes in development and disease. Cell 2013; 152:1298–307. https://doi.org/10.1016/j.cell.2013.02.012
crossref pmid pmc
12. Zhang C, Gao L, Xu EY. LncRNA, a new component of expanding RNA-protein regulatory network important for animal sperm development. Semin Cell Dev Biol 2016; 59:110–7. https://doi.org/10.1016/j.semcdb.2016.06.013
crossref pmid
13. Liu K, Li T, Ton H, Mao X, Chen Y. Advances of long noncoding RNAs-mediated regulation in reproduction. Chin Med J 2018; 131:226–34. https://doi.org/10.4103/0366-6999.222337
crossref pmid pmc
14. Liu Y, Sun Y, Li Y, et al. Analyses of long non-coding RNA and mRNA profiling using RNA sequencing in chicken testis with extreme sperm motility. Sci Rep 2017; 7:9055 https://doi.org/10.1038/s41598-017-08738-9
crossref pmid pmc
15. Zhang Y, Yang H, Han L, et al. Long noncoding RNA expression profile changes associated with dietary energy in the sheep testis during sexual maturation. Sci Rep 2017; 7:5180 https://doi.org/10.1038/s41598-017-05443-5
crossref pmid pmc
16. Weng B, Ran M, Chen B, He C, Dong L, Peng F. Genome-wide analysis of long non-coding RNAs and their role in postnatal porcine testis development. Genomics 2017; 109:446–56. https://doi.org/10.1016/j.ygeno.2017.07.001
crossref pmid
17. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods 2012; 9:357–9. https://doi.org/10.1038/nmeth.1923
crossref pmid pmc
18. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 2013; 14:R36 https://doi.org/10.1186/gb-2013-14-4-r36
crossref pmid pmc
19. Sun L, Luo H, Bu D, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res 2013; 41:e166 https://doi.org/10.1093/nar/gkt646
crossref pmid pmc
20. Kong L, Zhang Y, Ye Z, et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res 2007; 35:W345–9. https://doi.org/10.1093/nar/gkm391
crossref pmid pmc
21. El-Gebali S, Mistry J, Bateman A, et al. The Pfam protein families database in 2019. Nucleic Acids Res 2019; 47:D427–32. https://doi.org/10.1093/nar/gky995
crossref pmid
22. Wang L, Park H, Dasari S, Wang S, Kocher JP, Li W. CPAT: Coding-potential assessment tool using an alignment-free logistic regression model. Nucleic Acids Res 2013; 41:e74 https://doi.org/10.1093/nar/gkt006
crossref pmid pmc
23. Trapnell C, Williams B, Pertea G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 2010; 28:511–5. https://doi.org/10.1038/nbt.1621
crossref pmid pmc
24. Young M, Wakefield M, Smyth G, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol 2010; 11:R14 https://doi.org/10.1186/gb-2010-11-2-r14
crossref pmid pmc
25. Xie C, Mao X, Huang J, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res 2011; 39:W316–22. https://doi.org/10.1093/nar/gkr483
crossref pmid pmc
26. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 2008; 9:559 https://doi.org/10.1186/1471-2105-9-559
crossref pmid pmc
27. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 ΔΔ C T method. Methods 2001; 25:402–8. https://doi.org/10.1006/meth.2001.1262
crossref pmid
28. Veneziano D, Nigita G, Ferro A. Computational approaches for the analysis of ncRNA through deep sequencing techniques. Front Bioeng Biotechnol 2015; 3:77 https://doi.org/10.3389/fbioe.2015.00077
crossref pmid pmc
29. Bao J, Wu J, Schuster A, Hennig G, Yan W. Expression profiling reveals developmentally regulated lncRNA repertoire in the mouse male germline. Biol Reprod 2013; 89:107 https://doi.org/10.1095/biolreprod.113.113308
crossref pmid pmc
30. Joshi M, Rajender S. Long non-coding RNAs (lncRNAs) in spermatogenesis and male infertility. Reprod Biol Endocrinol 2020; 18:103 https://doi.org/10.1186/s12958-020-00660-6
crossref pmid pmc
31. Yang H, Wang F, Li F, et al. Comprehensive analysis of long noncoding RNA and mRNA expression patterns in sheep testicular maturation. Biol Reprod 2018; 99:650–61. https://doi.org/10.1093/biolre/ioy088
crossref pmid
32. Gao Y, Li S, Lai Z, et al. Analysis of long non-coding RNA and mRNA expression profiling in immature and mature bovine (Bos taurus) testes. Front Genet 2019; 10:646 https://doi.org/10.3389/fgene.2019.00646
crossref pmid pmc
33. Bo D, Jiang X, Liu G, Hu R, Chong Y. RNA-Seq implies divergent regulation patterns of lincRNA on spermatogenesis and testis growth in goats. Animals 2021; 11:625 https://doi.org/10.3390/ani11030625
crossref pmid pmc
34. Bai M, Sun L, Zhao J, et al. Histological analysis and identification of spermatogenesis-related genes in 2-, 6-, and 12- month-old sheep testes. Naturwissenschaften 2017; 104:84 https://doi.org/10.1007/s00114-017-1505-1
crossref pmid
35. Mayerhofer A. Human testicular peritubular cells: more than meets the eye. Reproduction (Cambridge, England) 2013; 145:R107–16. https://doi.org/10.1530/rep-12-0497
crossref pmid
36. Siu MKY, Cheng CY. Extracellular matrix: recent advances on its role in junction dynamics in the seminiferous epithelium during spermatogenesis. Biol Reprod 2004; 71:375–91. https://doi.org/10.1095/biolreprod.104.028225
crossref pmid
37. Tan KS, Kulkeaw K, Nakanishi Y, Sugiyama D. Expression of cytokine and extracellular matrix mRNAs in fetal hepatic stellate cells. Genes Cells 2017; 22:836–44. https://doi.org/10.1111/gtc.12517
crossref pmid
38. van Dam S, Võsa U, van der Graaf A, Franke L, de Magalhães JP. Gene co-expression analysis for functional classification and gene-disease predictions. Brief Bioinform 2018; 19:575–92. https://doi.org/10.1093/bib/bbw139
crossref pmid
39. Huang H, Zhao W, Tang Z, et al. Characterization of porcine MMP-2 and its association with immune traits. Gene 2009; 435:63–71. https://doi.org/10.1016/j.gene.2009.01.002
crossref pmid
40. Schrade A, Kyrönlahti A, Akinrinade O, et al. GATA4 regulates blood-testis barrier function and lactate metabolism in mouse sertoli cells. Endocrinology 2016; 157:2416–31. https://doi.org/10.1210/en.2015-1927
crossref pmid pmc
41. McMillen P, Holley SA. Integration of cell-cell and cell-ECM adhesion in vertebrate morphogenesis. Curr Opin Cell Biol 2015; 36:48–53. https://doi.org/10.1016/j.ceb.2015.07.002
crossref pmid pmc
42. Safian D, Morais R, Bogerd J, Schulz R. Igf binding proteins protect undifferentiated spermatogonia in the zebrafish testis against excessive differentiation. Endocrinology 2016; 157:4423–33. https://doi.org/10.1210/en.2016-1315
crossref pmid
43. Su YY, Nishimoto T, Hoffman S, et al. Insulin-like growth factor binding protein-4 exerts antifibrotic activity by reducing levels of connective tissue growth factor and the C-X-C chemokine receptor 4. FASEB bioAdvances 2019; 1:167–79. https://doi.org/10.1096/fba.2018-00015
crossref pmid pmc
44. Yao Q, Chen T. LINC01128 regulates the development of osteosarcoma by sponging miR-299-3p to mediate MMP2 expression and activating Wnt/β-catenin signalling pathway. J Cell Mol Med 2020; 24:14293–305. https://doi.org/10.1111/jcmm.16046
crossref pmid pmc
45. Chen H, Fok K, Yu S, et al. CD147 is required for matrix metalloproteinases-2 production and germ cell migration during spermatogenesis. Mol Hum Reprod 2011; 17:405–14. https://doi.org/10.1093/molehr/gar013
crossref pmid
46. Finkielstain GP, Forcinito P, Lui JCK, et al. An extensive genetic program occurring during postnatal growth in multiple tissues. Endocrinology 2009; 150:1791–800. https://doi.org/10.1210/en.2008-0868
crossref pmid
METRICS Graph View
  • 1 Crossref
  • 1 Scopus
  • 1,838 View
  • 137 Download
Related articles

Editorial Office
Asian-Australasian Association of Animal Production Societies(AAAP)
Room 708 Sammo Sporex, 23, Sillim-ro 59-gil, Gwanak-gu, Seoul 08776, Korea   
TEL : +82-2-888-6558    FAX : +82-2-888-6559   
E-mail : animbiosci@gmail.com               

Copyright © 2022 by Animal Bioscience.

Developed in M2PI

Close layer
prev next