Female tissues


The sow´s cycle is non-seasonal and polyestrous with an estrous cycle of 18–24 days and the uterus includes the uterine horn, see Figure 1.

The hormonal cycle of the sow can be described as two phases (luteal and follicular phase) separated by ovulation (Figure 2). The reproductive cycle is classified based on the sexual behavior, estrus - when the sow is sexually receptive or non-estrus.

Figure 1, Schematic drawing of the sow genitalia

In the pig transcriptomic analysis, ovary, fallopian tube, endometrium (uterine wall), cervix and breast are the samples included and used for quantitative RNA measurements, representing female tissues. The gene classification strategy highlights genes with an elevated level of expression in one or a group of tissues compared to all other tissues.

343 genes are classified as ovary elevated out of which 20 genes are highlighted as ovary enriched. Based on the expression in fallopian tube 432 genes are classified as elevated compared to other tissue types, out of which 4 genes are highlighted as tissue enriched in fallopian tube. RNA expression levels in endometrium resulted in 140 genes classified as elevated in endometrium compared to other tissue types, out of which 4 genes are highlighted as tissue enriched in endometrium. 169 genes are classified as cervix elevated out of which 4 genes are highlighted as cervix enriched. Based on the expression in breast 153 genes are classified as elevated compared to other tissue types, out of which 14 genes are highlighted as tissue enriched in breast.


Figure 2. Schematic drawing of the sow hormonal cycle, recreated from Lorenzen E et al. (2015)

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

Ovary

The ovary is a paired organ and is an important component of the sow genitalia system. The main function is to produce and release mature oocytes and produce hormones that are needed for reproduction. The released oocyte travels down the fallopian tube into the uterus. In the ovary, oocytes mature in fluid-filled follicles that changes during the maturation process. After rupture of the follicle and release of the mature oocyte, the follicle changes into a corpus luteum, secreting steroid hormones.

Gene expression in ovary 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, 14772 genes are detected above cut of (1NX) in pig ovary. The tissue distribution category highlights 5 genes only detected in ovary while 343 genes are classified as ovary elevated compared to other tissues. Table 1 shows the overlap for the ovary 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 ovary as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in ovary 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 ovary

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 2684 20
Group enriched 035193 57
Tissue enhanced 36611483 266
Total 510714190 343

Ovary enriched expression

In the list of genes classified as ovary enriched in pig, there are several genes with different expression profiles in human and pig. In fact, most of the top 10 enriched - genes in pig ovary is not expressed or expressed at much lower levels in the human ovary. For example, BMP15 show ovary enriched expression in pig but is not detected in the human ovary. TMEM56 is classified as ovary enriched in pig but liver enriched in human, SERPINE1 is highly expressed in the human placenta (not included in the pig dataset) and human ovary only show low expression, compared to the pig ovary expression enrichment. The prorelaxin precursor (RLN2) is only detected in the pig ovary, while in human show high expression in nucleus accumbens of the basal ganglia in brain.

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

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000005216 na prorelaxin precursor Detected in single 152.3 222
ENSSSCG00000040725 IL11 interleukin 11 Detected in some 130.1 36
ENSSSCG00000003599 na Detected in some 31.8 25
ENSSSCG00000012277 TIMP1 metalloproteinase inhibitor 1 precursor Detected in many 718.4 19
ENSSSCG00000005688 PTGES prostaglandin E synthase Detected in many 139.1 13
ENSSSCG00000025698 SERPINE1 plasminogen activator inhibitor 1 precursor Detected in many 397.7 12
ENSSSCG00000024689 HTR3C 5-hydroxytryptamine receptor 3C Detected in some 22.3 12
ENSSSCG00000039854 na Detected in many 138.1 11
ENSSSCG00000012310 BMP15 bone morphogenetic protein 15 Detected in single 7.5 11
ENSSSCG00000012583 ACSL4 long-chain-fatty-acid--CoA ligase 4 Detected in all 135.8 8

Fallopian tube

