A common sequence motif associated with recombination hot spots and genome instability in humans

Genetic analysis of genome-scale recombination rate evolution in house mice

The rate of meiotic recombination varies markedly between species and among individuals. Classical genetic experiments demonstrated a heritable component to population variation in recombination rate, and specific sequence variants that contribute to recombination rate differences between individuals have recently been identified. Despite these advances, the genetic basis of species divergence in recombination rate remains unexplored. Using a cytological assay that allows direct in situ imaging of recombination events in spermatocytes, we report a large (∼30%) difference in global recombination rate between males of two closely related house mouse subspecies (Mus musculus musculus and M. m. castaneus). To characterize the genetic basis of this recombination rate divergence, we generated an F2 panel of inter-subspecific hybrid males (n = 276) from an intercross between wild-derived inbred strains CAST/EiJ (M. m. castaneus) and PWD/PhJ (M. m. musculus). We uncover considerable heritable variation for recombination rate among males from this mapping population. Much of the F2 variance for recombination rate and a substantial portion of the difference in recombination rate between the parental strains is explained by eight moderate- to large-effect quantitative trait loci, including two transgressive loci on the X chromosome. In contrast to the rapid evolution observed in males, female CAST/EiJ and PWD/PhJ animals show minimal divergence in recombination rate (∼5%). The existence of loci on the X chromosome suggests a genetic mechanism to explain this male-biased evolution. Our results provide an initial map of the genetic changes underlying subspecies differences in genome-scale recombination rate and underscore the power of the house mouse system for understanding the evolution of this trait.

Homologous recombination is an indispensable feature of the mammalian meiotic program and an important mechanism for creating genetic diversity. Despite its central significance, recombination rates vary markedly between species and among individuals. Although recent studies have begun to unravel the genetic basis of recombination rate variation within populations, the genetic mechanisms of species divergence in recombination rate remain poorly characterized. In this study, we show that two closely related house mouse subspecies differ in their genomic recombination rates by ∼30%, providing an excellent model system for studying evolutionary divergence in this trait. Using quantitative genetic methods, we identify eight genomic regions that contribute to divergence in global recombination rate between these subspecies, including large effect loci and multiple loci on the X-chromosome. Our study uncovers novel genomic loci contributing to species divergence in global recombination rate and offers simple genetic explanations for rapid phenotypic divergence in this trait.

Mammalian recombination hot spots: properties, control and evolution

Recombination, together with mutation, generates the raw material of evolution, is essential for reproduction and lies at the heart of all genetic analysis. Recent advances in our ability to construct genome-scale, high-resolution recombination maps and new molecular techniques for analysing recombination products have substantially furthered our understanding of this important biological phenomenon in humans and mice: from describing the properties of recombination hot spots in male and female meiosis to the recombination landscape along chromosomes. This progress has been accompanied by the identification of trans-acting systems that regulate the location and relative activity of individual hot spots.

High recombination rates and hotspots in a Plasmodium falciparum genetic cross

Background

The human malaria parasite Plasmodium falciparum survives pressures from the host immune system and antimalarial drugs by modifying its genome. Genetic recombination and nucleotide substitution are the two major mechanisms that the parasite employs to generate genome diversity. A better understanding of these mechanisms may provide important information for studying parasite evolution, immune evasion and drug resistance.

Results

Here, we used a high-density tiling array to estimate the genetic recombination rate among 32 progeny of a P. falciparum genetic cross (7G8 × GB4). We detected 638 recombination events and constructed a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb per centimorgan and identified 54 candidate recombination hotspots. Similar to centromeres in other organisms, the sequences of P. falciparum centromeres are found in chromosome regions largely devoid of recombination activity. Motifs enriched in hotspots were also identified, including a 12-bp G/C-rich motif with 3-bp periodicity that may interact with a protein containing 11 predicted zinc finger arrays.

Conclusions

These results show that the P. falciparum genome has a high recombination rate, although it also follows the overall rule of meiosis in eukaryotes with an average of approximately one crossover per chromosome per meiosis. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination rate observed. The lack of recombination activity in centromeric regions is consistent with the observations of reduced recombination near the centromeres of other organisms.

Genome-wide analysis reveals novel molecular features of mouse recombination hotspots

Meiotic recombination predominantly occurs at discrete genomic loci called recombination hotspots, but the features defining these areas are still largely unknown (reviewed in refs 1-5). To allow a comprehensive analysis of hotspot-associated DNA and chromatin characteristics, we developed a direct molecular approach for mapping meiotic DNA double-strand breaks that initiate recombination. Here we present the genome-wide distribution of recombination initiation sites in the mouse genome. Hotspot centres are mapped with approximately 200-nucleotide precision, which allows analysis of the fine structural details of the preferred recombination sites. We determine that hotspots share a centrally distributed consensus motif, possess a nucleotide skew that changes polarity at the centres of hotspots and have an intrinsic preference to be occupied by a nucleosome. Furthermore, we find that the vast majority of recombination initiation sites in mouse males are associated with testis-specific trimethylation of lysine 4 on histone H3 that is distinct from histone H3 lysine 4 trimethylation marks associated with transcription. The recombination map presented here has been derived from a homogeneous mouse population with a defined genetic background and therefore lends itself to extensive future experimental exploration. We note that the mapping technique developed here does not depend on the availability of genetic markers and hence can be easily adapted to other species with complex genomes. Our findings uncover several fundamental features of mammalian recombination hotspots and underline the power of the new recombination map for future studies of genetic recombination, genome stability and evolution.

What are the genomic drivers of the rapid evolution of PRDM9?

Mammalian Prdm9 has been proposed to be a key determinant of the positioning of chromosome double-strand breaks during meiosis, a contributor to speciation processes, and the most rapidly evolving gene in human, and other animal, genomes. Prdm9 genes often exhibit substantial variation in their numbers of encoded zinc fingers (ZFs), not only between closely related species but also among individuals of a species. The near-identity of these ZF sequences appears to render them very unstable in copy number. The rare sequence differences, however, cluster within ZF sites that determine the DNA-binding specificity of PRDM9, and these substitutions are frequently positively selected. Here, possible drivers of the rapid evolution of Prdm9 are discussed, including selection for efficient pairing of homologous chromosomes or for recombination of deleterious linked alleles, and selection against depletion of recombination hotspots or against disease-associated genome rearrangement.

