AFTGAN: prediction of multi-type PPI based on attention free transformer and graph attention network.

MOTIVATION: Protein-protein interaction (PPI) networks and transcriptional regulatory networks are critical in regulating cells and their signaling. A thorough understanding of PPIs can provide more insights into cellular physiology at normal and disease states. Although numerous methods have been proposed to predict PPIs, it is still challenging for interaction prediction between unknown proteins. In this study, a novel neural network named AFTGAN was constructed to predict multi-type PPIs. Regarding feature input, ESM-1b embedding containing much biological information for proteins was added as a protein sequence feature besides amino acid co-occurrence similarity and one-hot coding. An ensemble network was also constructed based on a transformer encoder containing an AFT module (performing the weight operation on vital protein sequence feature information) and graph attention network (extracting the relational features of protein pairs) for the part of the network framework. RESULTS: The experimental results showed that the Micro-F1 of the AFTGAN based on three partitioning schemes (BFS, DFS and the random mode) on the SHS27K and SHS148K datasets was 0.685, 0.711 and 0.867, as well as 0.745, 0.819 and 0.920, respectively, all higher than that of other popular methods. In addition, the experimental comparisons confirmed the performance superiority of the proposed model for predicting PPIs of unknown proteins on the STRING dataset. AVAILABILITY AND IMPLEMENTATION: The source code is publicly available at https://github.com/1075793472/AFTGAN. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of chicken Relaxin3 gene: mRNA expression and response to reproductive hormone treatment in ovarian granulosa cells, and single nucleotide polymorphisms associated with egg laying traits in hens.

As a protein structurally similar to insulin, relaxin3 (RLN3) plays a role in promoting arousal, suppressing depressive or anxious behaviors. Two studies revealed the increase of RLN3 expression during chicken follicle selection. In this study, by real-time quantitative PCR and luciferase assay, mRNA expression and single nucleotide polymorphisms (SNPs) of chicken RLN3 were investigated. The mRNA expression of chicken RLN3 was higher in the granulosa cell of hierarchal follicles (Post-GCs) than that of pre-hierarchal follicles (Pre-GCs). In Pre-GCs, the mRNA expression of chicken RLN3 was stimulated by FSH and progesterone; in Post-GCs, it was stimulated by higher concentration of estrogen and FSH, however, was inhibited by progesterone. Four SNPs including g.-655G > C, g-592G > A, g.-372T > A and g.-282G > C were identified in the critical promoter region from -1291 bp to -207 bp of chicken RLN3, among which g.-655G > C, and g-592G > A were associated with age at first laying and clutch size, respectively, in Zaozhuang Sunzhi chickens. At g.-655G > C and g-592G > A, allele C and allele A had higher transcriptional activity, respectively. These data suggest that RLN3 plays an important role in chicken follicle development and SNPs in its promoter region are potential DNA markers for improving egg production traits.

Phosphoproteomic analysis on ovarian follicles reveals the involvement of LSD1 phosphorylation in Chicken follicle selection.