The fallopian tubes, also called oviducts, are the tubes in which the oocytes transports through from the ovary to the uterus. The inside of the fallopian tube is lined with simple columned epithelium with ciliated cells. Gene expression in fallopian tube 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, 15717 genes are detected above cut of (1NX) in pig fallopian tube. The tissue distribution category highlights 4 genes only detected in fallopian tube while 432 genes are classified as fallopian tube elevated compared to other tissues. Table 3 shows the overlap for the fallopian tube 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 fallopian tube as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in fallopian tube 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 fallopian tube

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 3100 4
Group enriched 029130 42
Tissue enhanced 118816631 386
Total 421817931 432

Fallopian tube enriched expression

The large portion of ciliated cells in the fallopian tube is evident when examining the genes classified as fallopian tube elevated. However, since several other tissue types also include ciliated cells, cilia associated genes are often shared and classified as group enriched or enhanced. The fallopian tube only resulted in 4 genes classified as tissue enriched, where the gene with highest tissue specificity is oviduct-specific glycoprotein precursor OVGP1. OVGP1 codes for a large glycoprotein and suggested to be secreted from non-ciliated cells in the fallopian tube.

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

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000006791 OVGP1 oviduct-specific glycoprotein precursor Detected in single 113.4 150
ENSSSCG00000039165 na Detected in single 10.9 21
ENSSSCG00000034869 na Detected in single 1.8 18
ENSSSCG00000036279 na Detected in some 21.5 15

Among the genes classified as group enriched in fallopian tube, the estrogen receptor alpha (ESR1) and progesterone receptor (PGR) are shared among female tissues while there are plenty of cilia-related proteins shared with tissues from the respiratory system, choroid plexus, testis and epididymis. Some genes are characterized and associated with cilia (CFAP65, BBOF1 and AXDND1) and some genes are less characterized, such as C2orf70, C1orf87 and ENSSSCG00000017814.

PACRG is an example of a protein classified as testis enriched and detected in the cilia of fallopian tube.


Immunohistochemically labeled PACRG using HPA066293

Endometrium

The endometrium refers to the inner epithelial layer and mucosal layer of the uterus. In pig the uterus body is rather small in comparison to human. The uterine body acts more lika a part of the birth canal in pig rather than location for the development of fetuses, which instead mainly takes place in the uterine horn. The samples used in this study was taken from the upper part of the uterine body. Gene expression in endometrium 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, 15439 genes are detected above cut of (1NX) in pig endometrium. The tissue distribution category highlights 2 genes only detected in adrenal gland while 140 genes are classified as endometrium elevated compared to other tissues. Table 5 shows the overlap for the endometrium 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 endometrium as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in endometrium 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 endometrium

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 1111 4
Group enriched 02540 29
Tissue enhanced 260378 107
Total 386429 140

Endometrium enriched expression

The endometrium tissue resulted only in 4 genes classified as tissue enriched. All without gene name, except SLC39A2, coding for a zink transporter, which in human is mainly expressed by male tissues.

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

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000034740 na Detected in some 23.2 10
ENSSSCG00000040485 na Detected in single 2.5 8
ENSSSCG00000039158 na Detected in all 43.3 5
ENSSSCG00000029754 SLC39A2 zinc transporter ZIP2 Detected in many 24.0 5

Examples of interesting findings among the genes classified as elevated in endometrium includes ADD2 shared between female tissues, brain and endocrine tissues. Both ACTA2 and AGTR2 are shared between male and female tissues.

Cervix

Gene expression in cervix 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, 14558 genes are detected above cut of (1NX) in pig cervix. The tissue distribution category highlights 2 genes only detected in cervix while 169 genes are classified as adrenal gland elevated compared to other tissues. Table 7 shows the overlap for the cervix 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 cervix as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in cervix 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 cervix

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 1021 4
Group enriched 01330 16
Tissue enhanced 1427828 149
Total 2558329 169

