Molecular Evolution of the Bactericidal/Permeability-Increasing Protein (BPIFA1) Regulating the Innate Immune Responses in Mammals

Bactericidal/permeability-increasing protein, a primary factor of the innate immune system of mammals, participates in natural immune protection against invading bacteria. BPIFA1 actively contributes to host defense via multiple mechanisms, such as antibacterial, surfactant, airway surface liquid control, and immunomodulatory activities. However, the evolutionary history and selection forces on the BPIFA1 gene in mammals during adaptive evolution are poorly understood. This study examined the BPIFA1 gene of humans compared with that of other mammalian species to estimate the selective pressure derived by adaptive evolution. To assess whether or not positive selection occurred, we employed several different possibility tests (M1 vs. M2 and M7 vs. M8). The proportions of positively selected sites were significant, with a likelihood log value of 93.63 for the BPIFA1 protein. The Selecton server was used on the same dataset to reconfirm positive selection for specific sites by employing the Mechanistic-Empirical Combination model, thus providing additional evidence supporting the findings of positive selection. There was convincing evidence for positive selection signals in the BPIFA1 genes of mammalian species, which was more significant for selection signs and creating signals. We performed probability tests comparing various models based on dN/dS ratios to recognize specific codons under positive selection pressure. We identified positively selected sites in the LBP-BPI domain of BPIFA1 proteins in the mammalian genome, including a lipid-binding domain with a very high degree of selectivity for DPPC. BPIFA1 activates the upper airway's innate immune system in response to numerous genetic signals in the mammalian genome. These findings highlight evolutionary advancements in immunoregulatory effects that play a significant role in the antibacterial and antiviral defenses of mammalian species.

Oviduct Transcriptomic Reveals the Regulation of mRNAs and lncRNAs Related to Goat Prolificacy in the Luteal Phase

Simple Summary

The kidding number is an important reproductive trait in domestic goats. The oviduct, as one of the most major organs, is directly involved in the reproductive process, providing nutrition and a location for early embryonic development. The current study provides genome-wide expression profiles of mRNA and long noncoding RNAs (lncRNAs) expression in Yunshang black goat, a new breed of meat goat bred in China with a high kidding number. During the luteal phases, oviduct mRNAs and lncRNAs associated with high- and low-fecundity Yunshang black goats were identified, and their potential biological functions were predicted using GO, KEGG, and GSEA enrichment analysis. These findings shed light on the oviduct-based prolificacy mechanism in goats.

Abstract

The oviduct is associated with embryo development and transportation and regulates the pregnancy success of mammals. Previous studies have indicated a molecular mechanism of lncRNAs in gene regulation and reproduction. However, little is known about the function of lncRNAs in the oviduct in modulating goat kidding numbers. Therefore, we combined RNA sequencing (RNA-seq) to map the expression profiles of the oviduct at the luteal phase from high- and low-fecundity goats. The results showed that 2023 differentially expressed mRNAs (DEGs) and 377 differentially expressed lncRNAs (DELs) transcripts were screened, and 2109 regulated lncRNA-mRNA pairs were identified. Subsequently, the genes related to reproduction (IGF1, FGFRL1, and CREB1) and those associated with embryonic development and maturation (DHX34, LHX6) were identified. KEGG analysis of the DEGs revealed that the GnRH- and prolactin-signaling pathways, progesterone-mediated oocyte maturation, and oocyte meiosis were related to reproduction. GSEA and KEGG analyses of the target genes of DELs demonstrated that several biological processes and pathways might interact with oviduct functions and the prolificacy of goats. Furthermore, the co-expression network analysis showed that XLOC_029185, XLOC_040647, and XLOC_090025 were the cis-regulatory elements of the DEGs MUC1, PPP1R9A, and ALDOB, respectively; these factors might be associated with the success of pregnancy and glucolipid metabolism. In addition, the GATA4, LAMA2, SLC39A5, and S100G were trans-regulated by lncRNAs, predominantly mediating oviductal transport to the embryo and energy metabolism. Our findings could pave the way for a better understanding of the roles of mRNAs and lncRNAs in fecundity-related oviduct function in goats.

TIMP1 may affect goat prolificacy by regulating biological function of granulosa cells