BACKGROUND: Follicle selection in chickens refers to the process of selecting one follicle from a group of small yellow follicles (SY, 6-8 mm in diameter) for development into 12-15 mm hierarchal follicles (usually F6 follicles), which is controlled by sex hormones including follicle-stimulating factor (FSH), estrogen and progesterone. Follicle selection is a critical process impacting egg production in chicken, therefore, is the focus of many studies. Phosphorylation is important for the proper function of proteins, thus, needs to be analyzed by proteomic level. RESULT: In this study, we compared the phosphoproteomes of SY and F6 follicles in laying hens and identified 2,386 phosphoproteins and 5,940 phosphosites, of which 4,235 sites of 1,963 phosphoproteins were quantified. From SY to F6 follicles, 190 phosphorylation sites of 149 proteins changed significantly, among which the phosphorylation level of lysine demethylase 1 A (LSD1) at the conserved 54th serine (LSD1Ser54p) was significantly upregulated in F6 follicles compared to SY follicles (p < 0.05); however, the expression of chicken LSD1 were not significantly different on both mRNA and protein levels. LSD1Ser54p is mainly located in the nucleus of both SY and F6 follicles, and was higher in F6 follicles than that of SY follicles revealed by both immunofluorescence and Western blotting. LSD1Ser54p level increased after treatment with 5 ng/mL and 10 ng/mL of FSH in the theca cells and the granulosa cells of pre-hierarchal follicles, and with 50 ng/mL in granulosa cells of hierarchal follicles. In the theca cells of hierarchal follicles, estrogen stimulated the level of LSD1Ser54p in a dosage-dependent manner, and in granulosa cells of pre-hierarchal follicles, 10 ng/mL of estrogen increased LSD1Ser54p expression. Treatment with 50 ng/mL of progesterone increased LSD1Ser54p expression in theca cells of pre-hierarchal follicles, and with 10 to 100 ng/ml enhanced LSD1Ser54p expression in the granulosa cells of hierarchal follicles. CONCLUSION: The expression dynamics of LSD1Ser54p in follicles from SY to F6 and its regulation by sex hormones suggest that it is involved in chicken follicle selection.

Identification of nucleotide polymorphisms in the key promoter region of chicken annexins A2 gene associatied with egg laying traits.

Annexin A2 (ANXA2) is a member of the A subfamily of a multifunctional calcium dependent membrane phospholipid binding protein family. The mRNA expression of ANXA2 is consistent with ovary function and egg laying in chickens. In this study, six nucleotide polymorphisms in the key promoter region of chicken ANXA2 gene (-2861 bp to -1394 bp), i.e.,: g.-2337 indel (GT), g.-2255 C > T, g. -2248 A > G, g.-2188 A > G, g.-2169 G > A, g.-2160 A > C, were identified. Their distributions in populations of Xinyang Brown, Recessive White Rock, Wenchang and Wenshang Barred chickens were analyzed. In the Recessive White Rock chicken population, CAA, CAG and TGG were three major haplotypes. Association analysis indicated that the individuals with diplotype TGG/TGG laid more eggs at 32 weeks, and the individual with diplotype CAG/TGG laid at the earlier age. Luciferase activity assay showed that mutation from C to T at -2255 increased trascriptional activity of chicken ANXA2, which is consistent with its effect on egg laying traits.

Identification of a promoter polymorphism affecting GPAT3 gene expression that is likely related to intramuscular fat content in pigs.

The intramuscular fat content (IMF) is an economically important trait in pigs and the Laiwu pig is famous for its excessively extremely high level of IMF. Our previous transcriptome study revealed that the dynamic expression of glycerol-phosphate acyltransferase 3 (GPAT3) is consistent with changes in the IMF of Laiwu pigs. In this study, we further analyzed the expression and polymorphism of GPAT3 in its promoter region. The results indicated that the expression of GPAT3 increased dramatically from 120 to 240 days and is consistent with changes in IMF deposition, and at both mRNA and protein levels, GPAT3 expression was markedly higher in the LD muscle of Laiwu pigs than that of Duroc × Landrace × Yorkshire pigs. Deletion from -1695 to -1187 of porcine GPAT3 greatly increased its transcription. Of the two SNPs identified, the transition from C to T at -1526 site increased the transcription of porcine GPAT3 and allele T mainly distributed in Laiwu pig population. These results collectively suggest that the SNP at -1526 site of GPAT3 may contribute to IMF deposition by affecting its expression in pigs.

Characterization of the chicken T cell receptor γ repertoire by high-throughput sequencing.

