The 2.7-Mb major pseudoautosomal region (PAR1) on the short arms of the human X and Y chromosomes plays a critical role in meiotic sex chromosome segregation and male fertility and has been regarded as evolutionarily stable. However, some European Y chromosomes belonging to Y haplogroups (Y-Hgs) R1b and I2a carry an ∼115-kb extension (ePAR [extended PAR]) arising from X-Y non-allelic homologous recombination (NAHR). To investigate the diversity, history, and dynamics of ePAR formation, we screened for its presence, and that of the predicted reciprocal X chromosome deletion, among ∼218,300 46,XY males of the UK Biobank (UKB), a cohort associated with longitudinal clinical data. The UKB incidence of ePAR is ∼0.77%, and that of the deletion is ∼0.02%. We found that Y-Hg I2a sub-lineages accounted for nearly 90% of ePAR cases but, by Y haplotyping and breakpoint sequencing, determined that, in total, there have been at least 18 independent ePAR origins, associated with nine different Y-Hgs. We found examples of ePAR linked to Y-Hg K among men of self-declared Pakistani ancestry and Y-Hg E1, typical of men with African ancestry, showing that ePAR is not restricted to Europeans. ePAR formation is likely random, with high frequencies in some Y-Hgs arising through drift and male-mediated expansions. Sequencing recombination junction fragments identified likely reciprocal events, and the heterogeneity of ePAR and X-deletion junctions highlighted the recurrent nature of the NAHR events. A phenome-wide association study revealed an association between ePAR and elevated levels of circulating IGF-1 as well as musculoskeletal phenotypes.

DNA methylation plays a critical role in the dynamics of gene expression regulation and the development of various disorders. Whole-genome bisulfite sequencing can provide single base resolution of CpG methylation levels and is the “gold standard” for DNA methylation quantification, but it also has a high cost. In contrast, targeted sequencing is optimal when focusing on specific candidate regions, while providing sufficient sequencing depth. Here, we present a targeted bisulfite sequencing approach to study the methylation status of regions of interest. We amplify selected regions from bisulfite-treated DNA and sequence them using Nanopore sequencing. In this work, we applied this workflow to candidate gene promoters for severe premature labor in a Latin American population.We successfully amplified fragments over 1 Kb in length using long PCR conditions for 12 genes that were barcoded per sample and pooled to be sequenced on MinION flow cells. This approach achieved high sequencing depths, ensuring reliable DNAm estimation. We found significant hypomethylation of the MIR155HG gene promoter in severe preterm birth samples, which is concordant with reported gene expression changes.We demonstrate that combining bisulfite DNA treatment with pooled long-read sequencing is a cost- and time-effective method to evaluate DNAm in several targeted regions and several samples in parallel. This study provides proof-of-concept for larger studies, demonstrating the applicability and high scalability of our assay to any locus of interest. Our experience suggests that this approach can be easily transferred to other diagnostic questions.### Competing Interest StatementThe authors have declared no competing interest.

Abstract Objectives The predominance of Bantu languages in sub‐Saharan Africa has sparked a large debate over the processes through which they came to disperse over time and space—the “Bantu expansion.” The overall genetic similarity shown by Bantu‐speaking populations indicates that movement of people occurred too, but the extent of the correlation between genetics, linguistics and geography has been a matter of debate among scholars of different disciplines. In this work, we aim to investigate how genetic, linguistic and geographic distances relate to each other in Bantu‐speaking populations. Methods We analyzed genome‐wide SNP array data from a set of 37 Bantu and non‐Bantu‐speaking populations together with related linguistic and geographic data. Due to the complex demographic relationships resulting from events of admixture in the history of these populations, we develop and implement a method for controlling the signatures of admixture. Results Genetic distances were only minimally correlated with linguistic and geographic distances, possibly as the result of gene flow from neighboring groups into Bantu‐speaking populations. When signatures of admixture are controlled for, the correlation of genetic data with linguistic and geographic distances significantly increases. Discussion The increase of the correlation between linguistic and genetic distances after the signatures of admixture are taken into account is in agreement with a scenario of spatial co‐dispersal of languages and people. Additional specific cultural and demographic dynamics have probably further affected the relationship between language and genetics, which will be necessary to take into account when integrating multidisciplinary data to reconstruct the history of populations.

