- Open Access
Integrated transcriptome and proteome revealed that the declined expression of cell cycle-related genes associated with follicular atresia in geese
BMC Genomics volume 24, Article number: 24 (2023)
Geese exhibit relatively low reproductive performance, and follicular atresia is an important factor that restricts the egg production of geese. Systematic analysis of the regulation of follicle atresia in geese through transcriptome and proteome levels could provide meaningful information on clarifying the mechanism of follicle atresia in poultry.
The granulosa cell layer was loose, disintegrated and showed apoptosis in atretic follicles and remained intact in normal follicles. The hormone levels of FSH and LH were significantly decreased in the atresia follicles compared to the normal follicles (P < 0.05). A total of 954 differentially expressed genes (DEGs, 315 increased and 639 decreased) and 161 differentially expressed proteins (DEPs, 61 increased and 100 decreased) were obtained in atresia follicles compared to normal follicles, of which, 15 genes were differentially expressed in both transcriptome and proteome. The DEGs were mainly enriched in sodium transmembrane transport, plasma membrane, and transmembrane transporter activity based on the GO enrichment analysis and in the cell cycle pathway based on the KEGG enrichment analysis. The DEPs were mainly enriched in localization, lysosome, and phospholipid-binding based on the GO enrichment analysis. Candidate genes Smad2/3, Smad4, Annexin A1 (ANXA1), Stromelysin-1 (MMP3), Serine/threonine-protein kinase (CHK1), DNA replication licensing factor (MCM3), Cyclin-A2 (CCNA2), mitotic spindle assembly checkpoint protein (MAD2), Cyclin-dependent kinase 1 (CDK1), fibroblast growth factor 12 (FGF12), and G1/S-specific cyclin-D1 (CCND1) were possibly responsible for the regulation of atresia.
The cell cycle is an important pathway for the regulation of follicular atresia. Sodium outflow and high expression of MMP3 and MMP9 could be responsible for structural destruction and apoptosis of follicular cells.
Follicular atresia is an important factor restricting egg production performance. In birds, less than 5% of the follicles could complete their development process for mature and ovulation, while the remaining 95% become atretic follicles and eventually get degraded . Goose, especially indigenous breeds, exhibits low egg production performance, which causes a relatively higher cost of goslings and lower marketing numbers as compared to broilers and meat ducks. The prevention of follicular atresia and promotion of the normal development of follicle could help increase egg production. Therefore, understanding the mechanism of follicular atresia in geese and screening key signal pathways or genes that regulate follicular atresia would be helpful for identifying molecular markers for selection of goose with high level of atresia.
Follicular atresia is a spontaneous process in which follicles naturally degenerate without an inflammatory response and is generally considered to be caused by defects in the oocyte itself or a lack of necessary growth factors and hormones . The TGFβ superfamily plays an important role in follicular development, in which, BMP-4 and BMP-7 are positive regulators in the transition from the primordial to the primary follicle, while the anti-Mullerian hormone (AMH) is a negative regulator for this transition . The bFGF inhibits follicular atresia through the PI3K-Akt pathway in chickens , while FGF10 promotes follicular atresia by decreasing the secretion of estradiol in cattle . In bovine, FGF18 inhibited the secretion of estradiol and progesterone, alters the progression of the cell cycle, and promotes follicular atresia by accelerating the apoptosis of granulosa cells . Hormones have been reported to be involved in follicular development, ovulation, atresia, and other processes. The follicle-stimulating hormone (FSH) could prevent follicular atresia through inhibiting the apoptosis of granulosa cells by inducing their differentiation, the formation of luteinizing hormone receptors, and the synthesis of follicular steroid hormones . It has been proven that atretic follicles exhibited a reduced production of estradiol in granulose cells of porcine .
Atretic follicle showed shrunken surface with hemorrhagic spots. The histological and ultrastructural analysis of atretic follicles in two fish breeds revealed the degradation of mitochondria and marked phagocytic activity of digestive vacuoles, myelin, and lipofuscin granules . Granulosa cells were considered to be the first ones to undergo apoptosis in the follicle, which then induced the overall apoptosis and atresia . The atresia of monkey follicles was characterized by three consecutive stages including morphological alterations, cessation of the proliferation and apoptosis, and fragmentation of DNA in granulosa cells that occurred only in late atresia . Furthermore, the follicular development or atresia is thought to be regulated by the survival and apoptosis of granulosa cells in mammals . Unlike the granulosa cells of mammals, which are scattered in the follicular fluid, follicles of poultry have distinct layers of granulosa cells, which make them useful experimental models for the study of follicular atresia.
Many reports have stated the mechanism about follicular atresia in mammals, fish, however, still less was conducted on the systematic research of the molecular mechanism of follicular atresia in poultry. In this study, the histology and reproductive hormones level were compared in normal and atretic follicles. The site of apoptosis of the atretic follicle was first determined in goose. Granulose cells were used for integrated transcriptome and proteome analyses to interpret the molecular mechanism of follicular atresia in geese.
Histology and hormone levels of atretic and normal follicles
The granulosa cell layer was loose and disintegrated in atretic follicles and remained intact in normal follicles (Fig. 1a). Apoptosis was observed in all parts of atretic follicles and was mostly exhibited in the collapsed granulosa cell layer, which was not observed in the normal follicles (Fig. 1a). The hormone levels of FSH and LH were significantly decreased in the atresia follicles compared to the normal follicles (P < 0.05). The levels of hormones E2, P4, and PRL showed no significant differences between atretic and healthy follicles (Fig. 1b).
De novo transcriptome analysis
A total of 141,051,229 raw reads and 135,590,434 clean reads were obtained. The error rates of sequencing were < 0.03, the Q20 and Q30 were > 90%, and the mapped rate was > 80% (Supplementary Table S1). After de novo assembly, 140,475 transcripts and 70,318 unigenes with a mean length of 2007 bp and 1436 bp, respectively, were obtained (Supplementary Table S2). The complete single-copied genes accounted for 83.3% of the all Benchmarking Universal Single-Copy Orthologs (BUSCO), indicating a qualified transcriptome assembled (Fig. 2b). A total of 46,696 unigenes were annotated, and among the seven databases, the Nr database had the highest annotations as compared to other six databases (Supplementary Table S3). According to the results of gene annotation, Anser cygnoides showed the highest similarity with the geese among the annotated species (Fig. 2a). The expression levels of all genes were expressed in FPKM values (Supplementary Table S4) and were mainly distributed in the range of 3.57–15 (Fig. 2c).
The principal component analysis (PCA) indicated a significant difference in mRNA expression between AF and NF groups (Fig. S1). A total of 954 DEGs were obtained, of which, 315 were upregulated, and 639 were downregulated in AF when compared to NF (Fig. 3a & b). A detailed analysis of the top 10 DEGs showed that matrix metallopeptidase-3 (MMP3), which was responsible for proteolysis, was 11.68 times upregulated in AF (Table 1), and matrix metallopeptidase-7 (MMP7), a disintegrin, metalloproteinase with thrombospondin motifs 18 (ADAMTS18) were the matrix metalloproteinases all present in the top 10 upregulated genes (Supplementary Table S5). The downregulated gene with the highest change was the Kazal-type serine protease inhibitor domain-containing protein 1-like gene (log2FC = − 6.19), which is responsible for the regulation of cell growth (Table 1). In addition, growth factors like fibroblast (FGF12) and platelet-derived growth factor subunit A (PDGFA) were also downregulated in AF (Supplementary Table S5).
The GO enrichment analysis was used for the functional analysis of DEGs. Among biological processes, DEGs were over-represented in the sodium ion transmembrane transport (GO: 0035725), pyruvate metabolic process (GO: 0006090), and arginine metabolic process (GO: 0006525). Regarding cellular components, DEGs were over-represented in the plasma membrane (GO: 0005886), cell periphery (GO: 0071944), andprotein-N(PI)-phosphohistidine-sugar phosphotransferase complex (GO: 0009357). In terms of the molecular functions, DEGs were over-represented in transmembrane transporter activity (GO: 0022857), oxaloacetate decarboxylase activity (GO: 0008948), and peptidase activity (GO: 0008233) (Supplementary Table S6, Fig. 4a). Based on the KEGG enrichment analysis, the DEGs were over-represented only in the cell cycle pathway (P < 0.05) in which Smad 2–3, Smad 4, cyclin D, cyclin A, cyclin-dependent kinase (CDK1, CDK2), mini-chromosome maintenance (MCM3), Myelin transcription factor 1 (Myt1), Serine/threonine-protein kinase (chk1, chk 2), mitotic spindle assembly checkpoint protein (mad2), and Securin (PTTG) were downregulated, while the ECREB-binding protein (p300) and cyclin-dependent kinase inhibitor (Ink4b) were upregulated (Supplementary Table S7, Fig. 4b&c).
The Protein−Protein Interaction Network (PPI) analysis was performed to further explore the key genes that regulate the occurrence of follicular atresia (Fig. 5a). In the PPI based on the transcriptome data, cyclic AMP-responsive element-binding protein 3-like protein 1 (CR3L1), Myb-related protein A (MYBA), Set1/Ash2 histone methyltransferase complex subunit (ASH2L), histone-lysine N-methyltransferase 2D (KMT2D), zinc finger protein 830 (ZN830), and retinoic acid receptor RXR-alpha (RXRA), which are involved in the regulation of cell cycle, upregulated the EP300 in atretic follicles. The downregulated lutropin-choriogonadotropin hormone receptor (LSHR), estradiol 17-beta-dehydrogenase 1 (DHB1), aromatase (CP19A), androgen receptor (ANDR), N-acetylglucosamine-6-sulfatase (GNS), endophilin-A1 (SH3G2), protein-patched homolog 1 (PTC1), sodium channel protein type 5 subunit alpha (SCN5A), and FGF12, and the upregulated low-density lipoprotein receptor-related protein 2 (LRP2) inhibited cell cycle in atretic follicles (Fig. 5b).
Mass spectrometry produced a total of 604,048 spectra, 61,114 matched spectra, 30,333 peptides, 4657 identified proteins, and 4654 quantifiable proteins (Supplementary Table S8). The length of most peptides varied between eight and 16 amino acids, protein coverage more than 10% was 54.07%, and the protein mass was mainly distributed in the range of 20 kDa–30 kDa and > 100 kDa (Supplementary Fig. S2). The principal component analysis (PCA) indicated a significant difference in protein expression between AF and NF groups (Fig. 6a). A low coefficient variance (CV) of 0.2 indicated good repeatability of each proteomic library (Fig. 6b).
A total of 161 DEPs was obtained, among which 61 were upregulated and 100 were downregulated in AF (Fig. 6c & d, Supplementary Table S9). The upregulated protein with the highest change was the MMP3 (log2FC = 1.41), responsible for proteolysis. The downregulated protein with the highest change was the U3 small nucleolar RNA-associated protein 18 homolog (UTP18, log2FC = − 1.28), responsible for protein processing (Table 2). The three proteins: vitellogenin-2-like (down), transgelin isoform X1 (up), and low-density lipoprotein receptor-related protein 1 isoform X1 (down) showed the highest expression change among DEPs (Supplementary Table S9).
Among the biological processes, DEPs were over-represented in localization (GO: 0051179), transport (GO: 0006810), and transport of organic substances (GO: 0071702). Regarding the cellular components, DEPs were over-represented in the lysosome (GO: 0005764), retromer complex (GO: 0030904), and the extracellular region (GO: 0044421). In the case of molecular function, the DEPs were over-represented in phospholipid binding (GO: 0005543), lipid binding (GO: 0008289), and transporter activity (GO: 0005215, Fig. 7A, Supplementary Table S10). The subcellular localization of DEPs was mainly in the nuclear protein (22.02%), cytoplasm protein (18.81%), and extracellular protein (16.06%, Fig. 7b).
The PPI of proteome showed that the upregulation of MMP3, MMP9, CD44, annexin A1 (ANXA1), ANXA5, and the downregulation of low-density lipoprotein receptor-related protein 1 (LRP1), disabled homolog 2 (DAB2), alpha-2-macroglobulin (A2M), vitronectin (VTN), and AMBP were involved in proteolysis (Fig. 8a & b). The upregulation of MMP3 and MMP9, along with the downregulation of A disintegrin and metalloproteinase with thrombospondin motifs 13 (ADAMTS13), sulfhydryl oxidase 1 (QSOX1), VTN, and LRP1 were involved in the disassembly of the extracellular matrix. The upregulated expression of MMP3 and MMP9 played a key role in regulating follicle atresia through the promotion of protein degradation (Fig. 8c). To test the accuracy of mass spectrometry, western blotting of MMP3 and MMP9 was performed, and the results revealed high expression of MMP3 and MMP9 in atretic follicles (Fig. 8d & e).
Integrated analysis of transcriptome and proteome
A weak correlation between the mRNA and protein expression was observed according to Pearson’s correlation coefficient (0.107) (Fig. 9a, Supplementary Table S11). Fifteen DEGs were identified both in proteome and transcriptome and eight showed the same expression trend (Fig. 9b, Table 3). Among them MMP3, neuron navigator 1 (NAV1), and ANXA1 were upregulated, while adrenodoxin-like (ADX), glucosamine (N-acetyl)-6-sulfatase (GNS), carboxypeptidase N subunit 1 (CPN1), C-factor (CSGA), and VTN were downregulated. The GO analysis indicated that the DEGs in both proteome and transcriptome were mainly enriched in the binding and extracellular region (Fig. 9c).
qRT-PCR analysis of the candidate genes
According to the transcriptome, twelve candidate genes were selected and verified. CHK1, MCM3, CCNA2, MAD2, CDK1, and CCND1 at the nodes and played important role in cell cycle pathway. Smad2–3 and Smad4 at the transcriptional factors of TGFβ for transmembrane signal transduction. ANXA1 and MMP3 that differentially expressed in both transcriptome and proteome were chosen for additional quantitative determination with qRT-PCR. The expression level of the selected genes followed a similar trend (R2 = 0.9764) in the relative expression levels between the log2FC (RNA-Seq) and log2FC (qPCR), indicating the reliability of our transcriptome data (Fig. 10).
Reproductive hormones play an important regulatory role in follicular development. The FSH and LH promote follicular maturation and ovulation and inhibit the apoptosis of granulosa cells in human . Ilha et al. suggested that the activation of LIF-STAT3 in follicular granulosa cells resulted in the decreased level of FSH and causes follicle atresia in cattle . The FSH and LH are glycoproteins that are secreted by the pituitary gland and need to be transported through the blood to the target site to be effective. A study by Feranil et al. demonstrated that the distribution of peripheral vasculature in atretic follicles was lesser than that of normal follicles in swamp buffalo , which possibly led to lower levels of FSH and LH. The in vitro inhibition of angiogenesis could directly induce the apoptosis of granulosa cells and thus caused follicular atresia in rats . Lower levels of FSH and LH were also detected in AF of Yangzhou geese in this study, which demonstrated that steroid hormones were involved in follicle atresia. The PDGFA, which regulated cell growth and angiogenesis by binding to its receptor, could increase follicle size and enhance the transition from primordial to primary follicles in the in vitro culture of sheep follicles , was downregulated in AF in this study. This suggested that decreased angiogenesis was associated with follicle atresia.
Granulosa cells represented the first site of the initiation of apoptosis in rat atretic follicles , and follicular growth or atresia were regulated by the survival or death of granulosa cells in mammalian ovaries . All sodium ion transmembrane transport-related genes were downregulated in the atretic follicles, especially in the plasma membrane, which caused the efflux of sodium ions and thus initiated the activation and shrinkage of apoptotic cells [18, 19], which may explain the obvious collapse in the atretic follicles of Yangzhou geese in this study. The caspase recruitment domain (CARD) 9 and 11 showed a significant increase in atretic follicles at the transcriptional level , while the mitochondrial proteins SmaC and HtrA2, which promoted cytochrome C-dependent caspase activation by the inhibition of IAP, were upregulated in atretic follicles . These studies demonstrated that apoptotic related genes were highly expressed to initiate the apoptosis of atretic follicles, which also need transcriptome and proteome levels changes to complete granulose cell apoptosis and follicle atresia.
Huet et al. believed that changes in ECM components could affect the apoptosis of follicular cells in sheep . Matrix metalloproteinase (MMP) participates in the disassembly of the extracellular matrix, and the upregulated proteins MMP3 and MMP9 could be involved in the collapse of the granulosa cell layer. As a secretion protein, MMP acted on growth factors and growth factor-binding proteins to affect various cellular functions, including cell invasion, cell migration, apoptosis, and angiogenesis . In the follicles of humans and rats, MMP promoted ovulation by destroying the follicle wall and the surrounding matrix through the stimulation of LH and HCG . After ovulation, the ovaries of poultry leave postovulatory follicles, which gradually degenerate and disappeared. Hrabia et al. reported that the increased activity of MMP2 and MMP9 promoted follicle regression in chicken postovulatory follicles . Zhu et al. also reported that the protein expression of MMP9 increased with the degree of regression of postovulatory follicles . The upregulation of MMP3 and MMP9 detected in this study in the atretic follicles of Yangzhou geese also suggested that these two genes were involved in the apoptosis of granulosa cells in AF, and might also promote the degeneration and disappearance of atretic follicles.
A complex composed of CDK and cyclin synergistically promoted the cell cycle, while the downregulation of CDK (CDK1 and CDK2) and cyclin (CyCD and CyCA) caused a decrease in the cell division activity in atretic follicles. The downregulation of Smad 2–3 and Smad 4 resulted in the interruption of cell growth or transformation. Tomic et al. demonstrated that the decreased expression of Smad3 resulted in the downregulation of CyCD2 and increased follicular atresia in mice . On the other hand, the increased expression of p300 and Ink4b caused a high level of acetylation of p53 and resulted in the increased expression of p21, which then bound to CDK and CyC to interrupt the cell cycle process. Regan et al. demonstrated that high LSHR was associated with a higher number of follicles, reduced apoptosis, and a higher ovulation rate in merino sheep . A study reported the reduction of ANDR in atretic porcine follicles , while another reported that the mRNA of ANDR could be regulated by FSH . FGF12, a growth factor, could promote the migration of endothelial cells, proliferation of smooth muscle cells, formation of new blood vessels, and the repair of damaged endothelial cells, which could act intracellularly to inhibit apoptosis [31, 32]. The downregulation of FGF12 could also be involved in the apoptosis of atretic follicles. The bFGF, FGF10, and FGF18 have been reported to be closely associated with follicular development and atresia [4,5,6]. In this experiment, the PPI of DEPs suggested that the downregulated genes LSHR, DHB1, CP19A, ANDR, and FGF12, which participated in the cell cycle pathway, could be the key factors regulating follicular atresia in Yangzhou geese.
In this study, there was weak weak correlation between transcriptome and proteome (R2 = 0.107). Zhao et al. showed that the low quantitative correlation between the transcriptome and proteome in wild-type and fads2-deletion zebrafish, with R2 values of 0.012 and 0.076, respectively . The low correlation between transcriptome and proteome results may be due to complex post-transcriptional regulation. Differentially expressed genes in both transcriptome and proteome mainly enriched in binding and extracellular region signal pathways. ANXA1 was a calcium-dependent phospholipid-binding protein that is expressed in almost all organs . It participated in membrane transport, exocytosis, the fusion of some membranes, signal transduction, cell proliferation, apoptosis, and ion channel formation . Zhu et al. reported that ANXA2 was involved in the angiogenesis of chicken ovarian follicles . The downregulation of VTN, an extracellular matrix protein, was the intermediate cause of decreased adhesion between granulosa cell layers. VTN could be degraded by MMPs , and its content was positively correlated with the follicular size in bovine . The decreased expression of VTN further suggested that it might be degraded in atretic follicle and caused the collapse of cell layer.
In this experiment, normal and atretic follicles from geese were compared using integrated transcriptome and proteome analysis. The cell cycle is an important pathway for the regulation of follicular atresia. Sodium outflow and high expression of MMP3 and MMP9 were important causes of structural destruction and apoptosis of follicular cells. The genes Smad2/3, Smad4, ANXA1, MMP3, ANXA1, CHK1, MCM3, CCNA2, MAD2, CDK1, FGF12, and CCND1, which were involved in the cell cycle, could be candidate genes for the regulation of follicular atresia.
Materials and methods
Experimental animals and sample collection
Six egg-laying geese (provided by Jiangsu Changzhou Four Seasons Poultry Industry Co., Ltd.) were selected for slaughter after being anesthetized with sodium pentobarbital. The follicles which displayed typical morphology of atresia were separated as the atretic follicles (AF). Normal developing follicles (NF) with similar size of AF were separated and set as the control. One part of the separated follicles was immediately frozen in liquid nitrogen and then stored at − 80 °C for hormone determination. In each follicle development stage, the large white follicles were selected as the object for transcriptome and proteome. From the second part, follicle granulose cell layer was immediately separated and stored for, RNA and protein extraction. The remaining part was stored in a tissue-fixative solution (4% paraformaldehyde) for the paraffin section.
ELISA and TUNEL assays
After 24 h of immersion in the tissue-fixative solution, the AF and NF were embedded in paraffin and sectioned into slices of 5 μm thickness. These slices were then subjected to TUNEL staining using the TUNEL Apoptosis Detection Kit (yeasen, China) according to the manufacturer’s instructions. The levels of FSH, LH, PRL, E2, and P4 hormones were measured using the Enzyme-Linked Immunosorbent assay (ELISA) according to the manufacturer’s instructions (Keshun tech, Shanghai, China). The microtiter plate was read at 450 nm using microtiter plater reader (SUNRISE, Astria). Twelve follicles in each of the NF and AF groups were used for the measurement of hormone levels. .
De novo transcriptome analysis
The total RNA was extracted from each follicle. Three extracted samples were equally mixed to form one replicate, and three replicates represented nine samples that were submitted for analysis. Total RNA was extracted using the Animal Tissue RNA Extraction Kit (Ambion, USA). The purity of the extracted RNA was assessed using the RNA Nano 6000 Assay Kit (Bioanalyzer 2100 system, Agilent Technologies, CA, USA). Purification and reverse transcription of the mRNA were performed using the NEBNext Ultra RNA Library Prep Kit (Illumina, NEB, USA). The cDNA fragment of 370 ~ 420 bp was selected for PCR and repurified using the AMPure XP system (Beckman Coulter, Beverly, USA). The quantity and quality of the final library was assessed using the Qubit 2.0 Fluorometer (Thermo Fisher, USA) and Agilent Bioanalyzer 2100 system (Agilent, USA), respectively.
The mRNA library was sequenced using Illumina NovaSeq 6000. Clean reads were obtained from raw data by removing reads containing the adapter, ploy-N, and low-quality. De novo assembly was performed using the software Trinity , and each cluster was defined as a “Transcript”. The “UniGenes” were obtained by removing redundant data from the “Transcripts” using the software Corset . The completeness of transcript assembly was assessed by the Benchmarking Universal Single-Copy Orthologs (BUSCO, http://busco.ezlab.org/).
Functional annotation of the UniGenes was performed using seven databases: NCBI non-redundant protein sequences (Nr), NCBI nucleotide sequences (Nt), Protein family (PFAM), eukaryotic ortholog groups (KOG), SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene ontology (GO). To calculate the expression of the UniGenes, clean reads were mapped to the the reference database assembled by Trinity using the software RSEM . The gene expression was quantified and normalized by the calculation method of “Reads per kilobase per million mapped reads” (RPKM) , and RPKM values > 0.3 represented that the gene was detected in follicle.
The total protein from each follicle was extracted using the acetone precipitation method. Protein samples were extracted from nine follicles each of normal and atretic follicles. Three samples were equallty mixed to form one replicate and three replicates were submitted for proteome analysis. The concentration and quality of the extracted protein were measured using the Bradford Protein Quantitative Kit (P0006, Beyotime, Shanghai, China) and 12% SDS-PAGE electrophoresis, respectively. After the protein samples were analyzed, digested, and desalinated, they were labeled using isobaric tags (iTRAQ) for relative and absolute quantification . Briefly, the protein samples were digested with trypsin overnight, and the pH was adjusted to around 3.0 using formic acid (FA). The samples were then centrifuged at 12000 g for 5 min and the supernatant was collected for desalting through the C18 desalting column with elution buffer of 0.1% FA and 70% acetonitrile. Each sample was reconstituted with 20 μL of 1 M TEAB buffer (Sigma, USA) and iTRAQ Reagent-8 plex kit (Sigma, USA) for 2 h, following which, 100 μL of 50 mM Tris-HCl (pH = 8) was added to terminate the reaction. The labeled samples were fractionated using a C18 Nano-Trap column (4.6 × 250 mm, 5 μm, Thermo Fisher, USA) on a Rigol L3000 UHPLC system and monitored at 214 nm. Finally, all fractions were dried under vacuum and reconstituted in 0.1% FA.
The liquid chromatography-tandem mass spectrometry (LC-MS/MS) was performed on an EASY-nLCTM 1200 UHPLC system (Thermo Fisher, Massachusetts, USA) coupled to a Q Exactive™ HF-X mass spectrometer (Thermo Fisher, Massachusetts, USA) and operated in a data-dependent acquisition mode . The fraction sample was injected into a C18 Nano-Trap column (4.5 cm × 75 μm, 3 μm, Thermo Fisher, USA). The peptides were separated on an analytical column (15 cm × 150 μm, 1.9 μm) using linear gradient elution for 60 min with a flow rate of 600 nL/min. The gradient elution included mobile phase A (0.1% FA) varying from 94 to 0% and mobile phase B (0.1% FA, 80% acetonitrile) varying from 6 to100%. A full MS scan ranging from 407 to 1500 m/z was conducted at 60,000 resolution with a target automatic gain control (AGC) value of 3 × 106 ions and a maximum ion injection time (IT) of 20 ms. The top 40 precursors of the highest abundance in the full scan were selected and fragmented using high-energy collisional dissociation (HCD). All MS/MS spectra were scanned using the following parameters: resolution = 15,000, AGC = 5 × 104 ions, maximum IT = 45 ms, dynamic exclusion duration = 20 s, normalized collision energy = 32%, and intensity threshold = 2.2 × 104.
The coding sequence (CDS) of the UniGenes from transcriptome was searched in the Nr and SwissProt protein databases. Once the CDS could not match to these two protein database, it would be predicted by the software ESTSCAN version 3.0.3. This combined database was the reference for the proteomic analysis. All spectra were searched against the transcriptome CDS prediction database using the software Proteome Discoverer 2.2 (PD 2.2, Thermo). The search parameters were set as follows: mass tolerance for precursor ion was 10 ppm, mass tolerance for production was 0.02 Da, the fixed modification was carbamidomethyl, dynamic modifications were Oxidation of methionine (M) and iTRAQ plex, N-Terminal modifications were acetylation and iTRAQ plex, and a maximum of two missed cleavage sites was allowed. To ensure the accuracy of the MS data analysis, the Peptide Spectrum Matches (PSMs), protein confidences, and false discovery rates (FDR) were set to > 99%, > 1.0, and ≤ 1.0%, respectively. The proteins were annotated using the program InterProscan against the GO, KEGG, IPR, and COG databases .
Total RNA was extracted from the goose follicle and reverse transcribed for cDNA synthesis using Hifair® II 1st Strand cDNA Synthesis Kit (Yeasen, China) according to the manufacturer’s protocol. Real-time quantitative PCR (qRT-PCR) was performed using the TransStart Green qPCR SuperMix Kit (TAKARA, China) on an ABI 7500 system (Thermo, USA). The primers were designed using the online platform Primer 3.0 (https://bioinfo.ut.ee/primer3–0.4.0/) (Table 4). The conditions of qRT-PCR were: 95 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 10s and annealing and extension at 60 °C for 20 s. The expression level was calculated by using the 2–ΔΔCT method.
Western blot analysis
Gel electrophoresis was performed using 12% SDS-PAGE, after which, the proteins were transferred onto PVDF membrane. The PVDF membranes were then incubated in 5% skim milk (diluted in TBST) for 2 h and then incubated in the diluted primary antibodies (rabbit anti-MMP3 with 1:1000 dilution, YT4465, ImmunoWay, Jiangsu, China; anti-MMP9 with 1:1000 dilution, YT1892, ImmunoWay, Jiangsu, China; anti-β-actin with 1:10000 dilution, ABclonal, Wuhan, China) at 4 °C overnight. After washing by tris-buffered saline three times, the PVDF membrane was incubated with HRP-labeled goat anti-rabbit IgG (1:30,000, BIOMIKY, Shanghai, China) for 60 min. Protein blots were visualized using chemiluminescence imaging system 398 (UVitec, Cambridge, United Kingdom) and quantified using Image J software (National Institutes of Health, USA).
Statistical and bioinformatics analysis
The differences in expression of genes and proteins between AF and NF were analyzed using the T-test in the statistical program R. P values were corrected using Benjamini and Hochberg. The genes with padj of ≤0.05 and fold change (FC) of > 2 were considered differentially expressed (DEGs). The proteins with padj of ≤0.05 and FC of > 1.3 were defined as differentially expressed proteins (DEPs). The relative mRNA expression through qRT-PCR was analyzed using student T-test in the statistical software SPSS version 22.0 (SPSS, Chicago, IL, USA). The histogram was made using the software GraphPad Prism version 8.0 (GraphPad, San Diego, USA). The GO and KEGG enrichment analyses were performed using the software GOseq and KOBAS [46, 47]. Protein-protein interaction (PPI) was predicted and visually edited using the software STRING DB and Cytoscape, respectively, with String database (http://string-db.org/) as the reference .
Availability of data and materials
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository with the dataset identifier PXD028836. The transcriptomics datasets generated during the current study are available in NCBI SRA (PRJNA767219, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/bioproject/PRJNA767219).
Enzyme linked immunosorbent assay
TDT-mediated dUTP Nick-end Labeling
Quantitative real-time PCR
isobaric tags for relative and absolute quantification
Principal component analysis
Protein−Protein Interaction Network
Tilly JL, Kowalski KI, Johnson AL, Hsueh AJ. Involvement of apoptosis in ovarian follicular atresia and postovulatory regression. Endocrinology. 1991;129(5):2799–801.
Tilly JL. Commuting the death sentence: how oocytes strive to survive. Nat Rev Mol Cell Biol. 2001;2(11):838–48.
Knight PG, Glister C. TGF-β superfamily members and ovarian follicle development. Reproduction. 2006;132(2):191–206.
Lin X, Ma Y, Qian T, Yao J, Mi Y, Zhang C. Basic fibroblast growth factor promotes prehierarchical follicle growth and yolk deposition in the chicken. Theriogenology. 2019;139:90–7.
Gasperin BG, Ferreira R, Rovani MT, Santos JT, Buratini J, Price CA, et al. FGF10 inhibits dominant follicle growth and estradiol secretion in vivo in cattle. Reproduction. 2012;143:815–23.
Portela VM, Dirandeh E, Guerrero-Netro HM, Zamberlam G, Barreta MH, Goetten AF, et al. The role of fibroblast growth factor-18 in follicular atresia in cattle. Biol Reprod. 2015;92(1):14, 11-18.
Hsueh AJ, Billig H, Tsafriri A. Ovarian follicle atresia: a hormonally controlled apoptotic process. Endocr Rev. 1994;15(6):707–24.
Guthrie H, Cooper B, Welch G, Zakaria A, Johnson L. Atresia in follicles grown after ovulation in the pig: measurement of increased apoptosis in granulosa cells and reduced follicular fluid estradiol-17β. Biol Reprod. 1995;52(4):920–7.
Miranda A, Bazzoli N, Rizzo E, Sato Y. Ovarian follicular atresia in two teleost species: a histological and ultrastructural study. Tissue Cell. 1999;31(5):480–8.
Palumbo A, Yeh J. In situ localization of apoptosis in the rat ovary during follicular atresia. Biol Reprod. 1994;51(5):888–95.
Giebel J, Hegele-Hartung C, Rune GM. Proliferation and apoptosis in follicles of the marmoset monkey (Callithrix jacchus) ovary. Ann Anat. 1997;179(5):413–9.
Matsuda F, Inoue N, Manabe N, Ohkura S. Follicular growth and atresia in mammalian ovaries: regulation by survival and death of granulosa cells. J Reprod Dev. 2012;58(1):44–50.
Suikkari A-M, Tulppala M, Tuuri T, Hovatta O, Barnes F. Luteal phase start of low-dose FSH priming of follicles results in an efficient recovery, maturation and fertilization of immature human oocytes. Hum Reprod. 2000;15(4):747–51.
Ilha GF, Rovani MT, Gasperin BG, Antoniazzi AQ, Goncalves P, Bordignon V, et al. Lack of FSH support enhances LIF-STAT3 signaling in granulosa cells of atretic follicles in cattle. Reproduction. 2015;150(4):395–403.
Feranil JB, Isobe N, Nakao T. Changes in the thecal vasculature during follicular atresia in the ovary of swamp buffalo. J Reprod Develop. 2004;50(3):315–21.
Garside SA, Harlow CR, Hillier SG, Fraser HM, Thomas FH. Thrombospondin-1 inhibits angiogenesis and promotes follicular atresia in a novel in vitro angiogenesis assay. Endocrinology. 2010;151(3):1280–9.
Smith PR, Quirke L, Juengel JL, Hurst PR. Expression of PDGFs and their receptors in the fetal and adult sheep ovary. Biol Reprod. 2010;83:651.
McCarthy JV, Cotter TG. Cell shrinkage and apoptosis: a role for potassium and sodium ion efflux. Cell Death Differ. 1997;4(8):756–70.
Bortner CD, Hughes FM, Cidlowski JA. A primary role for K+ and Na+ efflux in the activation of apoptosis. J Biol Chem. 1997;272(51):32436–42.
Bertin J, Guo Y, Wang L, Srinivasula SM, Jacobson MD, Poyet J-L, et al. CARD9 is a novel caspase recruitment domain-containing protein that interacts with BCL10/CLAP and activates NF-κB. J Biol Chem. 2000;275(52):41082–6.
Du C, Fang M, Li Y, Li L, Wang X. Smac, a mitochondrial protein that promotes cytochrome c–dependent caspase activation by eliminating IAP inhibition. Cell. 2000;102(1):33–42.
Huet C, Monget P, Pisselet C, Monniaux D. Changes in extracellular matrix components and steroidogenic enzymes during growth and atresia of antral ovarian follicles in the sheep. Biol Reprod. 1997;56(4):1025–34.
Puttabyatappa M, Jacot TA, Al-Alem LF, Rosewell KL, Duffy DM, Brännström M, et al. Ovarian membrane-type matrix metalloproteinases: induction of MMP14 and MMP16 during the periovulatory period in the rat, macaque, and human. Biol Reprod. 2014;91(2):34–12.
McCord LA, Li F, Rosewell KL, Brännström M, Curry TE. Ovarian expression and regulation of the stromelysins during the periovulatory period in the human and the rat. Biol Reprod. 2012;86(3):78–1.
Hrabia A, Socha JK, Sechman A. Involvement of matrix metalloproteinases (MMP-2,-7,-9) and their tissue inhibitors (TIMP-2,-3) in the regression of chicken postovulatory follicles. Gen Comp Endocrinol. 2018;260:32–40.
Zhu G, Kang L, Wei Q, Cui X, Wang S, Chen Y, et al. Expression and regulation of MMP1, MMP3, and MMP9 in the chicken ovary in response to gonadotropins, sex hormones, and TGFB1. Biol Reprod. 2014;90(3):57–1.
Tomic D, Miller KP, Kenny HA, Woodruff TK, Hoyer P, Flaws JA. Ovarian follicle development requires Smad3. Mol Endocrinol. 2004;18(9):2224–40.
Regan S, McFarlane JR, O’Shea T, Andronicos N, Arfuso F, Dharmarajan A, et al. Flow cytometric analysis of FSHR, BMRR1B, LHR and apoptosis in granulosa cells and ovulation rate in merino sheep. Reproduction. 2015;150(2):151–63.
Słomczyńska M, Tabarowski Z, Duda M, Burek M, Knapczyk K. Androgen receptor in early apoptotic follicles in the porcine ovary at pregnancy. Folia Histochem Cytobiol. 2006;44(3):185–8.
Weil S, Vendola K, Zhou J, Bondy CA. Androgen and follicle-stimulating hormone interactions in primate ovarian follicle development. J Clin Endocrinol Metab. 1999;84(8):2951–6.
Nakayama F, Müller K, Hagiwara A, Ridi R, Akashi M, Meineke V. Involvement of intracellular expression of FGF12 in radiation-induced apoptosis in mast cells. J Radiat Res. 2008;49(5):491–501.
Nakayama F, Yasuda T, Umeda S, Asada M, Imamura T, Meineke V, et al. Fibroblast growth factor-12 (FGF12) translocation into intestinal epithelial cells is dependent on a novel cell-penetrating peptide domain: involvement of internalization in the in vivo role of exogenous FGF12. J Biol Chem. 2011;286(29):25823–34.
Zhao Y, Yang G, Wu N, Cao X, Gao J. Integrated transcriptome and phosphoproteome analyses reveal that fads2 is critical for maintaining body LC-PUFA homeostasis. J Proteome. 2020;229:103967.
Patel DM, Ahmad SF, Weiss DG, Gerke V, Kuznetsov SA. Annexin A1 is a new functional linker between actin filaments and phagosomes during phagocytosis. J Cell Sci. 2011;124(4):578–88.
Iseki Y, Imoto A, Okazaki T, Harigae H, Takahashi S. Identification of annexin 1 as a PU. 1 target gene in leukemia cells. Leuk Res. 2009;33(12):1658–63.
Zhu G, Chen X, Mao Y, Kang L, Ma X, Jiang Y. Characterization of annexin A2 in chicken follicle development: evidence for its involvement in angiogenesis. Anim Reprod Sci. 2015;161:104–11.
Imai K, Shikata H, Okada Y. Degradation of vitronectin by matrix metalloproteinases-1,-2,-3,-7 and-9. FEBS Lett. 1995;369(2–3):249–51.
Peter A, Perrone M, Asem E. Bovine ovarian follicular fluid vitronectin content is influenced by the follicle size. Theriogenology. 1995;43(7):1239–47.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Davidson NM, Oshlack A. Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biol. 2014;15(7):410.
Dewey CN, Li B. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12(1):323.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Wu J, An Y, Pu H, Shan Y, Ren X, An M, et al. Enrichment of serum low-molecular-weight proteins using C18 absorbent under urea/dithiothreitol denatured environment. Anal Biochem. 2010;398(1):34–44.
Wu J, Xie X, Liu Y, He J, Benitez R, Buckanovich RJ, et al. Identification and confirmation of differentially expressed fucosylated glycoproteins in the serum of ovarian cancer patients using a lectin array and LC–MS/MS. J Proteome Res. 2012;11(9):4541–52.
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2012;41(D1):D808–15.
This work was financially supported by Graduate Innovation Fund of Anhui Agricultural University (2021yjs-19), National Innovative Training Program for College Student (XJDC2020576) and Science and Technology Program of Xinjiang Uygur Autonomous Region(2020B01004–2-1).
Ethics approval and consent to participate
The animal experiment was reviewed and approved by the Institutional Animal Care and Use Committee of the Anhui Agricultural University (no.SYDW-P20200600601). The experiments were performed according to the Regulations for the Administration of Affairs Concerning Experimental Animals and the Standards for the Administration of Experimental Practices, as well as the ARRIVE guidelines version 2.0.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Fig. S1.
Principal component analysis of DEGs in normal and atretic follicles.
Additional file 2: Fig. S2.
(a) The quality control of mass spectrometry, including peptide length, (b) precursor ion tolerance, (c) protein mass distribution, and (d) protein coverage.
Additional file 3: Fig. S3.
(a) Western blotting of MMP3 protein in normal and atretic follicles. The remaining three bands, also normal follicles, were cropped out of the manuscript picture. (b) Western blotting of MMP9 protein in normal and atretic follicles. (c) Western blotting of ACTIN protein. AF-1, AF-2 and AF-3 represent the three repeats of atreated follicles, while NF-1, NF-2 and NF-3 represent the three repeats of normal follicles. The red area is the cropped area in the manuscript. The PVDF membrane was cut before incubation with antibodies.
Additional file 4: Table S1.
Summary of transcriptome sequencing quality.
Additional file 5: Table S2.
Summary of Transcripts and Unigene.
Additional file 6: Table S3.
Summary of the functional annotation of assembled unigenes.
Additional file 7: Table S4.
The FPKM value of all genes in the goose follicles.
Additional file 8: Table S5.
The differentially expressed genes between normal and atretic follicles.
Additional file 9: Table S6.
GO enrichment of differentially expressed genes between normal and atretic follicles (molecular_function).
Additional file 10: Table S7.
KEGG enrichment of differentially expressed genes between normal and atretic follicles.
Additional file 11: Table S8.
The expression of all proteins in the goose follicles.
Additional file 12: Table S9.
The differentially expressed proteins between normal and atretic follicles.
Additional file 13: Table S10.
GO enrichment of differentially expressed proteins between normal and atretic follicles (p < 0.05).
Additional file 14: Table S11.
The expressions of mRNAs and proteins that corresponds to each other in the goose follicles.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yang, W., Chen, X., Liu, Z. et al. Integrated transcriptome and proteome revealed that the declined expression of cell cycle-related genes associated with follicular atresia in geese. BMC Genomics 24, 24 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-022-09088-1
- Atretic follicle
- Cell cycle
- Matrix metalloproteinase