BACKGROUND: As one of "γδ-high" species, chicken is an excellent model for the study of γδ T cells in non-mammalian animals. However, a comprehensive characterization of the TCRγδ repertoire is still missing in chicken. The objective of this study was to characterize the expressed TCRγ repertoire in chicken thymus using high-throughput sequencing. METHODS: In this study, we first obtained the detailed genomic organization of the TCRγ locus of chicken based on the latest assembly of the red jungle fowl genome sequences (GRCg6a) and then characterized the TCRγ repertoire in the thymus of four chickens by using 5' Rapid Amplification of cDNA Ends (5' RACE) along with high-throughput sequencing (HTS). RESULTS: The chicken TCRγ locus contains a single Cγ gene, three functional Jγ segments and 44 Vγ segments that could be classified into six subgroups, each containing six, nineteen, nine, four, three and three members. Dot-plot analysis of the chicken TCRγ locus against itself showed that almost all the entire zone containing Vγ segments had arisen through tandem duplication events, and the main homology unit, containing 9 or 10 Vγ gene segments, has tandemly duplicated for four times. For the analysis of chicken TCRγ repertoire, more than 100,000 unique Vγ-region nucleotide sequences were obtained from the thymus of each chicken. After alignment to the germline Vγ and Jγ segments identified above, we found that the four chickens had similar repertoire profile of TCRγ. In brief, four Vγ segments (including Vγ3.7, Vγ2.13, Vγ1.6 and Vγ1.3) and six Vγ-Jγ pairs (including Vγ3.7-Jγ3, Vγ2.13-Jγ1, Vγ2.13-Jγ3, Vγ1.6-Jγ3, Vγ3.7-Jγ1 and Vγ1.6-Jγ1) were preferentially utilized by all four individuals, and vast majority of the unique CDR3γ sequences encoded 4 to 22 amino acids with mean 12.90 amino acids, which exhibits a wider length distribution and/or a longer mean length than CDR3γ of human, mice and other animal species. CONCLUSIONS: In this study, we present the first in-depth characterization of the TCRγ repertoire in chicken thymus. We believe that these data will facilitate the studies of adaptive immunology in birds.

Transcriptomic and proteomic analyses of ovarian follicles reveal the role of VLDLR in chicken follicle selection.

BACKGROUND: Follicle selection in chickens refers to the process of selecting one follicle from a group of small yellow follicles (SY, 6-8 mm in diameter) for development into 12-15 mm hierarchical follicles (usually F6 follicles), which is an important process affecting laying performance in the poultry industry. Although transcriptomic analysis of chicken ovarian follicles has been reported, integrated analysis of chicken follicles for selection by using both transcriptomic and proteomic approaches is still rarely performed. In this study, we compared the proteomes and transcriptomes of SY and F6 follicles in laying hens and identified several genes involved in chicken follicle selection. RESULTS: Transcriptomic analysis revealed 855 differentially expressed genes (DEGs) between SY follicles and F6 follicles in laying hens, among which 202 were upregulated and 653 were downregulated. Proteomic analysis revealed 259 differentially expressed proteins (DEPs), including 175 upregulated and 84 downregulated proteins. Among the identified DEGs and DEPs, changes in the expression of seven genes, including VLDLR1, WIF1, NGFR, AMH, BMP15, GDF6 and MMP13, and nine proteins, including VLDLR, VTG1, VTG3, PSCA, APOB, APOV1, F10, ZP2 and ZP3L2, were validated. Further analysis indicated that the mRNA level of chicken VLDLR was higher in F6 follicles than in SY follicles and was also higher in granulosa cells (GCs) than in thecal cells (TCs), and it was stimulated by FSH in GCs. CONCLUSIONS: By comparing the proteomes and transcriptomes of SY and F6 follicles in laying hens, we identified several differentially expressed proteins/genes that might play certain roles in chicken follicle selection. These data may contribute to the identification of functional genes and proteins involved in chicken follicle selection.

Lin28 gene and mammalian puberty.

Lin28a and Lin28b, homologs of the Caenorhabditis elegans Lin28 gene, play important roles in cell pluripotency, reprogramming, and tumorigenicity. Recently, genome-wide association and transgenic studies showed that Lin28a and/or Lin28b gene were involved in the onset of mammalian puberty, the stage representing the attainment of reproduction capacity; however, the detailed mechanism of these genes in mammalian puberty remains largely unknown. The present paper reviews the research progress on the roles of Lin28a/b genes in the onset of mammalian puberty by analyzing the results coming from gene expression patterns, mutations, and transgenic studies, and put forward possible pathways for further studies on their roles in animal reproduction.