Birds of prey have suffered persecution for centuries through trapping, shooting, poisoning and theft from the wild to meet the demand from egg collectors and falconers; they were also amongst the earliest beneficiaries of DNA testing in wildlife forensics. Here we report the identification and characterisation of 14 novel tetramer, pentamer and hexamer short tandem repeat (STR) markers which can be typed either by capillary electrophoresis or massively parallel sequencing (MPS) and apply them to historical casework samples involving 49 peregrine falcons, 30 of which were claimed to be the captively bred offspring of nine pairs. The birds were initially tested in 1994 with a multilocus DNA fingerprinting probe, a sex test and eight single-locus minisatellite probes (SLPs) demonstrating that 23 birds were unrelated to the claimed parents. The multilocus and SLP approaches were highly discriminating but extremely time consuming and required microgram quantities of high molecular weight DNA and the use of radioisotopes. The STR markers displayed between 2 and 21 alleles per locus (mean = 7.6), lengths between 140 and 360 bp, and heterozygosities from 0.4 to 0.93. They produced wholly concordant conclusions with similar discrimination power but in a fraction of the time using a hundred-fold less DNA and with standard forensic equipment. Furthermore, eleven of these STRs were amplified in a single reaction and typed using MPS on the Illumina MiSeq platform revealing eight additional alleles (three with variant repeat structures and five solely due to flanking SNPs) across four loci. This approach gave a random match probability of < 1E-9, and a parental pair false inclusion probability of < 1E-5, with a further ten-fold reduction in the amount of DNA required (~3 ng) and the potential to analyse mixed samples. These STRs will be of value in monitoring wild populations of these key indicator species as well as for testing captive breeding claims and establishing a database of captive raptors. They have the potential to resolve complex cases involving trace, mixed and degraded samples from raptor persecution casework representing a significant advance over the previously applied methods.

The human genetic diversity of the Americas has been affected by several events of gene flow that have continued since the colonial era and the Atlantic slave trade. Moreover, multiple waves of migration followed by local admixture occurred in the last two centuries, the impact of which has been largely unexplored. Here, we compiled a genome-wide dataset of ∼12,000 individuals from twelve American countries and ∼6,000 individuals from worldwide populations and applied haplotype-based methods to investigate how historical movements from outside the New World affected (1) the genetic structure, (2) the admixture profile, (3) the demographic history, and (4) sex-biased gene-flow dynamics of the Americas. We revealed a high degree of complexity underlying the genetic contribution of European and African populations in North and South America, from both geographic and temporal perspectives, identifying previously unreported sources related to Italy, the Middle East, and to specific regions of Africa.

Novel raptor multiplex systems for forensics and conservation genetics purposes are developed in this pilot study for non-human forensics. The use of MPS multiplexes in wildlife contexts mirrors the developments made in human applications and leverages the additional information available from sequence data. Novel STRs developed for falcon and accipiter panels are characterised using MPS analysis tools. The newly characterised STRs have shown to be highly variable across and within a range of species.

