Variation in sequences and mRNA expression levels of inhibin subunits α (INHA) and βA (INHBA) genes between prolific and nonprolific goat breeds

CYP19A1 May Influence Lambing Traits in Goats by Regulating the Biological Function of Granulosa Cells

Simple Summary

Aromatase (CYP19A1), a member of the cytochrome family, is widely expressed in ovarian and granulosa cells and is primarily responsible for the conversion of androgens to estrogens. Increased expression of CYP19A1 in follicular granulosa cells has implications for cell proliferation, steroid hormone secretion, and the expression of related functional indicator genes. We hypothesize that CYP19A1 may indirectly influence lambing numbers in goats by regulating follicular cell growth and development, as well as ovarian ovulation.

Abstract

Abnormal expression of CYP19A1, a gene related to steroid hormone synthesis, causes steroid hormone disruption and leads to abnormal ovulation in granulosa cells. However, the exact mechanism of CYP19A1 regulation is unclear. In this study, we confirmed the localization of CYP19A1 in goat ovarian tissues using immunohistochemistry. Subsequently, we investigated the effects of CYP19A1 on granulosa cell proliferation, steroid hormone secretion, and expression of candidate genes for multiparous traits by overexpressing and silencing CYP19A1 in goat granulosa cells (GCs). The immunohistochemistry results showed that CYP19A1 was expressed in all types of follicular, luteal, and granulosa cells, with subcellular localization results revealing that CYP19A1 protein was mainly localized in the cytoplasm and nucleus. Overexpression of CYP19A1 significantly increased the mRNA levels of CYP19A1, FSHR, and INHBA, which are candidate genes for multiple birth traits in goats. It also promoted cell proliferation, PCNA and Cyclin E mRNA levels in granulosa cells, and secretion of estrogen and progesterone. However, it inhibited the mRNA levels of STAR, CYP11A1, and 3βSHD, which are genes related to steroid synthesis. Silencing CYP19A1 expression significantly reduced CYP19A1, FSHR, and INHBA mRNA levels in granulosa cells and inhibited granulosa cell proliferation and PCNA and Cyclin E mRNA levels. It also reduced estrogen and progesterone secretion but enhanced the mRNA levels of STAR, CYP11A1, and 3βSHD. CYP19A1 potentially influenced the lambing traits in goats by affecting granulosa cell proliferation, hormone secretion, and expression of candidate genes associated with traits for multiple births.

Population validation of reproductive gene mutation loci and association with the litter size in Nubian goat

Litter size is an important component trait of doe reproduction. By improving it, production efficiency and economic benefits can be significantly provided. Genetic marker-assisted selection (MAS) based on proven molecular indicators could enhance the efficacy of goat selection, as well as litter size trait. Many molecular markers have been identified that they can be used to improve litter size in different goat breeds. However, the presence and value of these markers vary among goat breeds. In the present study, we used the reported loci on other breeds of goat as candidate loci to detect whether these loci appear in this Nubian goat population; then we proceed to genotype and detect surrounding loci (50 bp) by multiplex PCR and sequencing technology. As a result, 69 mutations (59 SNPs and 10 indels) were screened out from 23 candidate genes in Nubian goat population, 12 loci were significantly associated with the litter size of first-parity individuals; 5 loci were significantly associated with the litter size of second-parity individuals; 3 loci were significantly associated with the litter size of third-parity individuals. In addition, five loci were significantly associated with the average litter size. The additive effect value of KITLG: g.18047318 G>A in first parity, KITLG: g.18152042G>A in third parity, KISS-1: g.1341674 C>G in first parity, and GHR: g.32134187G>A in second parity exceed more than 0.40, and the preponderant alleles are G, C, A and G, respectively. Further, linkage disequilibrium analysis of 21 mutation loci shows that 3 haplotype blocks are formed, and the litter size of combination type AACC in KISS-1 gene and AAGG in KITLG gene are significantly lower than that of other combinations genotype in first parity (P<0.05). These findings can provide effective candidate DNA markers for selecting superior individuals in Nubian goat breeding.

INHBA transfection regulates proliferation, apoptosis and hormone synthesis in sheep granulosa cells

Inhibin subunit beta A (INHBA) participates in the synthesis of inhibin A, activin A and activin AB. Here we investigated the effect and molecular mechanism of INHBA on proliferation, apoptosis and hormone synthesis in sheep granulosa cells (GCs) using in vitro transfection. We first noticed that INHBA expression increased with follicle diameter and was widely distributed in ovarian tissue. The proliferation rate of GCs was significantly increased and decreased with overexpression and silence of INHBA, respectively, compared with the negative controls. INHBA transfection affected GC proliferation and apoptosis, regulating the expression of many cell cycle-related and apoptosis-related genes. INHBA overexpression significantly decreased activin and estradiol secretion while increasing inhibin and progesterone secretion. The expression of follicle-stimulating hormone beta subunit was significantly decreased and increased with INHBA overexpression and knockdown, respectively. Notably, silence of INHBA inhibited the expression of many transforming growth factor beta-related genes. Overall, the functional molecule of INHBA gene may be associated with follicular development via regulating proliferation, apoptosis and folliculogenesis-related hormone secretion of sheep GCs. In addition, our findings may contribute to a better understanding of the law of follicular development and thus improve the reproductive performance of female animals.

A functional polymorphism of inhibin alpha subunit at miR-181b-1-3p-binding site regulates proliferation and apoptosis of chicken ovarian granular cells