Tissue inhibitor of metalloproteinase 1 (TIMP1) is associated with animal reproductive processes, such as follicular growth, ovulation, luteinization, and embryo development in mammals. The purposes of this study were to explore the expression and localization of TIMP1 in the ovarian tissues and determine the effect of TIMP1 on the function of granulosa cells and the association of TIMP1 with lambing-related genes of the goats. Immunohistochemical analysis showed that TIMP1 protein was strongly expressed by granulosa cells. Enzyme-linked immunosorbent assay (ELISA) results showed that TIMP1 overexpression promoted the secretion of estradiol of granulosa cells after 12, 24, and 48 h of transfection. Moreover, in vitro experiments indicated that TIMP1 had the ability to promote the cell proliferation and elevate the transcriptional levels of four genes associated with goat prolificacy, including BMPR-1B, BMP15, GDF9, and FSHB, in granulosa cells. In conclusion, TIMP1 could be an important molecule in regulating reproductive performance of the goats by affecting estrogen secretion and cell proliferation, as well as the expression of lambing-related genes of granulosa cells in the goats.

Genome-Wide Identification of mRNAs, lncRNAs, and Proteins, and Their Relationship With Sheep Fecundity

The exploration of multiple birth-related genes has always been a significant focus in sheep breeding. This study aimed to find more genes and proteins related to the litter size in sheep. Ovarian specimens of Small Tail Han sheep (multiple births) and Xinji Fine Wool sheep (singleton) were collected during the natural estrus cycle. Transcriptome and proteome of ovarian specimens were analyzed. The transcriptome results showed that "steroid hormone biosynthesis" and "ovarian steroidogenesis" were significantly enriched, in which HSD17B1 played an important role. The proteome data also confirmed that the differentially expressed proteins (DEPs) were enriched in the ovarian steroidogenesis pathway, and the CYP17A1 was the candidate DEP. Furthermore, lncRNA MSTRG.28645 was highly expressed in Small Tailed Han sheep but lowly expressed in Xinji fine wool sheep. In addition, MSTRG.28645, a hub gene in the co-expression network between mRNAs and lncRNAs, was selected as one of the candidate genes for subsequent verification. Expectedly, the overexpression and interference of HSD17B1 and MSTRG.28645 showed a significant effect on hormone secretion in granulosa cells. Therefore, this study confirmed that HSD17B1 and MSTRG.28645 might be potential genes related to the fecundity of sheep. It was concluded that both HSD17B1 and MSTRG.28645 were critical regulators in the secretion of hormones that affect the fecundity of the sheep.

Characterization of XR_311113.2 as a MicroRNA Sponge for Pre-ovulatory Ovarian Follicles of Goats via Long Noncoding RNA Profile and Bioinformatics Analysis

Long noncoding RNAs (lncRNAs) were identified recently as a large class of noncoding RNAs (ncRNAs) with a length ≥200 base pairs (bp). The function and mechanism of lncRNAs have been reported in a growing number of species and tissues. In contrast, the regulatory mechanism of lncRNAs in the goat reproductive system has rarely been reported. In the present study, we sequenced and analyzed the lncRNAs using bioinformatics to identify their expression profiles. As a result, 895 lncRNAs were predicted in the pre-ovulatory ovarian follicles of goats. Eighty-eight lncRNAs were differentially expressed in the Macheng black goat when compared with Boer goat. In addition, the lncRNA XR_311113.2 acted as a sponge of chi-miR-424-5p, as assessed via a luciferase activity assay. Taken together, our findings demonstrate that lncRNAs have potential effects in the ovarian follicles of goats and may represent a promising new research field to understand follicular development.

Hypothalamic Transcriptome Analysis Reveals the Crucial MicroRNAs and mRNAs Affecting Litter Size in Goats