Ontogenetic de novo copy number variations (CNVs) as a source of genetic individuality: studies on two families with MZD twins for schizophrenia

Genetic individuality is the foundation of personalized medicine, yet its determinants are currently poorly understood. One issue is the difference between monozygotic twins that are assumed identical and have been extensively used in genetic studies for decades [1]. Here, we report genome-wide alterations in two nuclear families each with a pair of monozygotic twins discordant for schizophrenia evaluated by the Affymetrix 6.0 human SNP array. The data analysis includes characterization of copy number variations (CNVs) and single nucleotide polymorphism (SNPs). The results have identified genomic differences between twin pairs and a set of new provisional schizophrenia genes. Samples were found to have between 35 and 65 CNVs per individual. The majority of CNVs (∼80%) represented gains. In addition, ∼10% of the CNVs were de novo (not present in parents), of these, 30% arose during parental meiosis and 70% arose during developmental mitosis. We also observed SNPs in the twins that were absent from both parents. These constituted 0.12% of all SNPs seen in the twins. In 65% of cases these SNPs arose during meiosis compared to 35% during mitosis. The developmental mitotic origin of most CNVs that may lead to MZ twin discordance may also cause tissue differences within individuals during a single pregnancy and generate a high frequency of mosaics in the population. The results argue for enduring genome-wide changes during cellular transmission, often ignored in most genetic analyses.

Have a break: determinants of meiotic DNA double strand break (DSB) formation and processing in plants

Meiosis is an essential process for sexually reproducing organisms, leading to the formation of specialized generative cells. This review intends to highlight current knowledge of early events during meiosis derived from various model organisms, including plants. It will particularly focus on cis- and trans-requirements of meiotic DNA double strand break (DSB) formation, a hallmark event during meiosis and a prerequisite for recombination of genetic traits. Proteins involved in DSB formation in different organisms, emphasizing the known factors from plants, will be introduced and their functions outlined. Recent technical advances in DSB detection and meiotic recombination analysis will be reviewed, as these new tools now allow analysis of early meiotic recombination in plants with incredible accuracy. To anticipate future directions in plant meiosis research, unpublished results will be included wherever possible.

Copy number gain at Xp22.31 includes complex duplication rearrangements and recurrent triplications

Genomic instability is a feature of the human Xp22.31 region wherein deletions are associated with X-linked ichthyosis, mental retardation and attention deficit hyperactivity disorder. A putative homologous recombination hotspot motif is enriched in low copy repeats that mediate recurrent deletion at this locus. To date, few efforts have focused on copy number gain at Xp22.31. However, clinical testing revealed a high incidence of duplication of Xp22.31 in subjects ascertained and referred with neurobehavioral phenotypes. We systematically studied 61 unrelated subjects with rearrangements revealing gain in copy number, using multiple molecular assays. We detected not only the anticipated recurrent and simple nonrecurrent duplications, but also unexpectedly identified recurrent triplications and other complex rearrangements. Breakpoint analyses enabled us to surmise the mechanisms for many of these rearrangements. The clinical significance of the recurrent duplications and triplications were assessed using different approaches. We cannot find any evidence to support pathogenicity of the Xp22.31 duplication. However, our data suggest that the Xp22.31 duplication may serve as a risk factor for abnormal phenotypes. Our findings highlight the need for more robust Xp22.31 triplication detection in that such further gain may be more penetrant than the duplications. Our findings reveal the distribution of different mechanisms for genomic duplication rearrangements at a given locus, and provide insights into aspects of strand exchange events between paralogous sequences in the human genome.

The Drosophila zinc finger protein trade embargo is required for double strand break formation in meiosis

Homologous recombination in meiosis is initiated by the programmed induction of double strand breaks (DSBs). Although the Drosophila Spo11 ortholog Mei-W68 is required for the induction of DSBs during meiotic prophase, only one other protein (Mei-P22) has been shown to be required for Mei-W68 to exert this function. We show here that the chromatin-associated protein Trade Embargo (Trem), a C2H2 zinc finger protein, is required to localize Mei-P22 to discrete foci on meiotic chromosomes, and thus to promote the formation of DSBs, making Trem the earliest known function in the process of DSB formation in Drosophila oocytes. We speculate that Trem may act by either directing the binding of Mei-P22 to preferred sites of DSB formation or by altering chromatin structure in a manner that allows Mei-P22 to form foci.

Substitution patterns are under different influences in primates and rodents

There are large-scale variations of the GC-content along mammalian chromosomes that have been called isochore structures. Primates and rodents have different isochore structures, which suggests that these lineages exhibit different modes of GC-content evolution. It has been shown that, in the human lineage, GC-biased gene conversion (gBGC), a neutral process associated with meiotic recombination, acts on GC-content evolution by influencing A or T to G or C substitution rates. We computed genome-wide substitution patterns in the mouse lineage from multiple alignments and compared them with substitution patterns in the human lineage. We found that in the mouse lineage, gBGC is active but weaker than in the human lineage and that male-specific recombination better predicts GC-content evolution than female-specific recombination. Furthermore, we were able to show that G or C to A or T substitution rates are predicted by a combination of different factors in both lineages. A or T to G or C substitution rates are most strongly predicted by meiotic recombination in the human lineage but by CpG odds ratio (the observed CpG frequency normalized by the expected CpG frequency) in the mouse lineage, suggesting that substitution patterns are under different influences in primates and rodents.

The choice in meiosis – defining the factors that influence crossover or non-crossover formation

