Nuclear lncRNA NORSF reduces E2 release in granulosa cells by sponging the endogenous small activating RNA miR-339.

BACKGROUND: Functioning as a competing endogenous RNA (ceRNA) is the main action mechanism of most cytoplasmic lncRNAs. However, it is not known whether this mechanism of action also exists in the nucleus. RESULTS: We identified four nuclear lncRNAs that are presented in granulosa cells (GCs) and were differentially expressed during sow follicular atresia. Notably, similar to cytoplasmic lncRNAs, these nuclear lncRNAs also sponge miRNAs in the nucleus of GCs through direct interactions. Furthermore, NORSF (non-coding RNA involved in sow fertility), one of the nuclear lncRNA acts as a ceRNA of miR-339. Thereby, it relieves the regulatory effect of miR-339 on CYP19A1 encoding P450arom, a rate-limiting enzyme for E2 synthesis in GCs. Interestingly, miR-339 acts as a saRNA that activates CYP19A1 transcription and enhances E2 release by GCs through altering histone modifications in the promoter by directly binding to the CYP19A1 promoter. Functionally, NORSF inhibited E2 release by GCs via the miR-339 and CYP19A1 axis. CONCLUSIONS: Our findings highlight an unappreciated mechanism of nuclear lncRNAs and show it acts as a ceRNA, which may be a common lncRNA function in the cytoplasm and nucleus. We also identified a potential endogenous saRNA for improving female fertility and treating female infertility.

Two single nucleotide variants in the miR-23a promoter affect granulosa cell apoptosis.

microRNAs (miRNAs) are well known to be important in mammalian female fertility. However, the genetic regulation of miRNAs associated with female fertility remains largely unknown. Here, we report that two single-nucleotide variants (SNVs) in the miR-23a promoter strongly influence miR-23a transcription and function in granulosa cell (GC) apoptosis. Two novel SNVs, g.-283G>C and g.-271C>T, were detected in the porcine miR-23a promoter by pooled-DNA sequencing. Furthermore, SNVs in the promoter region influenced miR-23a transcription in porcine GCs by altering its promoter activity. Functionally, SNVs in the promoter strongly influenced miR-23a regulation of early apoptosis in porcine GCs cultured in vitro. In addition, a preliminary association analysis showed that the combined genotypes of the two SNVs, rather than a single SNV, were tentatively associated with sow fertility traits in a Large White population. Overall, our findings suggest that the SNVs g.-283G>C and g.-271C>T in the miR-23a promoter are causal variants affecting GC apoptosis and miR-23a may be a potential small-molecule nonhormonal drug for regulating female fertility.

MEIS1 Is a Common Transcription Repressor of the miR-23a and NORHA Axis in Granulosa Cells.

MicroRNA-23a (miR-23a) is an endogenous small activating RNA (saRNA) involved in ovarian granulosa cell (GC) apoptosis and sow fertility by activating lncRNA NORHA transcription. Here, we reported that both miR-23a and NORHA were repressed by a common transcription factor MEIS1, which forms a small network regulating sow GC apoptosis. We characterized the pig miR-23a core promoter, and the putative binding sites of 26 common transcription factors were detected in the core promoters of both miR-23a and NORHA. Of them, transcription factor MEIS1 expression was the highest in the ovary, and widely distributed in various ovarian cells, including GCs. Functionally, MEIS1 is involved in follicular atresia by inhibiting GC apoptosis. Luciferase reporter and ChIP assays showed that transcription factor MEIS1 represses the transcription activity of miR-23a and NORHA through direct binding to their core promoters. Furthermore, MEIS1 represses miR-23a and NORHA expression in GCs. Additionally, MEIS1 inhibits the expression of FoxO1, a downstream of the miR-23a/NORHA axis, and GC apoptosis by repressing the miR-23a/NORHA axis. Overall, our findings point to MEIS1 as a common transcription repressor of miR-23a and NORHA, and develop the miR-23a/NORHA axis into a small regulatory network regulating GC apoptosis and female fertility.

TGFBR2 is a novel substrate and indirect transcription target of deubiquitylase USP9X in granulosa cells.