The hypothalamus was the coordination center of the endocrine system, which played an important role in goat reproduction. However, the molecular mechanism of hypothalamus regulating litter size in goats was still poorly understood. This study aims to investigate the key functional genes associated with prolificacy by hypothalamus transcriptome analysis of goats. In this research, an integrated analysis of microRNAs (miRNAs)-mRNA was conducted using the hypothalamic tissue of Yunshang black goats in the follicular stage. A total of 72,220 transcripts were detected in RNA-seq. Besides, 1,836 differentially expressed genes (DEGs) were identified between high fecundity goats at the follicular phase (FP-HY) and low fecundity goats at the follicular phase (FP-LY). DEGs were significantly enriched in 71 Gene Ontology (GO) terms and 8 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The transcriptome data suggested that DEGs such as BMPR1B, FGFR1, IGF1 and CREB1 are directly or indirectly involved in many processes like hypothalamic gonadal hormone secretion. The miRNA-seq identified 1,837 miRNAs, of which 28 differentially expressed miRNAs (DEMs). These DEMs may affect the nerve cells survival of goat hypothalamic regulating the function of target genes and further affect the hormone secretion activities related to reproduction. They were enriched in prolactin signaling pathway, Jak-STAT signaling pathway and GnRH signaling pathway, as well as various metabolic pathways. Integrated analysis of DEMs and DEGs showed that 87 DEGs were potential target genes of 28 DEMs. After constructing a miRNA-mRNA pathway network, we identified several mRNA-miRNAs pairs by functional enrichment analysis, which was involved in hypothalamic nerve apoptosis. For example, NTRK3 was co-regulated by Novel-1187 and Novel-566, as well as another target PPP1R13L regulated by Novel-566. These results indicated that these key genes and miRNAs may play an important role in the development of goat hypothalamus and represent candidate targets for further research. This study provides a basis for further explanation of the basic molecular mechanism of hypothalamus, but also provides a new idea for a comprehensive understanding of prolificacy characteristics in Yunshang black goats.

Identification and Comparative Analysis of Long Non-coding RNAs in High- and Low-Fecundity Goat Ovaries During Estrus

The ovary is the most important reproductive organ in goats and directly affects the fecundity. Long non-coding RNAs (lncRNAs) are involved in the biological process of oocyte maturation. However, in the context of reproduction in goats, few studies have explored the regulation of lncRNAs. Therefore, we herein used the ovaries of high and low fecundity Leizhou black goats to identify differentially expressed lncRNAs (DElncRNAs) by high-throughput RNA sequencing; moreover, we analyzed the target genes of lncRNAs by functional annotation to explore the role of DElncRNAs in ovarian development. Twenty DElncRNAs were identified, of which six were significantly upregulated and 14 were significantly downregulated in high fecundity goats. Gene Ontology analyses suggested that MSTRG.3782 positively influences the expression of the corresponding gene API5, exerting regulative effects on the development of follicles, through which litter size might show variations. The target gene KRR1 of ENSCHIT00000001883 is significantly enriched in cell components, and ENSCHIT00000001883 may regulate cell growth and thus affect follicular development. Further, as per Kyoto Encyclopedia of Genes and Genomes pathway analyses, MSTRG.2938 was found to be significantly enriched, and we speculate that MSTRG.2938 could regulate ribosomal biogenesis in the pre-snoRNP complex as well as cell transformation in eukaryotes. Quantitative real-time PCR results were consistent with sequencing data. To conclude, our research results indicate that some lncRNAs play a key role in regulating follicle development and cell growth during goat’ s ovarian development.

A Preliminary Study on the Characteristics of microRNAs in Ovarian Stroma and Follicles of Chuanzhong Black Goat during Estrus

microRNAs (miRNAs) play a significant role in ovarian follicular maturity, but miRNA expression patterns in ovarian stroma (OS), large follicles (LF), and small follicles (SF) have been rarely explored. We herein aimed to identify miRNAs, their target genes and signaling pathways, as well as their interaction networks in OS, LF, and SF of Chuanzhong black goats at the estrus phase using small RNA-sequencing. We found that the miRNA expression profiles of LF and SF were more similar than those of OS—32, 16, and 29 differentially expressed miRNAs were identified in OS vs. LF, OS vs. SF, and LF vs. SF, respectively. Analyses of functional enrichment and the miRNA-targeted gene interaction network suggested that miR-182 (SMC3), miR-122 (SGO1), and miR-206 (AURKA) were involved in ovarian organogenesis and hormone secretion by oocyte meiosis. Furthermore, miR-202-5p (EREG) and miR-485-3p (FLT3) were involved in follicular maturation through the MAPK signaling pathway, and miR-2404 (BMP7 and CDKN1C) played a key role in follicular development through the TGF-β signaling pathway and cell cycle; nevertheless, further research is warranted. To our knowledge, this is the first study to investigate miRNA expression patterns in OS, LF, and SF of Chuanzhong black goats during estrus. Our findings provide a theoretical basis to elucidate the role of miRNAs in follicular maturation. These key miRNAs might provide candidate biomarkers for the diagnosis of follicular maturation and will assist in developing new therapeutic targets for female goat infertility.

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.

Adaptation to host-specific bacterial pathogens drive rapid evolution of novel PhoP/PhoQ regulation pathway modulating the virulence