Meiotic crossovers are essential for ensuring correct chromosome segregation as well as for creating new combinations of alleles for natural selection to take place. During meiosis, excess meiotic double-strand breaks (DSBs) are generated; a subset of these breaks are repaired to form crossovers, whereas the remainder are repaired as non-crossovers. What determines where meiotic DSBs are created and whether a crossover or non-crossover will be formed at any particular DSB remains largely unclear. Nevertheless, several recent papers have revealed important insights into the factors that control the decision between crossover and non-crossover formation in meiosis, including DNA elements that determine the positioning of meiotic DSBs, and the generation and processing of recombination intermediates. In this review, we focus on the factors that influence DSB positioning, the proteins required for the formation of recombination intermediates and how the processing of these structures generates either a crossover or non-crossover in various organisms. A discussion of crossover interference, assurance and homeostasis, which influence crossing over on a chromosome-wide and genome-wide scale – in addition to current models for the generation of interference – is also included. This Commentary aims to highlight recent advances in our understanding of the factors that promote or prevent meiotic crossing over.

[What defines the genetic map? The specification of meiotic recombination sites]

During meiosis, homologous reciprocal recombination events or crossing-over determine the genetic map and are known not to be randomly distributed in the genome. Recent studies in yeasts and mammals reveal some key features of the molecular mechanism involved in this distribution. Through different molecular processes, specific histone post-translational modifications are induced at specific genomic sites, called hotspots, where initiation of meiotic recombination takes place. These sites are some transcription promoters in S. cerevisiae or binding sites for transcription factors in S. pombe, where chromatin modifiers are recruited. In mammals, the sites are DNA sequences recognized by the PRDM9 protein which has the ability both to bind DNA and to induce the trimethylation of the lysine 4 of histone H3. The properties of the chromatin at these sites, and potentially the binding of additional factors, allow the recruitment of proteins involved in the formation of DNA double strand breaks that initiate meiotic recombination.

A de novo unequal cross-over mutation between CYP11B1 and CYP11B2 genes causes familial hyperaldosteronism type I

Familial hyperaldosteronism type I (FH-I) is an autosomal dominant disorder caused by an unequal cross-over of the gene encoding steroid 11β-hydroxylase (CYP11B1) and aldosterone synthase (CYP11B2), giving rise to a chimeric CYP11B1/CYP11B2 gene that displays aldosterone synthase activity regulated by ACTH instead of angiotensin II.

AIM

To report an unprecedented case of a de novo unequal crossover mutation between CYP11B1 and CYP11B2 genes causing FH-I.

PATIENTS AND METHODS

The index case is a 45-yr-old Chilean male diagnosed with primary aldosteronism (PA). All family members were also studied: his biological parents, 1 brother, 6 sisters, 2 daughters, and 1 son. Plasma renin activity, serum aldosterone, and its ratio were measured in all patients. Genetic analyses were performed using long-extension PCR (XL-PCR), DNA sequencing and Southern blot methods.

RESULTS

PA was diagnosed for the index case, 1 of his daughters, his son but not for his parents or siblings. XLPCR and Southern blotting demonstrated the presence of the chimeric CYP11B1/CYP11B2 gene solely in PA-affected subjects, suggesting a case of a de novo mutation. Sequence analysis showed the unequal cross-over CYP11B1/CYP11B2 at intron 2 (c.2600-273 CYP11B2). We also identified a polymorphism at the same intron (c.2600-145C>A CYP11B2) in the genome of the index case's father.

CONCLUSION

We describe an unprecedented case of unequal cross-over mutation for the chimeric CYP11B1/CYP11B2 gene causing FH-I, which may be linked to a polymorphism in the index case's father germ line.

Intrachromosomal mitotic nonallelic homologous recombination is the major molecular mechanism underlying type-2 NF1 deletions

Nonallelic homologous recombination (NAHR) is responsible for the recurrent rearrangements that give rise to genomic disorders. Although meiotic NAHR has been investigated in multiple contexts, much less is known about mitotic NAHR despite its importance for tumorigenesis. Because type-2 NF1 microdeletions frequently result from mitotic NAHR, they represent a good model in which to investigate the features of mitotic NAHR. We have used microsatellite analysis and SNP arrays to distinguish between the various alternative recombinational possibilities, thereby ascertaining that 17 of 18 type-2 NF1 deletions, with breakpoints in the SUZ12 gene and its highly homologous pseudogene, originated via intrachromosomal recombination. This high proportion of intrachromosomal NAHR causing somatic type-2 NF1 deletions contrasts with the interchromosomal origin of germline type-1 NF1 microdeletions, whose breakpoints are located within the NF1-REPs (low-copy repeats located adjacent to the SUZ12 sequences). Further, meiotic NAHR causing type-1 NF1 deletions occurs within recombination hotspots characterized by high GC-content and DNA duplex stability, whereas the type-2 breakpoints associated with the mitotic NAHR events investigated here do not cluster within hotspots and are located within regions of significantly lower GC-content and DNA stability. Our findings therefore point to fundamental mechanistic differences between the determinants of mitotic and meiotic NAHR.

Characterization of the past and current duplication activities in the human 22q11.2 region

Background

Segmental duplications (SDs) on 22q11.2 (LCR22), serve as substrates for meiotic non-allelic homologous recombination (NAHR) events resulting in several clinically significant genomic disorders.

Results

To understand the duplication activity leading to the complicated SD structure of this region, we have applied the A-Bruijn graph algorithm to decompose the 22q11.2 SDs to 523 fundamental duplication sequences, termed subunits. Cross-species syntenic analysis of primate genomes demonstrates that many of these LCR22 subunits emerged very recently, especially those implicated in human genomic disorders. Some subunits have expanded more actively than others, and young Alu SINEs, are associated much more frequently with duplicated sequences that have undergone active expansion, confirming their role in mediating recombination events. Many copy number variations (CNVs) exist on 22q11.2, some flanked by SDs. Interestingly, two chromosome breakpoints for 13 CNVs (mean length 65 kb) are located in paralogous subunits, providing direct evidence that SD subunits could contribute to CNV formation. Sequence analysis of PACs or BACs identified extra CNVs, specifically, 10 insertions and 18 deletions within 22q11.2; four were more than 10 kb in size and most contained young AluYs at their breakpoints.

Conclusions

Our study indicates that AluYs are implicated in the past and current duplication events, and moreover suggests that DNA rearrangements in 22q11.2 genomic disorders perhaps do not occur randomly but involve both actively expanded duplication subunits and Alu elements.

Variable tandem repeats accelerate evolution of coding and regulatory sequences

