BioMed Research International

BioMed Research International / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 8572837 | 11 pages | https://doi.org/10.1155/2019/8572837

Analysis of Expression and Single Nucleotide Polymorphisms of INHA Gene Associated with Reproductive Traits in Chickens

Academic Editor: Kui Li
Received26 Feb 2019
Revised02 Jul 2019
Accepted15 Jul 2019
Published08 Aug 2019

Abstract

Inhibin α (INHA) is a candidate gene controlling ovulation in poultry. As the functional center of inhibin, INHA is a molecular marker associated with egg-laying performance. The objective of the current study was to analyze the expression differences of INHA in reproductive system and single nucleotide polymorphisms (SNPs) associations with reproductive traits in chickens. A total of 260 LuHua chickens (barred-feather chicken) were adopted. Twelve SNPs were detected in INHA gene. Among the exonic SNPs, three (g. 22177991A>G, g. 22178249G>C, and g. 22178414G>A) were missense mutations, resulting in the amino acid substitutions Val→Ala, Ala→Gly, and Ala→Gly, respectively. Four SNPs in the 3 untranslated region of INHA were predicted to either disturb or create microRNA-target interactions. Five SNPs (g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A) were significantly associated with the number of eggs at 300 d of age (EN) (P < 0.05). Birds carrying GA genotype exhibited more EN than those with AA genotype (P < 0.01). In addition, quantitative real-time PCR revealed that INHA is mainly expressed in follicles on d 300 in chickens. Firstly, INHA expression increased and then decreased. The highest INHA mRNA abundance was found in the fifth largest preovulatory follicle (F5) (P < 0.01). In the prehierarchical follicles, INHA mRNA expression increased dramatically in small yellow follicles (SYF) (P < 0.01). Western blotting analysis showed that the INHA protein expression profile in the follicle was similar to its mRNA counterpart with greater expression in F5 and SYF follicles and lowest expression in F1 follicles (P < 0.05). These results suggest that INHA is a potential candidate gene improving reproductive traits in chickens.

1. Introduction

Inhibin, a member of the transforming growth factor-β (TGF-β) superfamily, plays an important role in modulating the reproductive axis and affects all reproductive events [1]. It is a gonadal glycoprotein hormone and is principally produced by the granulosa cells of ovarian follicles in females and by the sertoli cells of the testes in males [2]. Inhibin comprises multiple disulfide-linked dimers that shares a common α-subunit (encoded by INHA gene) and differs in β-subunit; when the α-subunit binds to βA-subunit, inhibin A (αβA) is formed, and when it binds to the βB-subunit, inhibin B (αβB) is formed. Both inhibin A and inhibin B have the ability to specifically suppress follicle-stimulating hormone (FSH) secretion in pituitary cells without affecting LH secretion [37]. In addition to its endocrine function, inhibin also shown to exert a variety of autocrine/paracrine in mammalian species. Many studies indicated a local role for these factors in modulating the growth of small follicles by regulating cell proliferation as well as expression of the inhibin subunits and gonadotropin receptors [8, 9].

Inhibin A is produced and secreted by the granulosa layer of the large pre-ovulatory follicle [10]. Immunizations against the INHA have resulted in increased ovulation in sheep, pigs, chickens, mice, and cows [1115]. Therefore, the INHA is thought to be the functional center of inhibin and a potential candidate gene increasing the ovulation rate in poultry. In a murine knockout model, the INHA was found to play a tumor suppressive role in gonadal tissue after gonadectomy for the adrenal cortex, with 99% of INHA mice developing adrenocortical steroid-secreting carcinoma after gonadectomy [16, 17]. Pathways involved in this effect include the differentiation into granulosa cell-like cells with the expression of fetal or gonadal markers such as LHR, FSHR, and Cyp17a1 [18]. Recently, Li et al. found that INHA gene polymorphisms are significantly associated with the presence of follicular cysts in Large White sows [19]. In poultry, Chen et al. reported that the downregulation of INHA expression has an antagonistic effect on granulosa cell apoptosis [20]. Although many studies have shown that inhibin was a critical regulator of gonadal function, whether its expression level in the developing follicles or polymorphisms are associated with reproductive traits in chicken is still unknown.

In characterizing candidate genes, it is important to detect single nucleotide polymorphisms (SNPs) and analyze the associations between these SNPs and reproductive traits, together with assessing the role of gene expression. Therefore, the aim of this study was to characterize the role of the chicken INHA gene in determining chicken reproductive traits through SNP association and expression analyses.

2. Materials and Methods

2.1. Chicken Populations

The LuHua chicken (LH) is a local breed in Shandong province, P. R. China. It was included in the Key Poultry Breeds Requiring Protection in 2002. A total of 260 birds after hatching were raised in cages with the same housing conditions and diets. Their performances were measured in Poultry Breeding Farm of Sichuan Agricultural University (Ya’an, China). Blood was collected, and genomic DNA was isolated using a standard phenol/chloroform method. The DNA purity and concentration were measured with a Nucleic Acid Protein Analyzer (Thermo Scientific, USA). Based on the original concentration, samples were diluted to the final concentration of 100 ng/μL with Tris-EDTA buffer. All DNA samples were stored at -20°C until use.