The presence of the PhoP-PhoQ system is usually different in various bacterial groups, suggesting that PhoP can control the expression of different genes in species. However, little is known about the evolution of the PhoP-PhoQ system among bacterial pathogens. Here, we study the evolution of PhoP and PhoQ regulation in 15 species of Enterobacteriaceae family. We have determined that the regulatory objectives adopted by PhoP and PhoQ are mainly different, due to the result of horizontal gene transfer events and even the change in the genetic content between closely related species. We have compared many possibilities tests (M1 vs. M2 and M7 with M8) to determine the positive selection. Estimating parameters at M1 and M2, with positive selection in M2 of the two proteins. The proportions of positive selection sites significant with ω = 4.53076 for PhoP and ω = 4.21041 PhQ. M8 was significant for PhoP and PhQ proteins. To further confirm the positive selection results, we used the Selecton server to confer positive selection on individual sites using the Mechanistic-Empirical Combination model, and we noticed that several sites had been identified under selection pressure during the evolution. There was a strong indication for the positive selection in bacterial genes of PhoP and PhoQ showed the results. By the use of REL and IFEL, the positive selection for PhoP was detected 14 and 11 sites respectively at different codon positions. The positively selected sites of amino acids such as Arginine, Alanine, Lysine, and Leucine are more important for the production of signals. Our results suggest that the positive selection of PhoP-PhoQ genes in host adaptation during evolution raises an intriguing possibility causes subtle variations in actions of PhoP-PhoQ and also increases the opportunities that cause modification in protein structure for the evolution of increasing pathogenicity in bacterial pathogens.

Cloning and expression analysis of the follicle-stimulating hormone receptor (FSHR) gene in the reproductive axis of female yaks (Bos grunniens)

Follicle-stimulating hormone receptor (FSHR) plays a central role in promoting follicle maturation through the follicle-stimulating hormone (FSH)-mediated cAMP pathway in animals. The objectives of the present study were to clone the FSHR gene of yaks (Bos grunniens) and compare differences in FSHR mRNA expression in the reproductive axis between yaks and cattle. Hypothalamus, anterior pituitary, oviduct, ovary, and uterus tissue samples were collected from adult female yaks (n = 5) and cattle (n = 5) during the follicular phase. Using reverse transcriptase-polymerase chain reaction (RT-PCR), we found that the FSHR coding region of the yak is 2088 bp and encodes 695 amino acids. Its amino acid sequence showed 99.38%-72.22% similarity to the homologous genes of cattle, goats, sheep, cats, donkeys, horses, humans, chickens, monkeys, mice, rats, and wild boar. Real-time PCR analysis revealed that the FSHR gene was expressed in all tissues examined. Expression of the FSHR gene in the yak was higher in the uterus than other tissues (P < 0.05) but, in cattle, was higher in the ovary than other tissues (P < 0.05). The FSHR gene expression level in the cattle ovary was significantly higher than that in the yak ovary (P < 0.01). These results indicate that the FSHR gene is relatively conserved in the course of animal evolution. The variation in sequence and expression level of FSHR between the two species might be associated with the difference in their reproduction.

Comparison of the sequences and expression levels of genes related to follicular development and atresia between prolific and nonprolific goat breeds

This study investigated the variations of the nucleotide sequences and ovarian expression levels of genes related to follicular development and atresia in prolific Jintang black goats and nonprolific Tibetan goats. Eight genes, FSHB, LHB, FSHR, LHCGR, ESR2, B4GANT2, BCL2 and BAX, were examined using reverse transcription‐polymerase chain reaction and quantitative real‐time PCR. The results showed that the nucleotide and deduced amino acid sequences of the LHB and BAX genes were not different, but there was one base change in the FSHR genes between the two breeds. There was one base change in the FSHB gene, which resulted in one amino acid substitution; there were nine base changes in the LHCGR gene, which resulted in five amino acid substitutions; and there were six base changes in the B4GANT2 gene, which resulted in four amino acid substitutions. The expression levels of the FSHR, LHCGR, ESR2, B4GANT2, BCL2 and BAX genes in the ovaries were not different between the two breeds. The plasma concentrations of FSH were not different, but the plasma concentrations of LH, P₄ and E₂ were lower in prolific Jintang black goats than in nonprolific Tibetan goats (P ˂ 0.05) at 40 hr after removal of the Controlled Internal Drug Release Devices. These results provide some foundations elucidating the endocrine and molecular mechanisms controlling ovulation rate in goats, but these need to be further verified.