The human X and Y chromosomes are heteromorphic but share a region of homology at the tips of their short arms, pseudoautosomal region 1 (PAR1), that supports obligate crossover in male meiosis. Although the boundary between pseudoautosomal and sex-specific DNA has traditionally been regarded as conserved among primates, it was recently discovered that the boundary position varies among human males, due to a translocation of ~110 kb from the X to the Y chromosome that creates an extended PAR1 (ePAR). This event has occurred at least twice in human evolution. So far, only limited evidence has been presented to suggest this extension is recombinationally active. Here, we sought direct proof by examining thousands of gametes from each of two ePAR-carrying men, for two subregions chosen on the basis of previously published male X-chromosomal meiotic double-strand break (DSB) maps. Crossover activity comparable to that seen at autosomal hotspots was observed between the X and the ePAR borne on the Y chromosome both at a distal and a proximal site within the 110-kb extension. Other hallmarks of classic recombination hotspots included evidence of transmission distortion and GC-biased gene conversion. We observed good correspondence between the male DSB clusters and historical recombination activity of this region in the X chromosomes of females, as ascertained from linkage disequilibrium analysis; this suggests that this region is similarly primed for crossover in both male and female germlines, although sex-specific differences may also exist. Extensive resequencing and inference of ePAR haplotypes, placed in the framework of the Y phylogeny as ascertained by both Y microsatellites and single nucleotide polymorphisms, allowed us to estimate a minimum rate of crossover over the entire ePAR region of 6-fold greater than genome average, comparable with pedigree estimates of PAR1 activity generally. We conclude ePAR very likely contributes to the critical crossover function of PAR1.

Recombination plays a fundamental role in meiosis. Non-exchange gene conversion (non-crossover, NCO) may facilitate homologue pairing, while reciprocal crossover (CO) physically connects homologues so they orientate appropriately on the meiotic spindle. In males, X-Y homologous pairing and exchange occurs within the two pseudoautosomal regions (PARs) together comprising <5% of the human sex chromosomes. Successful meiosis depends on an obligatory CO within PAR1, while the nature and role of exchange within PAR2 is unclear. Here, we describe the identification and characterization of a typical ~1 kb wide recombination hotspot within PAR2. We find that both COs and NCOs are strongly modulated in trans by the presumed chromatin remodelling protein PRDM9, and in cis by a single nucleotide polymorphism (SNP) located at the hotspot centre that appears to influence recombination initiation and which causes biased gene conversion in SNP heterozygotes. This, the largest survey to date of human NCOs reveals for the first time substantial inter-individual variation in the NCO:CO ratio. Although the extent of biased transmission at the central marker in COs is similar across men, it is highly variable among NCO recombinants. This suggests that cis-effects are mediated not only through recombination initiation frequencies varying between haplotypes but also through subsequent processing, with the potential to significantly intensify meiotic drive of hotspot-suppressing alleles. The NCO:CO ratio and extent of transmission distortion among NCOs appear to be inter-related, suggesting the existence of two NCO pathways in humans.

Alec Jeffreys and colleagues report that variation at PRDM9 in humans influences sperm recombination hotspot activity, independent of a consensus binding motif, as well as meiotic genome instability. PRDM9 has recently been identified as a likely trans regulator of meiotic recombination hot spots in humans and mice1,2,3. PRDM9 contains a zinc finger array that, in humans, can recognize a short sequence motif associated with hot spots4, with binding to this motif possibly triggering hot-spot activity via chromatin remodeling5. We now report that human genetic variation at the PRDM9 locus has a strong effect on sperm hot-spot activity, even at hot spots lacking the sequence motif. Subtle changes within the zinc finger array can create hot-spot nonactivating or enhancing variants and can even trigger the appearance of a new hot spot, suggesting that PRDM9 is a major global regulator of hot spots in humans. Variation at the PRDM9 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.

Traditional methods for surveying meiotic recombination in humans are limited to pedigree and linkage disequilibrium analyses. We have developed assays that allow the direct detection of crossover and gene conversion molecules in batches of sperm DNA. To date, we have characterized 26 recombination hotspots by allele-specific PCR and selectively amplified recombinant DNA molecules from these regions. These analyses have revealed that meiotic crossover hotspots in humans are highly localized and flanked by DNA segments where recombination is suppressed. The centers of crossover hotspots are also active in noncrossover recombination, displaying short conversion tracts.