The ubiquitin-specific peptidase 9 X-linked (USP9X) is one of the highly conserved members belonging to the ubiquitin-specific proteases (USPs) family, which has been reported to control substrates-mediated biological functions through deubiquitinating and stabilizing substrates. Here, we have found that TGFBR2, the type II receptor of the transforming growth factor beta (TGF-ß) signaling pathway, is a novel substrate and indirect transcription target of deubiquitylase USP9X in granulosa cells (GCs). Mechanically, USP9X positively influences the expression of TGFBR2 at different levels through two independent ways: (i) directly targets and deubiquitinates TGFBR2, which maintains the protein stability of TGFBR2 through avoiding degradation mediated by ubiquitin-proteasome system; (ii) indirectly maintains TGFBR2 messenger RNA (mRNA) expression via SMAD4/miR-143 axis. Specifically, SMAD4, another substrate of USP9X, acts as a transcription factor and suppresses miR-143 which inhibits the mRNA level of TGFBR2 by directly binding to its 3'-untranslated region. Functionally, the maintenance of TGFBR2 by USP9X activates the TGF-ß signaling pathway, which further represses GC apoptosis. Our study highlights a functional micro-regulatory network composed of deubiquitinase (USP9X), small noncoding RNA (miR-143) and the TGF-ß signaling pathway, which plays a crucial role in the regulation of GC apoptosis and female fertility.

Lnc2300 is a cis-acting long noncoding RNA of CYP11A1 in ovarian granulosa cells.

The high level of progesterone and 17ß-estradiol ratio (P4/E2) in follicular fluid has been considered as a biomarker of follicular atresia. CYP11A1, the crucial gene encoding the rate-limiting enzyme for steroid hormone synthesis, has been reported differently expressed in the ovary during follicular atresia. However, the regulation mechanism of CYP11A1 expression during follicular atresia still remains unclear. Here, we have demonstrated that lnc2300, a novel pig ovary-specific highly expressed cis-acting long noncoding RNA (lncRNA) transcribed from chromosome 7, has the ability to induce the expression of CYP11A1 and inhibit the apoptosis of porcine granulosa cells (GCs). Mechanistically, lnc2300, mainly located in the cytoplasm of porcine GCs, sponges and suppresses the expression of miR-365-3p through acting as a competing endogenous RNA (ceRNA), which further relieves the inhibitory effects of miR-365-3p on the expression of CYP11A1. Besides, CYP11A1 is validated as a direct functional target of miR-365-3p in porcine GCs. Functionally, lnc2300 is an antiapoptotic lncRNA that reduces porcine GC apoptosis by inhibiting the proapoptotic function of miR-365-3p. In summary, our findings reveal a cis-acting regulation mechanism of CYP11A1 through lncRNA, and define a novel signaling pathway, lnc2300/miR-365-3p/CYP11A1 axis, which is involved in the regulation of GC apoptosis and follicular atresia.

SMARCA2 is regulated by NORFA-miR-29c, a novel pathway that controls granulosa cell apoptosis and is related to female fertility.

SMARCA2, an evolutionarily conserved catalytic ATPase subunit of SWI/SNF complexes, has been implicated in development and diseases; however, its role in mammalian ovarian function and female fertility is unknown. Here, we identified and characterized the 3'-UTR of the porcine SMARCA2 gene and identified a novel adenylate number variation. Notably, this mutation was significantly associated with sow litter size traits and SMARCA2 levels, due to its influence on the stability of SMARCA2 mRNA in ovarian granulosa cells (GCs). Immunohistochemistry and functional analysis showed that SMARCA2 is involved in the regulation of follicular atresia by inhibiting GC apoptosis. In addition, miR-29c, a pro-apoptotic factor, was identified as a functional miRNA that targets SMARCA2 in GCs and mediates regulation of SMARCA2 expression via the NORFA-SMAD4 axis. Although a potential miR-29c-responsive element was identified within NORFA, negative regulation of miR-29c expression by NORFA was not due to activity as a competing endogenous RNA. In conclusion, our findings demonstrate that SMARCA2 is a candidate gene for sow litter size traits, because it regulates follicular atresia and GC apoptosis. Additionally, we have defined a novel candidate pathway for sow fertility, the NORFA-TGFBR2-SMAD4-miR-29c-SMARCA2 pathway.This article has an associated First Person interview with the first author of the paper.

circSLC41A1 Resists Porcine Granulosa Cell Apoptosis and Follicular Atresia by Promoting SRSF1 through miR-9820-5p Sponging.

