1
|
Chen SY, Luo Z, Jia X, Zhou J, Lai SJ. Evaluating genomic inbreeding of two Chinese yak (Bos grunniens) populations. BMC Genomics 2024; 25:712. [PMID: 39044139 PMCID: PMC11267844 DOI: 10.1186/s12864-024-10640-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2024] [Accepted: 07/19/2024] [Indexed: 07/25/2024] Open
Abstract
BACKGROUND Yaks are a vital livestock in the Qinghai-Tibetan Plateau area for providing food products, maintaining sustainable ecosystems, and promoting cultural heritage. Because of uncontrolled mating, it is impossible to estimate inbreeding level of yak populations using the pedigree-based approaches. With the aims to accurately evaluate inbreeding level of two Chinese yak populations (Maiwa and Jiulong), we obtained genome-wide single nucleotide polymorphisms (SNPs) by DNA sequencing and calculated five SNP-by-SNP estimators ([Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]), as well as two segment-based estimators of runs of homozygosity (ROH, [Formula: see text]) and homozygous-by-descent (HBD, [Formula: see text]). Functional implications were analyzed for the positional candidate genes located within the related genomic regions. RESULTS A total of 151,675 and 190,955 high-quality SNPs were obtained from 71 Maiwa and 30 Jiulong yaks, respectively. Jiulong had greater genetic diversity than Maiwa in terms of allele frequency and nucleotide diversity. The two populations could be genetically distinguished by principal component analysis, with the mean differentiation index (Fst) of 0.0054. The greater genomic inbreeding levels of Maiwa yaks were consistently supported by all five SNP-by-SNP estimators. Based on simple proportion of homozygous SNPs ([Formula: see text]), a lower inbreeding level was indicated by three successfully sequenced old leather samples that may represent historical Maiwa yaks about five generations ago. There were 3304 ROH detected among all samples, with mean and median length of 1.97 Mb and 1.0 Mb, respectively. A total of 94 HBD segments were found among all samples, whereas 92 of them belonged to the shortest class with the mean length of 10.9 Kb. Based on the estimates of [Formula: see text] and [Formula: see text], however, there was no difference in inbreeding level between Maiwa and Jiulong yaks. Within the genomic regions with the significant Fst or enriched by ROH, we found several candidate genes and pathways that have been reported to be related to diverse production traits in farm animals. CONCLUSIONS We successfully evaluated the genomic inbreeding level of two Chinese yak populations. Although different estimators resulted in inconsistent conclusions on their genomic inbreeding levels, our results may be helpful to implement the genetic conservation and utilization programs for the two yak populations.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 211# Huimin Road, Wenjiang, Sichuan, 611130, China.
| | - Zhihao Luo
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 211# Huimin Road, Wenjiang, Sichuan, 611130, China
| | - Xianbo Jia
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 211# Huimin Road, Wenjiang, Sichuan, 611130, China
| | - Junkun Zhou
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 211# Huimin Road, Wenjiang, Sichuan, 611130, China
| | - Song-Jia Lai
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 211# Huimin Road, Wenjiang, Sichuan, 611130, China.
| |
Collapse
|
2
|
Benjelloun B, Leempoel K, Boyer F, Stucki S, Streeter I, Orozco-terWengel P, Alberto FJ, Servin B, Biscarini F, Alberti A, Engelen S, Stella A, Colli L, Coissac E, Bruford MW, Ajmone-Marsan P, Negrini R, Clarke L, Flicek P, Chikhi A, Joost S, Taberlet P, Pompanon F. Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints. Mol Ecol 2023:e17257. [PMID: 38149334 DOI: 10.1111/mec.17257] [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: 09/15/2022] [Revised: 08/18/2023] [Accepted: 12/05/2023] [Indexed: 12/28/2023]
Abstract
The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.
Collapse
Affiliation(s)
- Badr Benjelloun
- Livestock Genomics Laboratory, Regional Center of Agricultural Research Tadla, National Institute of Agricultural Research INRA, Rabat, Morocco
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| | - Kevin Leempoel
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Frédéric Boyer
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| | - Sylvie Stucki
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Ian Streeter
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, Cambridge, UK
| | - Pablo Orozco-terWengel
- School of Biosciences, Cardiff University, Wales, UK
- Sustainable Places Research Institute, Cardiff University, Cardiff, UK
| | - Florian J Alberto
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| | - Bertrand Servin
- GenPhySE, Université de Toulouse, INRAE, INPT, ENVT, Castanet-Tolosan, France
| | - Filippo Biscarini
- Institute of Agricultural Biology and Biotechnology, Consiglio Nazionale delle Ricerche (CNR), Milan, Italy
| | - Adriana Alberti
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ. Evry, Université Paris-Saclay, Evry, France
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), Gif-sur-Yvette, France
| | - Stefan Engelen
- Genoscope, Institut de biologie François-Jacob, Commissariat à l'Energie Atomique CEA, Université Paris-Saclay, Evry, France
| | - Alessandra Stella
- Institute of Agricultural Biology and Biotechnology, Consiglio Nazionale delle Ricerche (CNR), Milan, Italy
| | - Licia Colli
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Facoltà di Scienze Agrarie, Alimentari e Ambientali, Università Cattolica del S. Cuore, Piacenza, Italy
- BioDNA - Centro di Ricerca sulla Biodiversità e sul DNA Antico, Facoltà di Scienze Agrarie, Alimentari e Ambientali, Università Cattolica del S. Cuore, Piacenza, Italy
| | - Eric Coissac
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| | - Michael W Bruford
- School of Biosciences, Cardiff University, Wales, UK
- Sustainable Places Research Institute, Cardiff University, Cardiff, UK
| | - Paolo Ajmone-Marsan
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Facoltà di Scienze Agrarie, Alimentari e Ambientali, Università Cattolica del S. Cuore, Piacenza, Italy
- BioDNA - Centro di Ricerca sulla Biodiversità e sul DNA Antico, Facoltà di Scienze Agrarie, Alimentari e Ambientali, Università Cattolica del S. Cuore, Piacenza, Italy
| | - Riccardo Negrini
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Facoltà di Scienze Agrarie, Alimentari e Ambientali, Università Cattolica del S. Cuore, Piacenza, Italy
- AIA Associazione Italiana Allevatori, Roma, Italy
| | - Laura Clarke
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, Cambridge, UK
| | - Paul Flicek
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, Cambridge, UK
| | - Abdelkader Chikhi
- Livestock Genomics Laboratory, Regional Center of Agricultural Research Tadla, National Institute of Agricultural Research INRA, Rabat, Morocco
| | - Stéphane Joost
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Pierre Taberlet
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| | - François Pompanon
- Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, LECA, Grenoble, France
| |
Collapse
|
3
|
Integrated Analysis Reveals a lncRNA-miRNA-mRNA Network Associated with Pigeon Skeletal Muscle Development. Genes (Basel) 2021; 12:genes12111787. [PMID: 34828393 PMCID: PMC8625974 DOI: 10.3390/genes12111787] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Revised: 11/06/2021] [Accepted: 11/10/2021] [Indexed: 12/13/2022] Open
Abstract
Growing evidence has demonstrated the emerging role of long non-coding RNA as competitive endogenous RNA (ceRNA) in regulating skeletal muscle development. However, the mechanism of ceRNA regulated by lncRNA in pigeon skeletal muscle development remains unclear. To reveal the function and regulatory mechanisms of lncRNA, we first analyzed the expression profiles of lncRNA, microRNA (miRNA), and mRNA during the development of pigeon skeletal muscle using high-throughput sequencing. We then constructed a lncRNA-miRNA-mRNA ceRNA network based on differentially expressed (DE) lncRNAs, miRNAs, and mRNAs according to the ceRNA hypothesis. Functional enrichment and short time-series expression miner (STEM) analysis were performed to explore the function of the ceRNA network. Hub lncRNA-miRNA-mRNA interactions were identified by connectivity degree and validated using dual-luciferase activity assay. The results showed that a total of 1625 DE lncRNAs, 11,311 DE mRNAs, and 573 DE miRNAs were identified. A ceRNA network containing 9120 lncRNA-miRNA-mRNA interactions was constructed. STEM analysis indicated that the function of the lncRNA-associated ceRNA network might be developmental specific. Functional enrichment analysis identified potential pathways regulating pigeon skeletal muscle development, such as cell cycle and MAPK signaling. Based on the connectivity degree, lncRNAs TCONS_00066712, TCONS_00026594, TCONS_00001557, TCONS_00001553, and TCONS_00003307 were identified as hub genes in the ceRNA network. lncRNA TCONS_00026594 might regulate the FSHD region gene 1 (FRG1)/ SRC proto-oncogene, non-receptor tyrosine kinase (SRC) by sponge adsorption of cli-miR-1a-3p to affect the development of pigeon skeletal muscle. Our findings provide a data basis for in-depth elucidation of the lncRNA-associated ceRNA mechanism underlying pigeon skeletal muscle development.
Collapse
|