Genotype-to-phenotype mapping commonly focuses on two major classes of mutations: single nucleotide polymorphisms (SNPs) and copy number variation (CNV). Here, we discuss an underestimated third class of genotypic variation: changes in microsatellite and minisatellite repeats. Such tandem repeats (TRs) are ubiquitous, unstable genomic elements that have historically been designated as nonfunctional "junk DNA" and are therefore mostly ignored in comparative genomics. However, as many as 10% to 20% of eukaryotic genes and promoters contain an unstable repeat tract. Mutations in these repeats often have fascinating phenotypic consequences. For example, changes in unstable repeats located in or near human genes can lead to neurodegenerative diseases such as Huntington disease. Apart from their role in disease, variable repeats also confer useful phenotypic variability, including cell surface variability, plasticity in skeletal morphology, and tuning of the circadian rhythm. As such, TRs combine characteristics of genetic and epigenetic changes that may facilitate organismal evolvability.

Structural variation of the human genome: mechanisms, assays, and role in male infertility

Genomic disorders are defined as diseases caused by rearrangements of the genome incited by a genomic architecture that conveys instability. Y-chromosome related dysfunctions such as male infertility are frequently associated with gross DNA rearrangements resulting from its peculiar genomic architecture. The Y-chromosome has evolved into a highly specialized chromosome to perform male functions, mainly spermatogenesis. Direct and inverted repeats, some of them palindromes with highly identical nucleotide sequences that can form DNA cruciform structures, characterize the genomic structure of the Y-chromosome long arm. Some particular Y chromosome genomic deletions can cause spermatogenic failure likely because of removal of one or more transcriptional units with a potential role in spermatogenesis. We describe mechanisms underlying the formation of human genomic rearrangements on autosomes and review Y-chromosome deletions associated with male infertility.

A human genome structural variation sequencing resource reveals insights into mutational mechanisms

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

Initiation of meiotic recombination in mammals

Meiotic recombination is initiated by the induction of programmed DNA double strand breaks (DSBs). DSB repair promotes homologous interactions and pairing and leads to the formation of crossovers (COs), which are required for the proper reductional segregation at the first meiotic division. In mammals, several hundred DSBs are generated at the beginning of meiotic prophase by the catalytic activity of SPO11. Currently it is not well understood how the frequency and timing of DSB formation and their localization are regulated. Several approaches in humans and mice have provided an extensive description of the localization of initiation events based on CO mapping, leading to the identification and characterization of preferred sites (hotspots) of initiation. This review presents the current knowledge about the proteins known to be involved in this process, the sites where initiation takes place, and the factors that control hotspot localization.

Variation in crossover rates across a 3-Mb contig of bread wheat (Triticum aestivum) reveals the presence of a meiotic recombination hotspot

In bread wheat (Triticum aestivum L.), initial studies using deletion lines indicated that crossover (CO) events occur mainly in the telomeric regions of the chromosomes with a possible correlation with the presence of genes. However, little is known about the distribution of COs at the sequence level. To investigate this, we studied in detail the pattern of COs along a contig of 3.110 Mb using two F2 segregating populations (Chinese Spring × Renan (F2-CsRe) and Chinese Spring × Courtot (F2-CsCt)) each containing ~2,000 individuals. The availability of the sequence of the contig from Cs enabled the development of 318 markers among which 23 co-dominant polymorphic markers (11 SSRs and 12 SNPs) were selected for CO distribution analyses. The distribution of CO events was not homogeneous throughout the contig, ranging from 0.05 to 2.77 cM/Mb, but was conserved between the two populations despite very different contig recombination rate averages (0.82 cM/Mb in F2-CsRe vs 0.35 cM/Mb in F2-CsCt). The CO frequency was correlated with the percentage of coding sequence in Cs and with the polymorphism rate between Cs and Re or Ct in both populations, indicating an impact of these two factors on CO distribution. At a finer scale, COs were found in a region covering 2.38 kb, spanning a gene coding for a glycosyl transferase (Hga3), suggesting the presence of a CO hotspot. A non-crossover event covering at least 453 bp was also identified in the same interval. From these results, we can conclude that gene content could be one of the factors driving recombination in bread wheat.

The Red Queen theory of recombination hotspots

Recombination hotspots are small chromosomal regions, where meiotic crossover events happen with high frequency. Recombination is initiated by a double-strand break (DSB) that requires the intervention of the molecular repair mechanism. The DSB repair mechanism may result in the exchange of homologous chromosomes (crossover) and the conversion of the allelic sequence that breaks into the one that does not break (biased gene conversion). Biased gene conversion results in a transmission advantage for the allele that does not break, thus preventing recombination and rendering recombination hotspots transient. How is it possible that recombination hotspots persist over evolutionary time (maintaining the average chromosomal crossover rate) when they are self-destructive? This fundamental question is known as the recombination hotspot paradox and has attracted much attention in recent years. Yet, that attention has not translated into a fully satisfactory answer. No existing model adequately explains all aspects of the recombination hotspot paradox. Here, we formulate an intragenomic conflict model resulting in Red Queen dynamics that fully accounts for all empirical observations regarding the molecular mechanisms of recombination hotspots, the nonrandom targeting of the recombination machinery to hotspots and the evolutionary dynamics of hotspot turnover.

A map of human genome variation from population-scale sequencing

The 1000 Genomes Project aims to provide a deep characterization of human genome sequence variation as a foundation for investigating the relationship between genotype and phenotype. Here we present results of the pilot phase of the project, designed to develop and compare different strategies for genome-wide sequencing with high-throughput platforms. We undertook three projects: low-coverage whole-genome sequencing of 179 individuals from four populations; high-coverage sequencing of two mother-father-child trios; and exon-targeted sequencing of 697 individuals from seven populations. We describe the location, allele frequency and local haplotype structure of approximately 15 million single nucleotide polymorphisms, 1 million short insertions and deletions, and 20,000 structural variants, most of which were previously undescribed. We show that, because we have catalogued the vast majority of common variation, over 95% of the currently accessible variants found in any individual are present in this data set. On average, each person is found to carry approximately 250 to 300 loss-of-function variants in annotated genes and 50 to 100 variants previously implicated in inherited disorders. We demonstrate how these results can be used to inform association and functional studies. From the two trios, we directly estimate the rate of de novo germline base substitution mutations to be approximately 10(-8) per base pair per generation. We explore the data with regard to signatures of natural selection, and identify a marked reduction of genetic variation in the neighbourhood of genes, due to selection at linked sites. These methods and public data will support the next phase of human genetic research.

