1
|
Nisa FU, Naqvi RZ, Arshad F, Ilyas I, Asif M, Amin I, Mrode R, Mansoor S, Mukhtar Z. Assessment of Genomic Diversity and Selective Pressures in Crossbred Dairy Cattle of Pakistan. Biochem Genet 2024; 62:4137-4156. [PMID: 38664326 DOI: 10.1007/s10528-024-10809-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Accepted: 04/08/2024] [Indexed: 09/28/2024]
Abstract
Improving the low productivity levels of native cattle breeds in smallholder farming systems is a pressing concern in Pakistan. Crossbreeding high milk-yielding holstein friesian (HF) breed with the adaptability and heat tolerance of Sahiwal cattle has resulted in offspring that are well-suited to local conditions and exhibit improved milk yield. The exploration of how desirable traits in crossbred dairy cattle are selected has not yet been investigated. This study aims to provide the first overview of the selective pressures on the genome of crossbred dairy cattle in Pakistan. A total of eighty-one crossbred, thirty-two HF and twenty-four Sahiwal cattle were genotyped, and additional SNP genotype data for HF and Sahiwal were collected from a public database to equate the sample size in each group. Within-breed selection signatures in crossbreds were investigated using the integrated haplotype score. Crossbreds were also compared to each of their parental breeds to discover between-population signatures of selection using two approaches: cross-population extended haplotype homozygosity and fixation index. We identified several overlapping genes associated with production, immunity, and adaptation traits, including U6, TMEM41B, B4GALT7, 5S_rRNA, RBM27, POU4F3, NSD1, PRELID1, RGS14, SLC34A1, TMED9, B4GALT7, OR2AK3, OR2T16, OR2T60, OR2L3, and CTNNA1. Our results suggest that regions responsible for milk traits have generally experienced stronger selective pressure than others.
Collapse
Affiliation(s)
- Fakhar Un Nisa
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
- Department of Animal Breeding and Genetics, University of Veterinary and Animal Sciences, Lahore, 54000, Pakistan
| | - Rubab Zahra Naqvi
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
| | - Fazeela Arshad
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
| | - Iram Ilyas
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
| | - Muhammad Asif
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
| | - Imran Amin
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
| | - Raphael Mrode
- Animal Biosciences, International Livestock Research Institute, Nairobi, Kenya
- Animal and Veterinary Sciences, Scotland's Rural College, Edinburgh, UK
| | - Shahid Mansoor
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan
- International Centre for Chemical and Biological Sciences, University of Karachi, Karachi, Pakistan
| | - Zahid Mukhtar
- Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering College (NIBGE-C), Faisalabad, 38000, Pakistan.
- Pakistan Institute of Engineering & Applied Sciences (PIEAS), Nilore, Islamabad, 45650, Pakistan.
| |
Collapse
|
2
|
Fonseca PAS, Alonso-García M, Pelayo R, Marina H, Esteban-Blanco C, Mateo J, Gutiérrez-Gil B, Arranz JJ, Suárez-Vega A. Integrated analyses of the methylome and transcriptome to unravel sex differences in the perirenal fat from suckling lambs. Front Genet 2022; 13:1035063. [PMID: 36386829 PMCID: PMC9663842 DOI: 10.3389/fgene.2022.1035063] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 10/18/2022] [Indexed: 11/25/2022] Open
Abstract
In sheep, differences were observed regarding fat accumulation and fatty acid (FA) composition between males and females, which may impact the quality and organoleptic characteristics of the meat. The integration of different omics technologies is a relevant approach for investigating biological and genetic mechanisms associated with complex traits. Here, the perirenal tissue of six male and six female Assaf suckling lambs was evaluated using RNA sequencing and whole-genome bisulfite sequencing (WGBS). A multiomic discriminant analysis using multiblock (s)PLS-DA allowed the identification of 314 genes and 627 differentially methylated regions (within these genes), which perfectly discriminate between males and females. These candidate genes overlapped with previously reported QTLs for carcass fat volume and percentage of different FAs in milk and meat from sheep. Additionally, differentially coexpressed (DcoExp) modules of genes between males (nine) and females (three) were identified that harbour 22 of these selected genes. Interestingly, these DcoExp were significantly correlated with fat percentage in different deposits (renal, pelvic, subcutaneous and intramuscular) and were associated with relevant biological processes for adipogenesis, adipocyte differentiation, fat volume and FA composition. Consequently, these genes may potentially impact adiposity and meat quality traits in a sex-specific manner, such as juiciness, tenderness and flavour.
Collapse
Affiliation(s)
- Pablo A. S. Fonseca
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - María Alonso-García
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Rocio Pelayo
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Hector Marina
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Cristina Esteban-Blanco
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Javier Mateo
- Departamento de Higiene y Tecnología de Los Alimentos, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Beatriz Gutiérrez-Gil
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| | - Juan-José Arranz
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain,*Correspondence: Juan-José Arranz,
| | - Aroa Suárez-Vega
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, León, Spain
| |
Collapse
|
3
|
Pan C, Yang C, Wang S, Ma Y. Identifying Key Genes and Functionally Enriched Pathways of Diverse Adipose Tissue Types in Cattle. Front Genet 2022; 13:790690. [PMID: 35237299 PMCID: PMC8884536 DOI: 10.3389/fgene.2022.790690] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 01/26/2022] [Indexed: 12/26/2022] Open
Abstract
Background: Fat is a tissue that not just stores energy and plays a protective role; it is also a vital endocrine organ that generates and integrates signals to influence metabolism. Meanwhile, the excessive accumulation of lipids in adipose tissue can lead to metabolic disturbance and diseases. To date, the complicated molecular mechanisms of bovine adipose tissue are still unknown. This study aimed to identify key genes and functionally enriched pathways in various adipose tissue types. Results: The RNAseq data of 264 samples were downloaded from Gene Expression Omnibus (GEO) and analyzed by weighted gene co-expression network analysis (WGCNA). We identified 19 modules that significantly associated with at least one adipose tissue type. The brown module from GSE39618 was most closely associated with intramuscular fat tissue, which contained 550 genes. These genes were significantly enriched in pathways that related to inflammation and disease, such as TNF signaling pathway, IL-17 signaling pathway, and NF-kappa B signaling pathway. The pink module (GSE39618) that contained 58 genes was most closely associated with omental fat tissue. The turquoise (GSE39618), blue (GSE116775), and yellow (GSE65125) module were most closely associated with subcutaneous fat tissue. Genes in these modules were significantly enriched in pathways related to fat metabolism, such as the PPAR signaling pathway, fatty acid metabolism and PI3K-Akt signaling pathway. At last, key genes for intramuscular fat (PTGS2 and IL6), omental fat (ARHGEF5 and WT1), and subcutaneous fat (KIT, QR6Q1, PKD2L1, etc.) were obtained and verified. In addition, it was found that IL10 and VCAM1 might be potential genes to distinguish adipose and muscle. Conclusion: The study applied WGCNA to generate a landscape of adipose tissue and provide a basis for identifying potential pathways and hub genes of different adipose tissue types.
Collapse
Affiliation(s)
- Cuili Pan
- School of Agriculture, Ningxia University, Yinchuan, China
- Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia Hui Autonomous Region, Ningxia University, Yinchuan, China
| | - Chaoyun Yang
- School of Agriculture, Ningxia University, Yinchuan, China
- Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia Hui Autonomous Region, Ningxia University, Yinchuan, China
| | - Shuzhe Wang
- School of Agriculture, Ningxia University, Yinchuan, China
- Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia Hui Autonomous Region, Ningxia University, Yinchuan, China
| | - Yun Ma
- School of Agriculture, Ningxia University, Yinchuan, China
- Key Laboratory of Ruminant Molecular and Cellular Breeding, Ningxia Hui Autonomous Region, Ningxia University, Yinchuan, China
- *Correspondence: Yun Ma,
| |
Collapse
|
4
|
Multipopulational transcriptome analysis of post-weaned beef cattle at arrival further validates candidate biomarkers for predicting clinical bovine respiratory disease. Sci Rep 2021; 11:23877. [PMID: 34903778 PMCID: PMC8669006 DOI: 10.1038/s41598-021-03355-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Accepted: 11/24/2021] [Indexed: 12/13/2022] Open
Abstract
Bovine respiratory disease (BRD) remains the leading infectious disease in post-weaned beef cattle. The objective of this investigation was to contrast the at-arrival blood transcriptomes from cattle derived from two distinct populations that developed BRD in the 28 days following arrival versus cattle that did not. Forty-eight blood samples from two populations were selected for mRNA sequencing based on even distribution of development (n = 24) or lack of (n = 24) clinical BRD within 28 days following arrival; cattle which developed BRD were further stratified into BRD severity cohorts based on frequency of antimicrobial treatment: treated once (treated_1) or treated twice or more and/or died (treated_2+). Sequenced reads (~ 50 M/sample, 150 bp paired-end) were aligned to the ARS-UCD1.2 bovine genome assembly. One hundred and thirty-two unique differentially expressed genes (DEGs) were identified between groups stratified by disease severity (healthy, n = 24; treated_1, n = 13; treated_2+, n = 11) with edgeR (FDR ≤ 0.05). Differentially expressed genes in treated_1 relative to both healthy and treated_2+ were predicted to increase neutrophil activation, cellular cornification/keratinization, and antimicrobial peptide production. Differentially expressed genes in treated_2+ relative to both healthy and treated_1 were predicted to increase alternative complement activation, decrease leukocyte activity, and increase nitric oxide production. Receiver operating characteristic (ROC) curves generated from expression data for six DEGs identified in our current and previous studies (MARCO, CFB, MCF2L, ALOX15, LOC100335828 (aka CD200R1), and SLC18A2) demonstrated good-to-excellent (AUC: 0.800–0.899; ≥ 0.900) predictability for classifying disease occurrence and severity. This investigation identifies candidate biomarkers and functional mechanisms in at arrival blood that predicted development and severity of BRD.
Collapse
|
5
|
Doyle JL, Purfield DC, Moore T, Carthy TR, Walsh SW, Veerkamp RF, Evans RD, Berry DP. Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle. J Anim Sci 2021; 99:6157006. [PMID: 33677555 DOI: 10.1093/jas/skab070] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2020] [Accepted: 02/26/2021] [Indexed: 12/11/2022] Open
Abstract
Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.
Collapse
Affiliation(s)
- Jennifer L Doyle
- Teagasc, Animal and Grassland Research and Innovation Centre, Moorepark, Fermoy, Co. Cork, Ireland.,Department of Science, Waterford Institute of Technology, Cork Road, Co. Waterford, Ireland
| | - Deirdre C Purfield
- Department of Biological Sciences, Cork Institute of Technology, Bishopstown, Co. Cork, Ireland
| | - Tom Moore
- School of Biochemistry and Cell Biology, University College Cork, Western Gateway Building, Western Road, Cork, Ireland
| | - Tara R Carthy
- Teagasc, Animal and Grassland Research and Innovation Centre, Moorepark, Fermoy, Co. Cork, Ireland
| | - Siobhan W Walsh
- Department of Science, Waterford Institute of Technology, Cork Road, Co. Waterford, Ireland
| | - Roel F Veerkamp
- Animal Breeding and Genomics Centre, Wageningen University and Research Centre, Livestock Research, Wageningen, the Netherlands
| | - Ross D Evans
- Irish Cattle Breeding Federation, Bandon, Co. Cork, Ireland
| | - Donagh P Berry
- Teagasc, Animal and Grassland Research and Innovation Centre, Moorepark, Fermoy, Co. Cork, Ireland
| |
Collapse
|
6
|
Chen Q, Shen J, Hanif Q, Chen N, Huang Y, Dang R, Lan X, Chen H, Lei C. Whole genome analyses revealed genomic difference between European taurine and East Asian taurine. J Anim Breed Genet 2020; 138:56-68. [PMID: 32770713 DOI: 10.1111/jbg.12501] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2020] [Revised: 06/29/2020] [Accepted: 07/10/2020] [Indexed: 12/19/2022]
Abstract
European taurine and East Asian taurine are two main clades in Bos taurus, but their genomic differences are not clearly elucidated. Here, we sequenced 16 Mongolian cattle genomes and compared them to the 92 genomes of 10 representative breeds worldwide. We found the highest LD level in Mishima cattle and the fastest LD decay in European taurine. Phylogenetic analysis revealed that Mongolian, Hanwoo and Mishima cattle were clustered into East Asian taurine. From selective sweep, gene annotation, functional enrichment and differential expression analysis, we identified selective signals including genes and/or pathways related to rapid growth and large body size in European taurine, and superior meat quality in East Asian taurine. Our findings will help us understand the evolutionary history and formation process of the breeds and provide theoretical materials regarding the genetic mechanism underlying breed characteristics and molecular breeding programmes of the taurine clades in the future.
Collapse
Affiliation(s)
- Qiuming Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Jiafei Shen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Quratulain Hanif
- National Institute for Biotechnology and Genetic Engineering, Pakistan Institute of Engineering and Applied Sciences, Faisalabad, Pakistan.,Pakistan Institute of Engineering and Applied Sciences, Nilore, Islamabad, Pakistan
| | - Ningbo Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Yongzhen Huang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Ruihua Dang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Xianyong Lan
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Hong Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| |
Collapse
|
7
|
Wang Y, Jiang H, Yang L. Transcriptome Analysis of Zebrafish Olfactory Epithelium Reveal Sexual Differences in Odorant Detection. Genes (Basel) 2020; 11:genes11060592. [PMID: 32471067 PMCID: PMC7349279 DOI: 10.3390/genes11060592] [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: 03/03/2020] [Revised: 05/21/2020] [Accepted: 05/25/2020] [Indexed: 01/04/2023] Open
Abstract
Animals have evolved a large number of olfactory receptor genes in their genome to detect numerous odorants in their surrounding environments. However, we still know little about whether males and females possess the same abilities to sense odorants, especially in fish. In this study, we used deep RNA sequencing to examine the difference of transcriptome between male and female zebrafish olfactory epithelia. We found that the olfactory transcriptomes between males and females are highly similar. We also found evidence of some genes showing differential expression or alternative splicing, which may be associated with odorant-sensing between sexes. Most chemosensory receptor genes showed evidence of expression in the zebrafish olfactory epithelium, with a higher expression level in males than in females. Taken together, our results provide a comprehensive catalog of the genes mediating olfactory perception and pheromone-evoked behavior in fishes.
Collapse
Affiliation(s)
- Ying Wang
- Hubei Engineering Research Center for Protection and Utilization of Special Biological Resources in the Hanjiang River Basin, School of Life Sciences, Jianghan University, Wuhan 430056, Hubei, China;
| | - Haifeng Jiang
- The Key Laboratory of Aquatic Biodiversity and Conservation of Chinese Academy of Sciences, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, Hubei, China;
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Liandong Yang
- The Key Laboratory of Aquatic Biodiversity and Conservation of Chinese Academy of Sciences, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, Hubei, China;
- Correspondence: ; Tel.: +86-27-6878-0281
| |
Collapse
|
8
|
Kassam I, Wu Y, Yang J, Visscher PM, McRae AF. Tissue-specific sex differences in human gene expression. Hum Mol Genet 2020; 28:2976-2986. [PMID: 31044242 DOI: 10.1093/hmg/ddz090] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Revised: 04/12/2019] [Accepted: 04/24/2019] [Indexed: 02/07/2023] Open
Abstract
Despite extensive sex differences in human complex traits and disease, the male and female genomes differ only in the sex chromosomes. This implies that most sex-differentiated traits are the result of differences in the expression of genes that are common to both sexes. While sex differences in gene expression have been observed in a range of different tissues, the biological mechanisms for tissue-specific sex differences (TSSDs) in gene expression are not well understood. A total of 30 640 autosomal and 1021 X-linked transcripts were tested for heterogeneity in sex difference effect sizes in n = 617 individuals across 40 tissue types in Genotype-Tissue Expression (GTEx). This identified 65 autosomal and 66 X-linked TSSD transcripts (corresponding to unique genes) at a stringent significance threshold. Results for X-linked TSSD transcripts showed mainly concordant direction of sex differences across tissues and replicate previous findings. Autosomal TSSD transcripts had mainly discordant direction of sex differences across tissues. The top cis-expression quantitative trait loci (eQTLs) across tissues for autosomal TSSD transcripts are located a similar distance away from the nearest androgen and estrogen binding motifs and the nearest enhancer, as compared to cis-eQTLs for transcripts with stable sex differences in gene expression across tissue types. Enhancer regions that overlap top cis-eQTLs for TSSD transcripts, however, were found to be more dispersed across tissues. These observations suggest that androgen and estrogen regulatory elements in a cis region may play a common role in sex differences in gene expression, but TSSD in gene expression may additionally be due to causal variants located in tissue-specific enhancer regions.
Collapse
Affiliation(s)
| | - Yang Wu
- Institute for Molecular Bioscience
| | - Jian Yang
- Institute for Molecular Bioscience.,Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | - Peter M Visscher
- Institute for Molecular Bioscience.,Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | - Allan F McRae
- Institute for Molecular Bioscience.,Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| |
Collapse
|
9
|
Naqvi S, Godfrey AK, Hughes JF, Goodheart ML, Mitchell RN, Page DC. Conservation, acquisition, and functional impact of sex-biased gene expression in mammals. Science 2019; 365:eaaw7317. [PMID: 31320509 PMCID: PMC6896219 DOI: 10.1126/science.aaw7317] [Citation(s) in RCA: 121] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2019] [Accepted: 06/12/2019] [Indexed: 12/15/2022]
Abstract
Sex differences abound in human health and disease, as they do in other mammals used as models. The extent to which sex differences are conserved at the molecular level across species and tissues is unknown. We surveyed sex differences in gene expression in human, macaque, mouse, rat, and dog, across 12 tissues. In each tissue, we identified hundreds of genes with conserved sex-biased expression-findings that, combined with genomic analyses of human height, explain ~12% of the difference in height between females and males. We surmise that conserved sex biases in expression of genes otherwise operating equivalently in females and males contribute to sex differences in traits. However, most sex-biased expression arose during the mammalian radiation, which suggests that careful attention to interspecies divergence is needed when modeling human sex differences.
Collapse
Affiliation(s)
- Sahin Naqvi
- Whitehead Institute, Cambridge, MA 02142, USA
- Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Alexander K Godfrey
- Whitehead Institute, Cambridge, MA 02142, USA
- Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | | | - Mary L Goodheart
- Whitehead Institute, Cambridge, MA 02142, USA
- Howard Hughes Medical Institute, Whitehead Institute, Cambridge, MA 02142, USA
| | - Richard N Mitchell
- Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA
| | - David C Page
- Whitehead Institute, Cambridge, MA 02142, USA.
- Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Howard Hughes Medical Institute, Whitehead Institute, Cambridge, MA 02142, USA
| |
Collapse
|
10
|
Saitou M, Gokcumen O. Resolving the Insertion Sites of Polymorphic Duplications Reveals a HERC2 Haplotype under Selection. Genome Biol Evol 2019; 11:1679-1690. [PMID: 31124564 PMCID: PMC6587411 DOI: 10.1093/gbe/evz107] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/19/2019] [Indexed: 12/18/2022] Open
Abstract
Polymorphic duplications in humans have been shown to contribute to phenotypic diversity. However, the evolutionary forces that maintain variable duplications across the human genome are largely unexplored. We developed a linkage-disequilibrium based method to detect insertion sites of polymorphic duplications not represented in reference genomes. This method also allows resolution of haplotypes harboring the duplications. Using this approach, we conducted genome-wide analyses and identified the insertion sites of 22 common polymorphic duplications. We found that the majority of these duplications is intrachromosomal and only one of them is an interchromosomal insertion. Further characterization of these duplications revealed significant associations to blood and skin phenotypes. On the basis of population genetics analyses, we found that the duplication of a well-characterized pigmentation-related region, including the HERC2 gene, may be selected against in European populations. We further demonstrated that the haplotype harboring this duplication significantly affects the expression of the HERC2P9 gene in multiple tissues. Our study sheds light onto the evolutionary impact of understudied polymorphic duplications in human populations and presents methodological insights for future studies.
Collapse
Affiliation(s)
- Marie Saitou
- Department of Biological Sciences, SUNY at Buffalo
| | | |
Collapse
|
11
|
Jahuey-Martínez FJ, Parra-Bracamonte GM, Sifuentes-Rincón AM, Moreno-Medina VR. Signatures of selection in Charolais beef cattle identified by genome-wide analysis. J Anim Breed Genet 2019; 136:378-389. [PMID: 31020734 DOI: 10.1111/jbg.12399] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Revised: 03/23/2019] [Accepted: 03/25/2019] [Indexed: 11/29/2022]
Abstract
Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8-59.3 Mb on BTA 4; 67.5-69.3 on BTA 5; 39.5-41.0 Mb on BTA 6; and 26.4-29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.
Collapse
|
12
|
Pareek CS, Sachajko M, Jaskowski JM, Herudzinska M, Skowronski M, Domagalski K, Szczepanek J, Czarnik U, Sobiech P, Wysocka D, Pierzchala M, Polawska E, Stepanow K, Ogłuszka M, Juszczuk-Kubiak E, Feng Y, Kumar D. Comparative Analysis of the Liver Transcriptome among Cattle Breeds Using RNA-seq. Vet Sci 2019; 6:vetsci6020036. [PMID: 30934933 PMCID: PMC6631511 DOI: 10.3390/vetsci6020036] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Revised: 03/16/2019] [Accepted: 03/26/2019] [Indexed: 12/26/2022] Open
Abstract
Global gene expression in liver transcriptome varies among cattle breeds. The present investigation was aimed to identify the differentially expressed genes (DEGs), metabolic gene networks and metabolic pathways in bovine liver transcriptome of young bulls. In this study, we comparatively analyzed the bovine liver transcriptome of dairy (Polish Holstein Friesian (HF); n = 6), beef (Hereford; n = 6), and dual purpose (Polish-Red; n = 6) cattle breeds. This study identified 895, 338, and 571 significant (p < 0.01) differentially expressed (DE) gene-transcripts represented as 745, 265, and 498 hepatic DE genes through the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-HF versus Polish-Red breeds comparisons, respectively. By combining all breeds comparisons, 75 hepatic DE genes (p < 0.01) were identified as commonly shared among all the three breed comparisons; 70, 160, and 38 hepatic DE genes were commonly shared between the following comparisons: (i) Polish-Red versus Hereford and Polish-HF versus Hereford; (ii) Polish-Red versus Hereford and Polish-HF versus Polish-Red; and (iii) Polish-HF versus Hereford and Polish-HF versus Polish-Red, respectively. A total of 440, 82, and 225 hepatic DE genes were uniquely observed for the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-Red versus Polish-HF comparisons, respectively. Gene ontology (GO) analysis identified top-ranked enriched GO terms (p < 0.01) including 17, 16, and 31 functional groups and 151, 61, and 140 gene functions that were DE in all three breed liver transcriptome comparisons. Gene network analysis identified several potential metabolic pathways involved in glutamine family amino-acid, triglyceride synthesis, gluconeogenesis, p38MAPK cascade regulation, cholesterol biosynthesis (Polish-Red versus Hereford); IGF-receptor signaling, catecholamine transport, lipoprotein lipase, tyrosine kinase binding receptor (Polish-HF versus Hereford), and PGF-receptor binding, (Polish-HF versus Polish-Red). Validation results showed that the relative expression values were consistent to those obtained by RNA-seq, and significantly correlated between the quantitative reverse transcription PCR (RT-qPCR) and RNA-seq (Pearson’s r > 0.90). Our results provide new insights on bovine liver gene expressions among dairy versus dual versus beef breeds by identifying the large numbers of DEGs markers submitted to NCBI gene expression omnibus (GEO) accession number GSE114233, which can serve as useful genetic tools to develop the gene assays for trait-associated studies as well as, to effectively implement in genomics selection (GS) cattle breeding programs in Poland.
Collapse
Affiliation(s)
- Chandra Shekhar Pareek
- Centre for Modern Interdisciplinary Technologies, Nicolaus Copernicus University, 87-100 Torun, Poland.
- Centre of Veterinary Sciences, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Mateusz Sachajko
- Centre for Modern Interdisciplinary Technologies, Nicolaus Copernicus University, 87-100 Torun, Poland.
- Centre of Veterinary Sciences, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Jedrzej M Jaskowski
- Centre of Veterinary Sciences, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Magdalena Herudzinska
- Centre of Veterinary Sciences, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Mariusz Skowronski
- Centre of Veterinary Sciences, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Krzysztof Domagalski
- Centre for Modern Interdisciplinary Technologies, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Joanna Szczepanek
- Centre for Modern Interdisciplinary Technologies, Nicolaus Copernicus University, 87-100 Torun, Poland.
| | - Urszula Czarnik
- Faculty of Animal Bio-engineering, University of Warmia and Mazury in Olsztyn, 10-719 Olsztyn, Poland.
| | - Przymeslaw Sobiech
- Faculty of Veterinary Medicine, University of Warmia and Mazury in Olsztyn, 10719 Olsztyn, Poland.
| | - Dominika Wysocka
- Faculty of Veterinary Medicine, University of Warmia and Mazury in Olsztyn, 10719 Olsztyn, Poland.
| | - Mariusz Pierzchala
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, 05-552 Jastrzębiec, Poland.
| | - Ewa Polawska
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, 05-552 Jastrzębiec, Poland.
| | - Kamila Stepanow
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, 05-552 Jastrzębiec, Poland.
| | - Magdalena Ogłuszka
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, 05-552 Jastrzębiec, Poland.
| | - Edyta Juszczuk-Kubiak
- Faculty of Veterinary Medicine, Warsaw University of Life Sciences, 02-787 Warsaw, Poland.
| | - Yaping Feng
- Waksman Institute of Microbiology, Rutgers, The State University of New Jersey, Piscataway, NJ 08 854, USA.
| | - Dibyendu Kumar
- Waksman Institute of Microbiology, Rutgers, The State University of New Jersey, Piscataway, NJ 08 854, USA.
| |
Collapse
|
13
|
Piles M, Fernandez-Lozano C, Velasco-Galilea M, González-Rodríguez O, Sánchez JP, Torrallardona D, Ballester M, Quintanilla R. Machine learning applied to transcriptomic data to identify genes associated with feed efficiency in pigs. Genet Sel Evol 2019; 51:10. [PMID: 30866799 PMCID: PMC6417084 DOI: 10.1186/s12711-019-0453-y] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Accepted: 03/04/2019] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND To date, the molecular mechanisms that underlie residual feed intake (RFI) in pigs are unknown. Results from different genome-wide association studies and gene expression analyses are not always consistent. The aim of this research was to use machine learning to identify genes associated with feed efficiency (FE) using transcriptomic (RNA-Seq) data from pigs that are phenotypically extreme for RFI. METHODS RFI was computed by considering within-sex regression on mean metabolic body weight, average daily gain, and average backfat gain. RNA-Seq analyses were performed on liver and duodenum tissue from 32 high and 33 low RFI pigs collected at 153 d of age. Machine-learning algorithms were used to predict RFI class based on gene expression levels in liver and duodenum after adjusting for batch effects. Genes were ranked according to their contribution to the classification using the permutation accuracy importance score in an unbiased random forest (RF) algorithm based on conditional inference. Support vector machine, RF, elastic net (ENET) and nearest shrunken centroid algorithms were tested using different subsets of the top rank genes. Nested resampling for hyperparameter tuning was implemented with tenfold cross-validation in the outer and inner loops. RESULTS The best classification was obtained with ENET using the expression of 200 genes in liver [area under the receiver operating characteristic curve (AUROC): 0.85; accuracy: 0.78] and 100 genes in duodenum (AUROC: 0.76; accuracy: 0.69). Canonical pathways and candidate genes that were previously reported as associated with FE in several species were identified. The most remarkable pathways and genes identified were NRF2-mediated oxidative stress response and aldosterone signalling in epithelial cells, the DNAJC6, DNAJC1, MAPK8, PRKD3 genes in duodenum, and melatonin degradation II, PPARα/RXRα activation, and GPCR-mediated nutrient sensing in enteroendocrine cells and SMOX, IL4I1, PRKAR2B, CLOCK and CCK genes in liver. CONCLUSIONS ML algorithms and RNA-Seq expression data were found to provide good performance for classifying pigs into high or low RFI groups. Classification was better with gene expression data from liver than from duodenum. Genes associated with FE in liver and duodenum tissue that can be used as predictive biomarkers for this trait were identified.
Collapse
Affiliation(s)
- Miriam Piles
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| | - Carlos Fernandez-Lozano
- Computer Science Department, University of A Coruña, Campus Elviña s/n, 15071 A Coruña, Spain
| | - María Velasco-Galilea
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| | - Olga González-Rodríguez
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| | - Juan Pablo Sánchez
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| | - David Torrallardona
- Animal Nutrition Program, Institute of Agriculture and Food Research and Technology (IRTA), Mas de Bover, 43120 Constantí, Spain
| | - Maria Ballester
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| | - Raquel Quintanilla
- Animal Breeding and Genetics Program, Institute of Agriculture and Food Research and Technology (IRTA), Torre Marimon s/n, 08140 Caldes de Montbui, Barcelona, Spain
| |
Collapse
|
14
|
Yoon J, Kim H. Multi-tissue observation of the long non-coding RNA effects on sexually biased gene expression in cattle. ASIAN-AUSTRALASIAN JOURNAL OF ANIMAL SCIENCES 2019; 32:1044-1051. [PMID: 30744377 PMCID: PMC6603329 DOI: 10.5713/ajas.18.0516] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/06/2018] [Accepted: 11/14/2018] [Indexed: 11/27/2022]
Abstract
OBJECTIVE Recent studies have implied that gene expression has high tissue-specificity, and therefore it is essential to investigate gene expression in a variety of tissues when performing the transcriptomic analysis. In addition, the gradual increase of long non-coding RNA (lncRNA) annotation database has increased the importance and proportion of mapped reads accordingly. METHODS We employed simple statistical models to detect the sexually biased/dimorphic genes and their conjugate lncRNAs in 40 RNA-seq samples across two factors: sex and tissue. We employed two quantification pipeline: mRNA annotation only and mRNA+lncRNA annotation. RESULTS As a result, the tissue-specific sexually dimorphic genes are affected by the addition of lncRNA annotation at a non-negligible level. In addition, many lncRNAs are expressed in a more tissue-specific fashion and with greater variation between tissues compared to protein-coding genes. Due to the genic region lncRNAs, the differentially expressed gene list changes, which results in certain sexually biased genes to become ambiguous across the tissues. CONCLUSION In a past study, it has been reported that tissue-specific patterns can be seen throughout the differentially expressed genes between sexes in cattle. Using the same dataset, this study used a more recent reference, and the addition of conjugate lncRNA information, which revealed alterations of differentially expressed gene lists that result in an apparent distinction in the downstream analysis and interpretation. We firmly believe such misquantification of genic lncRNAs can be vital in both future and past studies.
Collapse
Affiliation(s)
- Joon Yoon
- Department of Natural Science, Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul 08826, Korea
| | - Heebal Kim
- Department of Natural Science, Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul 08826, Korea.,Department of Agricultural Biotechnology, Animal Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul 08826, Korea
| |
Collapse
|
15
|
Duan JE, Shi W, Jue NK, Jiang Z, Kuo L, O'Neill R, Wolf E, Dong H, Zheng X, Chen J, Tian XC. Dosage Compensation of the X Chromosomes in Bovine Germline, Early Embryos, and Somatic Tissues. Genome Biol Evol 2019; 11:242-252. [PMID: 30566637 PMCID: PMC6354180 DOI: 10.1093/gbe/evy270] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/12/2018] [Indexed: 12/15/2022] Open
Abstract
Dosage compensation of the mammalian X chromosome (X) was proposed by Susumu Ohno as a mechanism wherein the inactivation of one X in females would lead to doubling the expression of the other. This would resolve the dosage imbalance between eutherian females (XX) versus male (XY) and between a single active X versus autosome pairs (A). Expression ratio of X- and A-linked genes has been relatively well studied in humans and mice, despite controversial results over the existence of upregulation of X-linked genes. Here we report the first comprehensive test of Ohno’s hypothesis in bovine preattachment embryos, germline, and somatic tissues. Overall an incomplete dosage compensation (0.5 < X:A < 1) of expressed genes and an excess X dosage compensation (X:A > 1) of ubiquitously expressed “dosage-sensitive” genes were seen. No significant differences in X:A ratios were observed between bovine female and male somatic tissues, further supporting Ohno’s hypothesis. Interestingly, preimplantation embryos manifested a unique pattern of X dosage compensation dynamics. Specifically, X dosage decreased after fertilization, indicating that the sperm brings in an inactive X to the matured oocyte. Subsequently, the activation of the bovine embryonic genome enhanced expression of X-linked genes and increased the X dosage. As a result, an excess compensation was exhibited from the 8-cell stage to the compact morula stage. The X dosage peaked at the 16-cell stage and stabilized after the blastocyst stage. Together, our findings confirm Ohno’s hypothesis of X dosage compensation in the bovine and extend it by showing incomplete and over-compensation for expressed and “dosage-sensitive” genes, respectively.
Collapse
Affiliation(s)
| | - Wei Shi
- Department of Statistics, University of Connecticut, Storrs, CT
| | - Nathaniel K Jue
- School of Natural Sciences, California State University, Monterey Bay, CA
| | - Zongliang Jiang
- School of Animal Science, Louisiana State University, Agricultural Center, Baton Rouge, LA
| | - Lynn Kuo
- Department of Statistics, University of Connecticut, Storrs, CT
| | - Rachel O'Neill
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT
| | - Eckhard Wolf
- Gene Center and Department of Biochemistry, Ludwig-Maximilians-Universität Muünchen, Germany
| | - Hong Dong
- Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi, Xinjiang, P.R. China
| | - Xinbao Zheng
- Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi, Xinjiang, P.R. China
| | - Jingbo Chen
- Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi, Xinjiang, P.R. China
| | | |
Collapse
|
16
|
Sabino M, Carmelo VAO, Mazzoni G, Cappelli K, Capomaccio S, Ajmone-Marsan P, Verini-Supplizi A, Trabalza-Marinucci M, Kadarmideen HN. Gene co-expression networks in liver and muscle transcriptome reveal sex-specific gene expression in lambs fed with a mix of essential oils. BMC Genomics 2018; 19:236. [PMID: 29618337 PMCID: PMC5885410 DOI: 10.1186/s12864-018-4632-y] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2017] [Accepted: 03/27/2018] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND Essential oil (EO) dietary supplementation is a new strategy to improve animal health. EO compounds have antiparasitic, antimicrobial, antiviral, antimycotic, antioxidant and anti-inflammatory proprieties. Nutrigenomics investigations represent innovative approaches in understanding the relation between diet effect and gene expression related to the animal performance. Few nutrigenomics studies have used a high-throughput RNA-Sequencing (RNA-Seq) approach, despite great potential of RNA-Seq data in gene expression quantification and in co-expression network analyses. Our aim is to use the potential of RNA-Sequencing data in order to evaluate the effect of an EO supplementary diet on gene expression in both lamb liver and muscle. RESULTS Using a treatment and sex interaction model, 13 and 4 differentially expressed genes were identified in liver and muscle respectively. Sex-specific differentially expressed (DE) genes were identified in both sexes. Using network based analysis, different clusters of co-expressed genes that were highly correlated to the diet were detected in males vs. females, in agreement with DE analysis. A total of five regulatory genes in liver tissue associated to EO diet were identified: DNAJB9, MANF, UFM1, CTNNLA1 and NFX1. Our study reveals a sex-dependent effect of EO diet in both tissues, and an influence on the expression of genes mainly involved in immune, inflammatory and stress pathway. CONCLUSION Our analysis suggests a sex-dependent effect of the EO dietary supplementation on the expression profile of both liver and muscle tissues. We hypothesize that the presence of EOs could have beneficial effects on wellness of male lamb and further analyses are needed to understand the biological mechanisms behind the different effect of EO metabolites based on sex. Using lamb as a model for nutrigenomics studies, it could be interesting to investigate the effects of EO diets in other species and in humans.
Collapse
Affiliation(s)
- Marcella Sabino
- Dipartimento di Medicina Veterinaria, University of Perugia, Perugia, Italy
| | | | - Gianluca Mazzoni
- Department of Bio and Health Informatics, Technical University of Denmark, Kongens Lyngby, Copenhagen, Denmark
| | - Katia Cappelli
- Dipartimento di Medicina Veterinaria, University of Perugia, Perugia, Italy
| | - Stefano Capomaccio
- Dipartimento di Medicina Veterinaria, University of Perugia, Perugia, Italy
| | - Paolo Ajmone-Marsan
- Istituto di Zootecnica, Catholic University of the Sacred Heart, Piacenza, Italy
| | | | | | - Haja N Kadarmideen
- Department of Bio and Health Informatics, Technical University of Denmark, Kongens Lyngby, Copenhagen, Denmark.
| |
Collapse
|
17
|
Kassam I, Lloyd-Jones L, Holloway A, Small KS, Zeng B, Bakshi A, Metspalu A, Gibson G, Spector TD, Esko T, Montgomery GW, Powell JE, Yang J, Visscher PM, McRae AF. Autosomal genetic control of human gene expression does not differ across the sexes. Genome Biol 2016; 17:248. [PMID: 27908293 PMCID: PMC5134098 DOI: 10.1186/s13059-016-1111-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2016] [Accepted: 11/18/2016] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND Despite their nearly identical genomes, males and females differ in risk, incidence, prevalence, severity and age-at-onset of many diseases. Sexual dimorphism is also seen in human autosomal gene expression, and has largely been explored by examining the contribution of genotype-by-sex interactions to variation in gene expression. RESULTS In this study, we use data from a mixture of pedigree and unrelated individuals with verified European ancestry to investigate the sex-specific genetic architecture of gene expression measured in whole blood across n=1048 males and n=1005 females by treating gene expression intensities in the sexes as two distinct traits and estimating the genetic correlation (r G) between them. These correlations measure the similarity of the combined additive genetic effects of all single-nucleotide polymorphisms across the autosomal chromosomes, and thus the level of common genetic control of gene expression across the sexes. Genetic correlations are estimated across the sexes for the expression levels of 12,528 autosomal gene expression probes using bivariate GREML, and tested for differences in autosomal genetic control of gene expression across the sexes. Overall, no deviation of the distribution of test statistics is observed from that expected under the null hypothesis of a common autosomal genetic architecture for gene expression across the sexes. CONCLUSIONS These results suggest that males and females share the same common genetic control of gene expression.
Collapse
Affiliation(s)
- Irfahan Kassam
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia.
| | - Luke Lloyd-Jones
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | - Alexander Holloway
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | - Kerrin S Small
- Department of Twin Research and Genetic Epidemiology, King's College London, London, UK
| | - Biao Zeng
- School of Biology and Centre for Integrative Genomics, Georgia Institute of Technology, Atlanta, USA
| | - Andrew Bakshi
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | | | - Greg Gibson
- School of Biology and Centre for Integrative Genomics, Georgia Institute of Technology, Atlanta, USA
| | - Tim D Spector
- Department of Twin Research and Genetic Epidemiology, King's College London, London, UK
| | - Tonu Esko
- Estonian Genome Centre, University of Tartu, Tartu, Estonia
| | - Grant W Montgomery
- Institute for Molecular Bioscience, University of Queensland, Brisbane, Australia.,QIMR Berghofer Medical Research Institute, Brisbane, Australia
| | - Joseph E Powell
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia.,Institute for Molecular Bioscience, University of Queensland, Brisbane, Australia
| | - Jian Yang
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| | - Peter M Visscher
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia.,University of Queensland Diamantina Institute, Translational Research Institute, The University of Queensland, Brisbane, Australia
| | - Allan F McRae
- Queensland Brain Institute, The University of Queensland, Brisbane, Australia
| |
Collapse
|
18
|
Guebel DV, Torres NV. Sexual Dimorphism and Aging in the Human Hyppocampus: Identification, Validation, and Impact of Differentially Expressed Genes by Factorial Microarray and Network Analysis. Front Aging Neurosci 2016; 8:229. [PMID: 27761111 PMCID: PMC5050216 DOI: 10.3389/fnagi.2016.00229] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2016] [Accepted: 09/14/2016] [Indexed: 01/09/2023] Open
Abstract
Motivation: In the brain of elderly-healthy individuals, the effects of sexual dimorphism and those due to normal aging appear overlapped. Discrimination of these two dimensions would powerfully contribute to a better understanding of the etiology of some neurodegenerative diseases, such as “sporadic” Alzheimer. Methods: Following a system biology approach, top-down and bottom-up strategies were combined. First, public transcriptome data corresponding to the transition from adulthood to the aging stage in normal, human hippocampus were analyzed through an optimized microarray post-processing (Q-GDEMAR method) together with a proper experimental design (full factorial analysis). Second, the identified genes were placed in context by building compatible networks. The subsequent ontology analyses carried out on these networks clarify the main functionalities involved. Results: Noticeably we could identify large sets of genes according to three groups: those that exclusively depend on the sex, those that exclusively depend on the age, and those that depend on the particular combinations of sex and age (interaction). The genes identified were validated against three independent sources (a proteomic study of aging, a senescence database, and a mitochondrial genetic database). We arrived to several new inferences about the biological functions compromised during aging in two ways: by taking into account the sex-independent effects of aging, and considering the interaction between age and sex where pertinent. In particular, we discuss the impact of our findings on the functions of mitochondria, autophagy, mitophagia, and microRNAs. Conclusions: The evidence obtained herein supports the occurrence of significant neurobiological differences in the hippocampus, not only between adult and elderly individuals, but between old-healthy women and old-healthy men. Hence, to obtain realistic results in further analysis of the transition from the normal aging to incipient Alzheimer, the features derived from the sexual dimorphism in hippocampus should be explicitly considered.
Collapse
Affiliation(s)
- Daniel V Guebel
- Biotechnology Counselling ServicesBuenos Aires, Argentina; Systems Biology and Mathematical Modelling Group, Departamento de Bioquímica, Microbiología, Biología Celular y Genética, Facultad de Ciencias, Universidad de La LagunaSan Cristóbal de La Laguna, España
| | - Néstor V Torres
- Systems Biology and Mathematical Modelling Group, Departamento de Bioquímica, Microbiología, Biología Celular y Genética, Facultad de Ciencias, Universidad de La Laguna San Cristóbal de La Laguna, España
| |
Collapse
|
19
|
Ka S, Ahn H, Seo M, Kim H, Kim JN, Lee HJ. Status of dosage compensation of X chromosome in bovine genome. Genetica 2016; 144:435-44. [PMID: 27376899 DOI: 10.1007/s10709-016-9912-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Accepted: 06/27/2016] [Indexed: 01/18/2023]
Abstract
Dosage compensation system with X chromosome upregulation and inactivation have evolved to overcome the genetic imbalance between sex chromosomes in both male and female of mammals. Although recent development of chromosome-wide technologies has allowed us to test X upregulation, discrete data processing and analysis methods draw disparate conclusions. A series of expression studies revealed status of dosage compensation in some species belonging to monotremes, marsupials, rodents and primates. However, X upregulation in the Artiodactyla order including cattle have not been studied yet. In this study, we surveyed the genome-wide transcriptional upregulation in X chromosome in cattle RNA-seq data using different gene filtration methods. Overall examination of RNA-seq data revealed that X chromosome in the pituitary gland expressed more genes than in other peripheral tissues, which was consistent with the previous results observed in human and mouse. When analyzed with globally expressed genes, a median X:A expression ratio was 0.94. The ratio of 1-to-1 ortholog genes between chicken and mammals, however, showed considerable reduction to 0.68. These results indicate that status of dosage compensation for cattle is not deviated from those found in rodents and primate, and this is consistent with the evolutionary history of cattle.
Collapse
Affiliation(s)
- Sojeong Ka
- Department of Agricultural Biotechnology, College of Agriculture and Life Sciences, Seoul National University, Seoul, 08826, Republic of Korea
| | - Hyeonju Ahn
- Department of Agricultural Biotechnology, College of Agriculture and Life Sciences, Seoul National University, Seoul, 08826, Republic of Korea
| | - Minseok Seo
- Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, 08826, Republic of Korea
- C&K Genomics Inc., Seoul National University Research Park, Seoul, 08826, Republic of Korea
| | - Heebal Kim
- Department of Agricultural Biotechnology, College of Agriculture and Life Sciences, Seoul National University, Seoul, 08826, Republic of Korea
- Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, 08826, Republic of Korea
| | - Jin Nam Kim
- C&K Genomics Inc., Seoul National University Research Park, Seoul, 08826, Republic of Korea.
| | - Hyun-Jeong Lee
- Animal Nutritional Physiology Team, National Institute of Animal Science, Jeonju, Jeollabuk-Do, 55365, Republic of Korea.
| |
Collapse
|
20
|
RNA-seq analysis for detecting quantitative trait-associated genes. Sci Rep 2016; 6:24375. [PMID: 27071914 PMCID: PMC4829873 DOI: 10.1038/srep24375] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Accepted: 03/22/2016] [Indexed: 02/06/2023] Open
Abstract
Many recent RNA-seq studies were focused mainly on detecting the differentially expressed genes (DEGs) between two or more conditions. In contrast, only a few attempts have been made to detect genes associated with quantitative traits, such as obesity index and milk yield, on RNA-seq experiment with large number of biological replicates. This study illustrates the linear model application on trait associated genes (TAGs) detection in two real RNA-seq datasets: 89 replicated human obesity related data and 21 replicated Holsteins’ milk production related RNA-seq data. Based on these two datasets, the performance between suggesting methods, such as ordinary regression and robust regression, and existing methods: DESeq2 and Voom, were compared. The results indicate that suggesting methods have much lower false discoveries compared to the precedent two group comparisons based approaches in our simulation study and qRT-PCR experiment. In particular, the robust regression outperforms existing DEG finding method as well as ordinary regression in terms of precision. Given the current trend in RNA-seq pricing, we expect our methods to be successfully applied in various RNA-seq studies with numerous biological replicates that handle continuous response traits.
Collapse
|