Mutation at most human minisatellites is driven by complex interallelic processes that give rise to a high degree of length polymorphism and internal structural variation. MSY1, the only highly variable minisatellite on the non-recombining region of the Y chromosome, is constitutively haploid and therefore precluded from interallelic interactions, yet maintains high diversity in both length and structure. To investigate the basis of its mutation processes, an unbiased structural analysis of >500 single-molecule MSY1 PCR products from matched sperm and blood samples from a single donor was undertaken. The overall mutation frequencies in sperm and blood DNAs were not significantly different, at 2.68% and 1.88%, respectively. Sperm DNA showed significantly more length mutants than blood DNA, with mutants in both tissues involving small-scale (1–3 repeat units in a 77 repeat progenitor allele) increases or decreases in repeat block lengths, with no gain or loss bias. Isometric mutations altering structure but not length were found in both tissues, and involved either the apparent shift of a boundary between repeat unit blocks (a ‘boundary switch’) or the conversion of a repeat within a block to a different repeat type (‘modular structure’ mutant). There was a significant excess of boundary switch mutants and deficit of modular structure mutants in sperm. A comparison of mutant structures with phylogenetically matched alleles in population samples showed that alleles with structures resembling the blood mutants were unlikely to arise in populations. Mutation seems likely to involve gene conversion via synthesis-dependent strand annealing, and the blood-sperm differences may reflect more relaxed constraint on sister chromatid alignment in blood.

Nat. Genet. 36, 151–156 (2004). In the version of this paper initially published online, several lines were missing from the bottom of each page of the Supplementary Note. This mistake has been corrected online.

Meiotic recombination plays a key role in the maintenance of sequence diversity in the human genome. However, little is known about the fine-scale distribution and processes of recombination in human chromosomes, or how these impact on patterns of human diversity. We have therefore developed sperm typing systems that allow human recombination to be analysed at very high resolution. The emerging picture is that human crossovers are far from randomly distributed but instead are targeted into very narrow hot spots that can profoundly influence patterns of haplotype diversity in the human genome. These hot spots provide fundamental information on processes of human crossover and gene conversion, as well as evidence that they can violate basic rules of Mendelian inheritance.

Detecting rare sequence variants in genomic DNA is central to the analysis of de novo mutation and recombination events and the detection of rare pathological mutations in mixed cell populations. Current PCR techniques suffer from noise that limits detection to variants present at a frequency of at least 10(-4)-10(-5) per cell. We now describe an alternative approach that recovers genomic DNA molecules containing a known single-nucleotide variant by hybridization selection using a biotinylated allele-specific oligonucleotide, followed by hybrid capture on streptavidin-coated paramagnetic beads and subsequent analysis by PCR. This technique of DNA enrichment by allele-specific hybridization (DEASH) is fast, effective for all tested single-nucleotide polymorphisms (SNPs), and can recover large (>10 kb) single-stranded molecules. A single round of DEASH is effective in separating haplotypes from genomic DNA and can not only readily detect and validate DNA molecules containing a single base change at a frequency of 10(-5) per cell, but can also place these changes within the context of an extended haplotype. This technique offers a new approach to the analysis of mutation and recombination, and has the potential to detect very rare de novo base substitutions.

Crossover between the human sex chromosomes during male meiosis is restricted to the terminal pseudoautosomal pairing regions. An obligatory exchange occurs in PAR1, an Xp/Yp pseudoautosomal region of 2.6 Mb, which creates a male-specific recombination 'hot domain' with a recombination rate that is about 20 times higher than the genome average1,2,3. Low-resolution analysis of PAR1 suggests that crossovers are distributed fairly randomly4. By contrast, linkage disequilibrium (LD)5 and sperm crossover analyses6 indicate that crossovers in autosomal regions tend to cluster into 'hot spots' of 1–2 kb that lie between islands of disequilibrium of tens to hundreds of kilobases7. To determine whether at high resolution this autosomal pattern also applies to PAR1, we have examined linkage disequilibrium over an interval of 43 kb around the gene SHOX8. Here we show that in northern European populations, disequilibrium decays rapidly with physical distance, which is consistent with this interval of PAR1 being recombinationally active in male meiosis. Analysis of a subregion of 9.9 kb in sperm shows, however, that crossovers are not distributed randomly, but instead cluster into an intense recombination hot spot that is very similar in morphology to autosomal hot spots. Thus, PAR1 crossover activity may be influenced by male-specific hot spots that are highly suitable for characterization by sperm DNA analysis.