Ovarian granulosa cell (GC) apoptosis is the major cause of follicular atresia. Regulation of non-coding RNAs (ncRNAs) was proved to be involved in regulatory mechanisms of GC apoptosis. circRNAs have been recognized to play important roles in cellular activity. However, the regulatory network of circRNAs in follicular atresia has not been fully validated. In this study, we report a new circRNA, circSLC41A1, which has higher expression in healthy follicles compared to atretic follicles, and confirm its circular structure using RNase R treatment. The resistant function of circSLC41A1 during GC apoptosis was detected by si-RNA transfection and the competitive binding of miR-9820-5p by circSLC41A1 and SRSF1 was detected with a dual-luciferase reporter assay and co-transfection of their inhibitors or siRNA. Additionally, we predicted the protein-coding potential of circSLC41A1 and analyzed the structure of circSLC41A1-134aa. Our study revealed that circSLC41A1 enhanced SRSF1 expression through competitive binding of miR-9820-5p and demonstrated a circSLC41A1-miR-9820-5p-SRSF1 regulatory axis in follicular GC apoptosis. The study adds to knowledge of the post-transcriptional regulation of follicular atresia and provides insight into the protein-coding function of circRNA.

TRIM32 Deficiency Impairs Synaptic Plasticity by Excitatory-Inhibitory Imbalance via Notch Pathway.

Synaptic plasticity is the neural basis of physiological processes involved in learning and memory. Tripartite motif-containing 32 (TRIM32) has been found to play many important roles in the brain such as neural stem cell proliferation, neurogenesis, inhibition of nerve proliferation, and apoptosis. TRIM32 has been linked to several nervous system diseases including autism spectrum disorder, depression, anxiety, and Alzheimer's disease. However, the role of TRIM32 in regulating the mechanism of synaptic plasticity is still unknown. Our electrophysiological studies using hippocampal slices revealed that long-term potentiation of CA1 synapses was impaired in TRIM32 deficient (KO) mice. Further research found that dendritic spines density, AMPA receptors, and synaptic plasticity-related proteins were also reduced. NMDA receptors were upregulated whereas GABA receptors were downregulated in TRIM32 deficient mice, explaining the imbalance in excitatory and inhibitory neurotransmission. This caused overexcitation leading to decreased neuronal numbers in the hippocampus and cortex. In summary, this study provides this maiden evidence on the synaptic plasticity changes of TRIM32 deficiency in the brain and proposes that TRIM32 relates the notch signaling pathway and its related mechanisms contribute to this deficit.

SMAD4 Feedback Activates the Canonical TGF-ß Family Signaling Pathways.

TGF-ß family signaling pathways, including TGF-ß and BMP pathways, are widely involved in the regulation of health and diseases through downstream SMADs, which are also regulated by multiple validated mechanisms, such as genetic regulation, epigenetic regulation, and feedback regulation. However, it is still unclear whether R-SMADs or Co-SMAD can feedback regulate the TGF-ß family signaling pathways in granulosa cells (GCs). In this study, we report a novel mechanism underlying the feedback regulation of TGF-ß family signaling pathways, i.e., SMAD4, the only Co-SMAD, positive feedback activates the TGF-ß family signaling pathways in GCs with a basal level of TGF-ß ligands by interacting with the core promoters of its upstream receptors. Mechanistically, SMAD4 acts as a transcription factor, and feedback activates the transcription of its upstream receptors, including ACVR1B, BMPR2, and TGFBR2, of the canonical TGF-ß signaling pathways by interacting with three coactivators (c-JUN, CREB1, and SP1), respectively. Notably, three different interaction modes between SMAD4 and coactivators were identified in SMAD4-mediated feedback regulation of upstream receptors through reciprocal ChIP assays. Our findings in the present study indicate for the first time that SMAD4 feedback activates the canonical TGF-ß family signaling pathways in GCs, which improves and expands the regulatory mechanism, especially the feedback regulation modes of TGF-ß family signaling pathways in ovarian GCs.

TGF-ß/SMAD4 signaling pathway activates the HAS2-HA system to regulate granulosa cell state.