Evolutionary Analysis of Makorin Ring Finger Protein 3 Reveals Positive Selection in Mammals

Makorin ring finger proteins (MKRNs) are part the of ubiquitin-proteasome system; a complex system important for cell functions. Ubiquitin fate through proteolytic, non-proteolytic pathways varies, depending on covalent linkage between ubiquitin and protein substrates. Makorin ring finger protein 3 is an integral part of covalent linkage of ubiquitin to protein substrates. Similar to others imprinted genes, MKRN3 also evolve under positive selection; however, which codons are specifically selected in MKRN3 during evolution are needed to be explored. Different maximum-likelihood (ML) codon-based methodologies were used to ascertain positive selection signatures in 22 mammalian sequences of MKRN3 to probe an individual codon for positive selection signatures. By applying the HyPhy software package implemented in the Data Monkey Web Server and CODEML implemented in PAML, evolutionary analysis based on two Ml frameworks were conducted. The analysis was executed by comparing M1a against M2a, M7 against M8, and PAML models and 2∆Lnl (LRT) was resulted by likelihood logs. M1a contributed ω1 (dN/dS) with LRT value (∆Lnl) 12.01, and positive selection was found in M2a with ω3 = 2.23603. To further improve selection test, M8 was compared to M7 with 2∆Lnl (LRT) 30.17, and M8 showed positive selection with ω = 1.55759. The data were fit to M8 than M7, which suggests that M8 was the most significant model of selection. M8 was judged encouraging for this analysis and used to establish a positive selection of MKRN3 proteins. We found Gly312 as a positively selected amino acid in a zinc finger motif/Really Interesting New Gene (RING) finger motif; the former ones’ region is involved in RNA binding and the later ones in ubiquitin ligase activity of the protein, vital for protein function. Selection analyses of MKRNs might advance the developments in unique approaches that could lead to genetic progress over the selection of superior individuals with the breeding values higher for certain traits as ancestries to get the next generation.

Expression Analysis of the Prolific Candidate Genes, BMPR1B, BMP15, and GDF9 in Small Tail Han Ewes with Three Fecundity (FecB Gene) Genotypes

Simple Summary

As important prolific candidate genes, BMPR1B, BMP15, and GDF9 may affect the lambing performance of sheep. Therefore, regarding the three FecB genotypes of Small Tail Han (STH) sheep (FecB BB, FecB B+, and FecB ++), this study explored the gene expression characteristics of different tissues using reverse transcription PCR (RT-PCR) and real-time quantitative PCR (qPCR). The results showed that BMPR1B, BMP15, and GDF9 expression differed between the selected tissues, with all being highly expressed in the ovaries. Further analysis indicated that there was no significant difference in BMPR1B expression among the three FecB genotypes, but both GDF9 and BMP15 had the highest expression in FecB B+. As for other non-ovarian tissues, expression also varied. This study is relevant to understanding the high prolificacy of the STH breed.

Abstract

The expression characteristics of the prolific candidate genes, BMPR1B, BMP15, and GDF9, in the major visceral organs and hypothalamic–pituitary–gonadal (HPG) axis tissues of three FecB genotypes (FecB BB, FecB B+, and FecB ++) were explored in STH ewes using RT-PCR and qPCR. The results were as follows, BMPR1B was expressed in all FecB BB genotype (Han BB) tissues, and GDF9 was expressed in all selected tissues, but BMP15 was specifically expressed in the ovaries. Further study of ovarian expression indicated that there was no difference in BMPR1B expression between genotypes, but the FecB B+ genotype (Han B+) had greater expression of GDF9 and BMP15 than Han BB and FecB ++ genotype (Han ++) (p < 0.05, p < 0.01). BMP15 expression was lower in the ovaries of Han BB than in Han ++ sheep, but the reverse was shown for GDF9. The gene expression in non-ovarian tissues was also different between genotypes. Therefore, we consider that the three genes have an important function in ovine follicular development and maturation. This is the first systematic analysis of the tissue expression pattern of BMPR1B, BMP15, and GDF9 genes in STH sheep of the three FecB genotypes. These results contribute to the understanding of the molecular regulatory mechanism for ovine reproduction.

Molecular characterization, expression profile of the FSHRgene and its association with egg production traits in muscovy duck