Patterns of recombination activity on mouse chromosome 11 revealed by high resolution mapping

The success of high resolution genetic mapping of disease predisposition and quantitative trait loci in humans and experimental animals depends on the positions of key crossover events around the gene of interest. In mammals, the majority of recombination occurs at highly delimited 1-2 kb long sites known as recombination hotspots, whose locations and activities are distributed unevenly along the chromosomes and are tightly regulated in a sex specific manner. The factors determining the location of hotspots started to emerge with the finding of PRDM9 as a major hotspot regulator in mammals, however, additional factors modulating hotspot activity and sex specificity are yet to be defined. To address this limitation, we have collected and mapped the locations of 4829 crossover events occurring on mouse chromosome 11 in 5858 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. This chromosome was chosen for its medium size and high gene density and provided a comparison with our previous analysis of recombination on the longest mouse chromosome 1. Crossovers were mapped to an average resolution of 127 kb, and thirteen hotspots were mapped to <8 kb. Most crossovers occurred in a small number of the most active hotspots. Females had higher recombination rate than males as a consequence of differences in crossover interference and regional variation of sex specific rates along the chromosome. Comparison with chromosome 1 showed that recombination events tend to be positioned in similar fashion along the centromere-telomere axis but independently of the local gene density. It appears that mammalian recombination is regulated on at least three levels, chromosome-wide, regional, and at individual hotspots, and these regulation levels are influenced by sex and genetic background but not by gene content.

xnd-1 regulates the global recombination landscape in Caenorhabditis elegans

Meiotic crossover (CO) recombination establishes physical linkages between homologous chromosomes that are required for their proper segregation into developing gametes, and promotes genetic diversity by shuffling genetic material between parental chromosomes. COs require the formation of double strand breaks (DSBs) to create the substrate for strand exchange. DSBs occur in small intervals called hotspots and significant variation in hotspot usage exists between and among individuals. This variation is thought to reflect differences in sequence identity and chromatin structure, DNA topology and/ or chromosome domain organization. Chromosomes show different frequencies of nondisjunction (NDJ), reflecting inherent differences in meiotic crossover control, yet the underlying basis of these differences remains elusive. Here we show that a novel chromatin factor, X non-disjunction factor 1 (xnd-1), is responsible for the global distribution of COs in C. elegans. xnd-1 is also required for formation of double-strand breaks (DSBs) on the X, but surprisingly XND-1 protein is autosomally enriched. We show that xnd-1 functions independently of genes required for X chromosome-specific gene silencing, revealing a novel pathway that distinguishes the X from autosomes in the germ line, and further show that xnd-1 exerts its effects on COs, at least in part, by modulating levels of H2A lysine 5 acetylation.

Important characteristics of sequence-specific recombination hotspots in Schizosaccharomyces pombe

In many organisms, meiotic recombination occurs preferentially at a limited number of sites in the genome known as hotspots. In the fission yeast Schizosaccharomyces pombe, simple sequence motifs determine the location of at least some, and possibly most or all, hotspots. Recently, we showed that a large number of different sequences can create hotspots. Among those sequences we identified some recurring motifs that fell into at least five distinct families, including the well-characterized CRE family of hotspots. Here we report the essential sequence for activity of two of the novel hotspots, the oligo-C and CCAAT hotspots, and identify associated trans-acting factors required for hotspot activity. The oligo-C hotspot requires a unique 8-bp sequence, CCCCGCAC, though hotspot activity is also significantly affected by adjacent nucleotides. The CCAAT hotspot requires a more complex and degenerate sequence, including the originally identified seven nucleotide CCAATCA sequence at its core. We identified transcription factors, the CCAAT-binding factor (CBF) and Rst2, which are required specifically for activity of the CCAAT hotspots and oligo-C hotspots, respectively. Each of these factors binds to its respective motifs in vitro. However, unlike CRE, the sequence required for hotspot activity is larger than the sequence required for binding, suggesting the involvement of additional factors.

Informativeness of a novel multiallelic marker-set comprising an F8 intron 21 and three tightly linked loci for haemophilia A carriership analysis

The extraordinary heterogeneous nature of the deleterious mutations in the F8 gene that lead to functional deficiency of clotting factor VIII in haemophilia A makes routine direct mutation profiling difficult. When direct mutation analysis cannot be performed or a causative/candidate mutation is not found, a second-line approach to track the defective F8 gene within at-risk families is linkage genetic analysis with, tried-and-tested, F8-intragenic and/or extragenic non-recombining multiallelic short tandem repeats (STR). Although several typing STR loci within and around F8 have been described, there is need for improving assessment, because the combined informativeness of available assays rarely reaches 100%. Here, we characterized a newly identified 0.28 cM-resolution marker-set, consisting of a dinucleotide STR located on F8 intron 21 (F8Int21; [AC](n)) and three extragenic tetranucleotide STR located on GAB3 intron 1 (GAB3Int1; [TAAA](n)) and TMLHE intron 1 (TMLHEInt1.1; [GAAA](n) and TMLHEInt1.3; [ATTC](n)). Heterozygosity rates determined in 100 unrelated females ranged from 0.25 (GAB3Int1) to 0.63 (F8Int21). The set rendered a combined informativeness of 0.91 for at least one marker and 0.60 for a minimum of two loci, with at least one F8-intragenic. Multiallelic interlocus non-random association analysis revealed that GAB3Int1 is not in significant gametic disequilibrium (GD) with F8Int21, F8Int9.2, TMLHEInt1.3 or TMLHEInt1.1. Gametic disequilibrium breakdown attests historical recombination between GAB3Int1 and the F8 gene. Through computational analysis of reference assembly sequence data, we note in the GD breakdown region and in the F8 gene a higher than average density of the 13-mer CCNCCNTNNCCNC consensus motif, commonly associated with recombination hotspots.