Variation of brood sex ratio was studied in a Finnish population of Eurasian Kestrels Falco tinnunculus breeding in an unpredictably variable environment. From those young that survived until 2-4 weeks of age, blood was collected and their sex determined from polymorphic DNA profiles produced by hybridisation with a human minisatellite probe. The sex ratio was male-biased during a year of food (vole) scarcity. Furthermore, in broods without mortality, contrasting seasonal trends in sex ratios emerged. In this subsample, the proportion of males increased with later laying date during years of low and moderate food supply, whereas the opposite was true in a year of relatively high food supply. These trends may indicate circumstances that Favour the raising of different sex. The proportion of males in the brood was also negatively correlated with body condition of both male and female parents, reflecting an adaptive condition-dependent sex-ratio adjustment, or alternatively the inability of the parents to meet the requirements of the more energetically expensive female offspring. We discuss the limitations that unpredictable conditions during brood raising can impose on adaptive sex-ratio manipulation, particularly in species with sexual size dimorphism and consequent differences in the cost of raising the two sexes.

We report the results of a collaborative study aimed at developing reliable, direct assays for mutation in human cells. The project used common lymphoblastoid cell lines, both with and without mutagen treatment, as a shared resource to validate the development of new molecular methods for the detection of low-level mutations in the presence of a large excess of normal alleles. As the "gold standard, " hprt mutation frequencies were also measured on the same samples. The methods under development included i) the restriction site mutation (RSM) assay, in which mutations lead to the destruction of a restriction site; ii) minisatellite length-change mutation, in which mutations lead to alleles containing new numbers of tandem repeat units; iii) loss of heterozygosity for HLA epitopes, in which antibodies can be used to direct selection for mutant cells; iv) multiple fluorescence-based long linker arm nucleotides assay (mf-LLA) technology, for the detection of substitutional mutations; v) detection of alterations in the TP53 locus using a (CA) array as the target for the screening; and vi) PCR analysis of lymphocytes for the presence of the BCL2 t(14:18) translocation. The relative merits of these molecular methods are discussed, and a comparison made with more "traditional" methods.

We have examined compound haplotypes from a highly informative region of human chromosome 16, in which information from the rapid evolution of a highly unstable minisatellite is integrated with data on the longer-term evolution of this segment from 10 flanking substitutional polymorphisms. Combined with sequence data from non-human primates, analysis of relationships between these compound haplotypes allows the reconstruction of a rooted network of the evolutionary pathways between them. Most relationships can be explained via simple substitutional mutations, although the origins of some haplotypes involve recurrent events at a hotspot for substitutional mutation and/or gene conversion. For compound haplotypes including the minisatellite array, the network found in a range of world-wide populations constitutes a highly informative data set for the analysis of population history (437 different compound haplotypes were discriminated among 658 studied). Since the mutation rates and processes of the minisatellite array are known from direct studies, ages for individual lineages have been estimated using associated minisatellite diversity. These analyses suggest that the higher information content and sampling depth of these compound haplotypes may allow more precise calibration of lineage ages than is possible using coalescent analysis of DNA sequence. Using this method we have dated the oldest Eurasian lineage as 52,000-66,000 years and the oldest European specific lineage as 37,600-56,200 years.

Screening pedigrees for inherited minisatellite length changes provides an efficient means of monitoring repeat DNA instability but has given rise to apparently contradictory results regarding the effects of radiation on the human germline. To explore this further in individuals with known radiation doses and to potentially gain information on the timing of mutation induction, we have used an extremely sensitive single molecule approach to quantify the frequencies of mutation at the hypervariable minisatellites B6.7 and CEB1 in the sperm of three seminoma patients following hemipelvic radiotherapy. Scattered radiation doses to the testicles were monitored and pre-treatment sperm DNA was compared with sperm derived from irradiated pre-meiotic, meiotic and post-meiotic cells. We show no evidence for mutation induction in any of the patients and discuss this finding in the context of previous population studies using minisatellites as reporter systems, one of which provided evidence for radiation-induced germline mutation.