Male tissues


Figure 1. Schematic drawing of the male reproductive tissues of pig.

In the pig transcriptomic analysis, testis, epididymis, ductus deferens, penis and male glands (seminal vesicle, prostate, bulbourethral gland and urethral gland) were sampled and used for quantitative RNA measurements, representing male reproductive tissues. Sperm cells are produced in the testis, stored and matured in the epididymis, transported further via the ductus deferens and finally exit the body via the penis after the male glands have added various secreted products. The gene classification strategy highlights genes with an elevated level of expression in one or a group of tissues compared to all other tissues.

2718 genes are classified as testis elevated out of which 1004 genes are highlighted as testis enriched. Based on the expression in epididymis 685 genes are classified as elevated compared to other tissue types, out of which 115 genes are highlighted as tissue enriched in epididymis. RNA expression levels in ductus deferens resulted in 152 genes classified as elevated in ductus deferens compared to other tissue types, out of which 5 genes are highlighted as tissue enriched in ductus deferens. 117 genes are classified as penis elevated out of which 0 genes are highlighted as penis enriched. Based on the expression in male glands 1018 genes are classified as elevated compared to other tissue types, out of which 51 genes are highlighted as tissue enriched in male glands.

Histological image of the pig tissues used in the analysis can be found in the pig tissue dictionary.

Testis

The testis is essential for reproduction, responsible for producing sperm cells through spermatogenesis. Testis also produces hormones important for the development of male sex characteristics. The two major parts of the testis is the seminiferous ducts, where spermatogenesis takes place, and the Leydig cells located outside the seminiferous ducts producing hormones.

Gene expression in testis is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 2 summarizes the number of genes in respective category. In total, 16867 genes are detected above cut of (1NX) in pig testis. The tissue distribution category highlights 454 genes only detected in testis while 2718 genes are classified as testis elevated compared to other tissues. Table 1 shows the overlap for the testis elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 2, (A) The distribution of all genes across the five categories based on transcript specificity in testis as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in testis as well as in all other tissues. The combination of the two categories is shown in table 1.

Table 1, Number of genes in the subdivided categories of elevated expression and tissue distribution in testis

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 39741212768 1004
Group enriched 03104515 370
Tissue enhanced 57621386280 1344
Total 4541343558363 2718

Testis enriched expression

Cells found in the testis are cells not found elsewhere in the body, and therefore testis is the tissue type with the most uniquely expressed genes (expression only detected in testis and not detected above cut of in other tissues types). Among the genes classified as testis enriched there is a very good overlap to the corresponding human analysis of testis expression. TSSK6, PRSS37 and CRISP2 are all highlighted as testis enriched in pig as well as human.

Table 2, The 10 genes with the highest level of enriched expression in testis. "mRNA (tissue)" shows the transcript level in testis as NX values. "Tissue specificity score (TS)" corresponds to the fold-change between the expression level in testis and the tissue with second highest expression level.

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000017413 HSPB9 heat shock protein family B (small) member 9 Detected in single 169.4 1694
ENSSSCG00000013937 TSSK6 testis-specific serine/threonine-protein kinase 6 Detected in single 84.2 843
ENSSSCG00000005121 na Detected in single 77.9 780
ENSSSCG00000029612 PRSS37 serine protease 37 Detected in single 138.3 701
ENSSSCG00000034728 RNF151 ring finger protein 151 Detected in single 66.9 670
ENSSSCG00000001736 CRISP2 cysteine-rich secretory protein 2 precursor Detected in single 267.5 664
ENSSSCG00000030671 DAZL deleted in azoospermia like Detected in single 67.6 601
ENSSSCG00000039238 CABS1 calcium binding protein, spermatid associated 1 Detected in single 107.2 597
ENSSSCG00000004681 TRIM69 tripartite motif containing 69 Detected in single 55.7 557
ENSSSCG00000031522 ASB17 ankyrin repeat and SOCS box containing 17 Detected in single 55.4 555