Six reproductive traits were measured: body weight at first egg (BWAFE), first egg weight (FEW), age at first egg (AFE), total egg number at 300 days of age (EN), body weight at 300 days of age (BWTA), and egg weight at 300 days of age (EWTA). Traits determination was conducted in accordance with the Committee on Experimental Animal Management of Sichuan Agricultural University and carried out strictly according to the Regulations for the Administration of Affairs Concerning Experimental Animals of the State Council of China.

2.2. Sequencing of the INHA Gene

Primers for amplifying and sequencing the chicken INHA gene (Table 1) were designed with Primer Premier 5.0 software based on the complete DNA sequences of Gallus gallus INHA genes (ENSGALG00000054770). A DNA pool containing 100 ng DNA from each of the 60 least closely related chickens was constructed. PCR was carried out using a Gene Amp PCR System 9700 (Bio-Rad, Hercules, CA, USA) thermal cycler in a final volume of 25 μL containing 8.6 μL distilled H2O, 15 μL 2×Taq PCR Master Mix (including Mg2+, dNTPs, and Taq DNA polymerase; Beijing Tian Wei Biology Technique Corporation, Beijing, China), 0.3 μL each primer (10 nmol/L), and 0.8 μL chicken genomic DNA template. PCR was performed with the following cycling conditions: denaturation at 94°C for 5 min; 35 cycles of 94°C for 30 s, 56-60°Cfor 35 s (Table 1), and 72°C for 40 s, followed by a final extension at 72°C for 5 min. PCR products were purified and sequenced via an ABI 377 DNA sequencer (Shanghai Sangon Biological Engineering Technology, Shanghai, China). All sequences were edited, assembled, and aligned with DNASTAR and Codon Code Aligner software (http://www.codoncode.com/aligner). SNPs were identified by the presence of multiple peaks at the same base by direct sequencing.


NameSequences (5'-3')Binding regionsProduct size (bp)Tm (°C)

P1F-TATTCACGGCGAGCAGACAExon 156258.7
R-CCTCAGCCCTCCCCATCT
P2F-ATCCCACAGCCCCAAGACCExon 278361.0
R-GGCAGTAGTGGAAAACGAAGC
P3F-CACTGGACCGTGTTTGACTTCGExon 2 Partial58951.5
R-GGGATGGGCTCATCATCTGG
P4F-CGGGAACTGTGCCGAAGG5'UTR74958.7
R-GCACCCCGCTGCATAACC
P5F-TGTCCCAAACTCTGTCCAATG3'UTR70756.5
R-CTCAAATGCTCCAGCACCC

UTR: untranslated region.
2.3. Genotypes Identification

Based on the DNA pools sequence results, single nucleotide polymorphisms were found in the amplification products of five primer pairs P1, P2, P3, P4, and P5. Further polymorphism genotyping was carried out one by one using DNA samples of 260 chickens. PCR and amplified products were carried out as describe as above.

2.4. Statistical Analyses

SAS 9.4 (Statistical Analysis Systems Institute Inc., Cary, NC, USA) was used to determine the relationship between genotypes and reproductive traits in chicken. The model used was: = μ + + , where is a reproductive trait measured in the chickens, μ is the population mean of the trait, is a fixed effect genotype, and is the residual error. Significance was determined using Duncan’s multiple range tests. Pearson’s chi-square test was used to assess the Hardy-Weinberg equilibrium of the four SNPs discovered in the present study. The linkage disequilibrium (LD), D and r2 values of the SNPs were estimated by Haploview [21]. The polymorphism information content (PIC) was determined following a previous method [22]. In addition, haplotypes were constructed using the PHASE program v. 2.0. Haplotypes were analyzed using the model for the single marker association test by considering birds with 0, 1, or 2 copies of the haplotype in question. The PROC REG procedure of SAS (version 9.4, SAS Institute Inc.) was used to perform the analysis. Values were significant at P < 0.05.

2.5. Expression Abundances of INHA mRNA

To determine the reproductive system tissue-specific expression patterns of INHA mRNA, total RNA was extracted from the tissues of five chickens at 25 weeks, including the fallopian tubes (infundibulum, isthmus, and uterus), hierarchical follicles (F1, F2, F3, F4, and F5), and prehierarchical follicles [small white follicles (SWF), small yellow follicles (SYF), and large white follicles (LWF)], using TRIzol reagent (TakaRa Biotech Co. Ltd., Dalian, China) and was dissolved in RNase-free H2O (Tiangen Biotech Co., Ltd, Beijing, China). The integrity of the RNA was evaluated via electrophoresis on 1 % agarose gels, and the concentration and purity were analyzed with a NanoDrop 2000 by determining the absorbance ratio of 260/280 nm (Thermo Scientific, Waltham, MA, USA). The cDNA was synthesized by reverse-transcription PCR using 1 μg total RNA, 2 μL 5 × RT buffer, 0.5 μL RT enzyme mix, 0.5 μL primer mix, and 6 μL nuclease-free water (Toyobo Life Science Department, Shanghai, China). The reverse transcription reaction was maintained at 37°C for 15 min, followed by incubation at 98°C for 5 min. The cDNA samples were stored at -20°C. Gene-specific primers (Fw: 5-ACTACTGCCACGGGAACTGT-3, Rv: 5-GGAGTAGCCACCATCAGAGG-3) for qRT-PCR were designed using Primer 5 software according to the coding sequence of the chicken INHA gene (GenBank accession no. NM_001031257). qRT-PCR was conducted via a CFX96 Real-time System (Bio-Rad, Hercules, CA, USA) with the following conditions: 98°C for 2 min; 39 cycles of 98°C for 5 s and 55.7°C for 10 s; followed by a final extension at 95°C for 10 s. Each PCR reaction contained 5 μL Ssofast Evagreen supermix (Bio- Rad, Hercules, CA, USA), 1 μL of cDNA (50 ng/μL), 0.8 μL each primer (10 μM), and 2.4 μL ddH2O to a final reaction volume of 10 μL. All samples were triplicated. The housekeeping gene β-actin (GenBank accession no. AF047874; Fw: 5-GAGAAATTGTGCGTGACATCA-3, Rv: 5-CCTGAACCTCTCATTGCCA-3) was used for normalization of target gene expression [23]. The relative gene expression levels of INHA were calculated using the comparative method [24], where ΔCt = Ct target gene - Ct housekeeping gene. Differences in INHA mRNA expression were examined by one-way analysis of variance (ANOVA). Data are presented as the mean ± standard error means (SEM), and significances were determined at P < 0.05.

2.6. Expression Abundances of INHA Protein

INHA protein expression levels in different follicles were detected by Western Blotting. Total protein was extracted from the tissues of five chickens at 25 weeks, including the hierarchical follicles (F1, F2, F3, F4, and F5) and prehierarchical follicles (SWF, SYF, and LWF) using a Protein Extraction Kit (BestBio Biotech Co. Ltd., Shanghai, China). The concentration and purity of the protein samples were determined using a BCA Protein Assay Kit (BestBio Biotech Co. Ltd., Shanghai, China); standard curves were drawn to calculate the protein concentration. Samples were triplicated. The β-actin protein was used as a reference. A total of 25 μg protein was resolved by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene fluoride (PVDF) membranes. After blocking with 5 % non-fat milk in 1×Tris-buffered saline with Tween (TBST) buffer for 2 h at room temperature, membranes were incubated with rabbit anti-chicken INHA monoclonal antibody (Abcam, Cambridge, UK; 1:1000) and rabbit anti-chicken β-actin monoclonal antibody (Abcam, Cambridge, UK; 1:1000) overnight at 4°C. Blots were then washed in 1 × TBST buffer and probed with goat anti-rabbit horseradish peroxidase (HRP)-conjugated IgG secondary antibody (diluted 1:2000 in 1 × TBST; Abcam, Cambridge, UK) for 1 h at room temperature. Binding was visualized with enhanced chemiluminescence (ECL) reagent (Beyotime Institute of Biotechnology, Jiangsu, China) using a ChemiDoc XRS instrument (Bio-Rad, Hercules, CA, USA). Quantity One Software (Bio-Rad, Hercules, CA, USA) was used for densitometric analysis [25].

3. Results

3.1. Allele and Genotype Frequencies of Chicken INHA Gene SNPs

All exons and part of the untranslated region (UTR) of the chicken INHA gene were amplified and sequenced. A total of 12 SNPs were detected (Figure 1). Among them, one was located in the promoter region, three were in exon 1, three were in exon 2, and five were in the 3 UTR (Table 2). Among the exonic SNPs, g. 22177991 A>G, g. 22178249 G>C, and g. 22178414 G>A were missense mutations, which resulted in the amino acid substitutions Val→Ala, Ala→Gly, and Ala→Gly, respectively. SNPs in the 3 UTR were predicted to either disturb or create microRNA- (miRNA-) target interactions. Ten miRNA interactions weakened and four miRNA interactions strengthened (Table 3).


Primer pairs no.VariationsChr. positionGene regionFunction

P1g. 22176805 C>T22176805Promoter region-
P1g. 22176870T>C22176870Exon 1Synonymous
P2g. 22177100T>C22177100Exon 1Synonymous
P2g. 22177149T>C22177149Exon 1Synonymous
P3g. 22177991A>G22177991Exon 2Missense
P3g. 22178249G>C22178249Exon 2Missense
P4g. 22178414G>A22178414Exon 2Missense
P4g. 22178618T>C221786183' UTR-
P5g. 22178681A>T221786813' UTR-
P5g. 22178728C>T221787283' UTR-
P5g. 22178835A>G221788353' UTR-
P5g. 22178975G>A221789753' UTR-

UTR: Untranslated region. N =260.

VariationsDisturbed miRNAsCreated miRNAsBinding area

g. 22178618T>CmiR-128-3p-22178613-22178619
miR-6630-3p22178616-22178621

g. 22178728C>T-miR-125b-5p22178727-22178733
miR-34b-5p22178723-22178728
miR-6673-3p22178727-22178732

g. 22178835A>GmiR-17-3p-22178831-22178837
miR-6647-5p22178834-22178839
miR-6669-3p22178830-22178835
miR-7456-5p22178830-22178837
miR-7482-3p22178831-22178837

g. 22178975G>AmiR-181b-1-3p-22178972-22178978
miR-6613-3p22178972-22178978
miR-6677-5p22178973-22178978
-miR-460b-5p22178974-22178980

The genotype and allele frequencies and related genetic information for the 12 SNPs in the chicken INHA gene are summarized in Table 4. For all SNPs, their genotype and allele frequencies were above 5 %, which indicated that it was appropriate to conduct INHA gene analysis. For SNP1, the CT genotype frequency (0.6121) was higher than those of CC (0.3103) and TT (0.0776), and the allele frequency of C (0.61635) was higher than that of T (0.38365). The heterozygous genotype frequency was higher than the homozygous genotype frequencies for all 12 SNPs. The original allele frequency was higher than the mutant allele frequency for SNP1, SNP2, SNP3, SNP4, SNP7, SNP8, and SNP9, respectively, while the original allele frequency was lower than the mutant allele frequency for SNP5, SNP6, SNP10, SNP11, and SNP12. The PIC test results indicated that all SNPs could be considered intermediate polymorphisms, making them good genetic markers. All SNPs except SNP1 were in Hardy-Weinberg equilibrium (P > 0.05).


SNPsGenotype frequency (%)Allele frequency (%)PICP-value

SNP1CCCTTTCT
C22176805 T31.0361.217.7661.63538.3650.3611P=0.0068
SNP2TTTCCCTC
T22176870C35.3555.179.4862.93537.0650.3577P=0.0821
SNP3TTTCCCTC
T22177100C37.0754.318.6264.22535.7750.3539P=0.0835
SNP4TTTCCCTC
T22177149C34.4956.039.4862.50537.4950.3589P=0.0835
SNP5AAAGGGAG
A22177991G10.3452.5937.0736.63563.3650.3565P=0.235
SNP6GGGCCCGC
G22178249C13.8055.1731.0341.38558.6150.3675P=0.2132
SNP7GGGAAAGA
G22178414A30.1756.0313.8058.18541.8150.3682P=0.1611
SNP8TTTCCCTC
T22178618C31.0355.1713.8058.61541.3850.3675P=0.2132
SNP9AAATTTAT
A22178681T31.0355.1713.8058.61541.3850.3675P=0.2132
SNP10CCCTTTCT
C22178728T14.6655.1730.1742.24557.7550.3689P=0.2403
SNP11AAAGGGAG
A22178835G12.9456.0331.0340.95559.0450.3667P=0.2132
SNP12GGGAAAGA
G22178975A19.8357.7622.4148.7151.290.3748P=0.1467

P-value is the result of test of Hardy-Weinberg equilibrium; PIC < 0.25 indicated low polymorphism, 0.25 < PIC < 0.50 indicated intermediate polymorphism, and PIC > 0.50 indicated high polymorphism.
3.2. Relationships between Genotypes and Reproductive Traits

The results of the association analyses showed that g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A were significantly associated with EN (P < 0.05, Table 5). No significant associations were found between other SNPs and the reproductive traits. Therefore, subsequent genotype association analyses were performed for g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A (Table 6). The results showed that genotypes at the g. 22176870T>C, g. 22177100T>C, g. 22177149T>C SNPs were significantly associated with EN (P < 0.05). Birds with the TC genotype at these SNPs had significantly higher ENs than those with the CC genotype. At the INHA exon 2 SNP g. 22177991A>G was significantly associated with EN (P = 0.025), birds with the AG genotype had higher EN values. SNP g. 22178975G>A genotypes in INHA gene was significantly associated with EN (P < 0.05), with genotype GA had more EN than AA genotype (P < 0.01). Conversely, chickens with genotype AA had heavier FEW than those with GA genotype.


PolymorphismTraits (P value of significant test)
AFE(days) BWFE(g)FEW(g)BWTA(g)EWTA(g)EN(count)

g. 22176805 C>T0.8380.4070.9620.3080.1920.353
g. 22176870T>C0.2140.4060.9980.4570.8760.017
g. 22177100T>C0.2960.2590.9540.2220.8460.022
g. 22177149T>C0.5160.2030.9670.2110.8200.033
g. 22177991A>G0.2790.1790.9900.3210.7680.025
g. 22178249G>C0.5690.2400.8120.5040.3840.150
g. 22178414G>A0.7070.3930.7570.6180.3800.259
g. 22178618T>C0.5690.1850.8300.3440.5730.124
g. 22178681A>T0.5690.1850.8300.3440.5730.124
g. 22178728C>T0.6280.2420.8590.4520.4160.117
g. 22178835A>G0.5690.1850.8300.3440.5730.124
g. 22178975G>A0.2820.1260.1060.1330.1170.018

P ≤ 0.05; AFE=age at first egg; BWFE=body weight at first egg; FEW= first egg weight; BWTA=body weight at 300 days of age; EWTA=egg weight at 300 days of age; EN=total number of eggs at 300 days of age.

SNPGenotypesP-value

g. 22176870T>CTTTCCC
111.512±1.928b117.516±1.544a109.000±3.723b0.017
g. 22177100T>CTTTCCC
111.927±1.939b117.250±1.552a109.000±3.744b0.022
g. 22177149T>CTTTCCC
112.070±1.887b117.365±1.559a107.900±3.913b0.033
g. 22177991A>GAAAGGG
111.083±3.576b117.607±1.586a111.279±1.889b0.025
g. 22178975G>AGGGAAA
108.769±2.422b114.130±2.575b117.000±1.509a0.018

Results are expressed as mean ± standard errors. Values in the table were the total number of eggs at 300 days of age for different genotypes. Different letters indicate significant differences (P < 0.05).
3.3. Construction of Haplotypes and Their Associations with Chicken Reproductive Traits

Analysis of LD between SNPs in the INHA gene was shown in Figure 2. Haplotypes were constructed based on the 12 SNPs identified in LH chickens using the Haploview program (Table 7). H1 was the most frequent haplotype at 53.7 %. Ten diplotypes were obtained based on these five haplotypes (Table 8). However, no significant associations were determined between the reproductive traits in the chicken population and the haplotypes according to regression coefficient analysis.


HaplotypeSNP1SNP2SNP3SNP4SNP5SNP6SNP7SNP8SNP9SNP10SNP11Frequency

H1CTTTGCGTATG0.54
H2TCCCAGACTCA0.33
H3CTTTGGACTCA0.05
H4TTTTGCGTATG0.02
H5CCCCAGACTCA0.02


DiplotypesH1H1H1H2H1H3H1H4H2H2H2H3H2H4H2H5H3H4H4H4

Frequency (%)24.1447.414.30.865.173.450.863.450.860.86

3.4. Expression of INHA in the Chicken Reproductive System

According to expression analysis, INHA is mainly expressed in the ovary and follicles (F1-F5, SYF, and LWF). During development of chicken follicles, the relative abundance of INHA mRNA initially increased and then decreased, with significantly higher expression in F5 than others (P < 0.05). In pre-grade follicles, the INHA mRNA expression increased abruptly in SWF follicle (P < 0.05). Expression levels in the F1, F2, and F3 follicles were low and did not significantly differ from each other (P > 0.05) (Figure 3(a)).

The protein expression levels of INHA in different chicken period follicles were detected by Western blotting, with β-actin protein used as a reference. The results showed higher expression of the INHA protein in F5 compared with other follicles (P < 0.01) (Figure 3(b)). In the prehierarchical follicles, INHA protein was highly expressed in SYF. These results suggested that INHA is a potential candidate gene for improving chicken reproductive traits.

4. Discussion

Inhibins play key roles in folliculogenesis, oocyte maturation, and embryo development [26]. The INHA gene encodes the functional center of inhibin and may exhibit potential for increasing the ovulation rate in poultry. For example, the downregulation of INHA gene expression in cultured goose granulosa cells resulted in significant increases in apoptosis and proliferation indexes, a reduced percentage of cells in the G1 phase, and a correspondingly elevated percentage of cells in the S phase [20]. Moreover, transgenic mice overexpressing the rat INHA gene exhibited a reduced litter size and longer intervals between pregnancies when compared with control mice [27]. In some studies, polymorphisms in the INHA gene have been found to be significantly associated with follicular cysts in humans [28], pigs [19, 29], and other mammals [3032]. However, INHA polymorphisms associated with reproductive traits in chicken were previously unknown. The purpose of this experiment was to improve the reproductive performance of LH chickens and provide molecular markers for the selective breeding of laying hens.

Ovulation rate is an important reproductive trait, and the important function of INHA makes it a strong candidate for improving the poultry ovulation rate. We hypothesized that INHA may be a major gene affecting chicken egg production. The AFE, as well as other egg production and egg weight parameters, are important traits used in the breeding of high-quality layers. Therefore, understanding the effects of the INHA gene on these reproductive traits is essential. However, information previously available in the literature regarding the relationships between INHA polymorphisms and reproductive traits was inconclusive.

Recently, several studies in poultry have been identified crucial genes and explored their relationship with phenotypes, such as Liu with 279 Dongxiang blue-shelled chickens [33], Yu with 188 female Muchuan black-boned chickens [34], and EI-Sabrout with 200 Lohmann Brown hens [35]. In this study, 12 SNPs were detected in INHA gene. All SNPs except for SNP1 were in Hardy-Weinberg equilibrium (P > 0.05), and their genotype and allele frequencies were above 5 %, which indicate that the population is not affected by artificial selection and its sample size is appropriate for association analysis [36].

Five SNPs (g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A) that were significantly associated with chicken EN. This result is similar to that studied in Boer goats, in which the c.651G>A mutation in the INHA gene affected the lambing rate [31]. Egg production is an important economic indicator in the poultry industry, and in poultry breeding programs, EN is considered the most valuable indicator of total egg production potential. Our findings therefore confirm the influence of INHA polymorphisms on the female reproductive traits of LH chickens. We therefore deduce that INHA is an important gene for improving the EN in chickens.

The SNP g. 22177991A>G results in the amino acid substitutions Val→Ala. Zi et al. found that sequence variation in INHA was associated with the prolificacy of goat breeds, suggesting that amino acid substitutions in this gene may affect reproductive traits in chicken [37]. While SNP3 results in a synonymous mutation that causes no amino acid change, this mutation was still associated with a reproductive trait. One explanation for this is that this mutation may affect INHA function by altering the stability of the mRNA and its translation efficiency. Therefore, this SNP may be associated with the mRNA expression level and increase in INHA concentration may lead to decreases in FSH concentrations [38].

Zhu et al. analyzed miRNA-related SNPs and found that two SNPs in the 3 UTRs of target genes were predicted to either disturb or create miRNA-target interactions [39]. In the present study, SNP g. 22178618T>C in the 3 UTR of INHA was predicted to disturb miR-128-3p combination. Yu et al. confirmed the suppression of CYP2C9 by the miRNA hsa-miR-128-3p in human liver cells and its association with hepatocellular carcinoma [40], while miR-128-3p was found to suppress hepatocellular carcinoma proliferation [41]. SNP g. 22178728C>T was predicted to create bindings sites for two miRNAs miR-125b-5p and miR-34b-5p. miR-125b-5p serves as a novel biomarker for HBV-positive hepatocellular carcinoma [42], and miR-34b-5p inhibits the expression of Bcl-2 in ovarian cancer cells [43]. This suggests that miRNA-related SNPs in the 3 UTR of INHA may affect poultry reproductive performance. Further studies are needed to confirm that miRNA-related SNPs regulate mRNA and protein expression levels.

Good laying performance of poultry depends on the growth and development of follicles in the ovary. In poultry, only a few of the large number of follicles mature, with only about 5 % developing into SYF. Improving the laying rate and reproductive performance of native chicken breeds has become an urgent concern for chicken breeding and production. The expression levels and forms of inhibin are closely related to physiological activities such as follicular development, recruitment, and dominance selection. Furthermore, as INHA is the functional center of inhibin, we hypothesis that the expression level of the chicken INHA gene plays a dominant role in follicle development and is related to chicken reproductive traits. Although many studies have shown INHA to be a critical regulator of gonadal function, little is known about its expression in chicken developing follicles and the association between its expression and reproductive traits. In this study, we found that the INHA gene only expressed in the developing follicles. Intriguingly, INHA mRNA expression was the lowest in F1 and highest in F5, consistent with previous results in humans [44]. In the prehierarchical follicles, INHA mRNA expression increased sharply in the SWF. Zi et al. previously found that the mRNA expression levels of INHA affected prolificacy in goats [37]. The development of LWF into SYF is an important process in the recruitment of follicles in poultry [45]. Additionally, the INHA protein expression profile in the follicle was similar to that of its mRNA counterpart with greater expression in F5 and SYF follicles. These results indicate that the expression of INHA is related to follicle development. It can therefore hypothesize that INHA is involved in the regulation of follicle development, playing a critical role in follicle recruitment. Further study is needed to illuminate the specific action mechanism of INHA.

5. Conclusion

Generally, 12 SNPs were identified in chicken INHA gene. Five of them were significant associated with egg numbers. Among the exonic SNPs, g. 22177991A>G, g. 22178249G>C, and g. 22178414G>A were missense mutations, which resulted in the amino acid substitutions Val→Ala, Ala→Gly, and Ala→Gly. In addition, INHA highly expressed in F5 and SYF follicles. Therefore, we conclude that INHA is a candidate gene affecting egg production, and it plays a critical role in the recruitment of follicles in chickens. INHA SNPs are possible molecular markers for the genetic selection of layers.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declared no conflicts of interest.

Authors’ Contributions

Zhifu Cui and Lingbin Liu contributed equally to this work.

Acknowledgments

The authors gratefully acknowledge Lin Ye for sample collection. This work was financially supported by scientific research fund of Sichuan Provincial education department (15ZA0025), the Thirteenth Five Year Plan for Breeding Program in Sichuan (2016NYZ0050), and the China Agriculture Research System (CARS-41).

References

  1. C. Welt, Y. Sidis, H. Keutmann, and A. Schneyer, “Activins, inhibins, and follistatins: from endocrinology to signaling. A paradigm for the new millennium,” Experimental Biology and Medicine, vol. 227, no. 9, pp. 724–752, 2002. View at: Publisher Site | Google Scholar
  2. S. Ying, “Inhibins and activins: chemical properties and biological activity,” Experimental Biology and Medicine, vol. 186, no. 3, pp. 253–264, 1987. View at: Publisher Site | Google Scholar
  3. T. K. Woodruff and J. P. Mather, “Inhibin, activin and the female reproductive axis,” Annual Review of Physiology, vol. 57, no. 1, pp. 219–244, 1995. View at: Publisher Site | Google Scholar
  4. P. Knight, “Roles of inhibins, activins, and follistatin in the female reproductive system,” Frontiers in Neuroendocrinology, vol. 17, no. 4, pp. 476–509, 1996. View at: Publisher Site | Google Scholar
  5. A. L. L. Durlinger, P. Kramer, B. Karels et al., “Control of primordial follicle recruitment by anti-mullerian hormone in the mouse ovary,” Endocrinology, vol. 140, no. 12, pp. 5789–5796, 1999. View at: Publisher Site | Google Scholar
  6. S. Yang and S. Yu, “The advance on apoptosis regulated by hormones in ovarian follicles,” Dongwu Yixue Jinzhan, vol. 23, no. 2, pp. 1–5, 2001. View at: Google Scholar
  7. X. Chen and X. Jiang, “The apoptosis and mechanism of ovary cells,” Heilongjiang Journal of Animal Reproduction, vol. 12, pp. 13–18, 2004. View at: Google Scholar
  8. J. Findlay, A. Drummond, K. Britt et al., “The roles of activins, inhibins and estrogen in early committed follicles,” Molecular and Cellular Endocrinology, vol. 163, no. 1-2, pp. 81–87, 2000. View at: Publisher Site | Google Scholar
  9. P. G. Knight and C. Glister, “Potential local regulatory functions of inhibins, activins and follistatin in the ovary,” Reproduction, vol. 121, no. 4, pp. 503–512, 2001. View at: Publisher Site | Google Scholar
  10. P. Johnson, C. Brooks, and A. Davis, “Pattern of secretion of immunoreactive inhibin/activin subunits by avian granulosa cells,” General and Comparative Endocrinology, vol. 141, no. 3, pp. 233–239, 2005. View at: Publisher Site | Google Scholar
  11. D. Tannetta, S. Feist, E. Bleach, N. Groome, L. Evans, and P. Knight, “Effects of active immunization of sheep against an amino terminal peptide of the inhibin alpha C subunit on intrafollicular levels of activin A, inhibin A and follistatin,” Journal of Endocrinology, vol. 157, no. 1, pp. 157–168, 1998. View at: Publisher Site | Google Scholar
  12. F. Shi, K. Mochida, O. Suzuki et al., “Development of embryos in superovulated guinea pigs following active immunization against the inhibin α-subunit,” Endocrine Journal, vol. 47, no. 4, pp. 451–459, 2000. View at: Publisher Site | Google Scholar
  13. T. M. Lovell, P. G. Knight, N. P. Groome, and R. T. Gladwell, “Changes in plasma inhibin A levels during sexual maturation in the female chicken and the effects of active immunization against inhibin α-subunit on reproductive hormone profiles and ovarian function,” Biology of Reproduction, vol. 64, no. 1, pp. 188–196, 2001. View at: Publisher Site | Google Scholar
  14. X. Jiang, L. Yang, G. Liu et al., “Effects of inhibin α (1–32) gene immunization on mouse reproductive traits,” Chinese Journal of Veterinary Science, vol. 22, pp. 368–370, 2002. View at: Google Scholar
  15. M. S. Medan, T. Takedom, Y. Aoyagi et al., “The effect of active immunization against inhibin on gonadotropin secretions and follicular dynamics during the estrous cycle in cows,” The Journal of Reproduction and Development, vol. 52, no. 1, pp. 107–113, 2006. View at: Publisher Site | Google Scholar
  16. M. M. Matzuk, M. J. Finegold, J.-G. J. Su, A. J. W. Hsueh, and A. Bradley, “Alpha-inhibin is a tumour-suppressor gene with gonadal specificity in mice,” Nature, vol. 360, no. 6402, pp. 313–319, 1992. View at: Publisher Site | Google Scholar
  17. M. M. Matzuk, M. J. Finegold, J. P. Mather, L. Krummen, H. Lu, and A. Bradley, “Development of cancer cachexia-like syndrome and adrenal tumors in inhibin-deficient mice,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 91, no. 19, pp. 8817–8821, 1994. View at: Publisher Site | Google Scholar
  18. B. D. Looyenga and G. D. Hammer, “Origin and identity of adrenocortical tumors in inhibin knockout mice: implications for cellular plasticity in the adrenal cortex,” Molecular Endocrinology, vol. 20, no. 11, pp. 2848–2863, 2006. View at: Publisher Site | Google Scholar
  19. W. Li, L. Sun, S. Chen et al., “Association of inhibin-α gene polymorphisms with follicular cysts in large white sows,” Theriogenology, vol. 84, no. 6, pp. 839–845, 2015. View at: Publisher Site | Google Scholar
  20. F. Chen, X. Jiang, X. Chen, G. Liu, and J. Ding, “Effects of downregulation of inhibin α gene expression on apoptosis and proliferation of goose granulosa cells,” Journal of Genetics and Genomics, vol. 34, no. 12, pp. 1106–1113, 2007. View at: Publisher Site | Google Scholar
  21. J. C. Barrett, B. Fry, J. Maller, and M. J. Daly, “Haploview: analysis and visualization of LD and haplotype maps,” Bioinformatics, vol. 21, no. 2, pp. 263–265, 2005. View at: Publisher Site | Google Scholar
  22. D. Botstein, R. L. White, M. Skolnick, and R. W. Davis, “Construction of a genetic linkage map in man using restriction fragment length polymorphisms,” American Journal of Human Genetics, vol. 32, no. 3, pp. 314–331, 1980. View at: Google Scholar
  23. H.-Y. Xu, Y. Wang, Y.-P. Liu, J.-W. Wang, and Q. Zhu, “Polymorphisms and expression of the chicken POU1F1 gene associated with carcass traits,” Molecular Biology Reports, vol. 39, no. 8, pp. 8363–8371, 2012. View at: Publisher Site | Google Scholar
  24. K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method,” Methods, vol. 25, no. 4, pp. 402–408, 2001. View at: Publisher Site | Google Scholar
  25. L. Liu, D. Li, E. R. Gilbert et al., “Effect of monochromatic light on expression of estrogen receptor (ER) and progesterone receptor (PR) in ovarian follicles of Chicken,” PLoS ONE, vol. 10, no. 12, Article ID 0144102, 2015. View at: Publisher Site | Google Scholar
  26. A. V. Sirotkin, “Cytokines: signalling molecules controlling ovarian functions,” The International Journal of Biochemistry & Cell Biology, vol. 43, no. 6, pp. 857–861, 2011. View at: Publisher Site | Google Scholar
  27. B.-N. Cho, M. L. McMullen, L. Pei, C. J. Yates, and K. E. Mayo, “Reproductive deficiencies in transgenic mice expressing the rat inhibin α-subunit gene,” Endocrinology, vol. 142, no. 11, pp. 4994–5004, 2001. View at: Publisher Site | Google Scholar
  28. H. Rah, Y. J. Jeon, J. J. Ko et al., “Association of inhibin α gene promoter polymorphisms with risk of idiopathic primary ovarian insufficiency in Korean women,” Maturitas, vol. 77, no. 2, pp. 163–167, 2014. View at: Publisher Site | Google Scholar
  29. W. Li, S. Chen, H. Li et al., “A new insertion/deletion fragment polymorphism of inhibin-α gene associated with follicular cysts in Large White sows,” Journal of Veterinary Medical Science, vol. 78, no. 3, pp. 473–476, 2016. View at: Publisher Site | Google Scholar
  30. S. E. Harris, A. L. Chand, I. M. Winship et al., “INHA promoter polymorphisms are associated with premature ovarian failure,” MHR: Basic science of reproductive medicine, vol. 11, no. 11, pp. 779–784, 2005. View at: Publisher Site | Google Scholar
  31. W. Wu, G. Hua, L. Yang et al., “Association analysis of the INHA gene with litter size in Boer goats,” Small Ruminant Research, vol. 82, no. 2-3, pp. 139–143, 2009. View at: Publisher Site | Google Scholar
  32. K. Tang, S. Li, W. Yang et al., “An MspI polymorphism in the inhibin alpha gene and its associations with superovulation traits in Chinese Holstein cows,” Molecular Biology Reports, vol. 38, no. 1, pp. 17–21, 2011. View at: Publisher Site | Google Scholar
  33. L. Liu, Z. Cui, Q. Xiao et al., “Polymorphisms in the chicken growth differentiation factor 9 gene associated with reproductive traits,” BioMed Research International, vol. 2018, Article ID 9345473, 11 pages, 2018. View at: Publisher Site | Google Scholar
  34. S. Yu, J. Liao, M. Tang et al., “A functional single nucleotide polymorphism in the tyrosinase gene promoter affects skin color and transcription activity in the black-boned chicken,” Poultry Science, vol. 96, no. 11, pp. 4061–4067, 2017. View at: Publisher Site | Google Scholar
  35. K. El-Sabrout and S. Aggag, “Association of single nucleotide polymorphism in melanocortin receptor gene with egg production traits in Lohmann Brown chickens,” Kafkas Universitesi Veteriner Fakultesi Dergisi, vol. 24, no. 1, pp. 145–148, 2018. View at: Google Scholar
  36. G. Salanti and J. P. A. Ioannidis, “Hardy-Weinberg equilibrium,” Encyclopedia of Ecology, vol. 17, no. 1, pp. 1844–1846, 2008. View at: Google Scholar
  37. X.-D. Zi, H.-W. Xu, and Y. Wang, “Variation in sequences and mRNA expression levels of inhibin subunits α (INHA) and βA (INHBA) genes between prolific and nonprolific goat breeds,” Molecular Reproduction and Development, vol. 79, no. 4, p. 238, 2012. View at: Publisher Site | Google Scholar
  38. F. Capon, M. H. Allen, M. Ameen et al., “A synonymous SNP of the corneodesmosin gene leads to increased mRNA stability and demonstrates association with psoriasis across diverse ethnic groups,” Human Molecular Genetics, vol. 13, no. 20, pp. 2361–2368, 2004. View at: Publisher Site | Google Scholar
  39. Y. Zhu, W. Xue, J. Wang et al., “Identification of common carp (Cyprinus carpio) microRNAs and microRNA-related SNPs,” BMC Genomics, vol. 13, no. 1, p. 413, 2012. View at: Publisher Site | Google Scholar
  40. D. Yu, B. Green, A. Marrone et al., “Suppression of CYP2C9 by MicroRNA hsa-miR-128-3p in Human Liver Cells and Association with Hepatocellular Carcinoma,” Scientific Reports, vol. 5, no. 1, pp. 8534–8534, 2015. View at: Publisher Site | Google Scholar
  41. C. Huang, X. Huang, J. Zhu et al., “miR-128-3p suppresses hepatocellular carcinoma proliferation by regulating PIK3R1 and is correlated with the prognosis of HCC patients,” Oncology Reports, vol. 33, no. 6, pp. 2889–2898, 2015. View at: Publisher Site | Google Scholar
  42. B. G. Giray, G. Emekdas, S. Tezcan et al., “Profiles of serum microRNAs; miR-125b-5p and miR223-3p serve as novel biomarkers for HBV-positive hepatocellular carcinoma,” Molecular Biology Reports, vol. 41, no. 7, pp. 4513–4519, 2014. View at: Publisher Site | Google Scholar
  43. Y. L. Xie, Y. J. Yang, C. Tang et al., “Estrogen combined with progesterone decreases cell proliferation and inhibits the expression of Bcl-2 via microRNA let-7a and miR-34b in ovarian cancer cells,” Clinical and Translational Oncology, vol. 16, no. 10, pp. 898–905, 2014. View at: Publisher Site | Google Scholar
  44. C. K. Welt and A. L. Schneyer, “Differential regulation of inhibin B and inhibin A by follicle-stimulating hormone and local growth factors in human granulosa cells from small antral follicles,” The Journal of Clinical Endocrinology & Metabolism, vol. 86, no. 1, pp. 330–336, 2001. View at: Publisher Site | Google Scholar
  45. M. Safi, N. Buys, O. M. Onagbesan, B. Vleugels, and E. Decuypere, “Quantification of inhibin/activin α and β(A) subunit messenger ribonucleic acid by competitive reverse transcription-polymerase chain reaction in chicken granulosa cells during follicular development,” Biology of Reproduction, vol. 59, no. 5, pp. 1047–1054, 1998. View at: Publisher Site | Google Scholar

Copyright © 2019 Zhifu Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

731 Views | 295 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.