Suppression of lymphocyte apoptosis in spleen by CXCL13 after porcine circovirus type 2 infection and regulatory mechanism of CXCL13 expression in pigs.

Porcine circovirus-associated disease (PCVAD) is one of the most serious infectious diseases in pigs worldwide. The primary causative agent of PCVAD is porcine circovirus type 2 (PCV2), which can cause lymphoid depletion and immunosuppression in pigs. Our previous study demonstrated that Laiwu (LW) pigs, a Chinese indigenous pig breed, have stronger resistance to PCV2 infection than Yorkshire × Landrace (YL) pigs. In this study, we found that the YL pigs showed more severe lymphocyte apoptosis and higher viral load in the spleen tissue than LW pigs. To illustrate the differential gene expression between healthy and infected spleens, transcriptome profiling of spleen tissues from PCV2-infected and control YL pigs was compared by RNA sequencing. A total of 90 differentially expressed genes (DEGs) was identified, including CD207, RSAD2, OAS1, OAS2, MX2, ADRB3, CXCL13, CCR1, and ADRA2C, which were significantly enriched in gene ontology (GO) terms related to the defense response to virus and cell-cell signaling, and another nine DEGs, KLF11, HGF, PTGES3, MAP3K11, XDH, CYCS, ACTC1, HSPH1, and RYR2, which were enriched in GO terms related to regulation of cell proliferation or apoptosis. Among these DEGs, the CXCL13 gene, which can suppress lymphocyte apoptosis during PCV2 infection, was significantly down-regulated in response to PCV2 infection in YL but not in LW pigs. By analysis of the regulatory elements in the promoter and 3'-untranslated region (3'-UTR) of porcine CXCL13, we found that the single nucleotide polymorphism (SNP) -1014 G (LW) > A (YL) and the Sus scrofa microRNA-296-5p (ssc-miR-296-5p) participated in regulating CXCL13 expression during the response to PCV2 infection.

Establishment and evaluation of a PRRSV-sensitive porcine endometrial epithelial cell line by transfecting SV40 large T antigen.

BACKGROUND: PRRSV is an infectious illness causing lung injury and abortion in sows. Cells apoptosis in the interface between the endometrium and fetal placenta is a crucial factor causing abortion. Previous study confirmed PRRSV could cause apoptosis of macrophages but rarely produced an obvious change in porcine endometrial epithelial cells (PECs). Recently, PRRSV-induced abortion was attributed to fetal placental and endometrium epithelial cells (Sn+ and CD163+) apoptosis. However, the mechanism of abortion is still unrevealed because of the limit of porcine endometrium epithelial cells (PEC). The aim of this study was to establish a stable immortalized PECs lines and use it to reveal the abortion mechanism. RESULTS: In this study, highly purified primary PECs were harvested through differential digestion, and their characteristics were confirmed by CK18, ERɑ and PR staining. Cells were then immortalized by transfecting a lentiviral vector that expressed SV40 large T antigen. PECs lines were obtained after puromycin screening. Proliferation of cell line was evaluated by cell growth curve and cell cycle assays. Cell lines exhibited faster proliferation capacity than primary cells. Biological characteristics of cell line were assessed by Western blot, karyotype analysis and staining, which confirmed that the cell line retained the endometrium characteristics. Finally, PRRSV sensitivity was assessed; expression of Sn and CD163 indicated that primary PECs and cell lines were all potentially sensitive to PRRSV. PRRSV infection tests showed an obvious increase in apoptotic rate in the infected PEC cell line, which suggested its susceptibility. CONCLUSION: The newly constructed cell line is a useful tool for studying the mechanism of abortion caused by PRRSV.

Expression dynamics of gonadotropin-releasing hormone-I and its mutual regulation with luteinizing hormone in chicken ovary and follicles.