Immunohistochemical labeling of testis, targeting TNP1 using HPA044387 and CYP19A1 using HPA051194.

A notable differences between species is that the pig testis includes a larger number of Leydig cells compared to human testis, therefore resulting in several genes expressed by Leydig cells classified as testis elevated in pig but not human, for example genes that in human are classified as adrenal gland enriched and low expression in testis shows higher expression in pig testis and therefore classified as group enriched, for example CYP11A1, FDX1 and STAR.

Epididymis

After the sperm production in testis, spermatozoa are stored and protected in the epididymis, where an interactive and dynamic microenvironment in the fluid enable further maturation of the sperm cells and function as a protected storage of the cells until their passage further to the ductus deferens.

Gene expression in epididymis is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 3 summarizes the number of genes in respective category. In total, 15311 genes are detected above cut of (1NX) in pig epididymis. The tissue distribution category highlights 39 genes only detected in epididymis while 685 genes are classified as epididymis elevated compared to other tissues. Table 3 shows the overlap for the epididymis elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 3, (A) The distribution of all genes across the five categories based on transcript specificity in epididymis as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in epididymis as well as in all other tissues. The combination of the two categories is shown in table 3.

Table 3, Number of genes in the subdivided categories of elevated expression and tissue distribution in epididymis

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 24611713 115
Group enriched 093258 126
Tissue enhanced 1514920278 444
Total 3930324499 685

Epididymis enriched expression

The first line of defense against invading pathogens is the antimicrobial proteins. A vast majority of these proteins are produced and secreted by immune cells, phagocytic and dendritic cells, but also by all cavity-lining epithelial cells in the body. The beta defensin gene family consists of diverse, antimicrobial and antiviral peptides that have evolved through repeated gene duplication, and are found in most vertebrates. The epididymis is highly exposed to infections by pathogens due to its anatomical position close to the urethra, and the expression of several beta defensins are entirely restricted to the epididymis, therefore many of the genes classified as enriched in epididymis are defensins, such as DEFB125, DEFB129 and DEFB128.

Table 4, The 10 genes with the highest level of enriched expression in epididymis. "mRNA (tissue)" shows the transcript level in epididymis as NX values. "Tissue specificity score (TS)" corresponds to the fold-change between the expression level in epididymis and the tissue with second highest expression level.

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000007217 DEFB125 beta-defensin 125 precursor Detected in single 135.8 1022
ENSSSCG00000040952 PATE1 prostate and testis expressed 1 Detected in single 99.3 799
ENSSSCG00000007218 DEFB129 beta-defensin 129 precursor Detected in single 327.3 679
ENSSSCG00000027178 MAN2B2 epididymis-specific alpha-mannosidase precursor Detected in single 415.8 489
ENSSSCG00000014544 na Detected in single 55.3 164
ENSSSCG00000034764 na Detected in some 170.5 154
ENSSSCG00000007122 CST11 cystatin 11 Detected in some 311.5 126
ENSSSCG00000029743 SPAG11 sperm associated antigen 11 Detected in some 234.6 124
ENSSSCG00000030743 na Detected in some 146.3 117
ENSSSCG00000034080 DEFB128 defensin beta 128 Detected in some 252.3 100




Immunohistochemical labelling of ADAM7 in pig epididymis, using HPA008879.

Ductus deferens

The main purpose of the ductus deferens is to transport sperm cells from the epididymis to the ampulla close to the seminal vesicle and out to the urethra. The ductus deferens consists of a tube lined with ciliated pseudo stratified columnar epithelium, surrounded by smooth muscle.

Gene expression in ductus deferens is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 4 summarizes the number of genes in respective category. In total, 15527 genes are detected above cut of (1NX) in pig ductus deferens. The tissue distribution category highlights 2 genes only detected in ductus deferens while 152 genes are classified as ductus deferens elevated compared to other tissues. Table 5 shows the overlap for the ductus deferens elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 4, (A) The distribution of all genes across the five categories based on transcript specificity in ductus deferens as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in ductus deferens as well as in all other tissues. The combination of the two categories is shown in table 5.