Detecting sequence polymorphisms associated with meiotic recombination hotspots in the human genome

Background

Meiotic recombination events tend to cluster into narrow spans of a few kilobases long, called recombination hotspots. Such hotspots are not conserved between human and chimpanzee and vary between different human ethnic groups. At the same time, recombination hotspots are heritable. Previous studies showed instances where differences in recombination rate could be associated with sequence polymorphisms.

Results

In this work we developed a novel computational approach, LDsplit, to perform a large-scale association study of recombination hotspots with genetic polymorphisms. LDsplit was able to correctly predict the association between the FG11 SNP and the DNA2 hotspot observed by sperm typing. Extensive simulation demonstrated the accuracy of LDsplit under various conditions. Applying LDsplit to human chromosome 6, we found that for a significant fraction of hotspots, there is an association between variations in intensity of historical recombination and sequence polymorphisms. From flanking regions of the SNPs output by LDsplit we identified a conserved 11-mer motif GGNGGNAGGGG, whose complement partially matches 13-mer CCNCCNTNNCCNC, a critical motif for the regulation of recombination hotspots.

Conclusions

Our result suggests that computational approaches based on historical recombination events are likely to be more powerful than previously anticipated. The putative associations we identified may be a promising step toward uncovering the mechanisms of recombination hotspots.

Structures and molecular mechanisms for common 15q13.3 microduplications involving CHRNA7: benign or pathological?

We have investigated four approximately 1.6-Mb microduplications and 55 smaller 350-680-kb microduplications at 15q13.2-q13.3 involving the CHRNA7 gene that were detected by clinical microarray analysis. Applying high-resolution array-CGH, we mapped all 118 chromosomal breakpoints of these microduplications. We also sequenced 26 small microduplication breakpoints that were clustering at hotspots of nonallelic homologous recombination (NAHR). All four large microduplications likely arose by NAHR between BP4 and BP5 LCRs, and 54 small microduplications arose by NAHR between two CHRNA7-LCR copies. We identified two classes of approximately 1.6-Mb microduplications and five classes of small microduplications differing in duplication size, and show that they duplicate the entire CHRNA7. We propose that size differences among small microduplications result from preexisting heterogeneity of the common BP4-BP5 inversion. Clinical data and family histories of 11 patients with small microduplications involving CHRNA7 suggest that these microduplications might be associated with developmental delay/mental retardation, muscular hypotonia, and a variety of neuropsychiatric disorders. However, we conclude that these microduplications and their associated potential for increased dosage of the CHRNA7-encoded alpha 7 subunit of nicotinic acetylcholine receptors are of uncertain clinical significance at present. Nevertheless, if they prove to have a pathological effects, their high frequency could make them a common risk factor for many neurobehavioral disorders.

Genetic and evolutionary correlates of fine-scale recombination rate variation in Drosophila persimilis

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed "recombinational neighborhoods," where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5-7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.

A novel third type of recurrent NF1 microdeletion mediated by nonallelic homologous recombination between LRRC37B-containing low-copy repeats in 17q11.2

Large microdeletions encompassing the neurofibromatosis type-1 (NF1) gene and its flanking regions at 17q11.2 belong to the group of genomic disorders caused by aberrant recombination between segmental duplications. The most common NF1 microdeletions (type-1) span 1.4-Mb and have breakpoints located within NF1-REPs A and C, low-copy repeats (LCRs) containing LRRC37-core duplicons. We have identified a novel type of recurrent NF1 deletion mediated by nonallelic homologous recombination (NAHR) between the highly homologous NF1-REPs B and C. The breakpoints of these approximately 1.0-Mb ("type-3") NF1 deletions were characterized at the DNA sequence level in three unrelated patients. Recombination regions, spanning 275, 180, and 109-bp, respectively, were identified within the LRRC37B-P paralogues of NF1-REPs B and C, and were found to contain sequences capable of non-B DNA formation. Both LCRs contain LRRC37-core duplicons, abundant and highly dynamic sequences in the human genome. NAHR between LRRC37-containing LCRs at 17q21.31 is known to have mediated the 970-kb polymorphic inversions of the MAPT-locus that occurred independently in different primate species, but also underlies the syndromes associated with recurrent 17q21.31 microdeletions and reciprocal microduplications. The novel NF1 microdeletions reported here provide further evidence for the unusually high recombinogenic potential of LRRC37-containing LCRs in the human genome.

Efficient population assignment and outlier detection in human populations using biallelic markers chosen by principal component-based rankings

Whole-genome studies of genetic variation are now performed routinely and have accelerated the identification of disease-associated allelic variants, positive selection, recombination, and structural variation. However, these studies are sensitive to the presence of outlier data from individuals of different ancestry than the rest of the sample. Currently, the most common method of excluding outlier individuals is to collect a population sample and exclude outliers after genome-wide data have been collected. Here we show that a small collection of 20-27 polymorphic Alu insertions, selected using a principal component-based method with genetic ancestry estimates, may be used to easily assign Africans, East Asians, and Europeans to their population of origin. In addition, we show that samples from a geographically and genetically intermediate population (in our study, samples from India) can be identified within the original sample of Africans, East Asians, and Europeans. Finally, we show that outlier individuals from neighboring geographic regions (in our study, Yemen and sub-Saharan Africa) can be identified. These results will be of value in preselection of samples for more in-depth analysis as well as customized identification of maximally informative polymorphic markers for regional studies.

Meiosis: making a break for it

The perpetuation of most eukaryotic species requires differentiation of pluripotent progenitors into egg and sperm and subsequent fusion of these gametes to form a new zygote. Meiosis is a distinguishing feature of gamete formation as it leads to the twofold reduction in chromosome number thereby maintaining ploidy across generations. This process increases offspring diversity through the random segregation of chromosomes and the exchange of genetic material between homologous parental chromosomes, known as meiotic crossover recombination. These exchanges require the establishment of unique and dynamic chromatin configurations that facilitate cohesion, homolog pairing, synapsis, double strand break formation and repair. The precise orchestration of these events is critical for gamete survival as demonstrated by the majority of human aneuploidies that can be traced to defects in the first meiotic division (Hassold T, Hall H, Hunt P: The origin of human aneuploidy: where we have been, where we are going. Hum Mol Genet 2007, 16 Spec No. 2:R203-R208.). This review will focus on recent advances in our understanding of key meiotic events and how coordination of these events is occurring.

