MRLR: unraveling high-resolution meiotic recombination by linked reads

MOTIVATION Meiotic recombination facilitates the transmission of exchanged genetic material between homologous chromosomes and plays a crucial role in increasing the genetic variations in eukaryotic organisms. In humans, thousands of crossover events have been identified by genotyping related family members. However, most of these crossover regions span tens to hundreds of kb, which is not sufficient resolution to accurately identify the crossover breakpoints in a typical trio family. RESULTS We have developed MRLR, a software using 10X linked reads to identify crossover events at a high resolution. By reconstructing the gamete genome, MRLR only requires a trio family dataset and can efficiently discover the crossover events. Using MRLR, we revealed a fine-scale pattern of crossover regions in six human families. From the two closest heterozygous alleles around the crossovers, we determined that MRLR achieved a median resolution 4.5kb. This method can delineate a genome-wide landscape of crossover events at a precise scale, which is important for both functional and genomic features analysis of meiotic recombination. AVAILABILITY MRLR is freely available at https://github.com/ChongLab/MRLR, implemented in Perl. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Bernat Gel,et al.  regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests , 2015, Bioinform..

[2]  Victor Guryev,et al.  Dense and accurate whole-chromosome haplotyping of individual genomes , 2017, Nature Communications.

[3]  P. Flicek,et al.  The Ensembl Regulatory Build , 2015, Genome Biology.

[4]  D. Gudbjartsson,et al.  A high-resolution recombination map of the human genome , 2002, Nature Genetics.

[5]  Yongjun Zhao,et al.  DNA template strand sequencing of single-cells maps genomic rearrangements at high resolution , 2012, Nature Methods.

[6]  Li Ding,et al.  Multi-platform discovery of haplotype-resolved structural variation in human genomes , 2018, Nature Communications.

[7]  P. Patel,et al.  Recombination hot spots and human disease. , 1997, Genome research.

[8]  Terrence S. Furey,et al.  The UCSC Table Browser data retrieval tool , 2004, Nucleic Acids Res..

[9]  G. Coop,et al.  High-Resolution Mapping of Crossovers Reveals Extensive Variation in Fine-Scale Recombination Patterns Among Humans , 2008, Science.

[10]  A. Jeffreys,et al.  Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex , 2001, Nature Genetics.

[11]  Ronald W Davis,et al.  Meiotic recombination generates rich diversity in NK cell receptor genes, alleles, and haplotypes. , 2009, Genome research.

[12]  T. Hassold,et al.  Variation in human meiotic recombination. , 2004, Annual review of genomics and human genetics.

[13]  H. Kajiho,et al.  Characterization of RIN3 as a Guanine Nucleotide Exchange Factor for the Rab5 Subfamily GTPase Rab31* , 2011, The Journal of Biological Chemistry.

[14]  C. Walsh,et al.  NFIA Haploinsufficiency Is Associated with a CNS Malformation Syndrome and Urinary Tract Defects , 2007, PLoS genetics.

[15]  C. L. Baker,et al.  PRDM9 binding organizes hotspot nucleosomes and limits Holliday junction migration , 2014, Genome research.

[16]  Gabor T. Marth,et al.  Haplotype-based variant detection from short-read sequencing , 2012, 1207.3907.

[17]  Zhaohui S. Qin,et al.  A second generation human haplotype map of over 3.1 million SNPs , 2007, Nature.

[18]  Hanlee P. Ji,et al.  Chromosome-scale mega-haplotypes enable digital karyotyping of cancer aneuploidy , 2017, Nucleic acids research.

[19]  Peter Donnelly,et al.  A common sequence motif associated with recombination hot spots and genome instability in humans , 2008, Nature Genetics.

[20]  Michael C. Schatz,et al.  LRSim: A Linked-Reads Simulator Generating Insights for Better Genome Partitioning , 2017, Computational and structural biotechnology journal.

[21]  P. Sung,et al.  Mechanism of eukaryotic homologous recombination. , 2008, Annual review of biochemistry.

[22]  Hanlee P. Ji,et al.  Haplotyping germline and cancer genomes using high-throughput linked-read sequencing , 2015, Nature Biotechnology.

[23]  Chunlin Xiao,et al.  Reproducible integration of multiple sequencing datasets to form high-confidence SNP, indel, and reference calls for five human genome reference materials , 2018, bioRxiv.

[24]  R. Camerini-Otero,et al.  Recombination initiation maps of individual human genomes , 2014, Science.

[25]  Mark Hills,et al.  Single-cell template strand sequencing by Strand-seq enables the characterization of individual homologs , 2017, Nature Protocols.

[26]  P. Hunt,et al.  To err (meiotically) is human: the genesis of human aneuploidy , 2001, Nature Reviews Genetics.

[27]  Victor Guryev,et al.  Direct chromosome-length haplotyping by single-cell sequencing , 2016, Genome research.

[28]  S. Ralston,et al.  Targeted sequencing of the Paget's disease associated 14q32 locus identifies several missense coding variants in RIN3 that predispose to Paget's disease of bone , 2015, Human molecular genetics.

[29]  Yun-Shien Lee,et al.  Analysis of human meiotic recombination events with a parent-sibling tracing approach , 2011, BMC Genomics.

[30]  G. Coop,et al.  An evolutionary view of human recombination , 2007, Nature Reviews Genetics.

[31]  Michael G. Sehorn,et al.  Human meiotic recombinase Dmc1 promotes ATP-dependent homologous DNA strand exchange , 2004, Nature.