Table 5, Number of genes in the subdivided categories of elevated expression and tissue distribution in ductus deferens

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 0410 5
Group enriched 024114 39
Tissue enhanced 237636 108
Total 2657510 152

Ductus deferens enriched expression

Table 6, The 5 genes classified as enriched in ductus deferens. "mRNA (tissue)" shows the transcript level in ductus deferens as NX values. "Tissue specificity score (TS)" corresponds to the fold-change between the expression level in ductus deferens and the tissue with second highest expression level.

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000033173 SPINK14 serine peptidase inhibitor, Kazal type 14 (putative) Detected in some 155.6 24
ENSSSCG00000033906 na Detected in some 90.7 9
ENSSSCG00000039674 SPINK6 serine peptidase inhibitor, Kazal type 6 Detected in many 23.1 8
ENSSSCG00000000943 OTOGL otogelin like Detected in some 41.3 6
ENSSSCG00000012308 DGKK diacylglycerol kinase kappa Detected in some 25.4 5

Penis

The corkscrew-shaped penis of the pig was sampled at the middle of the penis body. Cross section of consecutive tissue can be found at the pig tissue dictionary. The penis tissue includes the urethra, lined with epithelia, no glandular tissue and is mainly composed by connective and erectile tissue.

Gene expression in penis is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 5 summarizes the number of genes in respective category. In total, 15066 genes are detected above cut of (1NX) in pig apenis. The tissue distribution category highlights 5 genes only detected in penis while 117 genes are classified as penis elevated compared to other tissues. Table 7 shows the overlap for the penis elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 5, (A) The distribution of all genes across the five categories based on transcript specificity in penis as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in penis as well as in all other tissues. The combination of the two categories is shown in table 7.

Table 7, Number of genes in the subdivided categories of elevated expression and tissue distribution in penis

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 0000 0
Group enriched 01320 15
Tissue enhanced 541524 102
Total 554544 117

In global expression analysis, penis cluster separate from the other male tissues and insterad show a global expression relation to ear cartilage, synovial tissue and other soft tissue rich tissue samples. Penis and ear cluster share the expression of elastin ELN, although much higher in the elastic cartilage of the ear, and the expresison in penis is one second place. Extracellular matrix related protines seem to play an important part in the expression profile of penis tissue, TGFB3 regulating extracellular matrix proteins is highly expressed in penis. The extracellular matrix protein NOV, is an example of gene with speces difference, in human NOV is higlh expressed in the adrenal gland and the pig adrenal gland show very low expression level while pig penis show high expression.

Other examples of genes highly expressed by penis are GREB1 shared with epididymis, prostate and female tissues and HOXD12 that is shared with male glands and urinary bladder.

Male glands

The tissues grouped as male glands represents the male accessory glands (seminal vesicle, prostate, bulbourethral gland and urethral glands) differ slightly between human and pig, both anatomically and histologically. The seminal vesicle is much larger in pig compared to human, the prostate is positioned a little differently. Pig has a large paired tubuloalveolar gland (bulbourethral gland) and urethral gland is surrounded by bulbourethral muscle. The glands are all secretory and the products both activate and nourish the spermatozoa in the semen. The pig accounts for the largest ejaculate volume among domestic animals, and the pig semen includes a large amout of gel.

The seminal fluild secreted from seminal vesicle function as a vehicle for the spermatozoa transport, while the secreted product of prostate formes a vaginal plug and also promotes movement of the spermatozoa. Bulbourethral gland produces mucus for lubrication and for clearing the urethra from urine, and finally The highest expression level among the four tissues was used as representative expression for male glands when comparing to other tissue types.

Gene expression in male glands is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 6 summarizes the number of genes in respective category. In total, 16395 genes are detected above cut of (1NX) in pig male glands. The tissue distribution category highlights 9 genes only detected in male glands while 1018 genes are classified as male glands elevated compared to other tissues. Table 9 shows the overlap for the male glands elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 6, (A) The distribution of all genes across the five categories based on transcript specificity in male glands as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in male glands as well as in all other tissues. The combination of the two categories is shown in table 9.

