Thesis Title: Using transcriptomic and proteomic analysis to characterize the biological basis of puberty in Brahman heifers
Loan To Nguyen
Bachelor of Biotechnology
Master of Biotechnology
A thesis submitted for the degree of Doctor of Philosophy at The University of Queensland in 2018 School of Chemistry and Molecular Biosciences
Abstract: The beef industry in Northern Australia is facing increased demands for efficient herd productivity. Brahman cattle, a breed of the Bos indicus sub-species, can withstand hostile conditions in northern production systems (Burrow 2012). However, there is a large variability in the age at onset of puberty in Bos indicus and Bos taurus crossbred heifers (Warnick 1963; Wythe 1970; Reynolds 1973). Brahman heifers are older and heavier at puberty in comparison with Hereford, Angus, Holstein and Jersey heifers, which are all Bos taurus breeds (Baker et al. 1989). Late onset of puberty in Bos indicus can lead to a decrease of lifetime productivity (Lesmeister et al. 1973; Johnston et al. 2013). Therefore, reducing the age at puberty to increase cattle productivity is a major aim for Bos indicus breeders.
In heifers, puberty is initiated when the hypothalamus-pituitary-gonad axis loses its sensitivity to negative feedback effects of steroid hormones, allowing an increase in gonadotropin-releasing hormone (GnRH) secretion from GnRH neurons in the hypothalamus. The increase then stimulates the secretion of gonadotropins luteinizing hormones (LH) and follicle stimulating hormone (FSH). FSH and LH then regulate gonadal development (Day et al., 1984; Day et al., 1987; Schillo et al., 1992; Day, 1998; Gasser et al., 2006). However, the process is controlled by multiple factors and the complex interactions between environmental and genetic factors regulating the process is only now coming to light. Therefore, the aim of this PhD project was to provide a transcriptomic and proteomic analyses of important tissues related to heifer puberty. In this work, I targeted the hypothalamus, the pituitary gland, the liver and the adipose tissue from Brahman heifers that were either pre- or post-pubertal.
The transcriptomic analysis using an RNA sequencing technique provided reliable sequence reads. An average of 60 million paired-end reads was produced from each sample, of which 60 – 70% were properly mapped to the Bos taurus reference genome assembly UMD3.1. Sequence data of these transcriptome data sets are available through the Functional Annotation of Animal Genomes project (http://data.faang.org/home). The study identified 392, 284, 452 and 525 differentially expressed (DE) genes among pre- and post-pubertal heifers in the hypothalamus, pituitary gland, liver and adipose tissue, respectively (P-value < 0.05). A regulatory impact factor metric identified key regulators (transcription factors, TF) for each studied tissue. A list of TF that were also DE genes associated with Brahman heifer puberty was revealed: BDNF, STAT6, PBX2, PBRM1, E2F8, NFAT5, SIX5, ZBTB38, ZNF605, HES7, HEY2, RORC, PPARG, PITX1, NFE2, HOXB5, and zinc finger family genes. The DE and key TF genes were mainly assigned to three pathways: Wnt, TGF-beta and PPAR signalling (P-
I
value < 0.05, corrected). These results support our knowledge about the complex nature of puberty, a trait influenced by numerous genes organized in functional pathways.
The project also investigated the protein abundance differences between pre- and post- pubertal heifers using the same tissues. The hypothalamus and pituitary gland appeared as tissues having the highest number of differentially abundant (DA) proteins after puberty (adjusted P-value < 0.01). These valuable proteome datasets are publicly available via The PRoteomics IDEntifications (PRIDE) database (https://www.ebi.ac.uk/pride/archive/login). The functional analysis of DA proteins showed numerous enriched pathways (FDR < 0.05) that might regulate the occurrence of puberty, such as GnRH signalling, estrogen signalling pathway, steroid hormone biosynthesis, axon guidance, glutamatergic synapse, GABAergic synapse, oxytocin signalling pathway, calcium signalling pathway, PI3K-AKT signalling pathway, PPAR signalling pathway, ribosome and specific amino acid metabolisms. The candidate proteins revealed in this study will be useful for understanding the molecular mechanisms of puberty in Bos indicus breeds.
When compared, proteins and their corresponding transcripts expression levels in the four studied tissues the Person correlations were weak: 0.05 in the Hypothalamus (P-value = 0.16), 0.006 in pituitary gland (P-value = 0.8), 0.13 in liver (P-value = 0.001) and -0.06 in the adipose tissue (P-value = 0.2). The non-significant global correlation between mRNA levels and protein abundance presented in this Thesis confirmed that the transcriptome data is not sufficient to predict or speculate regarding the corresponding proteome data when samples represent tissues. It is possible that in single cell analyses, the transcriptome would match the proteome better, but this data was not available. In other words, changes in mRNA expression might provide a limited understanding of changes in the expression of proteins in complex systems.
In summary, these are the first studies that performed transcriptomic and proteomic analyses of reproductive tissues (hypothalamus and pituitary) and metabolic-related tissues (liver and adipose tissue) in pre- and post-pubertal Brahman heifers. The results presented in this Thesis provide an important step forward in our efforts to understand the molecular biology underlying puberty. This knowledge can bring the research one step closer towards the goal of enhancing genetic improvement predictions for puberty in Bos indicus cattle.
II
Declaration by author
This thesis is composed of my original work, and contains no material previously published or written by another person except where due reference has been made in the text. I have clearly stated the contribution by others to jointly-authored works that I have included in my thesis.
I have clearly stated the contribution of others to my thesis as a whole, including statistical assistance, survey design, data analysis, significant technical procedures, professional editorial advice, and any other original research work used or reported in my thesis. The content of my thesis is the result of work I have carried out since the commencement of my research higher degree candidature and does not include a substantial part of work that has been submitted to qualify for the award of any other degree or diploma in any university or other tertiary institution. I have clearly stated which parts of my thesis, if any, have been submitted to qualify for another award.
I acknowledge that an electronic copy of my thesis must be lodged with the University Library and, subject to the policy and procedures of The University of Queensland, the thesis be made available for research and study in accordance with the Copyright Act 1968 unless a period of embargo has been approved by the Dean of the Graduate School.
I acknowledge that copyright of all material contained in my thesis resides with the copyright holder(s) of that material. Where appropriate I have obtained copyright permission from the copyright holder to reproduce material in this thesis.
III
Publications during candidature
Peer-reviewed papers M. R. S. Fortes, L. T. Nguyen, L.R. Porto Neto, A. Reverter, S. S. Moore, S. A. Lehnert, M. G. Thomas, (2016) Polymorphism and genes associated with puberty in heifers. Theriogenology 86: 333-339 M. R. S. Fortes, L. T. Nguyen, Weller, M.D.C.A., Cánovas. A., Islas-Trejo, A., Porto- Neto, L.R., Reverter, A., Lehnert, S.A., Boe-Hansen, G.B., Thomas, M.G., Medrano, J.F. and Moore, S.S. (2016) Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post-versus pre-pubertal Brahman heifers. Journal of Animal Science, 94(9):3693-3702 L.T. Nguyen, Reverter. A., Cánovas. A., Venus. B., Islas-Trejo. A., Porto-Neto. L.R., Lehnert. S.A., Medrano. J.F., Moore. S.S. and M.R. S. Fortes, (2017) Global differential gene expression in the pituitary gland and the ovaries of pre- and postpubertal Brahman heifers. Journal of Animal Science, 95:599–615 M.M. Dias, A. Cánovas, C. Mantilla-Rojas, D.G. Riley, P. Luna-Nevarez, S.J. Coleman, S.E. Speidel, R.M. Enns, A. Islas-Trejo, J.F. Medrano,S.S. Moore, M.R.S. Fortes, L.T. Nguyen, B. Venus, I.S.D.P. Diaz,F.R.P. Souza, L.F.S. Fonseca, F. Baldi, L.G. Albuquerque, M.G. Thomas and H.N. Oliveira, (2017) SNP detection using RNA- sequences of candidate genes associated with puberty in cattle. Genetics and molecular research: GMR 16(1)