Gonadotropin-releasing hormone-I (GnRH-I) has been identified in the ovaries of vertebrate species, and this decapeptide is a key regulator of reproductive functions. However, its biological action and regulatory mechanism in the chicken ovary remain to be characterized. In this study, the expression of GnRH-I gene in chicken hypothalamus and ovaries at different developmental stages and different sizes of follicles was investigated, and the effect of GnRH-I mRNA on chicken follicular cells was analyzed in vitro. The results showed that the expression of GnRH-I was dramatically decreased in the hen ovary compared to that in the hypothalamus after sexual maturation. In the mature ovarian follicles, GnRH-I mRNA levels were significantly higher in theca cells than that in granulosa cells. Overexpression of GnRH-I decreased the expression of luteinizing hormone receptor (LHR) mRNA in theca cells from preovulatory follicles but had no effect on granulosa cells. Treatment of theca cells with different concentrations of luteinizing hormone (LH) significantly increased GnRH-I mRNA expression at low doses (50 ng/ml) but significantly decreased it at higher doses (200 ng/ml). Furthermore, GnRH-I inhibited LH-induced LHR expression at the lower dose of LH (50 ng/ml). These findings provide strong evidence indicating that GnRH-I is an important regulator in the chicken ovary.

Genome-wide survey reveals dynamic effects of folate supplement on DNA methylation and gene expression during C2C12 differentiation.

Folic acid supplements taken during pregnancy can prevent neural tube defects and other developmental abnormalities. Here, we explored the effects of folate supplementation on gene expression and DNA methylation during C2C12 differentiation. Based on the folic acid concentration, this study comprised three groups: low folate (L), normal folate (N), and high-folate supplement (H). Our analyses revealed that differentiation and the mRNA expression of the gene myogenin in C2C12 cell were enhanced by folic acid; however, the overall methylation percentage in myogenin promoter between different treatment groups was not significantly different ( P > 0.05). The results of MeDIP-chip showed that hundreds of differentially methylated regions (DMRs) were identified between every two groups in both promoter and CpG islands, respectively. Genes with DMRs between N and L groups were mainly enriched in the processes of cell differentiation and cell development, whereas those with DMRs between H and N groups were frequently enriched in cellular process/cycle and cell metabolic processes. In addition, correlation analysis between methylation profile and expression profile revealed that some genes were regulated by methylation status directly. Together, these analyses suggest that folate deficiency and supplementation can influence the differentiation, genome-wide DNA methylation level and the expression of myogenesis-related genes including myogenin in the C2C12 cell line.

Identification and characterization of microRNA in the lung tissue of pigs with different susceptibilities to PCV2 infection.

Porcine circovirus type 2 (PCV2) is the primary cause of post-weaning multisystemic wasting syndrome (PMWS) and other PCV-associated diseases. According to our previous RNA-sequencing analysis, the differences in the susceptibility to PCV2 infection depended on the genetic differences between the Laiwu (LW) and Yorkshire × Landrace crossbred (YL) pigs, but the cellular microRNA (miRNA) that are differentially expressed between the LW and YL pigs before and after PCV2 infection remain to be determined. In this study, high-throughput sequencing was performed to determine the abundance and differential expression of miRNA in lung tissues from PCV2-infected and PCV2-uninfected LW and YL pigs. In total, 295 known and 95 novel miRNA were identified, and 23 known and 25 novel miRNA were significantly differentially expressed in the PCV2-infected vs. PCV2-uninfected LW pigs and/or the PCV2-infected vs. PCV2-uninfected YL pigs. The expression levels of ssc-miR-122, ssc-miR-192, ssc-miR-451, ssc-miR-486, and ssc-miR-504 were confirmed by quantitative real-time PCR (qRT-PCR). Analysis of the potential targets of the four up-regulated miRNA (i.e., ssc-miR-122, ssc-miR-192, ssc-miR-451 and ssc-miR-486) identified pathways and genes that may be important for disease resistance. Among the up-regulated miRNA, ssc-miR-122 can repress the protein expression and viral DNA replication of PCV2 and down-regulate the expression of the nuclear factor of activated T-cells 5 (NFAT5) and aminopeptidase puromycin sensitive (NPEPPS) by binding to their 3' untranslated region (3'UTR) in PK15 cells. Therefore, ssc-miR-122 may indirectly suppress PCV2 infection by targeting genes related to the host immune system, such as NFAT5 and NPEPPS.

