1
|
Yuan H, Wei W, Zhang Y, Li C, Zhao S, Chao Z, Xia C, Quan J, Gao C. Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing. Int J Mol Sci 2024; 25:5843. [PMID: 38892031 PMCID: PMC11172908 DOI: 10.3390/ijms25115843] [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: 04/10/2024] [Revised: 05/22/2024] [Accepted: 05/25/2024] [Indexed: 06/21/2024] Open
Abstract
Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.
Collapse
Affiliation(s)
- Haonan Yuan
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Wenjing Wei
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Yue Zhang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Changwen Li
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| | - Shengguo Zhao
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Zhe Chao
- Institute of Animal Science and Veterinary Medicine, Hainan Academy of Agricultural Sciences, Key Laboratory of Tropical Animal Breeding and Disease Research, Haikou 571100, China;
| | - Changyou Xia
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| | - Jinqiang Quan
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Caixia Gao
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| |
Collapse
|
2
|
Smith T, Olagunju T, Rosen B, Neibergs H, Becker G, Davenport K, Elsik C, Hadfield T, Koren S, Kuhn K, Rhie A, Shira K, Skibiel A, Stegemiller M, Thorne J, Villamediana P, Cockett N, Murdoch B. The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content. RESEARCH SQUARE 2024:rs.3.rs-4033388. [PMID: 38712074 PMCID: PMC11071540 DOI: 10.21203/rs.3.rs-4033388/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.
Collapse
Affiliation(s)
- Timothy Smith
- USDA, ARS, U.S. Meat Animal Research Center (USMARC)
| | | | | | | | | | | | | | | | - Sergey Koren
- Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health
| | | | - Arang Rhie
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | | | | | | | | | | | | | | |
Collapse
|
3
|
Greshnova A, Pál K, Martinez JFI, Canzar S, Makova KD. Transcript Isoform Diversity of Y Chromosome Ampliconic Genes of Great Apes Uncovered Using Long Reads and Telomere-to-Telomere Reference Genome Assemblies. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.02.587783. [PMID: 38617276 PMCID: PMC11014635 DOI: 10.1101/2024.04.02.587783] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
Y chromosomes of great apes harbor Ampliconic Genes (YAGs)-multi-copy gene families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) that encode proteins important for spermatogenesis. Previous work assembled YAG transcripts based on their targeted sequencing but not using reference genome assemblies, potentially resulting in an incomplete transcript repertoire. Here we used the recently produced gapless telomere-to-telomere (T2T) Y chromosome assemblies of great ape species (bonobo, chimpanzee, human, gorilla, Bornean orangutan, and Sumatran orangutan) and analyzed RNA data from whole-testis samples for the same species. We generated hybrid transcriptome assemblies by combining targeted long reads (Pacific Biosciences), untargeted long reads (Pacific Biosciences) and untargeted short reads (Illumina)and mapping them to the T2T reference genomes. Compared to the results from the reference-free approach, average transcript length was more than two times higher, and the total number of transcripts decreased three times, improving the quality of the assembled transcriptome. The reference-based transcriptome assemblies allowed us to differentiate transcripts originating from different Y chromosome gene copies and from their non-Y chromosome homologs. We identified two sources of transcriptome diversity-alternative splicing and gene duplication with subsequent diversification of gene copies. For each gene family, we detected transcribed pseudogenes along with protein-coding gene copies. We revealed previously unannotated gene copies of YAGs as compared to currently available NCBI annotations, as well as novel isoforms for annotated gene copies. This analysis paves the way for better understanding Y chromosome gene functions, which is important given their role in spermatogenesis.
Collapse
Affiliation(s)
- Aleksandra Greshnova
- Department of Biology, Penn State University, University Park, PA, USA
- Current address: Max Planck Institute for Evolutionary Biology, Plön, Germany
| | - Karol Pál
- Department of Biology, Penn State University, University Park, PA, USA
| | - Juan Francisco Iturralde Martinez
- Department of Computer Science and Engineering, The Pennsylvania State University, University Park, PA 16802, United States
- Huck Institutes of the Life Sciences. Pennsylvania State University, University Park, PA 16802, USA
| | - Stefan Canzar
- Faculty of Informatics and Data Science, University of Regensburg, Regensburg, Germany
- Department of Computer Science and Engineering, The Pennsylvania State University, University Park, PA 16802, United States
| | - Kateryna D Makova
- Department of Biology, Penn State University, University Park, PA, USA
| |
Collapse
|
4
|
Tomaszkiewicz M, Sahlin K, Medvedev P, Makova KD. Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes. Genome Biol Evol 2023; 15:evad205. [PMID: 37967251 PMCID: PMC10673640 DOI: 10.1093/gbe/evad205] [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: 04/14/2023] [Revised: 10/20/2023] [Accepted: 11/03/2023] [Indexed: 11/17/2023] Open
Abstract
Y chromosomal ampliconic genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been studied in great apes; however, the diversity of splicing variants remains unexplored. Here, we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this data set resulted in several findings. First, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Second, our results suggest that BPY2 transcripts and proteins originate from separate genomic regions in bonobo versus human, which is possibly facilitated by acquiring new promoters. Third, our analysis indicates that the PRY gene family, having the highest representation of noncoding transcripts, has been undergoing pseudogenization. Fourth, we have not detected signatures of selection in the five YAG families shared among great apes, even though we identified many species-specific protein-coding transcripts. Fifth, we predicted consensus disorder regions across most gene families and species, which could be used for future investigations of male infertility. Overall, our work illuminates the YAG isoform landscape and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.
Collapse
Affiliation(s)
- Marta Tomaszkiewicz
- Department of Biology, The Pennsylvania State University, University Park, PA 16802, USA
| | - Kristoffer Sahlin
- Department of Mathematics, Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Paul Medvedev
- Department of Computer Science and Engineering, The Pennsylvania State University, University Park, PA 16802, USA
- Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Medical Genomics, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Computational Biology and Bioinformatics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Kateryna D Makova
- Department of Biology, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Medical Genomics, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Computational Biology and Bioinformatics, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
5
|
Tian D, Sun D, Ren Q, Zhang P, Zhang Z, Zhang W, Luo H, Li X, Han B, Liu D, Zhao K. Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep. Sci Rep 2023; 13:17283. [PMID: 37828092 PMCID: PMC10570297 DOI: 10.1038/s41598-023-44402-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/28/2023] [Accepted: 10/07/2023] [Indexed: 10/14/2023] Open
Abstract
Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.
Collapse
Affiliation(s)
- Dehong Tian
- Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, 810001, Qinghai, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - De Sun
- Animal Husbandry and Veterinary Station of Huzhu County of Qinghai Province, Huzhu, 810500, Qinghai, China
| | - Qianben Ren
- Qinghai Sheep Breeding and Promotion Service Center, Gangcha, 812300, Qinghai, China
| | - Pei Zhang
- Qinghai Animal and Plant Quarantine Station, Xining, 810000, Qinghai, China
| | - Zian Zhang
- Qinghai Sheep Breeding and Promotion Service Center, Gangcha, 812300, Qinghai, China
| | - Wenkui Zhang
- Qinghai Sheep Breeding and Promotion Service Center, Gangcha, 812300, Qinghai, China
| | - Haizhou Luo
- Qinghai Sheep Breeding and Promotion Service Center, Gangcha, 812300, Qinghai, China
| | - Xue Li
- Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, 810001, Qinghai, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Buying Han
- Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, 810001, Qinghai, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Dehui Liu
- Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, 810001, Qinghai, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Kai Zhao
- Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, 810001, Qinghai, China.
- University of Chinese Academy of Sciences, Beijing, 100049, China.
| |
Collapse
|
6
|
Rhie A, Nurk S, Cechova M, Hoyt SJ, Taylor DJ, Altemose N, Hook PW, Koren S, Rautiainen M, Alexandrov IA, Allen J, Asri M, Bzikadze AV, Chen NC, Chin CS, Diekhans M, Flicek P, Formenti G, Fungtammasan A, Garcia Giron C, Garrison E, Gershman A, Gerton JL, Grady PGS, Guarracino A, Haggerty L, Halabian R, Hansen NF, Harris R, Hartley GA, Harvey WT, Haukness M, Heinz J, Hourlier T, Hubley RM, Hunt SE, Hwang S, Jain M, Kesharwani RK, Lewis AP, Li H, Logsdon GA, Lucas JK, Makalowski W, Markovic C, Martin FJ, Mc Cartney AM, McCoy RC, McDaniel J, McNulty BM, Medvedev P, Mikheenko A, Munson KM, Murphy TD, Olsen HE, Olson ND, Paulin LF, Porubsky D, Potapova T, Ryabov F, Salzberg SL, Sauria MEG, Sedlazeck FJ, Shafin K, Shepelev VA, Shumate A, Storer JM, Surapaneni L, Taravella Oill AM, Thibaud-Nissen F, Timp W, Tomaszkiewicz M, Vollger MR, Walenz BP, Watwood AC, Weissensteiner MH, Wenger AM, Wilson MA, Zarate S, Zhu Y, Zook JM, Eichler EE, O'Neill RJ, Schatz MC, Miga KH, Makova KD, Phillippy AM. The complete sequence of a human Y chromosome. Nature 2023; 621:344-354. [PMID: 37612512 PMCID: PMC10752217 DOI: 10.1038/s41586-023-06457-y] [Citation(s) in RCA: 71] [Impact Index Per Article: 71.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 07/19/2023] [Indexed: 08/25/2023]
Abstract
The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.
Collapse
Affiliation(s)
- Arang Rhie
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Sergey Nurk
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
- Oxford Nanopore Technologies Inc., Oxford, UK
| | - Monika Cechova
- Faculty of Informatics, Masaryk University, Brno, Czech Republic
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Savannah J Hoyt
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
| | - Dylan J Taylor
- Department of Biology, Johns Hopkins University, Baltimore, MD, USA
| | - Nicolas Altemose
- Department of Molecular and Cell Biology, University of California, Berkeley, CA, USA
| | - Paul W Hook
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Sergey Koren
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Mikko Rautiainen
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Ivan A Alexandrov
- Federal Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
- Center for Algorithmic Biotechnology, Saint Petersburg State University, St Petersburg, Russia
- Department of Anatomy and Anthropology and Department of Human Molecular Genetics and Biochemistry, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv-Yafo, Israel
| | - Jamie Allen
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Mobin Asri
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Andrey V Bzikadze
- Graduate Program in Bioinformatics and Systems Biology, University of California, San Diego, CA, USA
| | - Nae-Chyun Chen
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
| | - Chen-Shan Chin
- GeneDX Holdings Corp, Stamford, CT, USA
- Foundation of Biological Data Science, Belmont, CA, USA
| | - Mark Diekhans
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Paul Flicek
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
- Department of Genetics, University of Cambridge, Cambridge, UK
| | | | | | - Carlos Garcia Giron
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Erik Garrison
- Department of Genetics, Genomics and Informatics, University of Tennessee Health Science Center, Memphis, TN, USA
| | - Ariel Gershman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Jennifer L Gerton
- Stowers Institute for Medical Research, Kansas City, MO, USA
- University of Kansas Medical Center, Kansas City, MO, USA
| | - Patrick G S Grady
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
| | - Andrea Guarracino
- Department of Genetics, Genomics and Informatics, University of Tennessee Health Science Center, Memphis, TN, USA
- Genomics Research Centre, Human Technopole, Milan, Italy
| | - Leanne Haggerty
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Reza Halabian
- Institute of Bioinformatics, Faculty of Medicine, University of Münster, Münster, Germany
| | - Nancy F Hansen
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Robert Harris
- Department of Biology, Pennsylvania State University, University Park, PA, USA
| | - Gabrielle A Hartley
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
| | - William T Harvey
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Marina Haukness
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Jakob Heinz
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Thibaut Hourlier
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | | | - Sarah E Hunt
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Stephen Hwang
- XDBio Program, Johns Hopkins University, Baltimore, MD, USA
| | - Miten Jain
- Department of Bioengineering, Department of Physics, Northeastern University, Boston, MA, USA
| | - Rupesh K Kesharwani
- Human Genome Sequencing Center, Baylor College of Medicine, One Baylor Plaza, Houston, TX, USA
| | - Alexandra P Lewis
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Heng Li
- Department of Data Sciences, Dana-Farber Cancer Institute, Boston, MA, USA
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA
| | - Glennis A Logsdon
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Julian K Lucas
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Wojciech Makalowski
- Institute of Bioinformatics, Faculty of Medicine, University of Münster, Münster, Germany
| | - Christopher Markovic
- Genome Technology Access Center at the McDonnell Genome Institute, Washington University, St. Louis, MO, USA
| | - Fergal J Martin
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Ann M Mc Cartney
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Rajiv C McCoy
- Department of Biology, Johns Hopkins University, Baltimore, MD, USA
| | - Jennifer McDaniel
- Biosystems and Biomaterials Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | - Brandy M McNulty
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Paul Medvedev
- Department of Computer Science and Engineering, Pennsylvania State University, University Park, PA, USA
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, PA, USA
- Center for Computational Biology and Bioinformatics, Pennsylvania State University, University Park, PA, USA
| | - Alla Mikheenko
- Center for Algorithmic Biotechnology, Saint Petersburg State University, St Petersburg, Russia
- UCL Queen Square Institute of Neurology, UCL, London, UK
| | - Katherine M Munson
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Terence D Murphy
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
| | - Hugh E Olsen
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Nathan D Olson
- Biosystems and Biomaterials Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | - Luis F Paulin
- Human Genome Sequencing Center, Baylor College of Medicine, One Baylor Plaza, Houston, TX, USA
| | - David Porubsky
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Tamara Potapova
- Stowers Institute for Medical Research, Kansas City, MO, USA
| | - Fedor Ryabov
- Masters Program in National Research University Higher School of Economics, Moscow, Russia
| | - Steven L Salzberg
- Departments of Biomedical Engineering, Computer Science, and Biostatistics, Johns Hopkins University, Baltimore, MD, USA
| | | | - Fritz J Sedlazeck
- Human Genome Sequencing Center, Baylor College of Medicine, One Baylor Plaza, Houston, TX, USA
- Department of Computer Science, Rice University, Houston, TX, USA
| | | | | | - Alaina Shumate
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | | | - Likhitha Surapaneni
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Angela M Taravella Oill
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Françoise Thibaud-Nissen
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA
| | - Winston Timp
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Marta Tomaszkiewicz
- Department of Biology, Pennsylvania State University, University Park, PA, USA
- Department of Biomedical Engineering, Pennsylvania State University, State College, PA, USA
| | - Mitchell R Vollger
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
| | - Brian P Walenz
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Allison C Watwood
- Department of Biology, Pennsylvania State University, University Park, PA, USA
| | | | | | - Melissa A Wilson
- Center for Evolution and Medicine, School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Samantha Zarate
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
| | - Yiming Zhu
- Human Genome Sequencing Center, Baylor College of Medicine, One Baylor Plaza, Houston, TX, USA
| | - Justin M Zook
- Biosystems and Biomaterials Division, National Institute of Standards and Technology, Gaithersburg, MD, USA
| | - Evan E Eichler
- Department of Genome Sciences, University of Washington School of Medicine, Seattle, WA, USA
- Investigator, Howard Hughes Medical Institute, University of Washington, Seattle, WA, USA
| | - Rachel J O'Neill
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
- Institute for Systems Genomics, University of Connecticut, Storrs, CT, USA
- Department of Genetics and Genome Sciences, UConn Health, Farmington, CT, USA
| | - Michael C Schatz
- Department of Biology, Johns Hopkins University, Baltimore, MD, USA
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA
| | - Karen H Miga
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA, USA
- UC Santa Cruz Genomics Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Kateryna D Makova
- Department of Biology, Pennsylvania State University, University Park, PA, USA
| | - Adam M Phillippy
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
7
|
Tomaszkiewicz M, Sahlin K, Medvedev P, Makova KD. Transcript Isoform Diversity of Ampliconic Genes on the Y Chromosome of Great Apes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.02.530874. [PMID: 36993458 PMCID: PMC10054944 DOI: 10.1101/2023.03.02.530874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Y-chromosomal Ampliconic Genes (YAGs) are important for male fertility, as they encode proteins functioning in spermatogenesis. The variation in copy number and expression levels of these multicopy gene families has been recently studied in great apes, however, the diversity of splicing variants remains unexplored. Here we deciphered the sequences of polyadenylated transcripts of all nine YAG families (BPY2, CDY, DAZ, HSFY, PRY, RBMY, TSPY, VCY, and XKRY) from testis samples of six great ape species (human, chimpanzee, bonobo, gorilla, Bornean orangutan, and Sumatran orangutan). To achieve this, we enriched YAG transcripts with capture-probe hybridization and sequenced them with long (Pacific Biosciences) reads. Our analysis of this dataset resulted in several findings. First, we uncovered a high diversity of YAG transcripts across great apes. Second, we observed evolutionarily conserved alternative splicing patterns for most YAG families except for BPY2 and PRY. Our results suggest that BPY2 transcripts and predicted proteins in several great ape species (bonobo and the two orangutans) have independent evolutionary origins and are not homologous to human reference transcripts and proteins. In contrast, our results suggest that the PRY gene family, having the highest representation of transcripts without open reading frames, has been undergoing pseudogenization. Third, even though we have identified many species-specific protein-coding YAG transcripts, we have not detected any signatures of positive selection. Overall, our work illuminates the YAG isoform landscape and its evolutionary history, and provides a genomic resource for future functional studies focusing on infertility phenotypes in humans and critically endangered great apes.
Collapse
Affiliation(s)
- Marta Tomaszkiewicz
- Department of Biomedical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
| | - Kristoffer Sahlin
- Department of Mathematics, Science for Life Laboratory, Stockholm University, 106 91, Stockholm, Sweden
| | - Paul Medvedev
- Department of Computer Science and Engineering, The Pennsylvania State University
- Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Medical Genomics, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Computational Biology and Bioinformatics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Kateryna D Makova
- Center for Medical Genomics, The Pennsylvania State University, University Park, PA 16802, USA
- Center for Computational Biology and Bioinformatics, The Pennsylvania State University, University Park, PA 16802, USA
- Department of Biology, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
8
|
Rho NY, Mogas T, King WA, Favetta LA. Testis-Specific Protein Y-Encoded (TSPY) Is Required for Male Early Embryo Development in Bos taurus. Int J Mol Sci 2023; 24:ijms24043349. [PMID: 36834761 PMCID: PMC9959854 DOI: 10.3390/ijms24043349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2023] [Revised: 01/27/2023] [Accepted: 02/05/2023] [Indexed: 02/10/2023] Open
Abstract
TSPY is a highly conserved multi-copy gene with copy number variation (CNV) among species, populations, individuals and within families. TSPY has been shown to be involved in male development and fertility. However, information on TSPY in embryonic preimplantation stages is lacking. This study aims to determine whether TSPY CNV plays a role in male early development. Using sex-sorted semen from three different bulls, male embryo groups referred to as 1Y, 2Y and 3Y, were produced by in vitro fertilization (IVF). Developmental competency was assessed by cleavage and blastocyst rates. Embryos at different developmental stages were analyzed for TSPY CN, mRNA and protein levels. Furthermore, TSPY RNA knockdown was performed and embryos were assessed as per above. Development competency was only significantly different at the blastocyst stage, with 3Y being the highest. TSPY CNV and transcripts were detected in the range of 20-75 CN for 1Y, 20-65 CN for 2Y and 20-150 CN for 3Y, with corresponding averages of 30.2 ± 2.5, 33.0 ± 2.4 and 82.3 ± 3.6 copies, respectively. TSPY transcripts exhibited an inverse logarithmic pattern, with 3Y showing significantly higher TSPY. TSPY proteins, detected only in blastocysts, were not significantly different among groups. TSPY knockdown resulted in a significant TSPY depletion (p < 0.05), with no development observed after the eight-cell stage in male embryos, suggesting that TSPY is required for male embryo development.
Collapse
Affiliation(s)
- Na-Young Rho
- Reproductive Health and Biotechnology Lab, Department of Biomedical Sciences, Ontario Veterinary College, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - Teresa Mogas
- Department of Medicine and Animal Surgery, Autonomous University of Barcelona, Cerdanyola del Vallés, 08193 Barcelona, Spain
| | - W. Allan King
- Reproductive Health and Biotechnology Lab, Department of Biomedical Sciences, Ontario Veterinary College, University of Guelph, Guelph, ON N1G 2W1, Canada
- Karyotekk Inc., Ontario Veterinary College, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - Laura A. Favetta
- Reproductive Health and Biotechnology Lab, Department of Biomedical Sciences, Ontario Veterinary College, University of Guelph, Guelph, ON N1G 2W1, Canada
- Correspondence:
| |
Collapse
|
9
|
Castaneda C, Radović L, Felkel S, Juras R, Davis BW, Cothran EG, Wallner B, Raudsepp T. Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups. G3 (BETHESDA, MD.) 2022; 12:jkac278. [PMID: 36227030 PMCID: PMC9713435 DOI: 10.1093/g3journal/jkac278] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 10/08/2022] [Indexed: 11/03/2023]
Abstract
Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.
Collapse
Affiliation(s)
- Caitlin Castaneda
- Department of Veterinary Integrative Biosciences, School of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 7784-4458, USA
| | - Lara Radović
- Department of Biomedical Sciences, Institute of Animal Breeding and Genetics, University of Veterinary Medicine Vienna, Vienna 1210, Austria
- Vienna Graduate School of Population Genetics, University of Veterinary Medicine Vienna, Vienna 1210, Austria
| | - Sabine Felkel
- Department of Biomedical Sciences, Institute of Animal Breeding and Genetics, University of Veterinary Medicine Vienna, Vienna 1210, Austria
- Vienna Graduate School of Population Genetics, University of Veterinary Medicine Vienna, Vienna 1210, Austria
- Department of Biotechnology, Institute of Computational Biology, BOKU University of Life Sciences and Natural Resources, Vienna 1190, Austria
| | - Rytis Juras
- Department of Veterinary Integrative Biosciences, School of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 7784-4458, USA
| | - Brian W Davis
- Department of Veterinary Integrative Biosciences, School of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 7784-4458, USA
| | - Ernest Gus Cothran
- Department of Veterinary Integrative Biosciences, School of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 7784-4458, USA
| | - Barbara Wallner
- Department of Biomedical Sciences, Institute of Animal Breeding and Genetics, University of Veterinary Medicine Vienna, Vienna 1210, Austria
| | - Terje Raudsepp
- Department of Veterinary Integrative Biosciences, School of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 7784-4458, USA
| |
Collapse
|
10
|
PRAMEY: A Bovid-Specific Y-Chromosome Multicopy Gene Is Highly Related to Postnatal Testicular Growth in Hu Sheep. Animals (Basel) 2022; 12:ani12182380. [PMID: 36139240 PMCID: PMC9495132 DOI: 10.3390/ani12182380] [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: 07/28/2022] [Revised: 09/03/2022] [Accepted: 09/04/2022] [Indexed: 11/16/2022] Open
Abstract
PRAMEY (preferentially expressed antigen in melanoma, Y-linked) belongs to the cancer-testis antigens (CTAs) gene family and is predominantly expressed in testis, playing important roles in spermatogenesis and testicular development. This study cloned the full-length cDNA sequence of ovine PRAMEY using the rapid amplification of cDNA ends (RACE) method and analyzed the expression profile and copy number variation (CNV) of PRAMEY using quantitative real-time PCR (qPCR). The results revealed that the PRAMEY cDNA was 2099 bp in length with an open reading frame (ORF) of 1536 bp encoding 511 amino acids. PRAMEY was predominantly expressed in the testis and significantly upregulated during postnatal testicular development. The median copy number (MCN) of PRAMEY was 4, varying from 2 to 25 in 710 rams across eight sheep breeds. There was no significant correlation between the CNV of PRAMEY and testicular size, while a significant positive correlation was observed between the mRNA expression and testicular size in Hu sheep. The current study suggests that the expression levels of PRAMEY were closely associated with testicular size, indicating that PRAMEY may play an important role in testicular growth.
Collapse
|
11
|
Pei S, Xu H, Wang L, Li F, Li W, Yue X. Copy number variation of ZNF280BY across eight sheep breeds and its association with testicular size of Hu sheep. J Anim Sci 2022; 100:6624001. [PMID: 35775620 DOI: 10.1093/jas/skac232] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2022] [Accepted: 06/29/2022] [Indexed: 11/13/2022] Open
Abstract
ZNF280BY, a bovid-specific Y chromosome gene, was firstly found to be highly expressed in bovine testis, indicating it may play important roles in testicular development and male fertility. In this study, we firstly cloned the full-length cDNA of ovine ZNF280BY containing 1993 bp, and with a 1632 bp open reading frame. ZNF280BY was predominantly expressed in the testis, and its expression level was significantly higher in large testis than in small testis in Hu sheep at 6 months of age. In addition, the expression level of ZNF280BY significantly increased during testicular development, showing the highest expression level at 12 months of age. ZNF280BY showed copy number variation (CNV) in 723 rams from eight sheep breeds, ranging from 17 to 514 copies, with a median copy number of 188. Pearson correlation analysis showed that the CNV of ZNF280BY was negatively correlated with testis size in Hu sheep. Furthermore, its mRNA expression level in testis had no significant correlation with the CNV but was significantly correlated with testis size. This study concluded that the expression of ZNF280BY was closely related to testicular development, and the CNV of ZNF280BY could be used as an important genetic marker to evaluate the ram reproductive capacity at an early stage in Hu sheep.
Collapse
Affiliation(s)
- Shengwei Pei
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - Haiyue Xu
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - Li Wang
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - Fadi Li
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - Wanhong Li
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - Xiangpeng Yue
- State Key Laboratory of Grassland Agro-ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| |
Collapse
|
12
|
Zhang C, Zhao J, Guo Y, Xu Q, Liu M, Cheng M, Chao X, Schinckel AP, Zhou B. Genome-Wide Detection of Copy Number Variations and Evaluation of Candidate Copy Number Polymorphism Genes Associated With Complex Traits of Pigs. Front Vet Sci 2022; 9:909039. [PMID: 35847642 PMCID: PMC9280686 DOI: 10.3389/fvets.2022.909039] [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: 03/31/2022] [Accepted: 06/09/2022] [Indexed: 12/12/2022] Open
Abstract
Copy number variation (CNV) has been considered to be an important source of genetic variation for important phenotypic traits of livestock. In this study, we performed whole-genome CNV detection on Suhuai (SH) (n = 23), Chinese Min Zhu (MZ) (n = 11), and Large White (LW) (n = 12) pigs based on next-generation sequencing data. The copy number variation regions (CNVRs) were annotated and analyzed, and 10,885, 10,836, and 10,917 CNVRs were detected in LW, MZ, and SH pigs, respectively. Some CNVRs have been randomly selected for verification of the variation type by real-time PCR. We found that SH and LW pigs are closely related, while MZ pigs are distantly related to the SH and LW pigs by CNVR-based genetic structure, PCA, VST, and QTL analyses. A total of 14 known genes annotated in CNVRs were unique for LW pigs. Among them, the cyclin T2 (CCNT2) is involved in cell proliferation and the cell cycle. The FA Complementation Group M (FANCM) is involved in defective DNA repair and reproductive cell development. Ten known genes annotated in 47 CNVRs were unique for MZ pigs. The genes included glycerol-3-phosphate acyltransferase 3 (GPAT3) is involved in fat synthesis and is essential to forming the glycerol triphosphate. Glutathione S-transferase mu 4 (GSTM4) gene plays an important role in detoxification. Eleven known genes annotated in 23 CNVRs were unique for SH pigs. Neuroligin 4 X-linked (NLGN4X) and Neuroligin 4 Y-linked (NLGN4Y) are involved with nerve disorders and nerve signal transmission. IgLON family member 5 (IGLON5) is related to autoimmunity and neural activities. The unique characteristics of LW, MZ, and SH pigs are related to these genes with CNV polymorphisms. These findings provide important information for the identification of candidate genes in the molecular breeding of pigs.
Collapse
Affiliation(s)
- Chunlei Zhang
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Jing Zhao
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Yanli Guo
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Qinglei Xu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Mingzheng Liu
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Meng Cheng
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Xiaohuan Chao
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
| | - Allan P. Schinckel
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Bo Zhou
- College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, China
- *Correspondence: Bo Zhou
| |
Collapse
|
13
|
Larson EL, Kopania EEK, Hunnicutt KE, Vanderpool D, Keeble S, Good JM. Stage-specific disruption of X chromosome expression during spermatogenesis in sterile house mouse hybrids. G3 (BETHESDA, MD.) 2022; 12:jkab407. [PMID: 34864964 PMCID: PMC9210296 DOI: 10.1093/g3journal/jkab407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Accepted: 11/16/2021] [Indexed: 01/09/2023]
Abstract
Hybrid sterility is a complex phenotype that can result from the breakdown of spermatogenesis at multiple developmental stages. Here, we disentangle two proposed hybrid male sterility mechanisms in the house mice, Mus musculus domesticus and M. m. musculus, by comparing patterns of gene expression in sterile F1 hybrids from a reciprocal cross. We found that hybrid males from both cross directions showed disrupted X chromosome expression during prophase of meiosis I consistent with a loss of meiotic sex chromosome inactivation (MSCI) and Prdm9-associated sterility, but that the degree of disruption was greater in mice with an M. m. musculus X chromosome consistent with previous studies. During postmeiotic development, gene expression on the X chromosome was only disrupted in one cross direction, suggesting that misexpression at this later stage was genotype-specific and not a simple downstream consequence of MSCI disruption which was observed in both reciprocal crosses. Instead, disrupted postmeiotic expression may depend on the magnitude of earlier disrupted MSCI, or the disruption of particular X-linked genes or gene networks. Alternatively, only hybrids with a potential deficit of Sly copies, a Y-linked ampliconic gene family, showed overexpression in postmeiotic cells, consistent with a previously proposed model of antagonistic coevolution between the X- and Y-linked ampliconic genes contributing to disrupted expression late in spermatogenesis. The relative contributions of these two regulatory mechanisms and their impact on sterility phenotypes await further study. Our results further support the hypothesis that X-linked hybrid sterility in house mice has a variable genetic basis, and that genotype-specific disruption of gene regulation contributes to overexpression of the X chromosome at different stages of development. Overall, these findings underscore the critical role of epigenetic regulation of the X chromosome during spermatogenesis and suggest that these processes are prone to disruption in hybrids.
Collapse
Affiliation(s)
- Erica L Larson
- Department of Biological Sciences, University of Denver, Denver, CO 80208, USA
| | - Emily E K Kopania
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| | - Kelsie E Hunnicutt
- Department of Biological Sciences, University of Denver, Denver, CO 80208, USA
| | - Dan Vanderpool
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| | - Sara Keeble
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| | - Jeffrey M Good
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, USA
| |
Collapse
|
14
|
Vegesna R, Tomaszkiewicz M, Ryder OA, Campos-Sánchez R, Medvedev P, DeGiorgio M, Makova KD. Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels. Genome Biol Evol 2021; 12:842-859. [PMID: 32374870 PMCID: PMC7313670 DOI: 10.1093/gbe/evaa088] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/28/2020] [Indexed: 12/16/2022] Open
Abstract
Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.
Collapse
Affiliation(s)
- Rahulsimham Vegesna
- Bioinformatics and Genomics Graduate Program, The Huck Institutes for the Life Sciences, Pennsylvania State University, University Park
| | | | - Oliver A Ryder
- Institute for Conservation Research, San Diego Zoo Global, San Diego, California
| | | | - Paul Medvedev
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park.,Department of Computer Science and Engineering, Pennsylvania State University, University Park.,Center for Computational Biology and Bioinformatics, Pennsylvania State University, University Park.,Center for Medical Genomics, Pennsylvania State University, University Park
| | - Michael DeGiorgio
- Department of Biology, Pennsylvania State University, University Park.,Institute for Computational and Data Science, Pennsylvania State University, University Park
| | - Kateryna D Makova
- Department of Biology, Pennsylvania State University, University Park.,Center for Computational Biology and Bioinformatics, Pennsylvania State University, University Park.,Center for Medical Genomics, Pennsylvania State University, University Park
| |
Collapse
|
15
|
Vogt PH, Bender U, Deibel B, Kiesewetter F, Zimmer J, Strowitzki T. Human AZFb deletions cause distinct testicular pathologies depending on their extensions in Yq11 and the Y haplogroup: new cases and review of literature. Cell Biosci 2021; 11:60. [PMID: 33766143 PMCID: PMC7995748 DOI: 10.1186/s13578-021-00551-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 02/06/2021] [Indexed: 02/07/2023] Open
Abstract
Genomic AZFb deletions in Yq11 coined “classical” (i.e. length of Y DNA deletion: 6.23 Mb) are associated with meiotic arrest (MA) of patient spermatogenesis, i.e., absence of any postmeiotic germ cells. These AZFb deletions are caused by non-allelic homologous recombination (NAHR) events between identical sequence blocks located in the proximal arm of the P5 palindrome and within P1.2, a 92 kb long sequence block located in the P1 palindrome structure of AZFc in Yq11. This large genomic Y region includes deletion of 6 protein encoding Y genes, EIFA1Y, HSFY, PRY, RBMY1, RPS4Y, SMCY. Additionally, one copy of CDY2 and XKRY located in the proximal P5 palindrome and one copy of BPY1, two copies of DAZ located in the P2 palindrome, and one copy of CDY1 located proximal to P1.2 are included within this AZFb microdeletion. It overlaps thus distally along 2.3 Mb with the proximal part of the genomic AZFc deletion. However, AZFb deletions have been also reported with distinct break sites in the proximal and/or distal AZFb breakpoint intervals on the Y chromosome of infertile men. These so called “non-classical” AZFb deletions are associated with variable testicular pathologies, including meiotic arrest, cryptozoospermia, severe oligozoospermia, or oligoasthenoteratozoospermia (OAT syndrome), respectively. This raised the question whether there are any specific length(s) of the AZFb deletion interval along Yq11 required to cause meiotic arrest of the patient’s spermatogenesis, respectively, whether there is any single AZFb Y gene deletion also able to cause this “classical” AZFb testicular pathology? Review of the literature and more cases with “classical” and “non-classical” AZFb deletions analysed in our lab since the last 20 years suggests that the composition of the genomic Y sequence in AZFb is variable in men with distinct Y haplogroups especially in the distal AZFb region overlapping with the proximal AZFc deletion interval and that its extension can be “polymorphic” in the P3 palindrome. That means this AZFb subinterval can be rearranged or deleted also on the Y chromosome of fertile men. Any AZFb deletion observed in infertile men with azoospermia should therefore be confirmed as “de novo” mutation event, i.e., not present on the Y chromosome of the patient’s father or fertile brother before it is considered as causative agent for man’s infertility. Moreover, its molecular length in Yq11 should be comparable to that of the “classical” AZFb deletion, before meiotic arrest is prognosed as the patient’s testicular pathology.
Collapse
Affiliation(s)
- P H Vogt
- Division of Reproduction Genetics, Department of Gynaecol. Endocrinology & Infertility Disorders, Women Hospital, University of Heidelberg, Im Neuenheimer Feld 440, 69120, Heidelberg, Germany.
| | - U Bender
- Division of Reproduction Genetics, Department of Gynaecol. Endocrinology & Infertility Disorders, Women Hospital, University of Heidelberg, Im Neuenheimer Feld 440, 69120, Heidelberg, Germany
| | - B Deibel
- Division of Reproduction Genetics, Department of Gynaecol. Endocrinology & Infertility Disorders, Women Hospital, University of Heidelberg, Im Neuenheimer Feld 440, 69120, Heidelberg, Germany
| | - F Kiesewetter
- Department of Andrology, University Clinic of Dermatology, Erlangen, Germany
| | - J Zimmer
- Division of Reproduction Genetics, Department of Gynaecol. Endocrinology & Infertility Disorders, Women Hospital, University of Heidelberg, Im Neuenheimer Feld 440, 69120, Heidelberg, Germany
| | - T Strowitzki
- Department of Gynaecol. Endocrinology & Infertility Disorders, Women Hospital, University of Heidelberg, Heidelberg, Germany
| |
Collapse
|
16
|
Rogers TF, Pizzari T, Wright AE. Multi-Copy Gene Family Evolution on the Avian W Chromosome. J Hered 2021; 112:250-259. [PMID: 33758922 DOI: 10.1093/jhered/esab016] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 12/20/2020] [Accepted: 03/19/2021] [Indexed: 12/13/2022] Open
Abstract
The sex chromosomes often follow unusual evolutionary trajectories. In particular, the sex-limited chromosomes frequently exhibit a small but unusual gene content in numerous species, where many genes have undergone massive gene amplification. The reasons for this remain elusive with a number of recent studies implicating meiotic drive, sperm competition, genetic drift, and gene conversion in the expansion of gene families. However, our understanding is primarily based on Y chromosome studies as few studies have systematically tested for copy number variation on W chromosomes. Here, we conduct a comprehensive investigation into the abundance, variability, and evolution of ampliconic genes on the avian W. First, we quantified gene copy number and variability across the duck W chromosome. We find a limited number of gene families as well as conservation in W-linked gene copy number across duck breeds, indicating that gene amplification may not be such a general feature of sex chromosome evolution as Y studies would initially suggest. Next, we investigated the evolution of HINTW, a prominent ampliconic gene family hypothesized to play a role in female reproduction and oogenesis. In particular, we investigated the factors driving the expansion of HINTW using contrasts between modern chicken and duck breeds selected for different female-specific selection regimes and their wild ancestors. Although we find the potential for selection related to fecundity in explaining small-scale gene amplification of HINTW in the chicken, purifying selection seems to be the dominant mode of evolution in the duck. Together, this challenges the assumption that HINTW is key for female fecundity across the avian phylogeny.
Collapse
Affiliation(s)
- Thea F Rogers
- Department of Animal and Plant Sciences, University of Sheffield, UK
| | - Tommaso Pizzari
- Department of Animal and Plant Sciences, University of Sheffield, UK
| | - Alison E Wright
- Edward Grey Institute, Department of Zoology, University of Oxford, UK
| |
Collapse
|
17
|
Rogers MJ. Y chromosome copy number variation and its effects on fertility and other health factors: a review. Transl Androl Urol 2021; 10:1373-1382. [PMID: 33850773 PMCID: PMC8039628 DOI: 10.21037/tau.2020.04.06] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
The Y chromosome is essential for testis development and spermatogenesis. It is a chromosome with the lowest gene density owing to its medium size but paucity of coding genes. The Y chromosome is unique in that the majority of its structure is highly repetitive sequences, with the majority of these limited genes occurring in 9 amplionic sequences throughout the chromosome. The repetitive nature has its benefits as it can be protective against gene loss over many generations, but it can also predispose the Y chromosome to having wide variations of the number of gene copies present in these repeated sequences. This is known as copy number variation. Copy number variation is not unique to the Y chromosome but copy number variation is a well-known cause of male infertility and having effects on spermatogenesis. This is most commonly seen as deletions of the AZF sequences on the Y chromosome. However, there are other implications for copy number variation beyond just the AZF deletions that can affect spermatogenesis and potentially have other health implications. Copy number variations of TSPY1, DAZ, CDY1, RBMY1, the DYZ1 array, along with minor deletions of gr/gr, b1/b3, and b2/b3 have all be implicated in affecting spermatogenesis. UTY copy number variations have been implicated in risk for cardiovascular disease, and other deletions within gr/gr and the AZF sequences have been implicated in cancer and neuropsychiatric diseases. This review sets out to describe the Y chromosome and unique susceptibility to copy number variation and then to examine how this growing body of research impacts spermatogenesis and other health factors.
Collapse
Affiliation(s)
- Marc J Rogers
- Department of Urology, Medical University of South Carolina, Charleston, SC, USA
| |
Collapse
|
18
|
Martínez-Pacheco M, Tenorio M, Almonte L, Fajardo V, Godínez A, Fernández D, Cornejo-Páramo P, Díaz-Barba K, Halbert J, Liechti A, Székely T, Urrutia AO, Cortez D. Expression Evolution of Ancestral XY Gametologs across All Major Groups of Placental Mammals. Genome Biol Evol 2020; 12:2015-2028. [PMID: 32790864 PMCID: PMC7674692 DOI: 10.1093/gbe/evaa173] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/10/2020] [Indexed: 12/14/2022] Open
Abstract
Placental mammals present 180 million-year-old Y chromosomes that have retained a handful of dosage-sensitive genes. However, the expression evolution of Y-linked genes across placental groups has remained largely unexplored. Here, we expanded the number of Y gametolog sequences by analyzing ten additional species from previously unexplored groups. We detected seven remarkably conserved genes across 25 placental species with known Y repertoires. We then used RNA-seq data from 17 placental mammals to unveil the expression evolution of XY gametologs. We found that Y gametologs followed, on average, a 3-fold expression loss and that X gametologs also experienced some expression reduction, particularly in primates. Y gametologs gained testis specificity through an accelerated expression decay in somatic tissues. Moreover, despite the substantial expression decay of Y genes, the combined expression of XY gametologs in males is higher than that of both X gametologs in females. Finally, our work describes several features of the Y chromosome in the last common mammalian ancestor.
Collapse
Affiliation(s)
| | | | - Laura Almonte
- Center for Genome Sciences, UNAM, Cuernavaca, Mexico
| | | | - Alan Godínez
- Center for Genome Sciences, UNAM, Cuernavaca, Mexico
| | | | | | | | - Jean Halbert
- Center for Integrative Genomics, University of Lausanne, Switzerland
| | - Angelica Liechti
- Center for Integrative Genomics, University of Lausanne, Switzerland
| | - Tamas Székely
- Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Claverton Down, Bath, United Kingdom
| | - Araxi O Urrutia
- Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Claverton Down, Bath, United Kingdom.,Ecology Institute, UNAM, Mexico
| | - Diego Cortez
- Center for Genome Sciences, UNAM, Cuernavaca, Mexico
| |
Collapse
|
19
|
Abstract
The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered. The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.
Collapse
|