High-resolution Recombination Research Articles

A familial, telomere-to-telomere reference for human de novo mutation and recombination from a four-generation pedigree.

Using five complementary short- and long-read sequencing technologies, we phased and assembled >95% of each diploid human genome in a four-generation, 28-member family (CEPH 1463) allowing us to systematically assess de novo mutations (DNMs) and recombination. From this family, we estimate an average of 192 DNMs per generation, including 75.5 de novo single-nucleotide variants (SNVs), 7.4 non-tandem repeat indels, 79.6 de novo indels or structural variants (SVs) originating from tandem repeats, 7.7 centromeric de novo SVs and SNVs, and 12.4 de novo Y chromosome events per generation. STRs and VNTRs are the most mutable with 32 loci exhibiting recurrent mutation through the generations. We accurately assemble 288 centromeres and six Y chromosomes across the generations, documenting de novo SVs, and demonstrate that the DNM rate varies by an order of magnitude depending on repeat content, length, and sequence identity. We show a strong paternal bias (75-81%) for all forms of germline DNM, yet we estimate that 17% of de novo SNVs are postzygotic in origin with no paternal bias. We place all this variation in the context of a high-resolution recombination map (~3.5 kbp breakpoint resolution). We observe a strong maternal recombination bias (1.36 maternal:paternal ratio) with a consistent reduction in the number of crossovers with increasing paternal (r=0.85) and maternal (r=0.65) age. However, we observe no correlation between meiotic crossover locations and de novo SVs, arguing against non-allelic homologous recombination as a predominant mechanism. The use of multiple orthogonal technologies, near-telomere-to-telomere phased genome assemblies, and a multi-generation family to assess transmission has created the most comprehensive, publicly available "truth set" of all classes of genomic variants. The resource can be used to test and benchmark new algorithms and technologies to understand the most fundamental processes underlying human genetic variation.

In Brassica rapa meiotic recombination; chromosome axis remodeling is critical

The primary mechanism by which breeders generate biodiversity is meiotic recombination, which permits genetic rearrangement at every generation. It facilitates the accumulation of advantageous alleles and eliminates harmful mutations. With the development of one to very rarely more than three crossovers that are not randomly distributed, this mechanism is highly regulated. To ensure correct parental chromosomal segregation and the creation of unique allelic combinations, meiotic recombination is essential. Finding the variables affecting the rate and distribution of meiotic crossings (COs) is crucial because this process is strictly regulated, especially for plant breeding initiatives. Nevertheless, few high-resolution recombination maps exist for most crops, including the Brassica genus, and little is known regarding sex differences and intraspecific variation. Here, using progenies of cross populations from the crossed F1s population, we reveal fine-scale resolution recombination landscapes for two female and one male cross in Brassica rapa. Parents are from very different lines, which stands for higher quality and yield.

The fine-scale recombination rate variation and associations with genomic features in a butterfly.

Recombination is a key molecular mechanism that has profound implications on both micro- and macroevolutionary processes. However, the determinants of recombination rate variation in holocentric organisms are poorly understood, in particular in Lepidoptera (moths and butterflies). The wood white butterfly (Leptidea sinapis) shows considerable intraspecific variation in chromosome numbers and is a suitable system for studying regional recombination rate variation and its potential molecular underpinnings. Here, we developed a large whole-genome resequencing data set from a population of wood whites to obtain high-resolution recombination maps using linkage disequilibrium information. The analyses revealed that larger chromosomes had a bimodal recombination landscape, potentially caused by interference between simultaneous chiasmata. The recombination rate was significantly lower in subtelomeric regions, with exceptions associated with segregating chromosome rearrangements, showing that fissions and fusions can have considerable effects on the recombination landscape. There was no association between the inferred recombination rate and base composition, supporting a limited influence of GC-biased gene conversion in butterflies. We found significant but variable associations between the recombination rate and the density of different classes of transposable elements, most notably a significant enrichment of short interspersed nucleotide elements in genomic regions with higher recombination rate. Finally, the analyses unveiled significant enrichment of genes involved in farnesyltranstransferase activity in recombination coldspots, potentially indicating that expression of transferases can inhibit formation of chiasmata during meiotic division. Our results provide novel information about recombination rate variation in holocentric organisms and have particular implications for forthcoming research in population genetics, molecular/genome evolution, and speciation.

Fine mapping of meiotic crossovers in Brassica oleracea reveals patterns and variations depending on direction and combination of crosses.