Dynamic transcriptome and DNA methylome analyses on longissimus dorsi to identify genes underlying intramuscular fat content in pigs.

BACKGROUND: The intramuscular fat content (IMF) refers to the amount of fat within muscles, including the sum of phospholipids mainly found in cell membranes, triglycerides and cholesterol, and is determined both by hyperplasia and hypertrophy of adipocyte during the development of pigs. The IMF content is an important economic trait that is genetically controlled by multiple genes. The Laiwu pig is an indigenous fatty pig breed distributed in North China, characterized by excessively higher level of IMF content (9%~12%), therefore, is suitable for the identification of genes controlling IMF variations. To identify genes underlying IMF deposition, we performed genome-wide transcriptome and methylome analyses on longissimus dorsi (LD) muscle in Laiwu pigs across four developmental stages. RESULTS: A total of 22,524 expressed genes were detected and 1158 differentially expressed genes (DEGs) were hierarchically clustered in the LD muscle over four developmental stages from 60 d to 400 d. These genes were significantly clustered into four temporal expression profiles, and genes participating in fat cell differentiation and lipid biosynthesis processes were identified. From 120 d to 240 d, the period with the maximum IMF deposition rate, the lipid biosynthesis related genes (FOSL1, FAM213B and G0S2), transcription factors (TFs) (EGR1, KLF5, SREBF2, TP53 and TWIST1) and enriched pathways (steroid biosynthesis and fatty acid biosynthesis) were revealed; and fat biosynthesis relevant genes showing differences in DNA methylation in gene body or intergenic region were detected, such as FASN, PVALB, ID2, SH3PXD2B and EGR1. CONCLUSIONS: This study provides a comprehensive landscape of transcriptome of the LD muscle in Laiwu pigs ranging from 60 to 400 days old, and methylome of the LD muscle in 120 d and 240 d Laiwu pigs. A set of candidate genes and TFs involved in fat biosynthesis process were identified, which were probably responsible for IMF deposition. The results from this study would provide a reference for the identification of genes controlling IMF variation, and for exploring molecular mechanisms underlying IMF deposition in pigs.

Prospective stratification of patients at risk for emergency department revisit: resource utilization and population management strategy implications.

BACKGROUND: Estimating patient risk of future emergency department (ED) revisits can guide the allocation of resources, e.g. local primary care and/or specialty, to better manage ED high utilization patient populations and thereby improve patient life qualities. METHODS: We set to develop and validate a method to estimate patient ED revisit risk in the subsequent 6 months from an ED discharge date. An ensemble decision-tree-based model with Electronic Medical Record (EMR) encounter data from HealthInfoNet (HIN), Maine's Health Information Exchange (HIE), was developed and validated, assessing patient risk for a subsequent 6 month return ED visit based on the ED encounter-associated demographic and EMR clinical history data. A retrospective cohort of 293,461 ED encounters that occurred between January 1, 2012 and December 31, 2012, was assembled with the associated patients' 1-year clinical histories before the ED discharge date, for model training and calibration purposes. To validate, a prospective cohort of 193,886 ED encounters that occurred between January 1, 2013 and June 30, 2013 was constructed. RESULTS: Statistical learning that was utilized to construct the prediction model identified 152 variables that included the following data domains: demographics groups (12), different encounter history (104), care facilities (12), primary and secondary diagnoses (10), primary and secondary procedures (2), chronic disease condition (1), laboratory test results (2), and outpatient prescription medications (9). The c-statistics for the retrospective and prospective cohorts were 0.742 and 0.730 respectively. Total medical expense and ED utilization by risk score 6 months after the discharge were analyzed. Cluster analysis identified discrete subpopulations of high-risk patients with distinctive resource utilization patterns, suggesting the need for diversified care management strategies. CONCLUSIONS: Integration of our method into the HIN secure statewide data system in real time prospectively validated its performance. It promises to provide increased opportunity for high ED utilization identification, and optimized resource and population management.