INHA, the gene encoding the inhibin alpha subunit, was involved in folliculogenesis in mammals, but no study was reported for its working pathway in birds. Here we hypothesize that gene polymorphism in INHA 3'UTR might influence miRNAs binding efficiency and further affect the function of this gene. Thus, we investigated the association between the 3'UTR single-nucleotide polymorphisms (SNPs) in INHA and the laying performance in chickens and further explore their possible molecular cascades in granulosa cells (GC). Five SNPs were detected in Tianfu green-shell layers and g. 22,178,975 G > A was significantly associated with total egg numbers at the age of 300 days (EN, n = 286). Birds carrying the AA genotype laid more EN than those with GG (P < 0.05). The allele transition from G to A in the 3'UTR of INHA gene destroyed a binding site which was targeted by miR-181b-1-3p. The expression abundances of INHA mRNA increased firstly and then decreased with follicle growing, and reached the top in the sixth largest pre-ovulation follicle, whereas miR-181b-1-3p levels in chicken pre-hierarchical follicles had the contrary tendency. Further studies indicated that high levels of miR-181b-1-3p increased apoptosis and reduced GC proliferation while miR-181b-1-3p inhibitors decreased apoptosis and promoted GC proliferation. Additionally, depression of INHA increased apoptosis and reduced GC proliferation via a caspase-3-dependent mitochondrial pathway. Generally, the mutation in INHA 3'UTR was tightly correlated with egg production in chickens, and blocked a binding site of miR-181b-1-3p. miR-181b-1-3p inhibited GC proliferation and promoted apoptosis by targeting INHA.

Identification of Ovarian Circular RNAs and Differential Expression Analysis between MeiShan and Large White Pigs

MeiShan and Large White pigs differ in their female fecundity. However, the mechanisms behind the gene expression and regulation that cause these differences remain unclear. In this study, we profiled circRNAs and identified 5,879 circRNAs from the ovaries of MeiShan and Large White pigs. Eighty-five circRNAs were differentially expressed between the two pig breeds. Of these, 37 were up-regulated and 48 were down-regulated in MeiShan pigs. Gene ontology enrichment analysis suggested that the differentially expressed circRNA were involved in the hormone-mediated signaling pathway. We verified that circSCIN and its parent gene, scinderin (SCIN), were differentially expressed by reverse transcription and quantitative PCR (RT-qPCR). Luciferase assays demonstrated that circSCIN can target and sponge miR-133 and miR-148a/b. The identification of differentially expressed circRNAs (DECs) and their regulatory functions increased our understanding of the differences in reproductive efficiency between MeiShan and Large White pigs.

Comprehensive analysis of mRNAs and miRNAs in the ovarian follicles of uniparous and multiple goats at estrus phase

Background

Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique.

Results

The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells.

Conclusions

This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.

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

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.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA–miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm–egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

Genome-wide DNA methylation profiles of porcine ovaries in estrus and proestrus

DNA methylation is an important epigenetic modification involved in the estrous cycle and the regulation of reproduction. Here, we investigated the genome-wide profiles of DNA methylation in porcine ovaries in proestrus and estrus using methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in intergenic and intron regions. The methylation levels of coding regions were higher than those of the 5'- and 3'-flanking regions of genes. There were 4,813 differentially methylated regions (DMRs) of CpG islands in the estrus vs. proestrus ovarian genomes. Additionally, 3,651 differentially methylated genes (DMGs) were identified in pigs in estrus and proestrus. The DMGs were significantly enriched in biological processes and pathways related to reproduction and hormone regulation. We identified 90 DMGs associated with regulating reproduction in pigs. Our findings can serve as resources for DNA methylome research focused on porcine ovaries and further our understanding of epigenetically regulated reproduction in mammals.

Identification of Differentially Expressed Genes in Porcine Ovaries at Proestrus and Estrus Stages Using RNA-Seq Technique

Estrus is an important factor for the fecundity of sows, and it is involved in ovulation and hormone secretion in ovaries. To better understand the molecular mechanisms of porcine estrus, the expression patterns of ovarian mRNA at proestrus and estrus stages were analyzed using RNA sequencing technology. A total of 2,167 differentially expressed genes (DEGs) were identified (P ≤ 0.05, |log₂  Ratio| ≥ 1), of which 784 were upregulated and 1,383 were downregulated in the estrus compared with the proestrus group. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the cellular process, single-organism process, cell and cell part, and binding and metabolic process. In addition, a pathway analysis showed that these DEGs were significantly enriched in 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including cell adhesion molecules, ECM-receptor interaction, and cytokine-cytokine receptor interaction. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) confirmed the differential expression of 10 selected DEGs. Many of the novel candidate genes identified in this study will be valuable for understanding the molecular mechanisms of the sow estrous cycle.

Reference gene screening for analyzing gene expression across goat tissue

Real-time quantitative PCR (qRT-PCR) is one of the important methods for investigating the changes in mRNA expression levels in cells and tissues. Selection of the proper reference genes is very important when calibrating the results of real-time quantitative PCR. Studies on the selection of reference genes in goat tissues are limited, despite the economic importance of their meat and dairy products. We used real-time quantitative PCR to detect the expression levels of eight reference gene candidates (18S, TBP, HMBS, YWHAZ, ACTB, HPRT1, GAPDH and EEF1A2) in ten tissues types sourced from Boer goats. The optimal reference gene combination was selected according to the results determined by geNorm, NormFinder and Bestkeeper software packages. The analyses showed that tissue is an important variability factor in genes expression stability. When all tissues were considered, 18S, TBP and HMBS is the optimal reference combination for calibrating quantitative PCR analysis of gene expression from goat tissues. Dividing data set by tissues, ACTB was the most stable in stomach, small intestine and ovary, 18S in heart and spleen, HMBS in uterus and lung, TBP in liver, HPRT1 in kidney and GAPDH in muscle. Overall, this study provided valuable information about the goat reference genes that can be used in order to perform a proper normalisation when relative quantification by qRT-PCR studies is undertaken.