Meiotic recombination is crucial for assuring proper segregation of parental chromosomes and generation of novel allelic combinations. As this process is tightly regulated, identifying factors influencing rate, and distribution of meiotic crossovers (COs) is of major importance, notably for plant breeding programs. However, high-resolution recombination maps are sparse in most crops including the Brassica genus and knowledge about intraspecific variation and sex differences is lacking. Here, we report fine-scale resolution recombination landscapes for 10 female and 10 male crosses in Brassica oleracea, by analyzing progenies of five large four-way-cross populations from two reciprocally crossed F1s per population. Parents are highly diverse inbred lines representing major crops, including broccoli, cauliflower, cabbage, kohlrabi, and kale. We produced approximately 4.56T Illumina data from 1248 progenies and identified 15 353 CO across the 10 reciprocal crosses, 51.13% of which being mapped to <10kb. We revealed fairly similar Mb-scale recombination landscapes among all cross combinations and between the sexes, and provided evidence that these landscapes are largely independent of sequence divergence. We evidenced strong influence of gene density and large structural variations on CO formation in B. oleracea. Moreover, we found extensive variations in CO number depending on the direction and combination of the initial parents crossed with, for the first time, a striking interdependency between these factors. These data improve our current knowledge on meiotic recombination and are important for Brassica breeders.

Global dissection of the recombination landscape in soybean using a high-density 600K SoySNP array.

Recombination is crucial for crop breeding because it can break linkage drag and generate novel allele combinations. However, the high-resolution recombination landscape and its driving forces in soybean are largely unknown. Here, we constructed eight recombinant inbred line (RIL) populations and genotyped individual lines using the high-density 600K SoySNP array, which yielded a high-resolution recombination map with 5636 recombination sites at a resolution of 1.37 kb. The recombination rate was negatively correlated with transposable element density and GC content but positively correlated with gene density. Interestingly, we found that meiotic recombination was enriched at the promoters of active genes. Further investigations revealed that chromatin accessibility and active epigenetic modifications promoted recombination. Our findings provide important insights into the control of homologous recombination and thus will increase our ability to accelerate soybean breeding by manipulating meiotic recombination rate.

NanoCross: A pipeline that detecting recombinant crossover using ONT sequencing data