Both TGF-ß/SMAD4 signaling pathway and HAS2-HA system have been shown to control granulosa cell (GC) state in mammalian ovary. However, the regulatory relationship between TGF-ß/SMAD4 signaling pathway and HA system in GCs is not well known. Here, we report that the TGF-ß/SMAD4 signaling pathway activates the HAS2-HA system by binding directly to the HAS2 promoter, ultimately controlling the GC state via the CD44-Caspase3 axis. SMAD4-induced HAS2 expression, HAS2-driven HA secretion, and HAS2-mediated GC state (proliferation and apoptosis) by interacting directly with the promoter region of the HAS2 gene. The CD44-Caspase3 axis, located downstream of the HAS2-HA system, was also activated by SMAD4 and the TGF-ß/SMAD4 signaling pathway. However, there was no feedback regulation of the TGF-ß/SMAD4 signaling pathway by the HAS2-HA system in GCs. In addition, we found that miRNA-26b attenuated HAS2 expression via SMAD4-dependent and -independent mechanisms. Our findings provide compelling evidence that HAS2 is a direct transcriptional target of SMAD4. They also reveal a novel mechanism by which the TGF-ß/SMAD4 signaling pathway controls the GC state and alters the structural components of GCs in porcine ovaries.

GPR50 Distribution in the Mouse Cortex and Hippocampus.

G protein-coupled receptor 50 (GPR50) belongs to the G protein-coupled receptor which is highly homologous with the sequence of melatonin receptor MT1 and MT2. GPR50 expression has previously been reported in many brain regions, like cortex, midbrain, pons, amygdala. But, the distribution of GPR50 in the hippocampus and cortex and the cell types expressing GPR50 is not yet clear. In this study, we examined the distribution of GPR50 in adult male mice by immunofluorescence. Our results showed that GPR50 was localized in the CA1-3 pyramidal cells and the granule cells of the dentate gyrus. GPR50 was also expressed in excitatory and inhibitory neurons. As inhibitory neurons also contain many types, we found that GPR50 was localized in some interneurons in which it was co-expressed with the calcium-binding proteins calbindin, calretinin, and parvalbumin. Besides, similar results were seen in the cortex. The widespread expression of GPR50 in the hippocampus and cortex suggests that GPR50 may be associated with synaptic plasticity and cognitive function.

The distribution and expression of vascular endothelial growth factor A (VEGFA) during follicular development and atresia in the pig.

The involvement of vascular endothelial growth factor A (VEGFA) in ovarian physiological processes has been widely reported, but the location and role of VEGFA during follicular atresia remain unknown. This study investigated the distribution and expression of VEGFA during porcine follicular development and atresia. Pig ovaries were obtained, individual medium-sized (3-5mm in diameter) antral follicles were separated and classified into healthy, early atretic or progressively atretic groups. Immunobiology and quantitative techniques were used to investigate the varied follicular distribution of VEGFA at both the morphological and molecular level. The results indicated that VEGFA protein expression peaked in tertiary follicles, mostly distributed in the thecal and inner granulosa layers, during follicular development while VEGFA mRNA was mainly expressed in the inner granulosa layers. Additionally, healthy antral follicles showed a significantly higher expression of VEGFA than atretic follicles in both theca and granulosa cells. Knockdown of VEGFA using siRNA revealed an antiapoptosis effect of VEGFA in cultured pig granulosa cells. Our results increase the knowledge of VEGFA functions in follicles.

NF-κB "decoy" inhibits COX-2 expression in epileptic rat brain.

There is a need to investigate the role of nuclear factor kappa B in the regulation of cyclooxygenase-2 expression in the epileptic rat brain and cultured hippocampal neurons. Immunofluorescence and polymerase chain reaction was used to detect the expression of nuclear factor kappa B and cyclooxygenase-2. In cultured hippocampal neurons and rat brain: the control group compared with the normal group, nuclear factor kappa B expression in the hippocampal dentate gyrus, cerebral cortex, the piriform cortex brain regions were significantly increased (P < 0.01). This is accompanied by a significant increase in cyclooxygenase-2 protein and mRNA expressions in the hippocampus (P < 0.01). In the experimental group compared to the control group, the nuclear factor-kappa B expression in the hippocampal dentate gyrus, cerebral cortex, piriform cortex, and other brain regions was significantly lower (P < 0.01), with the accompanying decrease in cyclooxygenase-2 protein and mRNA expression (P < 0.01) in the hippocampus. In conclusion, κB-decoy can inhibit nuclear factor kappa B activation in epileptic rat brain and cyclooxygenase-2 overexpression.

A polymorphism in the transcriptional regulatory region strongly influences ovine FSHR mRNA decay.