PRDM9 variation strongly influences recombination hot-spot activity and meiotic instability in humans

PRDM9 has recently been identified as a likely trans-regulator of meiotic recombination hot spots in humans and mice1-3. The protein contains a zinc finger array that in humans can recognise a short sequence motif associated with hot spots4, with binding to this motif possibly triggering hot-spot activity via chromatin remodelling5. We now show that variation in the zinc finger array in humans has a profound effect on sperm hot-spot activity, even at hot spots lacking the sequence motif. Very subtle changes within the array can create hot-spot non-activating and enhancing alleles, and even trigger the appearance of a new hot spot. PRDM9 thus appears to be the preeminent global regulator of hot spots in humans. Variation at this locus also influences aspects of genome instability, specifically a megabase-scale rearrangement underlying two genomic disorders6 as well as minisatellite instability7, implicating PRDM9 as a risk factor for some pathological genome rearrangements.

Low-complexity regions in Plasmodium falciparum: missing links in the evolution of an extreme genome

Over the past decade, attempts to explain the unusual size and prevalence of low-complexity regions (LCRs) in the proteins of the human malaria parasite Plasmodium falciparum have used both neutral and adaptive models. This past research has offered conflicting explanations for LCR characteristics and their role in, and influence on, the evolution of genome structure. Here we show that P. falciparum LCRs (PfLCRs) are not a single phenomenon, but rather consist of at least three distinct types of sequence, and this heterogeneity is the source of the conflict in the literature. Using molecular and population genetics, we show that these families of PfLCRs are evolving by different mechanisms. One of these families, named here the HighGC family, is of particular interest because these LCRs act as recombination hotspots, both in genes under positive selection for high levels of diversity which can be created by recombination (antigens) and those likely to be evolving neutrally or under negative selection (metabolic enzymes). We discuss how the discovery of these distinct species of PfLCRs helps to resolve previous contradictory studies on LCRs in malaria and contributes to our understanding of the evolution of the of the parasite's unusual genome.

Genomic signatures of germline gene expression

Transcribed regions in the human genome differ from adjacent intergenic regions in transposable element density, crossover rates, and asymmetric substitution and sequence composition patterns. We tested whether these differences reflect selection or are instead a byproduct of germline transcription, using publicly available gene expression data from a variety of germline and somatic tissues. Crossover rate shows a strong negative correlation with gene expression in meiotic tissues, suggesting that crossover is inhibited by transcription. Strand-biased composition (G+T content) and A → G versus T → C substitution asymmetry are both positively correlated with germline gene expression. We find no evidence for a strand bias in allele frequency data, implying that the substitution asymmetry reflects a mutation rather than a fixation bias. The density of transposable elements is positively correlated with germline expression, suggesting that such elements preferentially insert into regions that are actively transcribed. For each of the features examined, our analyses favor a nonselective explanation for the observed trends and point to the role of germline gene expression in shaping the mammalian genome.

What drives recombination hotspots to repeat DNA in humans?

Recombination between homologous, but non-allelic, stretches of DNA such as gene families, segmental duplications and repeat elements is an important source of mutation. In humans, recent studies have identified short DNA motifs that both determine the location of 40 per cent of meiotic cross-over hotspots and are significantly enriched at the breakpoints of recurrent non-allelic homologous recombination (NAHR) syndromes. Unexpectedly, the most highly penetrant form of the motif occurs on the background of an inactive repeat element family (THE1 elements) and the motif also has strong recombinogenic activity on currently active element families including Alu and LINE2 elements. Analysis of genetic variation among members of these repeat families indicates an important role for NAHR in their evolution. Given the potential for double-strand breaks within repeat DNA to cause pathological rearrangement, the association between repeats and hotspots is surprising. Here we consider possible explanations for why selection acting against NAHR has not eliminated hotspots from repeat DNA including mechanistic constraints, possible benefits to repeat DNA from recruiting hotspots and rapid evolution of the recombination machinery. I suggest that rapid evolution of hotspot motifs may, surprisingly, tend to favour sequences present in repeat DNA and outline the data required to differentiate between hypotheses.

Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.

Nucleotide polymorphism of the TNF gene cluster in six Chinese populations

DNA variants in a 31-kb region of the human major histocompatibility complex, encompassing the tumor necrosis factor (TNF) gene cluster, were surveyed by direct sequencing of 283 unrelated individuals from six Chinese populations. A total of 273 polymorphic sites were identified, with nearly half of them novel. We observed an excess of rare variants and negative values of selection tests of the region, implying either that these populations experienced a historical expansion or that the surveyed region was subjected to natural selection. Different characteristics of the sequence variation in the six populations outline the genetic differentiation between Northern and Southern Chinese populations. The distributions of recombination rates are similar among all the populations, with variation in the magnitude and/or in the fine location of hot spots. Tag single-nucleotide polymorphisms (SNPs) selected from HapMap (Phase II) CHB data accounted for an average of 64% of common SNPs from the six Chinese populations. We also observed a limited transferability of tag SNPs between Chinese populations on the 31-kb region with an excess of untaggable SNPs and ragged linkage disequilibrium blocks. It suggested that the design and interpretation of future association studies should be more cautious, and that a resequencing approach may refine tag SNP selection on Chinese-specific disease mapping.

PRDM9 sticks its zinc fingers into recombination hotspots and between species

Meiotic recombination events typically cluster within narrow regions of the genome termed hotspots. A series of recent papers reveals that PRDM9, a C2H2-type zinc-finger protein with histone H3 lysine 4 methyltransferase activity, plays a major role in the specification of hotspots. The zinc fingers that contact DNA in a sequence-dependent manner evolve rapidly and are under positive selection, leading to differences in the location of recombination hotspots as well as hybrid sterility.

Mutation spectrum revealed by breakpoint sequencing of human germline CNVs