Meiotic recombination is crucial for eukaryotes but varies among taxonomic scales (between individuals, groups, species, etc.) and genome resolutions. Studying how and why recombination rates change can help us understand the molecular basis and mechanisms of genetics and evolution. We introduce a genome-wide identification script called NanoCross, which uses ONT sequences to detect pooled gamete DNA cross recombination events. NanoCross first reduced sequencing errors and then constructed individual haplotypes based on homopolymer-filtered ONT sequences. Then, each molecule read is used to estimate cross recombination. In the case of moderate heterozygous variation density and sequencing depth, simulations revealed that our technique offers a good level of sensitivity and specificity. We constructed a high-resolution recombination map of wild boar genomes using NanoCross and compared it to recombination maps of male breeding pig populations. NanoCross provides us with a method and scripts for constructing a high-resolution individual genome recombination map utilizing long-read sequencing, as well as a novel approach for examining the variation in individual recombination rate. The source code and data mechanism are hosted on GitHub (https://github.com/zuoquanchen/NanoCross).

Inferring recombination patterns in African populations.

Although several high-resolution recombination maps exist for European-descent populations, the recombination landscape of African populations remains relatively understudied. Given that there is high genetic divergence among groups in Africa, it is possible that recombination hotspots also diverge significantly. Both limitations and opportunities exist for developing recombination maps for these populations. In this review, we discuss various recombination inference methods, and the strengths and weaknesses of these methods in analyzing recombination in African-descent populations. Furthermore, we provide a decision tree and recommendations for which inference method to use in various research contexts. Establishing an appropriate methodology for recombination rate inference in a particular study will improve the accuracy of various downstream analyses including but not limited to local ancestry inference, haplotype phasing, fine-mapping of GWAS loci and genome assemblies.

High-resolution population-specific recombination rates and their effect on phasing and genotype imputation

Previous research has shown that using population-specific reference panels has a significant effect on downstream population genomic analyses like haplotype phasing, genotype imputation, and association, especially in the context of population isolates. Here, we developed a high-resolution recombination rate mapping at 10 and 50 kb scale using high-coverage (20–30×) whole-genome sequenced data of 55 family trios from Finland and compared it to recombination rates of non-Finnish Europeans (NFE). We tested the downstream effects of the population-specific recombination rates in statistical phasing and genotype imputation in Finns as compared to the same analyses performed by using the NFE-based recombination rates. We found that Finnish recombination rates have a moderately high correlation (Spearman’s ρ = 0.67–0.79) with NFE, although on average (across all autosomal chromosomes), Finnish rates (2.268 ± 0.4209 cM/Mb) are 12–14% lower than NFE (2.641 ± 0.5032 cM/Mb). Finnish recombination map was found to have no significant effect in haplotype phasing accuracy (switch error rates ~2%) and average imputation concordance rates (97–98% for common, 92–96% for low frequency and 78–90% for rare variants). Our results suggest that haplotype phasing and genotype imputation mostly depend on population-specific contexts like appropriate reference panels and their sample size, but not on population-specific recombination maps. Even though recombination rate estimates had some differences between the Finnish and NFE populations, haplotyping and imputation had not been noticeably affected by the recombination map used. Therefore, the currently available HapMap recombination maps seem robust for population-specific phasing and imputation pipelines, even in the context of relatively isolated populations like Finland.

Genetic dissection of an allotetraploid interspecific CSSLs guides interspecific genetics and breeding in cotton

BackgroundThe low genetic diversity of Upland cotton limits the potential for genetic improvement. Making full use of the genetic resources of Sea-island cotton will facilitate genetic improvement of widely cultivated Upland cotton varieties. The chromosome segments substitution lines (CSSLs) provide an ideal strategy for mapping quantitative trait loci (QTL) in interspecific hybridization.ResultsIn this study, a CSSL population was developed by PCR-based markers assisted selection (MAS), derived from the crossing and backcrossing of Gossypium hirsutum (Gh) and G. barbadense (Gb), firstly. Then, by whole genome re-sequencing, 11,653,661 high-quality single nucleotide polymorphisms (SNPs) were identified which ultimately constructed 1211 recombination chromosome introgression segments from Gb. The sequencing-based physical map provided more accurate introgressions than the PCR-based markers. By exploiting CSSLs with mutant morphological traits, the genes responding for leaf shape and fuzz-less mutation in the Gb were identified. Based on a high-resolution recombination bin map to uncover genetic loci determining the phenotypic variance between Gh and Gb, 64 QTLs were identified for 14 agronomic traits with an interval length of 158 kb to 27 Mb. Surprisingly, multiple alleles of Gb showed extremely high value in enhancing cottonseed oil content (SOC).ConclusionsThis study provides guidance for studying interspecific inheritance, especially breeding researchers, for future studies using the traditional PCR-based molecular markers and high-throughput re-sequencing technology in the study of CSSLs. Available resources include candidate position for controlling cotton quality and quantitative traits, and excellent breeding materials. Collectively, our results provide insights into the genetic effects of Gb alleles on the Gh, and provide guidance for the utilization of Gb alleles in interspecific breeding.

Residual Heterozygosity and Epistatic Interactions Underlie the Complex Genetic Architecture of Yield in Diploid Potato.

Deconvolution of the genetic architecture underlying yield is critical for understanding bases of genetic gain in species of agronomic importance. To dissect the genetic components of yield in potato, we adopted a reference-based recombination map composed of four segregating alleles from an interspecific pseudotestcross F1 potato population (n = 90). Approximately 1.5 million short nucleotide variants were utilized during map construction, resulting in unprecedented resolution for an F1 population, estimated by a median bin length of 146 kb and 11 genes per bin. Regression models uncovered 14 quantitative trait loci (QTL) underpinning yield, average tuber weight, and tubers produced per plant in a population exhibiting a striking 332% average midparent-value heterosis. Nearly 80% of yield-associated QTL were epistatic, and contained between 0 and 44 annotated genes. We found that approximately one-half of epistatic QTL overlap regions of residual heterozygosity identified in the inbred parental parent (M6). Genomic regions recalcitrant to inbreeding were associated with an increased density of genes, many of which demonstrated signatures of selection and floral tissue specificity. Dissection of the genome-wide additive and dominance values for yield and yield components indicated a widespread prevalence of dominance contributions in this population, enriched at QTL and regions of residual heterozygosity. Finally, the effects of short nucleotide variants and patterns of gene expression were determined for all genes underlying yield-associated QTL, exposing several promising candidate genes for future investigation.

Evolution of the Yeast Recombination Landscape.

Meiotic recombination comprises crossovers and noncrossovers. Recombination, crossover in particular, shuffles mutations and impacts both the level of genetic polymorphism and the speed of adaptation. In many species, the recombination rate varies across the genome with hot and cold spots. The hotspot paradox hypothesis asserts that recombination hotspots are evolutionarily unstable due to self-destruction. However, the genomic landscape of double-strand breaks (DSBs), which initiate recombination, is evolutionarily conserved among divergent yeast species, casting doubt on the hotspot paradox hypothesis. Nonetheless, because only a subset of DSBs are associated with crossovers, the evolutionary conservation of the crossover landscape could differ from that of DSBs. Here, we investigate this possibility by generating a high-resolution recombination map of the budding yeast Saccharomyces paradoxus through whole-genome sequencing of 50 meiotic tetrads and by comparing this recombination map with that of S. cerevisiae. We observe a 40% lower recombination rate in S. paradoxus than in S. cerevisiae. Compared with the DSB landscape, the crossover landscape is even more conserved. Further analyses indicate that the elevated conservation of the crossover landscape is explained by a near-subtelomeric crossover preference in both yeasts, which we find to be attributable at least in part to crossover interference. We conclude that the yeast crossover landscape is highly conserved and that the evolutionary conservation of this landscape can differ from that of the DSB landscape.

Genome wide analysis of meiotic recombination in yeast: For a few SNPs more.

Diploid organisms undergo meiosis to produce haploid germ cells. Crossover events during meiosis promote genetic diversity and facilitate accurate chromosome segregation. The baker's yeast Saccharomyces cerevisiae is extensively used as a model for analysis of meiotic recombination. Conventional methods for measuring recombination events in S. cerevisiae have been limited by the number and density of genetic markers. Next generation sequencing (NGS)‐based analysis of hybrid yeast genomes bearing thousands of heterozygous single nucleotide polymorphism (SNP) markers has revolutionized analysis of meiotic recombination. By facilitating analysis of marker segregation in the whole genome with unprecedented resolution, this method has resulted in the generation of high‐resolution recombination maps in wild‐type and meiotic mutants. These studies have provided novel insights into the mechanism of meiotic recombination. In this review, we discuss the methodology, challenges, insights and future prospects of using NGS‐based methods for whole genome analysis of meiotic recombination. The objective is to facilitate the use of these high through‐put sequencing methods for the analysis of meiotic recombination given their power to provide significant new insights into the process. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology, 70(8):743–752, 2018

Detection of major loci associated with the variation of 18 important agronomic traits between Solanum pimpinellifolium and cultivated tomatoes.

Wild species can be used to improve various agronomic traits in cultivars; however, a limited understanding of the genetic basis underlying the morphological differences between wild and cultivated species hinders the integration of beneficial traits from wild species. In the present study, we generated and sequenced recombinant inbred lines (RILs, 201 F10 lines) derived from a cross between Solanum pimpinellifolium and Solanum lycopersicum tomatoes. Based on a high-resolution recombination bin map to uncover major loci determining the phenotypic variance between wild and cultivated tomatoes, 104 significantly associated loci were identified for 18 agronomic traits. On average, these loci explained ~39% of the phenotypic variance of the RILs. We further generated near-isogenic lines (NILs) for four identified loci, and the lines exhibited significant differences for the associated traits. We found that two loci could improve the flower number and inflorescence architecture in the cultivar following introgression of the wild-species alleles. These findings allowed us to construct a trait-locus network to help explain the correlations among different traits based on the pleiotropic or linked loci. Our results provide insights into the morphological changes between wild and cultivated tomatoes, and will help to identify key genes governing important agronomic traits for the molecular selection of elite tomato varieties.

Low-coverage resequencing detects meiotic recombination pattern and features in tomato RILs

Traditional plant breeding relies on meiotic recombination for mixing of parental alleles to create novel allele combinations. Detailed analysis of recombination patterns in model organisms shows that recombination is tightly regulated within the genome, but frequencies vary extensively along chromosomes. Despite being a model organism for fruit developmental studies, high-resolution recombination patterns are lacking in tomato. In this study, we developed a novel methodology to use low-coverage resequencing to identify genome-wide recombination patterns and applied this methodology on 60 tomato Recombinant Inbred Lines (RILs). Our methodology identifies polymorphic markers from the low-coverage resequencing population data and utilizes the same data to locate the recombination breakpoints in individuals by using a variable sliding window. We identified 1,445 recombination sites comprising 112 recombination prone regions enriched for AT-rich DNA motifs. Furthermore, the recombination prone regions in tomato preferably occurred in gene promoters over intergenic regions, an observation consistent with Arabidopsis thaliana, Zea mays and Mimulus guttatus. Overall, our cost effective method and findings enhance the understanding of meiotic recombination in tomato and suggest evolutionarily conserved recombination associated genomic features.

Construction of high-resolution recombination maps in Asian seabass

BackgroundA high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.ResultsA high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.ConclusionsThese high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Dissecting meiotic recombination based on tetrad analysis by single-microspore sequencing in maize

Meiotic recombination drives eukaryotic sexual reproduction and the generation of genome diversity. Tetrad analysis, which examines the four chromatids resulting from a single meiosis, is an ideal method to study the mechanisms of homologous recombination. Here we develop a method to isolate the four microspores from a single tetrad in maize for the purpose of whole-genome sequencing. A high-resolution recombination map reveals that crossovers are unevenly distributed across the genome and are more likely to occur in the genic than intergenic regions, especially common in the 5′- and 3′-end regions of annotated genes. The direct detection of genomic exchanges suggests that conversions likely occur in most crossover tracts. Negative crossover interference and weak chromatid interference are observed at the population level. Overall, our findings further our understanding of meiotic recombination with implications for both basic and applied research.

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee.

BackgroundSocial hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.ResultsWe find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.ConclusionsOur data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.

The Drosophila early ovarian transcriptome provides insight to the molecular causes of recombination rate variation across genomes

BackgroundEvidence in yeast indicates that gene expression is correlated with recombination activity and double-strand break (DSB) formation in some hotspots. Studies of nucleosome occupancy in yeast and mice also suggest that open chromatin influences the formation of DSBs. In Drosophila melanogaster, high-resolution recombination maps show an excess of DSBs within annotated transcripts relative to intergenic sequences. The impact of active transcription on recombination landscapes, however, remains unexplored in a multicellular organism. We then investigated the transcription profile during early meiosis in D. melanogaster females to obtain a glimpse at the relevant transcriptional dynamics during DSB formation, and test the specific hypothesis that DSBs preferentially target transcriptionally active genomic regions.ResultsOur study of transcript profiles of early- and late-meiosis using mRNA-seq revealed, 1) significant differences in gene expression, 2) new genes and exons, 3) parent-of-origin effects on transcription in early-meiosis stages, and 4) a nonrandom genomic distribution of transcribed genes. Importantly, genomic regions that are more actively transcribed during early meiosis show higher rates of recombination, and we ruled out DSB preference for genic regions that are not transcribed.ConclusionsOur results provide evidence in a multicellular organism that transcription during the initial phases of meiosis increases the likelihood of DSB and give insight into the molecular determinants of recombination rate variation across the D. melanogaster genome. We propose that a model where variation in gene expression plays a role altering the recombination landscape across the genome could provide a molecular, heritable and plastic mechanism to observed patterns of recombination variation, from the high level of intra-specific variation to the known influence of environmental factors and stress conditions.