Identification of a short interspersed repetitive element insertion polymorphism in the porcine MX1 promoter associated with resistance to porcine reproductive and respiratory syndrome virus infection.

The myxovirus resistance (Mx) proteins belong to the dynamin superfamily and are important for innate host defence against RNA viruses. In this study, we demonstrate that positive elements are present in the two promoter regions of -2713 to -2565 and -688 to -431 in the porcine MX1 gene. Sequencing and alignment of the amplified porcine MX1 gene promoter region identified a short interspersed repetitive element (SINE) insertion of 275 bp at site -547. At this site, allele B (an insertion of 275 bp) is dominant in Chinese indigenous pig breeds but has a workable minor allele frequency in western lean-type pig breeds. Luciferase activity was compared between promoters with and without the insertion of the 275-bp fragment in transiently transfected MARC-145 cells. The insertion of the 275-bp fragment increased the luciferase activity significantly (P < 0.05) both prior to and post-porcine reproductive and respiratory syndrome (PRRS) virus inoculation. These results suggest that the SINE insertion polymorphism at site -547 of the MX1 gene promoter region is a potential DNA marker for PRRS resistance in pigs.

Expression and regulation of MMP1, MMP3, and MMP9 in the chicken ovary in response to gonadotropins, sex hormones, and TGFB1.

Matrix metalloproteinases (MMPs) are a specific class of proteolytic enzymes that play critical roles in follicular development and luteinization in mammals. However, the role of MMPs in avian ovary remains largely unknown. We found that three MMP genes (MMP1, MMP3, and MMP9) were significantly up-regulated in 23-wk-old (laying phase) chicken ovaries compared with 6-wk-old ovaries (prepubertal phase). In reproductively active chicken ovary, MMP1 expression (both mRNA and protein) remained low in prehierarchical and preovulatory follicles but increased in postovulatory follicles (POFs). Both MMP3 and MMP9 expression levels increased during follicular maturation. MMP3 reached maximal expression in the first largest follicle (F1), while MMP9 levels continued to rise in POF1 and POF2 after ovulation. Immunohistochemistry, Western blot analysis, and zymography experiments indicated that MMP1, MMP3, and MMP9 were synthesized and secreted by granulosa cells of different follicles in the chicken ovary. The mRNA expression of MMP1 and MMP3 in the granulosa cells was stimulated by follicle-stimulating hormone, luteinizing hormone, progesterone, and estrogen but not by transforming growth factor beta 1 (TGFB1). However, the mRNA of MMP9 was induced by TGFB1 but not follicle-stimulating hormone, luteinizing hormone, progesterone, or estrogen. Luciferase reporter and mutagenesis analysis indicated the AP1 and NFkappaB elements located in the promoter region from -1700 to -2400 bp were critical for both basal and TGFB1-induced MMP9 transcription. These data provide the first spatial-temporal expression analysis of MMP system in the chicken ovary.

Polymorphism, genetic effect and association with egg production traits of chicken matrix metalloproteinases 9 promoter.

Matrix metalloproteinases (MMP) are key enzymes involved in cell and tissue remodeling during ovarian follicle development and ovulation. The control of MMP9 transcription in ovarian follicles occurs through a core promoter region (-2,400 to -1,700 bp). The aim of this study was to screen genetic variations in the core promoter region and examine MMP9 transcription regulation and reproduction performance. A single cytosine deletion/insertion polymorphism was found at -1954 C(+)/C(-). Genetic association analysis indicated significant correlation between the deletion genotype (C(-)) with total egg numbers at 28 weeks (p = 0.031). Furthermore, luciferase-reporter assay showed the deletion genotype (C(-)) had significantly lower promoter activity than the insertion genotype (C(+)) in primary granulosa cells (p<0.01). Therefore, the identified polymorphism could be used for marker-assisted selection to improve chicken laying performance.