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

Abstract

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.

REFERENCES Ahmad, H. I., Liu, G., Jiang, X., Liu, C., Chong, Y., & Huarong, H. (2017). Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations. Ecology and Evolution, 7, 5170–5180. https://doi.org/10.1002/ece3.2919 Ahmad, H. I., Liu, G., Jiang, X., Liu, C., Fangzheng, X., Chong, Y., … Huarong, H. (2017). Adaptive selection at agouti gene inferred breed specific selection signature within the indigenous goat populations. Asian-Australasian Journal of Animal Sciences, 00, 1–7. https://doi.org/10.5713/ajas.16.0994 Asif, A. R., Awais, M., Qadri, S., Ahmad, H. I., & Du, X. (2017). Positive selection of IL-33 in adaptive immunity of domestic Chinese goats. Ecology and Evolution, 7, 1954–1963. Auclair, S., Rossetti, R., Meslin, C., Monestier, O., Di Pasquale, E., Pascal, G., … Fabre, S. (2013). Positive selection in bone morphogenetic protein 15 targets a natural mutation associated with primary ovarian insufficiency in human. PLoS ONE, 8, e78199. Bergström, T., & Gyllensten, U. (1995). Evolution of Mhc class II polymorphism: The rise and fall of class II gene function in primates. Immunological Reviews, 143, 13–31. Bielawski, J. P., & Yang, Z. (2003). Maximum likelihood methods for detecting adaptive evolution after gene duplication. Journal of Structural and Functional Genomics; 3 201–212. de Castro, F. C., Cruz, M. H. C., & Leal, C. L. V. (2016). Role of growth differentiation factor 9 and bone morphogenetic protein 15 in ovarian function and their importance in mammalian female fertility—A review. Asian-Australasian Journal of Animal Sciences, 29, 1065. Ceko, M. J., Hummitzsch, K., Hatzirodos, N., Bonner, W. M., Aitken, J. B., Russell, D. L., … Harris, H. H. (2014). X-Ray fluorescence imaging and other analyses identify selenium and GPX1 as important in female reproductive function. Metallomics, 7, 71–82. Chang, H., Brown, C. W., & Matzuk, M. M. (2002). Genetic analysis of the mammalian transforming growth factor-β superfamily. Endocrine Reviews, 23, 787–823. Cui, H. X., Zhao, S. M., Cheng, M. L., Guo, L., Ye, R. Q., Liu, W. Q., & Gao, S. Z. (2009). Cloning and expression levels of genes relating to the ovulation rate of the Yunling black goat. Biology of Reproduction, 80, 219–226. Dias, F. C. F., Khan, M. I. R., Adams, G. P., Sirard, M. A., & Singh, J. (2014). Granulosa cell function and oocyte competence: Super-follicles, super-moms and super-stimulation in cattle. Animal Reproduction Science, 149, 80–89. Franceschini, A., Szklarczyk, D., Frankild, S., Kuhn, M., Simonovic, M., Roth, A., … Jensen, L. J. (2012). STRING v9. 1: Protein–protein interaction networks, with increased coverage and integration. Nucleic Acids Research, 41, D808–D815. Gharib, W. H., & Robinson-Rechavi, M. (2013). The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC. Molecular Biology and Evolution, 30, 1675–1686. Glaser, F., Pupko, T., Paz, I., Bell, R. E., Bechor-Shental, D., Martz, E., & Ben-Tal, N. (2003). ConSurf: Identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics, 19, 163–164. Hatzirodos, N., Hummitzsch, K., Irving-Rodgers, H. F., Harland, M. L., Morris, S. E., & Rodgers, R. J. (2014). Transcriptome profiling of granulosa cells from bovine ovarian follicles during atresia. BMC Genomics, 15, 40. He, Y., Ma, X., Liu, X., Zhang, C., & Li, J. (2010). Candidate genes polymorphism and its association to prolificacy in Chinese goats. Journal of Agricultural Science, 2, 88. Kelley, L. A., & Sternberg, M. J. (2009). Protein structure prediction on the Web: A case study using the Phyre server. Nature Protocols, 4, 363–371. Li, H., Zhao, X., Wang, J., Zong, M., & Yang, H. (2017). Bioinformatics analysis of gene expression profile data to screen key genes involved in pulmonary sarcoidosis. Gene, 596, 98–104. Marini, N. J., Thomas, P. D., & Rine, J. (2010). The use of orthologous sequences to predict the impact of amino acid substitutions on protein function. PLoS Genetics, 6, e1000968. Meslin, C., Mugnier, S., Callebaut, I., Laurin, M., Pascal, G., Poupon, A., … Monget, P. (2012). Evolution of genes involved in gamete interaction: Evidence for positive selection, duplications and losses in vertebrates. PLoS ONE, 7, e44548. Peng, J., Li, Q., Wigglesworth, K., Rangarajan, A., Kattamuri, C., Peterson, R. T., … Matzuk, M. M. (2013). Growth differentiation factor 9: Bone morphogenetic protein 15 heterodimers are potent regulators of ovarian functions. Proceedings of the National Academy of Sciences, 110, E776–E785. Persani, L., Rossetti, R., Di Pasquale, E., Cacciatore, C., & Fabre, S. (2014). The fundamental role of bone morphogenetic protein 15 in ovarian function and its involvement in female fertility disorders. Human Reproduction Update, 20, 869–883. Poon, A. F., Frost, S. D., & Pond, S. L. K. (2009). Detecting signatures of selection from DNA sequences using Datamonkey. Bioinformatics for DNA Sequence Analysis, 537, 163–183. Pulkki, M. M., Myllymaa, S., Pasternack, A., Lun, S., Ludlow, H., Al-Qahtani, A., … Laitinen, M. (2011). The bioactivity of human bone morphogenetic protein-15 is sensitive to C-terminal modification: Characterization of the purified untagged processed mature region. Molecular and Cellular Endocrinology, 332, 106–115. Ratnakumar, A., Mousset, S., Glémin, S., Berglund, J., Galtier, N., Duret, L., & Webster, M. T. (2010). Detecting positive selection within genomes: The problem of biased gene conversion. Philosophical transactions of the Royal Society of London. Series B, Biological Sciences, 365, 2571–2580. Romiguier, J., Ranwez, V., Douzery, E. J., & Galtier, N. (2010). Contrasting GC-content dynamics across 33 mammalian genomes: Relationship with life-history traits and chromosome sizes. Genome Research, 20, 1001–1009. Sánchez, F., & Smitz, J. (2012). Molecular control of oogenesis. Biochimica et Biophysica Acta (BBA)-Molecular Basis of Disease, 1822, 1896–1912. Sengle, G., Ono, R. N., Sasaki, T., & Sakai, L. Y. (2011). Prodomains of transforming growth factor β (TGFβ) superfamily members specify different functions extracellular matrix interactions and growth factor bioavailability. Journal of Biological Chemistry, 286, 5087–5099. Tamura, K., Stecher, G., Peterson, D., Filipski, A., & Kumar, S. (2013). MEGA6: Molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution, 30, 2725–2729. Voight, B. F., Kudaravalli, S., Wen, X., & Pritchard, J. K. (2006). A map of recent positive selection in the human genome. PLoS Biology, 4, e72. Wigglesworth, K., Lee, K. B., O'Brien, M. J., Peng, J., Matzuk, M. M., & Eppig, J. J. (2013). Bidirectional communication between oocytes and ovarian follicular somatic cells is required for meiotic arrest of mammalian oocytes. Proceedings of the National Academy of Sciences, 110, E3723–E3729. Yang, Z. (2007). PAML 4: Phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution, 24, 1586–1591. Yang, J. R., Liao, B. Y., Zhuang, S. M., & Zhang, J. (2012). Protein misinteraction avoidance causes highly expressed proteins to evolve slowly. Proceedings of the National Academy of Sciences, 109, E831–E840.