Follicle-stimulating hormone (FSH) and its receptor play a key role in the follicular development and regulation of steroidogenesis in the ovary and spermatogenesis in the testis. The purpose of this study was to characterize themuscovy duck FSHR gene, identify SNPs and their association with egg production traits in muscovy ducks. Here, we cloned the complementary DNA (cDNA) sequence of FSHR, and examined the expression patterns of FSHR gene in adult female muscovy duck tissues. The cloned cDNA of the muscovy duck FSHR gene shared high similarity to those of pekin duck (Anas platyrhynchos) (95.7%) and chicken (93.2%). Three different muscovy duck FSHR transcripts were identified. Quantitative real-time PCR (RT-qPCR) results showed that the FSHR gene was expressed in all the 14 tested tissues, and the highest expression level was seen in the ovary. A total of 16 SNPs were identified, among which, four SNPs were located in the coding region of FSHR. The SNP C320T is significantly associated with egg production at 59 weeks of age (P < 0.05), whereas the SNP A227G is significantly associated with age at first egg stage (P < 0.05). These results suggest that the two SNPs (A227G and C320T) of FSHR gene are associated with egg production traits and could be potential markers that can be used for marker-assisted selection programmes to increase egg production in muscovy duck.

Positive selection drives the evolution of endocrine regulatory bone morphogenetic protein system in mammals

The rapid evolution of reproductive proteins might be driven by positive Darwinian selection. The bone morphogenetic protein family is the largest within the transforming growth factor (TGF) superfamily. A little have been known about the molecular evolution of bone morphogenetic proteins exhibiting potential role in mammalian reproduction. In this study we investigated mammalian bone morphogenetic proteins using maximum likelihood approaches of codon substitutions to identify positive Darwinian selection in various species. The proportion of positively selected sites was tested by different likelihood models for individual codon, and M8 were found to be the best model. The percentage of positively elected sites under M8 are 2.20% with ω = 1.089 for BMP2, 1.6% with ω = 1.61 for BMP 4 0.53% for BMP15 with ω = 1.56 and 0.78% for GDF9 with ω = 1.93. The percentage of estimated selection sites under M8 is strong statistical confirmation that divergence of bone morphogenetic proteins is driven by Darwinian selection. For the proteins, model M8 was found significant for all proteins with ω > 1. To further test positive selection on particular amino acids, the evolutionary conservation of amino acid were measured based on phylogenetic linkage among sequences. For exploring the impact of these somatic substitution mutations in the selection region on human cancer, we identified one pathogenic mutation in human BMP4 and one in BMP15, possibly causing prostate cancer and six neutral mutations in BMPs. The comprehensive map of selection results allows the researchers to perform systematic approaches to detect the evolutionary footprints of selection on specific gene in specific species.

The regulation of IMF deposition in pectoralis major of fast- and slow- growing chickens at hatching

Background

The lipid from egg yolk is largely consumed in supplying the energy for embryonic growth until hatching. The remaining lipid in the yolk sac is transported into the hatchling’s tissues. The gene expression profiles of fast- and slow-growing chickens, Arbor Acres (AA) and Beijing-You (BJY), were determined to identify global differentially expressed genes and enriched pathways related to lipid metabolism in the pectoralis major at hatching.

Results

Between these two breeds, the absolute and weight-specific amounts of total yolk energy (TYE) and intramuscular fat (IMF) content in pectoralis major of fast-growing chickens were significantly higher (P < 0.01, P < 0.01, P < 0.05, respectively) than those of the slow-growing breed. IMF content and u-TYE were significantly related (r = 0.9047, P < 0.01). Microarray analysis revealed that gene transcripts related to lipogenesis, including PPARG, RBP7, LPL, FABP4, THRSP, ACACA, ACSS1, DGAT2, and GK, were significantly more abundant in breast muscle of fast-growing chickens than in slow-growing chickens. Conversely, the abundance of transcripts of genes involved in fatty acid degradation and glycometabolism, including ACAT1, ACOX2, ACOX3, CPT1A, CPT2, DAK, APOO, FUT9, GCNT1, and B4GALT3, was significantly lower in fast-growing chickens. The results further indicated that the PPAR signaling pathway was directly involved in fat deposition in pectoralis major, and other upstream pathways (Hedgehog, TGF-beta, and cytokine–cytokine receptor interaction signaling pathways) play roles in its regulation of the expression of related genes.

Conclusions

Additional energy from the yolk sac is transported and deposited as IMF in the pectoralis major of chickens at hatching. Genes and pathways related to lipid metabolism (such as PPAR, Hedgehog, TGF-beta, and cytokine–cytokine receptor interaction signaling pathways) promote the deposition of IMF in the pectoralis major of fast-growing chickens compared with those that grow more slowly. These findings provide new insights into the molecular mechanisms underlying lipid metabolism and deposition in hatchling chickens.