Follicle-stimulating hormone receptor (FSHR) is an important G protein-coupled receptor, which is required for steroidogenesis, follicular development and female infertility. Here, we report a novel polymorphism in the 3'-UTR that strongly influences ovine FSHR mRNA decay. The partial 3'-UTR sequence of Hu sheep FSHR gene was isolated and characterized, and a polymorphism (c.2327A>G) was identified. Luciferase assay and qRT-PCR showed that c.2327A>G polymorphism in the 3'-UTR exerts a strong regulatory role in FSHR transcription. This regulatory role is achieved by affecting FSHR mRNA decay. Furthermore, the c.2327A>G mutation in the 3'-UTR influences ARE (AU-rich element, a cis-acting element promoting mRNA decay)-mediated mRNA decay of Hu sheep FSHR gene. Together, our study identified a novel polymorphism and elucidated a new mechanism underlying transcriptional regulation of FSHR in mammals.

Smad4 Feedback Enhances BMPR1B Transcription in Ovine Granulosa Cells.

BMPR1B is a type 1B receptor of the canonical bone morphogenetic protein (BMP)/Sma- and mad-related protein (Smad) signaling pathway and is well known as the first major gene associated with sheep prolificacy. However, little is known about the transcriptional regulation of the ovine BMPR1B gene. In this study, we identified the ovine BMPR1B gene promoter and demonstrated that its transcription was regulated by Smad4. In sheep ovarian follicles, three transcriptional variants of BMPR1B gene with distinct transcription start sites were identified using 5' RACE assay while variants II and III were more strongly expressed. Luciferase assay showed that the region -405 to -200 nt is the PII promoter region of variant II. Interestingly, two putative Smad4-binding elements (SBEs) were detected in this region. Luciferase and ChIP assay revealed that Smad4 enhances PII promoter activity of the ovine BMPR1B gene by directly interacting with SBE1 motif. Furthermore, in the ovine granulosa cells, Smad4 regulated BMPRIB expression, and BMPRIB-mediated granulosa cells apoptosis. Overall, our findings not only characterized the 5' regulatory region of the ovine BMPR1B gene, but also uncovered a feedback regulatory mechanism of the canonical BMP/Smad signaling pathway and provided an insight into the transcriptional regulation of BMPR1B gene and sheep prolificacy.

miR-181b-induced SMAD7 downregulation controls granulosa cell apoptosis through TGF-ß signaling by interacting with the TGFBR1 promoter.

SMAD7 disrupts the TGF-ß signaling pathway by influencing TGFBR1 stability and by blocking the binding of TGFBR1 to SMAD2/3. In this study, we showed that SMAD7 attenuated the TGF-ß signaling pathway in ovarian granulosa cells (GCs) by regulating TGFBR1 transcriptional activity. To function as a transcription factor, SMAD7 downregulated the mRNA levels of TGFBR1 via direct binding to the SMAD-binding elements (SBEs) within the promoter region of pig TGFBR1. We also showed that SMAD7 enhanced porcine GC apoptosis by interrupting TGFBR1 and the TGF-ß signaling pathway. Interestingly, miR-181b, a microRNA that is downregulated during porcine follicular atresia, was identified to be directly targeting SMAD7 at its 3'-UTR. By inhibiting SMAD7, miR-181b could inhibit GC apoptosis by activating the TGF-ß signaling pathway. Our findings provide new insights into the mechanisms underlying the regulation of the TGF-ß signaling pathway by SMAD7 and miR-181b.

Initiation of follicular atresia: gene networks during early atresia in pig ovaries.

In mammals, more than 99% of ovarian follicles undergo a degenerative process known as atresia. The molecular events involved in atresia initiation remain incompletely understood. The objective of this study was to analyze differential gene expression profiles of medium antral ovarian follicles during early atresia in pig. The transcriptome evaluation was performed on cDNA microarrays using healthy and early atretic follicle samples and was validated by quantitative PCR. Annotation analysis applying current database (Sus scrofa 11.1) revealed 450 significantly differential expressed genes between healthy and early atretic follicles. Among them, 142 were significantly upregulated in early atretic with respect to healthy group and 308 were downregulated. Similar expression trends were observed between microarray data and quantitative RT-PCR confirmation, which indicated the reliability of the microarray analysis. Further analysis of the differential expressed genes revealed the most significantly affected biological functions during early atresia including blood vessel development, regulation of DNA-templated transcription in response to stress and negative regulation of cell adhesion. The pathway and interaction analysis suggested that atresia initiation associates with (1) a crosstalk of cell apoptosis, autophagy and ferroptosis rather than change of typical apoptosis markers, (2) dramatic shift of steroidogenic enzymes, (3) deficient glutathione metabolism and (4) vascular degeneration. The novel gene candidates and pathways identified in the current study will lead to a comprehensive view of the molecular regulation of ovarian follicular atresia and a new understanding of atresia initiation.

