Summary The nonrandom distribution of meiotic recombination shapes patterns of inheritance and genome evolution, but chromosomal features governing this distribution are poorly understood. Formation of the DNA double-strand breaks (DSBs) that initiate recombination results in accumulation of Spo11 protein covalently bound to small DNA fragments. We show here that sequencing these fragments provides a genome-wide DSB map of unprecedented resolution and sensitivity. We use this map to explore the influence of large-scale chromosome structures, chromatin, transcription factors, and local sequence composition on DSB distributions. Our analysis supports the view that the recombination terrain is molded by combinatorial and hierarchical interaction of factors that work on widely different size scales. Mechanistic aspects of DSB formation and early processing steps are also uncovered. This map illuminates the occurrence of DSBs in repetitive DNA elements, repair of which can lead to chromosomal rearrangements. We discuss implications for evolutionary dynamics of recombination hotspots.

SUMMARY Heritability and genome stability are shaped by meiotic recombination, which is initiated via hundreds of DNA double-strand breaks (DSBs). The distribution of DSBs throughout the genome is not random, but mechanisms molding this landscape remain poorly understood. Here we exploit genome-wide maps of mouse DSBs at unprecedented nucleotide resolution to uncover previously invisible spatial features of recombination. At fine scale, we reveal a stereotyped hotspot structure—DSBs occur within narrow zones between methylated nucleosomes—and identify relationships between SPO11, chromatin, and the histone methyltransferase PRDM9. At large scale, DSB formation is suppressed on non-homologous portions of the sex chromosomes via the DSB-responsive kinase ATM, which also shapes the autosomal DSB landscape at multiple size scales. We also provide a genome-wide analysis of exonucleolytic DSB resection lengths and elucidate spatial relationships between DSBs and recombination products. Our results paint a comprehensive picture of features governing successive steps in mammalian meiotic recombination.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a “cosmopolitan” tagging approach to capture the genetic diversity across ∼2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Genome-wide association studies (GWASs) have identified many SNPs underlying variations in plasma-lipid levels. We explore whether additional loci associated with plasma-lipid phenotypes, such as high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), and triglycerides (TGs), can be identified by a dense gene-centric approach. Our meta-analysis of 32 studies in 66,240 individuals of European ancestry was based on the custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) covering ∼2,000 candidate genes. SNP-lipid associations were replicated either in a cohort comprising an additional 24,736 samples or within the Global Lipid Genetic Consortium. We identified four, six, ten, and four unreported SNPs in established lipid genes for HDL-C, LDL-C, TC, and TGs, respectively. We also identified several lipid-related SNPs in previously unreported genes: DGAT2, HCAR2, GPIHBP1, PPARG, and FTO for HDL-C; SOCS3, APOH, SPTY2D1, BRCA2, and VLDLR for LDL-C; SOCS3, UGT1A1, BRCA2, UBE3B, FCGR2A, CHUK, and INSIG2 for TC; and SERPINF2, C4B, GCK, GATA4, INSR, and LPAL2 for TGs. The proportion of explained phenotypic variance in the subset of studies providing individual-level data was 9.9% for HDL-C, 9.5% for LDL-C, 10.3% for TC, and 8.0% for TGs. This large meta-analysis of lipid phenotypes with the use of a dense gene-centric approach identified multiple SNPs not previously described in established lipid genes and several previously unknown loci. The explained phenotypic variance from this approach was comparable to that from a meta-analysis of GWAS data, suggesting that a focused genotyping approach can further increase the understanding of heritability of plasma lipids.

Chromatin state and oncogenic competence Although specific DNA mutations can lead to tumor generation, they are not transforming in all cellular contexts. This may be due to the intrinsic transcriptional program present in the cell of origin. Using zebrafish and human pluripotent stem cell cancer models, Baggiolini et al . report that neural crest cells and melanoblasts (precursors to melanocytes) are susceptible to specific mutation of the BRAF gene, whereas melanocytes are relatively resistant (see the Perspective by Vredevoogd and Peeper). The competent cells display higher levels of chromatin factors such as the protein ATAD2 compared with the less competent ones. ATAD2 forms a complex with the neural crest transcription factor SOX10 and establishes a chromatin state that makes them permissive to BRAF mutagenesis. These data indicate that developmental chromatin programs are a determinant of how cells respond to DNA mutations. —BAP