Classical Genetics Meets Next-Generation Sequencing: Uncovering a Genome-Wide Recombination Map in Drosophila melanogaster

Classical Genetics Meets Next-Generation Sequencing: Uncovering a Genome-Wide Recombination Map in Drosophila melanogaster

Fine-scale recombination rate differences between sexes, populations and individuals

Meiotic recombinations contribute to genetic diversity by yielding new combinations of alleles. Recently, high-resolution recombination maps were inferred from high-density single-nucleotide polymorphism (SNP) data using linkage disequilibrium (LD) patterns that capture historical recombination events. The use of these maps has been demonstrated by the identification of recombination hotspots and associated motifs, and the discovery that the PRDM9 gene affects the proportion of recombinations occurring at hotspots. However, these maps provide no information about individual or sex differences. Moreover, locus-specific demographic factors like natural selection can bias LD-based estimates of recombination rate. Existing genetic maps based on family data avoid these shortcomings, but their resolution is limited by relatively few meioses and a low density of markers. Here we used genome-wide SNP data from 15,257 parent-offspring pairs to construct the first recombination maps based on directly observed recombinations with a resolution that is effective down to 10 kilobases (kb). Comparing male and female maps reveals that about 15% of hotspots in one sex are specific to that sex. Although male recombinations result in more shuffling of exons within genes, female recombinations generate more new combinations of nearby genes. We discover novel associations between recombination characteristics of individuals and variants in the PRDM9 gene and we identify new recombination hotspots. Comparisons of our maps with two LD-based maps inferred from data of HapMap populations of Utah residents with ancestry from northern and western Europe (CEU) and Yoruba in Ibadan, Nigeria (YRI) reveal population differences previously masked by noise and map differences at regions previously described as targets of natural selection.