Home LiteratureArticle Details
PMID: 26934921 Published · ppublish English Journal Article Research Support, N.I.H., Extramural Research Support, Non-U.S. Gov't Research Support, U.S. Gov't, Non-P.H.S.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y.

Genome research ·Vol. 26 ·No. 4 ·2016-04-00 ·Pages 530-40

Tomaszkiewicz M, Rangavittal S, Cechova M, Campos Sanchez R, Fescemyer HW, Harris R, Ye D, O'Brien PC, Chikhi R, Ryder OA, Ferguson-Smith MA, Medvedev P, Makova KD

Abstract

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

MeSH Terms
Animals Computational Biology/methods Gene Rearrangement Genome Genomics Gorilla gorilla/genetics High-Throughput Nucleotide Sequencing Humans Inverted Repeat Sequences Male Mammals/genetics Microsatellite Repeats Pan troglodytes/genetics Repetitive Sequences, Nucleic Acid Sequence Analysis, DNA Y Chromosome
Authors & Affiliations
13 authors, click to expand affiliations / ORCID
Tomaszkiewicz Marta
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Rangavittal Samarth
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Cechova Monika
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Campos Sanchez Rebeca
Genetics Program, The Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Fescemyer Howard W
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Harris Robert
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Ye Danling
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
O'Brien Patricia C M
Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, United Kingdom;
Chikhi Rayan
University of Lille 1/CNRS 59655 Villeneuve d'Ascq, France; Department of Computer Science and Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, USA; The Genome Sciences Institute of the Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
Ryder Oliver A
San Diego Zoo Institute for Conservation Research, Escondido, California 92027, USA;
Ferguson-Smith Malcolm A
Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, United Kingdom;
Medvedev Paul
Department of Computer Science and Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, USA; The Genome Sciences Institute of the Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA; Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA.
Makova Kateryna D
Department of Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
References (70)
70 references, click to expand
  1. Comparative sex chromosome genomics in snakes: differentiation, evolutionary strata, and lack of global dosage compensation.
    PLoS Biol. 2013;11(8):e1001643 PMID: 24015111
  2. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler.
    Gigascience. 2012 Dec 27;1(1):18 PMID: 23587118
  3. REAPR: a universal tool for genome assembly evaluation.
    Genome Biol. 2013;14(5):R47 PMID: 23710727
  4. Reversal of an ancient sex chromosome to an autosome in Drosophila.
    Nature. 2013 Jul 18;499(7458):332-5 PMID: 23792562
  5. 28-way vertebrate alignment and conservation track in the UCSC Genome Browser.
    Genome Res. 2007 Dec;17(12):1797-808 PMID: 17984227
  6. High sensitivity detection and quantitation of DNA copy number and single nucleotide variants with single color droplet digital PCR.
    Anal Chem. 2014 Mar 4;86(5):2618-24 PMID: 24483992
  7. Increased resolution of Y chromosome haplogroup T defines relationships among populations of the Near East, Europe, and Africa.
    Hum Biol. 2011 Feb;83(1):39-53 PMID: 21453003
  8. Do variations in substitution rates and male mutation bias correlate with life-history traits? A study of 32 mammalian genomes.
    Evolution. 2011 Oct;65(10):2800-15 PMID: 21967423
  9. Multiplexed target detection using DNA-binding dye chemistry in droplet digital PCR.
    Anal Chem. 2013 Dec 3;85(23):11619-27 PMID: 24180464
  10. BLAT--the BLAST-like alignment tool.
    Genome Res. 2002 Apr;12(4):656-64 PMID: 11932250
  11. Copy number of DAZ genes in infertile men.
    Fertil Steril. 2005 Nov;84(5):1522-5 PMID: 16275261
  12. Reduced copy number of DAZ genes in subfertile and infertile men.
    Fertil Steril. 2002 Jan;77(1):68-75 PMID: 11779593
  13. Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes.
    Cell. 2014 Nov 6;159(4):800-13 PMID: 25417157
  14. Evolution of X-degenerate Y chromosome genes in greater apes: conservation of gene content in human and gorilla, but not chimpanzee.
    J Mol Evol. 2009 Feb;68(2):134-44 PMID: 19142680
  15. Assembling large genomes with single-molecule sequencing and locality-sensitive hashing.
    Nat Biotechnol. 2015 Jun;33(6):623-30 PMID: 26006009
  16. Y genetic variation and phenotypic diversity in health and disease.
    Biol Sex Differ. 2015 Mar 13;6:6 PMID: 25866616
  17. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data.
    Nat Methods. 2013 Jun;10(6):563-9 PMID: 23644548
  18. High-throughput droplet digital PCR system for absolute quantitation of DNA copy number.
    Anal Chem. 2011 Nov 15;83(22):8604-10 PMID: 22035192
  19. Sequencing genomes from single cells by polymerase cloning.
    Nat Biotechnol. 2006 Jun;24(6):680-6 PMID: 16732271
  20. The evolution of mammalian gene families.
    PLoS One. 2006;1:e85 PMID: 17183716
  21. Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content.
    Nature. 2010 Jan 28;463(7280):536-9 PMID: 20072128
  22. Comparative analysis of mammalian Y chromosomes illuminates ancestral structure and lineage-specific evolution.
    Genome Res. 2013 Sep;23(9):1486-95 PMID: 23788650
  23. Sex-biased dispersal in western lowland gorillas (Gorilla gorilla gorilla).
    Mol Ecol. 2007 Jun;16(11):2247-59 PMID: 17561888
  24. Sex linkage, sex-specific selection, and the role of recombination in the evolution of sexually dimorphic gene expression.
    Evolution. 2010 Dec;64(12):3417-42 PMID: 20874735
  25. A recent bottleneck of Y chromosome diversity coincides with a global change in culture.
    Genome Res. 2015 Apr;25(4):459-66 PMID: 25770088
  26. Rates and patterns of great ape retrotransposition.
    Proc Natl Acad Sci U S A. 2013 Aug 13;110(33):13457-62 PMID: 23884656
  27. Initial characterization of the large genome of the salamander Ambystoma mexicanum using shotgun and laser capture chromosome sequencing.
    Sci Rep. 2015;5:16413 PMID: 26553646
  28. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes.
    Nature. 2003 Jun 19;423(6942):825-37 PMID: 12815422
  29. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.
    Nat Protoc. 2013 Aug;8(8):1494-512 PMID: 23845962
  30. Full-length transcriptome assembly from RNA-Seq data without a reference genome.
    Nat Biotechnol. 2011 Jul;29(7):644-52 PMID: 21572440
  31. High-quality draft assemblies of mammalian genomes from massively parallel sequence data.
    Proc Natl Acad Sci U S A. 2011 Jan 25;108(4):1513-8 PMID: 21187386
  32. How Ebola impacts social dynamics in gorillas: a multistate modelling approach.
    J Anim Ecol. 2015 Jan;84(1):166-76 PMID: 24995485
  33. Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators.
    Nature. 2014 Apr 24;508(7497):494-9 PMID: 24759411
  34. SSPACE-LongRead: scaffolding bacterial draft genomes using long read sequence information.
    BMC Bioinformatics. 2014;15:211 PMID: 24950923
  35. TSPY1 copy number variation influences spermatogenesis and shows differences among Y lineages.
    J Clin Endocrinol Metab. 2009 Oct;94(10):4016-22 PMID: 19773397
  36. Efficient identification of Y chromosome sequences in the human and Drosophila genomes.
    Genome Res. 2013 Nov;23(11):1894-907 PMID: 23921660
  37. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.
    Genome Biol. 2013;14(4):R36 PMID: 23618408
  38. Mechanism of chimera formation during the Multiple Displacement Amplification reaction.
    BMC Biotechnol. 2007;7:19 PMID: 17430586
  39. The evolution of gene expression levels in mammalian organs.
    Nature. 2011 Oct 20;478(7369):343-8 PMID: 22012392
  40. Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes.
    Nature. 2012 Mar 1;483(7387):82-6 PMID: 22367542
  41. DSK: k-mer counting with very low memory usage.
    Bioinformatics. 2013 Mar 1;29(5):652-3 PMID: 23325618
  42. The development of chromosome microdissection and microcloning technique and its applications in genomic research.
    Curr Genomics. 2007 Mar;8(1):67-72 PMID: 18645627
  43. Distinct mutational behaviors differentiate short tandem repeats from microsatellites in the human genome.
    Genome Biol Evol. 2013;5(3):606-20 PMID: 23241442
  44. Insights into hominid evolution from the gorilla genome sequence.
    Nature. 2012 Mar 8;483(7388):169-75 PMID: 22398555
  45. The origin and evolution of human ampliconic gene families and ampliconic structure.
    Genome Res. 2007 Apr;17(4):441-50 PMID: 17185645
  46. Rate of molecular evolution of the seminal protein gene SEMG2 correlates with levels of female promiscuity.
    Nat Genet. 2004 Dec;36(12):1326-9 PMID: 15531881
  47. Y chromosome TSPY copy numbers and semen quality.
    Fertil Steril. 2010 Oct;94(5):1744-7 PMID: 19917507
  48. Mutation rate and specificity analysis of tetranucleotide microsatellite DNA alleles in somatic human cells.
    Mol Carcinog. 2002 Jul;34(3):140-50 PMID: 12112308
  49. The temporal dynamics of processes underlying Y chromosome degeneration.
    Genetics. 2008 Jul;179(3):1513-25 PMID: 18562655
  50. Review: Sex chromosome evolution and the expression of sex-specific genes in the placenta.
    Placenta. 2010 Mar;31 Suppl:S27-32 PMID: 20163856
  51. BESST--efficient scaffolding of large fragmented assemblies.
    BMC Bioinformatics. 2014;15:281 PMID: 25128196
  52. Structure of a hypervariable tandemly repeated DNA sequence on the short arm of the human Y chromosome.
    J Mol Biol. 1988 Oct 20;203(4):837-48 PMID: 3210241
  53. Fast gapped-read alignment with Bowtie 2.
    Nat Methods. 2012 Apr;9(4):357-9 PMID: 22388286
  54. Simian Y chromosomes: species-specific rearrangements of DAZ, RBM, and TSPY versus contiguity of PAR and SRY.
    Mamm Genome. 1998 Mar;9(3):226-31 PMID: 9501307
  55. Origins and functional evolution of Y chromosomes across mammals.
    Nature. 2014 Apr 24;508(7497):488-93 PMID: 24759410
  56. Mind the gap: upgrading genomes with Pacific Biosciences RS long-read sequencing technology.
    PLoS One. 2012;7(11):e47768 PMID: 23185243
  57. Copy number variation in Y chromosome multicopy genes is linked to a paternal parent-of-origin effect on CNS autoimmune disease in female offspring.
    Genome Biol. 2015;16:28 PMID: 25886764
  58. Fast and accurate short read alignment with Burrows-Wheeler transform.
    Bioinformatics. 2009 Jul 15;25(14):1754-60 PMID: 19451168
  59. Abundant gene conversion between arms of palindromes in human and ape Y chromosomes.
    Nature. 2003 Jun 19;423(6942):873-6 PMID: 12815433
  60. Sequencing of rhesus macaque Y chromosome clarifies origins and evolution of the DAZ (Deleted in AZoospermia) genes.
    Bioessays. 2012 Dec;34(12):1035-44 PMID: 23055411
  61. The pig X and Y Chromosomes: structure, sequence, and evolution.
    Genome Res. 2016 Jan;26(1):130-9 PMID: 26560630
  62. Scaffolding pre-assembled contigs using SSPACE.
    Bioinformatics. 2011 Feb 15;27(4):578-9 PMID: 21149342
  63. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.
    J Comput Biol. 2012 May;19(5):455-77 PMID: 22506599
  64. RBMY genes and AZFb deletions.
    J Endocrinol Invest. 2000 Nov;23(10):652-8 PMID: 11097429
  65. Mapping single molecule sequencing reads using basic local alignment with successive refinement (BLASR): application and theory.
    BMC Bioinformatics. 2012;13:238 PMID: 22988817
  66. The degeneration of Y chromosomes.
    Philos Trans R Soc Lond B Biol Sci. 2000 Nov 29;355(1403):1563-72 PMID: 11127901
  67. A comparative study of karyotypes of muntjacs by chromosome painting.
    Chromosoma. 1995 May;103(9):642-52 PMID: 7587587
  68. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets.
    Bioinformatics. 2003 Mar 22;19(5):651-2 PMID: 12651724
  69. Male-specific region of the bovine Y chromosome is gene rich with a high transcriptomic activity in testis development.
    Proc Natl Acad Sci U S A. 2013 Jul 23;110(30):12373-8 PMID: 23842086
  70. Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway.
    Nat Genet. 2006 Feb;38(2):158-67 PMID: 16388311
Article Info
Journal
Genome research
Abbr.
Genome Res
ISSN
1549-5469
Published
2016-04-00
Epub
2016-00-02
Pages
530-40
Language
English
Region
United States
NLM ID
9518021
PMCID
PMC4817776
Subset
IM
Grants
NIGMS NIH HHS · T32 GM102057 · United States
NIGMS NIH HHS · 1T32GM102057-0A1 · United States
Analysis Services
Analysis Services

Contact

No. 2 Wenbo Road, Zhangqiu District, Jinan, Shandong

Qilu Normal University · Genelibs Bioinformatics Lab

750 Shunhua Rd, Jinan

2F, Bldg F, University Science Park

Tel: 0531-88819269

WeChat Official Account

Follow our WeChat subscription account for real-time updates and the latest in medical and biological research.


Business Email

E-mail: product@genelibs.com