Height is a classic complex trait with common variants in a growing list of genes known to contribute to the phenotype. Using a genecentric genotyping array targeted toward cardiovascular-related loci, comprising 49,320 SNPs across approximately 2000 loci, we evaluated the association of common and uncommon SNPs with adult height in 114,223 individuals from 47 studies and six ethnicities. A total of 64 loci contained a SNP associated with height at array-wide significance (p < 2.4 × 10(-6)), with 42 loci surpassing the conventional genome-wide significance threshold (p < 5 × 10(-8)). Common variants with minor allele frequencies greater than 5% were observed to be associated with height in 37 previously reported loci. In individuals of European ancestry, uncommon SNPs in IL11 and SMAD3, which would not be genotyped with the use of standard genome-wide genotyping arrays, were strongly associated with height (p < 3 × 10(-11)). Conditional analysis within associated regions revealed five additional variants associated with height independent of lead SNPs within the locus, suggesting allelic heterogeneity. Although underpowered to replicate findings from individuals of European ancestry, the direction of effect of associated variants was largely consistent in African American, South Asian, and Hispanic populations. Overall, we show that dense coverage of genes for uncommon SNPs, coupled with large-scale meta-analysis, can successfully identify additional variants associated with a common complex trait.

Lineage plasticity is implicated in treatment resistance in multiple cancers. In lung adenocarcinomas (LUAD) amenable to targeted therapy, transformation to small cell lung cancer (SCLC) is a recognized resistance mechanism. Defining molecular mechanisms of neuroendocrine (NE) transformation in lung cancer has been limited by a paucity of pre/posttransformation clinical samples. Detailed genomic, epigenomic, transcriptomic, and protein characterization of combined LUAD/SCLC tumors, as well as pre/posttransformation samples, supports that NE transformation is primarily driven by transcriptional reprogramming rather than mutational events. We identify genomic contexts in which NE transformation is favored, including frequent loss of the 3p chromosome arm. We observed enhanced expression of genes involved in the PRC2 complex and PI3K/AKT and NOTCH pathways. Pharmacologic inhibition of the PI3K/AKT pathway delayed tumor growth and NE transformation in an EGFR-mutant patient-derived xenograft model. Our findings define a novel landscape of potential drivers and therapeutic vulnerabilities of NE transformation in lung cancer. Significance: The difficulty in collection of transformation samples has precluded the performance of molecular analyses, and thus little is known about the lineage plasticity mechanisms leading to LUAD-to-SCLC transformation. Here, we describe biological pathways dysregulated upon transformation and identify potential predictors and potential therapeutic vulnerabilities of NE transformation in the lung. See related commentary by Meador and Lovly, p. 2962. This article is highlighted in the In This Issue feature, p. 2945

Spo11, which makes DNA double-strand breaks (DSBs) essential for meiotic recombination, has long been recalcitrant to biochemical study. We provide molecular analysis of S. cerevisiae Spo11 purified with partners Rec102, Rec104 and Ski8. Rec102 and Rec104 jointly resemble the B subunit of archaeal Topoisomerase VI, with Rec104 occupying a position similar to the Top6B GHKL-type ATPase domain. Unexpectedly, the Spo11 complex is monomeric (1:1:1:1 stoichiometry), consistent with dimerization controlling DSB formation. Reconstitution of DNA binding reveals topoisomerase-like preferences for duplex-duplex junctions and bent DNA. Spo11 also binds noncovalently but with high affinity to DNA ends mimicking cleavage products, suggesting a mechanism to cap DSB ends. Mutations that reduce DNA binding in vitro attenuate DSB formation, alter DSB processing, and reshape the DSB landscape in vivo . Our data reveal structural and functional similarities between the Spo11 core complex and Topo VI, but also highlight differences reflecting their distinct biological roles.