1
|
Hoorn QA, Zayas GA, Rodriguez EE, Jensen LM, Mateescu RG, Hansen PJ. Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers. J Anim Sci Biotechnol 2023; 14:137. [PMID: 37932831 PMCID: PMC10629031 DOI: 10.1186/s40104-023-00940-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 09/10/2023] [Indexed: 11/08/2023] Open
Abstract
BACKGROUND In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production. RESULTS Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified. CONCLUSIONS In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.
Collapse
Affiliation(s)
- Quinn A Hoorn
- Department of Animal Sciences, Donald Henry Barron Reproductive and Perinatal Biology Research Program, and the Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Gabriel A Zayas
- Department of Animal Sciences, Donald Henry Barron Reproductive and Perinatal Biology Research Program, and the Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Eduardo E Rodriguez
- Department of Animal Sciences, Donald Henry Barron Reproductive and Perinatal Biology Research Program, and the Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Laura M Jensen
- Present address: School of Applied Systems Biology, La Trobe University, Bundoora, VIC, 3083, Australia
| | - Raluca G Mateescu
- Department of Animal Sciences, Donald Henry Barron Reproductive and Perinatal Biology Research Program, and the Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Peter J Hansen
- Department of Animal Sciences, Donald Henry Barron Reproductive and Perinatal Biology Research Program, and the Genetics Institute, University of Florida, Gainesville, FL, USA.
| |
Collapse
|
2
|
Marceau A, Wang J, Iqbal V, Jiang J, Liu GE, Ma L. Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods. Genes (Basel) 2023; 14:1789. [PMID: 37761929 PMCID: PMC10531232 DOI: 10.3390/genes14091789] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Revised: 09/01/2023] [Accepted: 09/08/2023] [Indexed: 09/29/2023] Open
Abstract
This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.
Collapse
Affiliation(s)
- Alexis Marceau
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA; (A.M.); (V.I.)
| | - Junjian Wang
- Department of Animal Science, North Carlonina State University, Raleigh, NC 27695, USA; (J.W.); (J.J.)
| | - Victoria Iqbal
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA; (A.M.); (V.I.)
| | - Jicai Jiang
- Department of Animal Science, North Carlonina State University, Raleigh, NC 27695, USA; (J.W.); (J.J.)
| | - George E. Liu
- Animal Genomics and Improvemennt Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA;
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA; (A.M.); (V.I.)
| |
Collapse
|
3
|
de Souza Fonseca PA, Suárez-Vega A, Cánovas A. Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs. Funct Integr Genomics 2022; 22:1361-1376. [PMID: 36001276 DOI: 10.1007/s10142-022-00893-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 07/16/2022] [Accepted: 08/09/2022] [Indexed: 11/25/2022]
Abstract
The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.
Collapse
Affiliation(s)
- Pablo Augusto de Souza Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Aroa Suárez-Vega
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
4
|
Nguyen LT, Lau LY, Fortes MRS. Proteomic Analysis of Hypothalamus and Pituitary Gland in Pre and Postpubertal Brahman Heifers. Front Genet 2022; 13:935433. [PMID: 35774501 PMCID: PMC9237413 DOI: 10.3389/fgene.2022.935433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 05/25/2022] [Indexed: 11/28/2022] Open
Abstract
The hypothalamus and the pituitary gland are directly involved in the complex systemic changes that drive the onset of puberty in cattle. Here, we applied integrated bioinformatics to elucidate the critical proteins underlying puberty and uncover potential molecular mechanisms from the hypothalamus and pituitary gland of prepubertal (n = 6) and postpubertal (n = 6) cattle. Proteomic analysis in the hypothalamus and pituitary gland revealed 275 and 186 differentially abundant (DA) proteins, respectively (adjusted p-value < 0.01). The proteome profiles found herein were integrated with previously acquired transcriptome profiles. These transcriptomic studies used the same tissues harvested from the same heifers at pre- and post-puberty. This comparison detected a small number of matched transcripts and protein changes at puberty in each tissue, suggesting the need for multiple omics analyses for interpreting complex biological systems. In the hypothalamus, upregulated DA proteins at post-puberty were enriched in pathways related to puberty, including GnRH, calcium and oxytocin signalling pathways, whereas downregulated proteins were observed in the estrogen signalling pathway, axon guidance and GABAergic synapse. Additionally, this study revealed that ribosomal pathway proteins in the pituitary were involved in the pubertal development of mammals. The reported molecules and derived protein-protein networks are a starting point for future experimental approaches that might dissect with more detail the role of each molecule to provide new insights into the mechanisms of puberty onset in cattle.
Collapse
Affiliation(s)
- Loan To Nguyen
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, Australia
- *Correspondence: Loan To Nguyen,
| | - Li Yieng Lau
- Agency of Science, Technology and Research, Singapore, Singapore
| | | |
Collapse
|
5
|
Mohamed AR, Naval-Sanchez M, Menzies M, Evans B, King H, Reverter A, Kijas JW. Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation. BMC Genomics 2022; 23:413. [PMID: 35650521 PMCID: PMC9158274 DOI: 10.1186/s12864-022-08514-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2021] [Accepted: 03/29/2022] [Indexed: 12/03/2022] Open
Abstract
Background Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation. Results Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation. Conclusion The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary – ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08514-8.
Collapse
|
6
|
Chen SY, Schenkel FS, Melo ALP, Oliveira HR, Pedrosa VB, Araujo AC, Melka MG, Brito LF. Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes. BMC Genomics 2022; 23:331. [PMID: 35484513 PMCID: PMC9052698 DOI: 10.1186/s12864-022-08555-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 04/12/2022] [Indexed: 02/06/2023] Open
Abstract
Background Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic. Results There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits. Conclusions These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08555-z.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 611130, Sichuan, China
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Ana L P Melo
- Department of Reproduction and Animal Evaluation, Rural Federal University of Rio de Janeiro, Seropédica, RJ, 23897-000, Brazil
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa, PR, 84030-900, Brazil
| | - Andre C Araujo
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA
| | - Melkaye G Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, River Falls, WI, 54022, USA
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA. .,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
7
|
Ross EM, Sanjana H, Nguyen LT, Cheng Y, Moore SS, Hayes BJ. Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle. Front Genet 2022; 13:784663. [PMID: 35401673 PMCID: PMC8990236 DOI: 10.3389/fgene.2022.784663] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 01/24/2022] [Indexed: 11/13/2022] Open
Abstract
Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.
Collapse
Affiliation(s)
- Elizabeth M. Ross
- Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, Australia
- *Correspondence: Elizabeth M. Ross,
| | - Hari Sanjana
- Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, Australia
| | - Loan T. Nguyen
- Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, Australia
| | - YuanYuan Cheng
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| | - Stephen S. Moore
- Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, Australia
| | - Ben J. Hayes
- Centre for Animal Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, Australia
| |
Collapse
|
8
|
Fu W, Wang R, Nanaei HA, Wang J, Hu D, Jiang Y. RGD v2.0: a major update of the ruminant functional and evolutionary genomics database. Nucleic Acids Res 2021; 50:D1091-D1099. [PMID: 34643708 PMCID: PMC8728256 DOI: 10.1093/nar/gkab887] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2021] [Revised: 09/10/2021] [Accepted: 10/05/2021] [Indexed: 12/26/2022] Open
Abstract
Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.
Collapse
Affiliation(s)
- Weiwei Fu
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| | - Rui Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| | - Hojjat Asadollahpour Nanaei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| | - Jinxin Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| | - Dexiang Hu
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| | - Yu Jiang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China.,Center for Ruminant Genetics and Evolution, Northwest A&F University, Yangling 712100, China
| |
Collapse
|
9
|
Alexandre PA, Naval-Sánchez M, Menzies M, Nguyen LT, Porto-Neto LR, Fortes MRS, Reverter A. Chromatin accessibility and regulatory vocabulary across indicine cattle tissues. Genome Biol 2021; 22:273. [PMID: 34548076 PMCID: PMC8454054 DOI: 10.1186/s13059-021-02489-7] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Accepted: 09/08/2021] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND Spatiotemporal changes in the chromatin accessibility landscape are essential to cell differentiation, development, health, and disease. The quest of identifying regulatory elements in open chromatin regions across different tissues and developmental stages is led by large international collaborative efforts mostly focusing on model organisms, such as ENCODE. Recently, the Functional Annotation of Animal Genomes (FAANG) has been established to unravel the regulatory elements in non-model organisms, including cattle. Now, we can transition from prediction to validation by experimentally identifying the regulatory elements in tropical indicine cattle. The identification of regulatory elements, their annotation and comparison with the taurine counterpart, holds high promise to link regulatory regions to adaptability traits and improve animal productivity and welfare. RESULTS We generate open chromatin profiles for liver, muscle, and hypothalamus of indicine cattle through ATAC-seq. Using robust methods for motif discovery, motif enrichment and transcription factor binding sites, we identify potential master regulators of the epigenomic profile in these three tissues, namely HNF4, MEF2, and SOX factors, respectively. Integration with transcriptomic data allows us to confirm some of their target genes. Finally, by comparing our results with Bos taurus data we identify potential indicine-specific open chromatin regions and overlaps with indicine selective sweeps. CONCLUSIONS Our findings provide insights into the identification and analysis of regulatory elements in non-model organisms, the evolution of regulatory elements within two cattle subspecies as well as having an immediate impact on the animal genetics community in particular for a relevant productive species such as tropical cattle.
Collapse
Affiliation(s)
- Pâmela A Alexandre
- CSIRO Agriculture & Food, 306 Carmody Rd., QLD, 4067, Brisbane, Australia.
| | - Marina Naval-Sánchez
- CSIRO Agriculture & Food, 306 Carmody Rd., QLD, 4067, Brisbane, Australia
- Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Moira Menzies
- CSIRO Agriculture & Food, 306 Carmody Rd., QLD, 4067, Brisbane, Australia
| | - Loan T Nguyen
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| | | | - Marina R S Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Antonio Reverter
- CSIRO Agriculture & Food, 306 Carmody Rd., QLD, 4067, Brisbane, Australia
| |
Collapse
|
10
|
Sweett H, Fonseca PAS, Suárez-Vega A, Livernois A, Miglior F, Cánovas A. Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle. Sci Rep 2020; 10:20102. [PMID: 33208801 PMCID: PMC7676258 DOI: 10.1038/s41598-020-75758-3] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Accepted: 10/16/2020] [Indexed: 12/20/2022] Open
Abstract
Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.
Collapse
Affiliation(s)
- H Sweett
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - P A S Fonseca
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Suárez-Vega
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Livernois
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Department of Pathobiology, Ontario Veterinary College, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - F Miglior
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Cánovas
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
11
|
Lau LY, Nguyen LT, Reverter A, Moore SS, Lynn A, McBride‐Kelly L, Phillips‐Rose L, Plath M, Macfarlane R, Vasudivan V, Morton L, Ardley R, Ye Y, Fortes MRS. Gene regulation could be attributed to TCF3 and other key transcription factors in the muscle of pubertal heifers. Vet Med Sci 2020; 6:695-710. [PMID: 32432381 PMCID: PMC7738712 DOI: 10.1002/vms3.278] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2019] [Revised: 03/13/2020] [Accepted: 04/09/2020] [Indexed: 01/17/2023] Open
Abstract
Puberty is a whole-body event, driven by the hypothalamic integration of peripheral signals such as leptin or IGF-1. In the process of puberty, reproductive development is simultaneous to growth, including muscle growth. To enhance our understanding of muscle function related to puberty, we performed transcriptome analyses of muscle samples from six pre- and six post-pubertal Brahman heifers (Bos indicus). Our aims were to perform differential expression analyses and co-expression analyses to derive a regulatory gene network associate with puberty. As a result, we identified 431 differentially expressed (DEx) transcripts (genes and non-coding RNAs) when comparing pre- to post-pubertal average gene expression. The DEx transcripts were compared with all expressed transcripts in our samples (over 14,000 transcripts) for functional enrichment analyses. The DEx transcripts were associated with "extracellular region," "inflammatory response" and "hormone activity" (adjusted p < .05). Inflammatory response for muscle regeneration is a necessary aspect of muscle growth, which is accelerated during puberty. The term "hormone activity" may signal genes that respond to progesterone signalling in the muscle, as the presence of this hormone is an important difference between pre- and post-pubertal heifers in our experimental design. The DEx transcript with the highest average expression difference was a mitochondrial gene, ENSBTAG00000043574 that might be another important link between energy metabolism and puberty. In the derived co-expression gene network, we identified six hub genes: CDC5L, MYC, TCF3, RUNX2, ATF2 and CREB1. In the same network, 48 key regulators of DEx transcripts were identified, using a regulatory impact factor metric. The hub gene TCF3 was also a key regulator. The majority of the key regulators (22 genes) are members of the zinc finger family, which has been implicated in bovine puberty in other tissues. In conclusion, we described how puberty may affect muscle gene expression in cattle.
Collapse
Affiliation(s)
- Li Yieng Lau
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Loan T. Nguyen
- Queensland Alliance for Agriculture and Food InnovationThe University of QueenslandBrisbaneQLDAustralia
| | - Antonio Reverter
- CSIRO Agriculture and FoodQueensland Biosciences PrecinctBrisbaneQLDAustralia
| | - Stephen S. Moore
- Queensland Alliance for Agriculture and Food InnovationThe University of QueenslandBrisbaneQLDAustralia
| | - Aaron Lynn
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Liam McBride‐Kelly
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Louis Phillips‐Rose
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Mackenzie Plath
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Rhys Macfarlane
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Vanisha Vasudivan
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Lachlan Morton
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Ryan Ardley
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Yunan Ye
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
| | - Marina R. S. Fortes
- School of Chemistry and Molecular BiologyThe University of QueenslandBrisbaneQLDAustralia
- Queensland Alliance for Agriculture and Food InnovationThe University of QueenslandBrisbaneQLDAustralia
| |
Collapse
|
12
|
Alexandre PA, Hudson NJ, Lehnert SA, Fortes MRS, Naval-Sánchez M, Nguyen LT, Porto-Neto LR, Reverter A. Genome-Wide Co-Expression Distributions as a Metric to Prioritize Genes of Functional Importance. Genes (Basel) 2020; 11:E1231. [PMID: 33092259 PMCID: PMC7593939 DOI: 10.3390/genes11101231] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Accepted: 10/15/2020] [Indexed: 12/17/2022] Open
Abstract
Genome-wide gene expression analysis are routinely used to gain a systems-level understanding of complex processes, including network connectivity. Network connectivity tends to be built on a small subset of extremely high co-expression signals that are deemed significant, but this overlooks the vast majority of pairwise signals. Here, we developed a computational pipeline to assign to every gene its pair-wise genome-wide co-expression distribution to one of 8 template distributions shapes varying between unimodal, bimodal, skewed, or symmetrical, representing different proportions of positive and negative correlations. We then used a hypergeometric test to determine if specific genes (regulators versus non-regulators) and properties (differentially expressed or not) are associated with a particular distribution shape. We applied our methodology to five publicly available RNA sequencing (RNA-seq) datasets from four organisms in different physiological conditions and tissues. Our results suggest that genes can be assigned consistently to pre-defined distribution shapes, regarding the enrichment of differential expression and regulatory genes, in situations involving contrasting phenotypes, time-series, or physiological baseline data. There is indeed a striking additional biological signal present in the genome-wide distribution of co-expression values which would be overlooked by currently adopted approaches. Our method can be applied to extract further information from transcriptomic data and help uncover the molecular mechanisms involved in the regulation of complex biological process and phenotypes.
Collapse
Affiliation(s)
- Pâmela A. Alexandre
- CSIRO Agriculture & Food, St Lucia, QLD 4067, Australia; (S.A.L.); (L.R.P.-N.); (A.R.)
| | - Nicholas J. Hudson
- School of Agriculture and Food Sciences, The University of Queensland, Brisbane, QLD 4072, Australia;
| | - Sigrid A. Lehnert
- CSIRO Agriculture & Food, St Lucia, QLD 4067, Australia; (S.A.L.); (L.R.P.-N.); (A.R.)
| | - Marina R. S. Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD 4072, Australia;
| | - Marina Naval-Sánchez
- Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia;
| | - Loan T. Nguyen
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD 4072, Australia;
| | - Laercio R. Porto-Neto
- CSIRO Agriculture & Food, St Lucia, QLD 4067, Australia; (S.A.L.); (L.R.P.-N.); (A.R.)
| | - Antonio Reverter
- CSIRO Agriculture & Food, St Lucia, QLD 4067, Australia; (S.A.L.); (L.R.P.-N.); (A.R.)
| |
Collapse
|
13
|
Fonseca PAS, Suárez-Vega A, Cánovas A. Weighted Gene Correlation Network Meta-Analysis Reveals Functional Candidate Genes Associated with High- and Sub-Fertile Reproductive Performance in Beef Cattle. Genes (Basel) 2020; 11:genes11050543. [PMID: 32408659 PMCID: PMC7290847 DOI: 10.3390/genes11050543] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 05/04/2020] [Accepted: 05/06/2020] [Indexed: 12/13/2022] Open
Abstract
Improved reproductive efficiency could lead to economic benefits for the beef industry, once the intensive selection pressure has led to a decreased fertility. However, several factors limit our understanding of fertility traits, including genetic differences between populations and statistical limitations. In the present study, the RNA-sequencing data from uterine samples of high-fertile (HF) and sub-fertile (SF) animals was integrated using co-expression network meta-analysis, weighted gene correlation network analysis, identification of upstream regulators, variant calling, and network topology approaches. Using this pipeline, top hub-genes harboring fixed variants (HF × SF) were identified in differentially co-expressed gene modules (DcoExp). The functional prioritization analysis identified the genes with highest potential to be key-regulators of the DcoExp modules between HF and SF animals. Consequently, 32 functional candidate genes (10 upstream regulators and 22 top hub-genes of DcoExp modules) were identified. These genes were associated with the regulation of relevant biological processes for fertility, such as embryonic development, germ cell proliferation, and ovarian hormone regulation. Additionally, 100 candidate variants (single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)) were identified within those genes. In the long-term, the results obtained here may help to reduce the frequency of subfertility in beef herds, reducing the associated economic losses caused by this condition.
Collapse
Affiliation(s)
- Pablo A. S. Fonseca
- Correspondence: (P.A.S.F.); (A.C.); Tel.: +1-519-824-4120 (ext. 56295) (A.C.)
| | | | - Angela Cánovas
- Correspondence: (P.A.S.F.); (A.C.); Tel.: +1-519-824-4120 (ext. 56295) (A.C.)
| |
Collapse
|
14
|
Nolte W, Weikard R, Brunner RM, Albrecht E, Hammon HM, Reverter A, Kühn C. Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls. Int J Mol Sci 2020; 21:ijms21093292. [PMID: 32384694 PMCID: PMC7247587 DOI: 10.3390/ijms21093292] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2020] [Revised: 05/02/2020] [Accepted: 05/04/2020] [Indexed: 12/18/2022] Open
Abstract
Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid β-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.
Collapse
Affiliation(s)
- Wietje Nolte
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany; (W.N.); (R.W.); (R.M.B.)
| | - Rosemarie Weikard
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany; (W.N.); (R.W.); (R.M.B.)
| | - Ronald M. Brunner
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany; (W.N.); (R.W.); (R.M.B.)
| | - Elke Albrecht
- Institute of Muscle Biology and Growth, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany;
| | - Harald M. Hammon
- Institute of Nutritional Physiology “Oskar Kellner”, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany;
| | - Antonio Reverter
- Commonwealth Scientific and Industrial Research Organisation (CSIRO) Agriculture and Food, Queensland Bioscience Precinct, St Lucia 4067 QLD, Australia;
| | - Christa Kühn
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), 18196 Dummerstorf, Germany; (W.N.); (R.W.); (R.M.B.)
- Faculty of Agricultural and Environmental Sciences, University Rostock, 18059 Rostock, Germany
- Correspondence:
| |
Collapse
|
15
|
Ehsani R, Drabløs F. Enhanced identification of significant regulators of gene expression. BMC Bioinformatics 2020; 21:134. [PMID: 32252623 PMCID: PMC7132893 DOI: 10.1186/s12859-020-3468-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Accepted: 03/24/2020] [Indexed: 12/29/2022] Open
Abstract
Background Diseases like cancer will lead to changes in gene expression, and it is relevant to identify key regulatory genes that can be linked directly to these changes. This can be done by computing a Regulatory Impact Factor (RIF) score for relevant regulators. However, this computation is based on estimating correlated patterns of gene expression, often Pearson correlation, and an assumption about a set of specific regulators, normally transcription factors. This study explores alternative measures of correlation, using the Fisher and Sobolev metrics, and an extended set of regulators, including epigenetic regulators and long non-coding RNAs (lncRNAs). Data on prostate cancer have been used to explore the effect of these modifications. Results A tool for computation of RIF scores with alternative correlation measures and extended sets of regulators was developed and tested on gene expression data for prostate cancer. The study showed that the Fisher and Sobolev metrics lead to improved identification of well-documented regulators of gene expression in prostate cancer, and the sets of identified key regulators showed improved overlap with previously defined gene sets of relevance to cancer. The extended set of regulators lead to identification of several interesting candidates for further studies, including lncRNAs. Several key processes were identified as important, including spindle assembly and the epithelial-mesenchymal transition (EMT). Conclusions The study has shown that using alternative metrics of correlation can improve the performance of tools based on correlation of gene expression in genomic data. The Fisher and Sobolev metrics should be considered also in other correlation-based applications.
Collapse
Affiliation(s)
- Rezvan Ehsani
- Department of Mathematics, University of Zabol, Zabol, Iran. .,Department of Bioinformatics, University of Zabol, Zabol, Iran.
| | - Finn Drabløs
- Department of Cancer Research and Molecular Medicine, NTNU - Norwegian University of Science and Technology, NO-7491, Trondheim, Norway.
| |
Collapse
|
16
|
Ghoreishifar SM, Moradi-Shahrbabak H, Fallahi MH, Jalil Sarghale A, Moradi-Shahrbabak M, Abdollahi-Arpanahi R, Khansefid M. Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis. BMC Genet 2020; 21:16. [PMID: 32041535 PMCID: PMC7011551 DOI: 10.1186/s12863-020-0824-y] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Accepted: 02/04/2020] [Indexed: 01/06/2023] Open
Abstract
Background Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (FROH), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI), and diagonal elements of the genomic relationship matrix (FGRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies. Results In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on FHOM, FUNI and FGRM were higher in AZ than they were in KHZ, which was in contrast to the FROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX). Conclusion The calculated FROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that FROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Hossein Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran.
| | - Mohammad Hossein Fallahi
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Ali Jalil Sarghale
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Mohammad Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Rostam Abdollahi-Arpanahi
- Departments of Animal and Poultry Science, College of Aburaihan, University of Tehran, Pakdasht, 33916-53755, Iran
| | - Majid Khansefid
- AgriBio Centre for AgriBioscience, Agriculture Victoria, Bundoora, VIC, 3083, Australia
| |
Collapse
|
17
|
Nolte W, Weikard R, Brunner RM, Albrecht E, Hammon HM, Reverter A, Kühn C. Biological Network Approach for the Identification of Regulatory Long Non-Coding RNAs Associated With Metabolic Efficiency in Cattle. Front Genet 2019; 10:1130. [PMID: 31824560 PMCID: PMC6883949 DOI: 10.3389/fgene.2019.01130] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2019] [Accepted: 10/17/2019] [Indexed: 12/17/2022] Open
Abstract
Background: Genomic regions associated with divergent livestock feed efficiency have been found predominantly outside protein coding sequences. Long non-coding RNAs (lncRNA) can modulate chromatin accessibility, gene expression and act as important metabolic regulators in mammals. By integrating phenotypic, transcriptomic, and metabolomic data with quantitative trait locus data in prioritizing co-expression network analyses, we aimed to identify and functionally characterize lncRNAs with a potential key regulatory role in metabolic efficiency in cattle. Materials and Methods: Crossbred animals (n = 48) of a Charolais x Holstein F2-population were allocated to groups of high or low metabolic efficiency based on residual feed intake in bulls, energy corrected milk in cows and intramuscular fat content in both genders. Tissue samples from jejunum, liver, skeletal muscle and rumen were subjected to global transcriptomic analysis via stranded total RNA sequencing (RNAseq) and blood plasma samples were used for profiling of 640 metabolites. To identify lncRNAs within the indicated tissues, a project-specific transcriptome annotation was established. Subsequently, novel transcripts were categorized for potential lncRNA status, yielding a total of 7,646 predicted lncRNA transcripts belonging to 3,287 loci. A regulatory impact factor approach highlighted 92, 55, 35, and 73 lncRNAs in jejunum, liver, muscle, and rumen, respectively. Their ensuing high regulatory impact factor scores indicated a potential regulatory key function in a gene set comprising loci displaying differential expression, tissue specificity and loci overlapping with quantitative trait locus regions for residual feed intake or milk production. These were subjected to a partial correlation and information theory analysis with the prioritized gene set. Results and Conclusions: Independent, significant and group-specific correlations (|r| > 0.8) were used to build a network for the high and the low metabolic efficiency group resulting in 1,522 and 1,732 nodes, respectively. Eight lncRNAs displayed a particularly high connectivity (>100 nodes). Metabolites and genes from the partial correlation and information theory networks, which each correlated significantly with the respective lncRNA, were included in an enrichment analysis indicating distinct affected pathways for the eight lncRNAs. LncRNAs associated with metabolic efficiency were classified to be functionally involved in hepatic amino acid metabolism and protein synthesis and in calcium signaling and neuronal nitric oxide synthase signaling in skeletal muscle cells.
Collapse
Affiliation(s)
- Wietje Nolte
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Rosemarie Weikard
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Ronald M Brunner
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Elke Albrecht
- Institute of Muscle Biology and Growth, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Harald M Hammon
- Institute of Nutritional Physiology "Oskar Kellner," Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Antonio Reverter
- Commonwealth Scientific and Industrial Research Organisation (CSIRO) Agriculture and Food, Queensland Bioscience Precinct, St Lucia, QLD, Australia
| | - Christa Kühn
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany.,Faculty of Agricultural and Environmental Sciences, University Rostock, Rostock, Germany
| |
Collapse
|
18
|
Li M, Lu X, Xia H, Zhang C, Wang X, Chen Z, Zhang H, Qu K, Huang B, Moore S, Mao Y, Yang Z. In-depth characterization of the pituitary transcriptome in Simmental and Chinese native cattle. Domest Anim Endocrinol 2019; 66:35-42. [PMID: 30391830 DOI: 10.1016/j.domaniend.2018.09.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Revised: 09/10/2018] [Accepted: 09/21/2018] [Indexed: 12/15/2022]
Abstract
The pituitary gland is a key endocrine organ responsible for growth and development. To get closer to understanding important molecular mechanisms at work in the bovine pituitary and identifying the core genes controlling growth, in the present study we have analyzed the transcriptome profiling of the pituitary glands of two cattle breeds (Wenshan and Simmental cattle) with extreme growth differences using high-throughput RNA sequencing. Our study revealed that the overall gene expression levels are quite similar between the two breeds. A total of 18,065 annotated genes were identified, which accounted for 85% of the annotated genes in cattle. The five most abundant hormone genes (GH, PRL, POMC, CGA, and LHB) were relatively stable in both breeds, indicating their pivotal roles in maintaining the basic functions of the pituitary. In addition, 105 genes were identified as differentially expressed between Wenshan and Simmental pituitary glands, including 83 known and 22 novel genes. Functional association analyses of the differentially expressed genes between the breeds revealed 60 enriched gene ontology terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Interestingly, the most enriched KEGG pathway, cell adhesion molecules, can modulate pituitary hormone secretion by cell-cell contact. Our findings demonstrated that SYTL2, SLC38A4, and NCAM2 are new candidates for crucial functions in the secretory pathways of the pituitary gland. These results will further understanding of the important molecular mechanisms at work in the bovine pituitary.
Collapse
Affiliation(s)
- Mingxun Li
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Xubin Lu
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Hailei Xia
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Chenglong Zhang
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Xiaolong Wang
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Zhi Chen
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Huimin Zhang
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Kaixing Qu
- Yunnan Academy of Grassland and Animal Science, Kunming, Yunnan 650212, China
| | - Bizhi Huang
- Yunnan Academy of Grassland and Animal Science, Kunming, Yunnan 650212, China
| | - Stephen Moore
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, Australia
| | - Yongjiang Mao
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China
| | - Zhangping Yang
- Key Laboratory of Animal Genetics, Breeding and Molecular Design of Jiangsu province, Yangzhou University, Yangzhou 225002, China; Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225002, China.
| |
Collapse
|
19
|
Fortes MRS, Zacchi LF, Nguyen LT, Raidan F, Weller MMDCA, Choo JJY, Reverter A, Rego JPA, Boe-Hansen GB, Porto-Neto LR, Lehnert SA, Cánovas A, Schulz BL, Islas-Trejo A, Medrano JF, Thomas MG, Moore SS. Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty. Anim Genet 2018; 49:539-549. [PMID: 30192028 DOI: 10.1111/age.12721] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/26/2018] [Indexed: 12/17/2022]
Abstract
Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.
Collapse
Affiliation(s)
- M R S Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia
| | - L F Zacchi
- School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia
| | - L T Nguyen
- School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia.,Faculty of Biotechnology, Vietnam National University of Agriculture, Gialam, Hanoi, Vietnam
| | - F Raidan
- Animal Science Department, Universidade Federal de Viçosa, Vicosa, Minas Gerais, 36570-900, Brazil
| | - M M D C A Weller
- CSIRO Agriculture and Food, Queensland Bioscience Precinct, Brisbane, QLD 4072, Australia
| | - J J Y Choo
- School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia
| | - A Reverter
- Animal Science Department, Universidade Federal de Viçosa, Vicosa, Minas Gerais, 36570-900, Brazil
| | - J P A Rego
- Instituto Federal de Educação, Ciência e Tecnologia do Ceara, Fortaleza, Ceará, 62930-000, Brazil
| | - G B Boe-Hansen
- School of Veterinary Sciences, The University of Queensland, Gatton, QLD 4343, Australia
| | - L R Porto-Neto
- Animal Science Department, Universidade Federal de Viçosa, Vicosa, Minas Gerais, 36570-900, Brazil
| | - S A Lehnert
- Animal Science Department, Universidade Federal de Viçosa, Vicosa, Minas Gerais, 36570-900, Brazil
| | - A Cánovas
- Department of Animal Biosciences, Centre of Genetic Improvement for Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - B L Schulz
- School of Chemistry and Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia
| | - A Islas-Trejo
- Department of Animal Science, University of California Davis, Davis, CA, 95616, USA
| | - J F Medrano
- Department of Animal Science, University of California Davis, Davis, CA, 95616, USA
| | - M G Thomas
- Department of Animal Science, Colorado State University, Fort Collins, CO, 80523, USA
| | - S S Moore
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD 4072, Australia
| |
Collapse
|