Table 9, Number of genes in the subdivided categories of elevated expression and tissue distribution in male glands

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 419919 51
Group enriched 078327 117
Tissue enhanced 5146403296 850
Total 9243444322 1018

Male glands enriched expression

Since the main function of the male glands is to secrete proteins, there is no suprise that most genes classified as male gland enriched are secretd proteins, and also highly abundant. Many of the highly enriched genes are due to the high expression in seminal vesicle, for example SPMI, PSP-II and BSP1. prostate also accounts for several highly enriched genes, such as MSMB.

*Table 10, The 10 genes with the highest level of enriched expression in male glands. "mRNA (tissue)" shows the transcript level in male glands as NX values. "Tissue specificity score (TS)" corresponds to the fold-change between the expression level in male glands and the tissue with second highest expression level.

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000021947 na carbohydrate-binding protein AQN-1 precursor Detected in single 463.6 1799
ENSSSCG00000002977 SPMI seminal plasma sperm motility inhibitor precursor Detected in some 2013.1 231
ENSSSCG00000002978 PSP-II major seminal plasma glycoprotein PSP-II precursor Detected in some 3633.2 178
ENSSSCG00000003053 BSP1 seminal plasma protein pB1 precursor Detected in some 499.5 133
ENSSSCG00000002980 na Detected in some 2217.0 60
ENSSSCG00000037112 na Detected in some 2205.1 53
ENSSSCG00000026852 NPPC C-type natriuretic peptide CNP-22 CNP-29 CNP-53 Detected in some 300.0 35
ENSSSCG00000029257 na sperm-associated acrosin inhibitor precursor Detected in some 566.2 29
ENSSSCG00000010400 MSMB beta-microseminoprotein precursor Detected in some 2458.8 25
ENSSSCG00000038157 na Detected in single 3.7 19

Expression comparison between the different male glands

The expression variation category enables further details related to expression heterogeneity within the tissues that include grouped samples. In this case, expression profile from the 4 male glands (seminal vesicle, prostate, bulbourethral gland and urethral gland) are compared and resulted in 3366 genes highlighted as variable.

An example of a protein with variable expression levels within the male glands is the PRR15, which is highly expressed in bulborethral, urathral and prostate while seminal vesicle show much lower expression.



Immunohistochemical detected of PRR15 in pig prostate and seminal vesicle using antibody HPA048121

Group enriched expression

The developmental origin and hormonal regulation of reproductive tissues can be observed by studyin the genes classified as group enriched. Several genes are shared among the male reproductive tissues, for example ABCB5 expressed by epididymis and ductus deferens, FOXR2specifically expressed by testis and epididymis, ADAM29 expressed by testis, epididymis and seminal vesicle. Several genes are shared between epididymis and ductus deferens, as well as epididymis and seminal vesicle. One eample of a protein shared between epididymis, ductus deferens and seminal veiscle is the transcription factor PAX2, important during development, which is also expressed in kidney and endometrium.





Immunohistochemical detection of PAX2 in pig epididymis, ductus deferens, seminal vesicle and endometrium using antibody HPA047704.

In order to illustrate the relation of male tissue to other tissue types, a network plot was generated, displaying the number of genes shared between different tissue types. Group enriched genes are defined as genes showing a 4-fold higher average level of mRNA expression in a group of 2-5 tissues, compared to all other tissues.

Figure 7. An interactive network plot of the male tissue enriched and group enriched genes connected to their respective enriched tissues (grey circles). Black circles shows tissues representing the male reproductive tissues. Red nodes represent the number of tissue enriched genes and orange nodes represent the number of genes that are group enriched. The sizes of the red and orange nodes are related to the number of genes displayed within the node. Each node is clickable and results in a list of all enriched genes connected to the highlighted edges. The network is limited to group enriched genes in combinations of up to 3 tissues, but the resulting lists show the complete set of group enriched genes in the particular tissue.