Electronic supplementary material

The online version of this article (10.1186/s40104-017-0207-z) contains supplementary material, which is available to authorized users.

Maximum-likelihood approaches reveal signatures of positive selection in BMP15 and GDF9 genes modulating ovarian function in mammalian female fertility

Bone morphogenetic proteins (BMPs) and the growth factors (GDFs) play an important role in ovarian folliculogenesis and essential regulator of processes of numerous granulosa cells. BMP15 gene variations linked to various ovarian phenotypic consequences subject to the species, from infertility to improved prolificacy in sheep, primary ovarian insufficiency in women or associated with minor subfertility in mouse. To study the evolving role of BMP15 and GDF9, a phylogenetic analysis was performed. To find out the candidate gene associated with prolificacy in mammals, the nucleotide sequence of BMP15 and GDF9 genes was recognized under positive selection in various mammalian species. Maximum‐likelihood approaches used on BMP15 and GDF9 genes exhibited a robust divergence and a prompted evolution as compared to other TGFβ family members. Furthermore, among 32 mammalian species, we identified positive selection signals in the hominidae clade resulting to 132D, 147E, 163Y, 191W, and 236P codon sites of BMP15 and 162F, 188K, 206R, 240A, 244L, 246H, 248S, 251D, 253L, 254F and other codon sites of GDF9. The positively selected amino acid sites such as Alanine, Lucien, Arginine, and lysine are important for signaling. In conclusion, this study evidences that GDF9 and BMP15 genes have rapid evolution than other TGFß family members and was subjected to positive selection in the mammalian clade. Selected sites under the positive selection are of remarkable significance for the particular functioning of the protein and consequently for female fertility.

Molecular characterization of the Chinese alligator follicle-stimulating hormone β subunit (FSHβ) and its expression during the female reproductive cycle

The Chinese alligator Alligator sinensis is an endangered species endemic to China, it has a highly specialized reproductive pattern with low fecundity. Up to date, little is known about the regulation of its female reproductive cycle. Follicle-stimulating hormone (FSH), a glycoprotein hormone, plays a key role in stimulating and regulating ovarian follicular development and egg production. In this study, the complete FSHβ cDNA from the ovary of the Chinese alligator was obtained for the first time, it consists of 843-bp nucleotides, including 120-bp nucleotides of the 5'-untranslated region (UTR), 396-bp of the open reading frame, and 3'-UTR of 327-bp nucleotides. It encodes a 131-amino acid precursor molecule of FSHβ with a signal peptide of 18 amino acids followed by a mature protein of 113 amino acids. Its deduced amino acid sequence shares high identities with the American alligator (100%) and birds (89-92%). Phylogenetic tree analysis of the FSHβ amino acid sequence indicated that alligators cluster into the bird branch. Tissue distribution analyses indicated that FSHβ mRNA is expressed in ovary, intestine and liver with the highest level in the ovary, while not in stomach, pancreas, heart, thymus and thyroid. Expression of FSHβ in ovary increases in May (breeding prophase) and peaks in July (breeding period), it is maintained at high levels through September, then decreases significantly in November (post-reproductive period) and remains relatively low from January to March (hibernating period). These temporal changes of FSHβ expression implicated that it might play an important role in promoting ovarian development during the female reproductive cycle.

Comparative messenger RNA expression of FSHβ, LHβ, FSHR, LHR, and ERβ in high and low prolific goat breeds

Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) have a central role in follicle growth and maturation, but no clear differences between breeds with different ovulation rates have been found. Therefore, this study investigated mRNA expression of FSHβ, LHβ, FSH receptor (FSHR), LH receptor (LHR), and estrogen receptor-β (ERβ) genes in prolific Lezhi black (LB) goats and nonprolific Tibetan (TB) goats by real-time PCR. Follicles and pituitaries were recovered from goats at 12-24 h after onset of estrus. Real-time PCR analysis revealed that the expression levels of FSHβ and LHβ mRNA were significantly higher (p < 0.01) in pituitary of LB than in TB does, but the expression levels of FSHR and LHR mRNA in follicle of TB were greater (p < 0.05). Expression level of follicular ER β was not different between the two breeds. Data provide evidence that the greater ovulation rate in the LB goat as compared to the TB breed is associated with a greater gonadotropin expression during follicular phase.