Precisely characterizing the breakpoints of copy number variants (CNVs) is crucial for assessing their functional impact. However, fewer than 10% of known germline CNVs have been mapped to the single-nucleotide level. We characterized the sequence breakpoints from a dataset of all CNVs detected in three unrelated individuals in previous array-based CNV discovery experiments. We used targeted hybridization-based DNA capture and 454 sequencing to sequence 324 CNV breakpoints, including 315 deletions. We observed two major breakpoint signatures: 70% of the deletion breakpoints have 1-30 bp of microhomology, whereas 33% of deletion breakpoints contain 1-367 bp of inserted sequence. The co-occurrence of microhomology and inserted sequence is low (10%), suggesting that there are at least two different mutational mechanisms. Approximately 5% of the breakpoints represent more complex rearrangements, including local microinversions, suggesting a replication-based strand switching mechanism. Despite a rich literature on DNA repair processes, reconstruction of the molecular events generating each of these mutations is not yet possible.

Discrete DNA sites regulate global distribution of meiotic recombination

Homologous recombination is induced to high levels in meiosis, is initiated by Spo11-catalyzed DNA double-strand breaks (DSBs) and is clustered at hotspots that regulate its positioning in the genome. Recombination is required for proper chromosome segregation in meiosis and defects in its frequency or positioning cause chromosome mis-segregation and, consequently, congenital birth defects such as Down's syndrome. Therefore, elucidating how meiotic recombination is positioned is of fundamental and biomedical interest. Our integration of historical and contemporary advances in the field, plus the re-analysis of published microarray data on the genome-wide distribution of recombination supports a unifying model for such regulation. We posit that discrete DNA sequence motifs position and regulate essentially all recombination across the genome, in much the same way that DNA sites position and regulate transcription. Moreover, we illustrate the use of overlapping mechanisms for the regulation of transcription and meiotic recombination. Bound transcription factors induce histone modifications that position recombination at hotspots.

The recombination landscape of the zebra finch Taeniopygia guttata genome

Understanding the causes and consequences of variation in the rate of recombination is essential since this parameter is considered to affect levels of genetic variability, the efficacy of selection, and the design of association and linkage mapping studies. However, there is limited knowledge about the factors governing recombination rate variation. We genotyped 1920 single nucleotide polymorphisms in a multigeneration pedigree of more than 1000 zebra finches (Taeniopygia guttata) to develop a genetic linkage map, and then we used these map data together with the recently available draft genome sequence of the zebra finch to estimate recombination rates in 1 Mb intervals across the genome. The average zebra finch recombination rate (1.5 cM/Mb) is higher than in humans, but significantly lower than in chicken. The local rates of recombination in chicken and zebra finch were only weakly correlated, demonstrating evolutionary turnover of the recombination landscape in birds. The distribution of recombination events was heavily biased toward ends of chromosomes, with a stronger telomere effect than so far seen in any organism. In fact, the recombination rate was as low as 0.1 cM/Mb in intervals up to 100 Mb long in the middle of the larger chromosomes. We found a positive correlation between recombination rate and GC content, as well as GC-rich sequence motifs. Levels of linkage disequilibrium (LD) were significantly higher in regions of low recombination, showing that heterogeneity in recombination rates have left a footprint on the genomic landscape of LD in zebra finch populations.

Copy number variation, chromosome rearrangement, and their association with recombination during avian evolution

Chromosomal rearrangements and copy number variants (CNVs) play key roles in genome evolution and genetic disease; however, the molecular mechanisms underlying these types of structural genomic variation are not fully understood. The availability of complete genome sequences for two bird species, the chicken and the zebra finch, provides, for the first time, an ideal opportunity to analyze the relationship between structural genomic variation (chromosomal and CNV) and recombination on a genome-wide level. The aims of this study were therefore threefold: (1) to combine bioinformatics, physical mapping to produce comprehensive comparative maps of the genomes of chicken and zebra finch. In so doing, this allowed the identification of evolutionary chromosomal rearrangements distinguishing them. The previously reported interchromosomal conservation of synteny was confirmed, but a larger than expected number of intrachromosomal rearrangements were reported; (2) to hybridize zebra finch genomic DNA to a chicken tiling path microarray and identify CNVs in the zebra finch genome relative to chicken; 32 interspecific CNVs were identified; and (3) to test the hypothesis that there is an association between CNV, chromosomal rearrangements, and recombination by correlating data from (1) and (2) with recombination rate data from a high-resolution genetic linkage map of the zebra finch. We found a highly significant association of both chromosomal rearrangements and CNVs with elevated recombination rates. The results thus provide support for the notion of recombination-based processes playing a major role in avian genome evolution.

Genome-wide mapping and assembly of structural variant breakpoints in the mouse genome

Structural variation (SV) is a rich source of genetic diversity in mammals, but due to the challenges associated with mapping SV in complex genomes, basic questions regarding their genomic distribution and mechanistic origins remain unanswered. We have developed an algorithm (HYDRA) to localize SV breakpoints by paired-end mapping, and a general approach for the genome-wide assembly and interpretation of breakpoint sequences. We applied these methods to two inbred mouse strains: C57BL/6J and DBA/2J. We demonstrate that HYDRA accurately maps diverse classes of SV, including those involving repetitive elements such as transposons and segmental duplications; however, our analysis of the C57BL/6J reference strain shows that incomplete reference genome assemblies are a major source of noise. We report 7196 SVs between the two strains, more than two-thirds of which are due to transposon insertions. Of the remainder, 59% are deletions (relative to the reference), 26% are insertions of unlinked DNA, 9% are tandem duplications, and 6% are inversions. To investigate the origins of SV, we characterized 3316 breakpoint sequences at single-nucleotide resolution. We find that approximately 16% of non-transposon SVs have complex breakpoint patterns consistent with template switching during DNA replication or repair, and that this process appears to preferentially generate certain classes of complex variants. Moreover, we find that SVs are significantly enriched in regions of segmental duplication, but that this effect is largely independent of DNA sequence homology and thus cannot be explained by non-allelic homologous recombination (NAHR) alone. This result suggests that the genetic instability of such regions is often the cause rather than the consequence of duplicated genomic architecture.