Cervix enriched expression

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

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000010535 HPSE2 heparanase 2 (inactive) Detected in many 86.3 6
ENSSSCG00000033216 HRK harakiri, BCL2 interacting protein Detected in many 26.7 5
ENSSSCG00000040375 na Detected in single 2.0 5
ENSSSCG00000005706 ABL1 ABL proto-oncogene 1, non-receptor tyrosine kinase Detected in all 131.8 4

Cervix includes a large proportion smooth muscle and a layer of epithelium, viewable in the tissue dictionary. Based on the tissue specificity classification cervix expression resulted in only 4 genes classified as enriched. HPSE2 overlaps with human expression profile, with highest expression in cervix, while HRK and ABL1 shows differences to human tissue expression. Cervix share several proteins with other female tissues, such as the estrogen receptor alpha (ESR1) and progesterone receptor (PGR), classified as group enriched and shared with endometrium and fallopian tube.

Surprisingly, there is a number of genes among the cervix elevated shared with the brain, for example GALNT9, NXPH3 and TRIM9.

Breast

In the sampling of the breast tissue, the fatty tissue close to the nipple was removed from the skin and half of each tissue sample was fixed and morphologically examined while the other half was used for RNA sequencing. Gene expression in breast is categorized based on two gene classification strategies, tissue detection and tissue specificity. Figure 7 summarizes the number of genes in respective category. In total, 15368 genes are detected above cut of (1NX) in pig breast. The tissue distribution category highlights 13 genes only detected in breast while 153 genes are classified as breast elevated compared to other tissues. Table 9 shows the overlap for the breast elevated genes and tissue distribution category.

A. Specificity

B. Distribution

Figure 7, (A) The distribution of all genes across the five categories based on transcript specificity in breast as well as in all other tissues. (B) The distribution of all genes across the six categories, based on transcript detection (NX≥1) in breast 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 breast

Distribution in the 44 tissues
Detected in singleDetected in someDetected in manyDetected in all Total
Specificity
Tissue enriched 4730 14
Group enriched 01810 19
Tissue enhanced 967395 120
Total 1392435 153

Breast enriched expression

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

Gene Gene name Description Tissue distribution mRNA (tissue) Tissue specificity score
ENSSSCG00000040042 PAEP Beta-lactoglobulin-1A/1C Detected in some 314.8 19
ENSSSCG00000025468 na ceruloplasmin (ferroxidase)-like precursor Detected in some 225.0 13
ENSSSCG00000032298 ABCC11 ATP binding cassette subfamily C member 11 Detected in some 43.3 13
ENSSSCG00000037616 na Detected in some 87.1 11
ENSSSCG00000009262 CSN1S1 Alpha-S1-casein Detected in some 55.1 10
ENSSSCG00000040209 na Detected in single 5.1 10
ENSSSCG00000011556 na Detected in many 51.4 6
ENSSSCG00000009261 CSN2 beta-casein precursor Detected in some 27.0 6
ENSSSCG00000039718 na Detected in single 3.0 6
ENSSSCG00000036704 na Detected in single 2.2 6

Beta-lactoglobulin- 1A/1C Beta-lactoglobulin-1A/1C PAEP is highly enriched in pig breast tissue compared to other tissues, however, the human orthologue is not enriched in breast tissue. Compared to ABCC11, CSN1S1 and CSN2, which all show breast enriched expression in both pig and human.

An other example of breast enriched protein that differ between pig and human is SCGN, a calcium binding protein with important neuroendocrine functions. In human SCGN is highly expressed in pituitary gland and very low expression in breast, whereas in pig highest expression levels is observed in breast and much lower in pitutiary gland.





Immunohistochemically labeled SCGN in pig tissues (from left, breast, olfactory bulb, esophagus, retina) using HPA006641.

Gene expression in female tissues compared to other tissues

In order to illustrate the relation of female 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 8. An interactive network plot of the female enriched and group enriched genes connected to their respective enriched tissues (grey circles). Black circles shows tissues representing the female 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.