A comprehensive transcriptomic view on the role of SMAD4 gene by RNAi-mediated knockdown in porcine follicular granulosa cells.

As a key mediator of the transforming growth factor-beta (TGF-ß) signaling pathway, which plays a pivotal role in regulating mammalian reproductive performance, Sma- and Mad-related protein 4 (SMAD4) is closely associated with the development of ovarian follicular. However, current knowledge of the genome-wide view on the role of SMAD4 gene in mammalian follicular granulosa cells (GCs) is still largely unknown. In the present study, RNA-Seq was performed to investigate the effects of SMAD4 knockdown by RNA interference (SMAD4-siRNA) in porcine follicular GCs. A total of 1025 differentially expressed genes (DEGs), including 530 upregulated genes and 495 downregulated genes, were identified in SMAD4-siRNA treated GCs compared with that treated with NC-siRNA. Furthermore, functional enrichment analysis indicated that upregulated DEGs in SMAD4-siRNA treated cells were mainly enriched in cell-cycle related processes, interferon signaling pathway, and immune system process, while downregulated DEGs in SMAD4-siRNA treated cells were mainly involved in extracellular matrix organization/disassembly, pathogenesis, and cell adhesion. In particular, cell cycle and TGF-ß signaling pathway were discovered as the canonical pathways changed under SMAD4-silencing. Taken together, our data reveals SMAD4 knockdown alters the expression of numerous genes involved in key biological processes of the development of follicular GCs and provides a novel global clue of the role of SMAD4 gene in porcine follicular GCs, thus improving our understanding of regulatory mechanisms of SMAD4 gene in follicular development.

Androgen receptor and miRNA-126* axis controls follicle-stimulating hormone receptor expression in porcine ovarian granulosa cells.

Androgen, which acts via the androgen receptor (AR), plays crucial roles in mammalian ovarian function. Recent studies showed that androgen/AR signaling regulates follicle-stimulating hormone receptor (FSHR) expression in follicles; however, the detailed mechanism underlying this regulation remained unknown. Here, we demonstrate that AR and miR-126* cooperate to inhibit FSHR expression and function in pig follicular granulosa cells (pGCs). In pGCs, overexpression of AR decreased, whereas knockdown increased, FSHR mRNA and protein expression; however, neither manipulation affected FSHR promoter activity. Using a dual-luciferase reporter assay, we found that the FSHR gene is a direct target of miR-126*, which inhibits FSHR expression and increases the rate of AR-induced apoptosis in pGCs. Collectively, our data show for the first time that the AR/miR-126* axis exerts synergetic effects in the regulation of FSHR expression and apoptosis in pGCs. Our findings thus define a novel pathway, AR/miR-126*/FSHR, that regulates mammalian GC functions.

RNA-seq analysis reveals new candidate genes for drip loss in a Pietrain × Duroc × Landrace × Yorkshire population.

Drip loss, one of the most important meat quality traits, is characterized by low heritability. To date, the genetic factors affecting the drip loss trait have not been clearly elucidated. The objective of this study was to identify critical candidate genes affecting drip loss. First, we generated a Pietrain × Duroc × Landrace × Yorkshire commercial pig population and obtained phenotypic values for the drip loss trait. Furthermore, we constructed two RNA libraries from pooled samples of longissimus dorsi muscles with the highest (H group) and lowest (L group) drip loss and identified the differentially expressed genes (DEGs) between these extreme phenotypes using RNA-seq technology. In total, 25 883 genes were detected in the H and L group libraries, and none was specifically expressed in only one library. Comparative analysis of gene expression levels found that 150 genes were differentially expressed, of which 127 were upregulated and 23 were downregulated in the H group relative to the L group. In addition, 68 drip loss quantitative trait loci (QTL) overlapping with 63 DEGs were identified, and these QTL were distributed mainly on chromosomes 1, 2, 5 and 6. Interestingly, the triadin (TRDN) gene, which is involved in muscle contraction and fat deposition, and the myostatin (MSTN) gene, which has a role in muscle growth, were localized to more than two drip loss QTL, suggesting that both are critical candidate genes responsible for drip loss.