Differential ovarian morphometry and follicular expression of BMP15, GDF9 and BMPR1B influence the prolificacy in goat

The variation in the kidding size of Black Bengal and Sirohi breed of goats makes them an interesting genetic material to study the underlying genetic mechanism of prolificacy. Accordingly, we studied the comparative ovarian morphometry including disparity in numbers of antral follicles of different sizes between these two breeds. Further, we evaluated the differential expression of the important candidate genes (viz., BMP15, GDF9 and BMPR1B) known to influence the ovulation rates and the prolificacy. The ovaries of Black Bengal (n = 20) goat were lighter (p < 0.01) in weight and smaller (p < 0.01) in diameter than those of Sirohi (n = 19) goats but possessed more numbers (p < 0.05) of corpus luteum (CL), large and small antral follicles. Quantitative real-time PCR (RT-qPCR) analysis revealed differential expression of mRNAs encoding for the BMP15 and GDF9. Small antral follicles of Black Bengal goats expressed 2.78-fold more (p < 0.05) of BMP 15 than those of Sirohi goat. Expression of BMP15 (p < 0.01) and GDF9 (p < 0.05) mRNAs was more abundant in the small than the large antral follicles of Black Bengal goat. The more numbers of antral follicles per unit of ovarian mass and differential expression of BMP15 and GDF9 may serve as an important clue for higher prolificacy.

FSH stimulates lipid biosynthesis in chicken adipose tissue by upregulating the expression of its receptor FSHR

Transcripts and protein for follicle-stimulating hormone receptor (FSHR) were demonstrated in abdominal adipose tissue of female chickens. There was no expression of the Fsh gene, but FSH and FSHR colocalized, suggesting that FSH was receptor bound. Partial correlations indicted that changes in abdominal fat (AF) content were most directly correlated with Fshr mRNA expression, and the latter was directly correlated with tissue FSH content. These relationships were consistent with FSH inducing Fshr mRNA expression and with the finding that FSH influenced the accumulation of AF in chickens, a novel role for the hormone. Chicken preadipocytes responded linearly to doubling concentrations of FSH in Fshr mRNA expression and quantities of FSHR and lipid, without discernable effect on proliferation. Cells exposed to FSH more rapidly acquired adipocyte morphology. Treatment of young chickens with chicken FSH (4 mIU/day, subcutaneous, days 7-13) did not significantly decrease live weight but increased AF weight by 54.61%, AF as a percentage of live weight by 55.45%, and FSHR transcripts in AF by 222.15% (2 h after injection). In cells stimulated by FSH, genes related to lipid metabolism, including Rdh10, Dci, RarB, Lpl, Acsl3, and Dgat2, were expressed differentially, compared with no FSH. Several pathways of retinal and fatty acid metabolism, and peroxisome proliferator-activated receptor (PPAR) signaling changed. In conclusion, FSH stimulates lipid biosynthesis by upregulating Fshr mRNA expression in abdominal adipose tissue of chickens. Several genes involved in fatty acid and retinal metabolism and the PPAR signaling pathway mediate this novel function of FSH.

Ovarian follicle development in Booroola sheep exhibiting impaired bone morphogenetic protein signalling pathway

The role of bone morphogenetic proteins (BMPs) in the regulation of ovarian function has been extensively investigated but the mechanism of regulation is not well understood. The aim of this study was to investigate the effect of mutation in the BMP receptor in Booroola sheep on the number of primordial follicles and rate of follicle recruitment in comparison with that in normal merino sheep in vivo. Whole sheep ovaries at the time of birth, 1.5 and 5 years old were collected and processed for the follicle quantification, using computerised stereological methods and statistical analyses. At birth, the total number of primordial follicles in Booroola sheep was significantly lower than in merino sheep. At 1.5 and 5 years, a reversed pattern in favour of Booroola ewes was seen with significantly more primordial follicles than merino. In parallel, the rate of primordial follicle recruitment to developing cohort was substantially lower in Booroola ewes with only 51 and 66% of primordial follicle consumption at 1.5 and 5 years respectively compared to 92 and 97% in merino ewes. On other hand, the mean numbers of developing primary follicles were smaller in Booroola sheep at the time of birth, yet, Booroola ewes possess more primary follicles than merino at 1.5 years. These findings suggest that attenuation of the intraovarian signalling pathway of BMPs may in fact be a successful means of rationalising follicle consumption, preventing unnecessary loss of follicles from the initial primordial follicle pool, hence increasing reproductive longevity and fertility.