Ross M. Fraser

A General Approach for Haplotype Phasing across the Full Spectrum of Relatedness

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally ‘unrelated’ individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Micrococcal Nuclease Does Not Substantially Bias Nucleosome Mapping

We have mapped sequence-directed nucleosome positioning on genomic DNA molecules using high-throughput sequencing. Chromatins, prepared by reconstitution with either chicken or frog histones, were separately digested to mononucleosomes using either micrococcal nuclease (MNase) or caspase-activated DNase (CAD). Both enzymes preferentially cleave internucleosomal (linker) DNA, although they do so by markedly different mechanisms. MNase has hitherto been very widely used to map nucleosomes, although concerns have been raised over its potential to introduce bias. Having identified the locations and quantified the strength of both the chicken or frog histone octamer binding sites on each DNA, the results obtained with the two enzymes were compared using a variety of criteria. Both enzymes displayed sequence specificity in their preferred cleavage sites, although the nature of this selectivity was distinct for the two enzymes. In addition, nucleosomes produced by CAD nuclease are 8–10 bp longer than those produced with MNase, with the CAD cleavage sites tending to be 4–5 bp further out from the nucleosomal dyad than the corresponding MNase cleavage sites. Despite these notable differences in cleavage behaviour, the two nucleases identified essentially equivalent patterns of nucleosome positioning sites on each of the DNAs tested, an observation that was independent of the histone type. These results indicate that biases in nucleosome positioning data collected using MNase are, under our conditions, not significant.

High-Resolution Mapping of Sequence-Directed Nucleosome Positioning on Genomic DNA

We have mapped in vitro nucleosome positioning on the sheep beta-lactoglobulin gene using high-throughput sequencing to characterise the DNA sequences recovered from reconstituted nucleosomes. This methodology surpasses previous approaches for coverage, accuracy and resolution and, most importantly, offers a simple yet rapid and relatively inexpensive method to characterise genomic DNA sequences in terms of nucleosome positioning capacity. We demonstrate an unambiguous correspondence between in vitro and in vivo nucleosome positioning around the promoter of the gene; identify discrete, sequence-specific nucleosomal structures above the level of the canonical core particle-a feature that has implications for regulatory protein access and higher-order chromatin packing; and reveal new insights into the involvement of periodically organised dinucleotide sequence motifs of the type GG and CC and not AA and TT, as determinants of nucleosome positioning-an observation that supports the idea that the core histone octamer can exploit different patterns of sequence organisation, or structural potential, in the DNA to bring about nucleosome positioning.

Resolving the ancestry of Austronesian-speaking populations

There are two very different interpretations of the prehistory of Island Southeast Asia (ISEA), with genetic evidence invoked in support of both. The “out-of-Taiwan” model proposes a major Late Holocene expansion of Neolithic Austronesian speakers from Taiwan. An alternative, proposing that Late Glacial/postglacial sea-level rises triggered largely autochthonous dispersals, accounts for some otherwise enigmatic genetic patterns, but fails to explain the Austronesian language dispersal. Combining mitochondrial DNA (mtDNA), Y-chromosome and genome-wide data, we performed the most comprehensive analysis of the region to date, obtaining highly consistent results across all three systems and allowing us to reconcile the models. We infer a primarily common ancestry for Taiwan/ISEA populations established before the Neolithic, but also detected clear signals of two minor Late Holocene migrations, probably representing Neolithic input from both Mainland Southeast Asia and South China, via Taiwan. This latter may therefore have mediated the Austronesian language dispersal, implying small-scale migration and language shift rather than large-scale expansion.

Nucleosome Positioning Signals in the DNA Sequence of the Human and Mouse H19 Imprinting Control Regions

We have investigated the sequences of the mouse and human H19 imprinting control regions (ICRs) to see whether they contain nucleosome positioning information pertinent to their function as a methylation-regulated chromatin boundary. Positioning signals were identified by an in vitro approach that employs reconstituted chromatin to comprehensively describe the contribution of the DNA to the most basic, underlying level of chromatin structure. Signals in the DNA sequence of both ICRs directed nucleosomes to flank and encompass the short conserved sequences that constitute the binding sites for the zinc finger protein CTCF, an essential mediator of insulator activity. The repeat structure of the human ICR presented a conserved array of strong positioning signals that would preferentially flank these CTCF binding sites with positioned nucleosomes, a chromatin structure that would tend to maintain their accessibility. Conversely, all four CTCF binding sites in the mouse sequence were located close to the centre of positioning signals that were stronger than those in their flanks; these binding sites might therefore be expected to be more readily incorporated into positioned nucleosomes. We found that CpG methylation did not effect widespread repositioning of nucleosomes on either ICR, indicating that allelic methylation patterns were unlikely to establish allele-specific chromatin structures for H19 by operating directly upon the underlying DNA-histone interactions; instead, epigenetic modulation of ICR chromatin structure is likely to be mediated principally at higher levels of control. DNA methylation did, however, both promote and inhibit nucleosome positioning at several sites in both ICRs and substantially negated one of the strongest nucleosome positioning signals in the human sequence, observations that underline the fact that this epigenetic modification can, nevertheless, directly and decisively modulate core histone-DNA interactions within the nucleosome.

The Image of the Barbarian in Early India

The concept of the barbarian in early India arises out of the curious situation of the arrival of Indo-Aryan-speaking nomadic pastoralists in northern India who came into contact with the indigenous population (possibly the remnants of the urban civilization of the Indus) and regarded them as barbarians. The earliest distinction made by the Aryan speakers was a linguistic distinction and, to a smaller extent, a physical distinction. The Indo-Aryan speakers spoke Sanskrit whereas the indigenous peoples probably spoke Dravidian and Munda. However the distinction was not one of binary opposition—in fact it admitted to many nuances and degrees of variation, hence the complication of trying to trace the history of the concept. The distinction was rarely clearly manifest and based either on language, ethnic origins or culture. Political status, ritual status and economic power, all tended to blur the contours of the distinction. Added to this has been the confusion introduced by those who tend to identify language with race and who thereby see all speakers of Sanskrit as members of that nineteenth-century myth, the Aryan race.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown—as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

Local Exome Sequences Facilitate Imputation of Less Common Variants and Increase Power of Genome Wide Association Studies

The analysis of less common variants in genome-wide association studies promises to elucidate complex trait genetics but is hampered by low power to reliably detect association. We show that addition of population-specific exome sequence data to global reference data allows more accurate imputation, particularly of less common SNPs (minor allele frequency 1–10%) in two very different European populations. The imputation improvement corresponds to an increase in effective sample size of 28–38%, for SNPs with a minor allele frequency in the range 1–3%.

A Comparison of In Vitro Nucleosome Positioning Mapped with Chicken, Frog and a Variety of Yeast Core Histones

Using high-throughput sequencing, we have mapped sequence-directed nucleosome positioning in vitro on four plasmid DNAs containing DNA fragments derived from the genomes of sheep, drosophila, human and yeast. Chromatins were prepared by reconstitution using chicken, frog and yeast core histones. We also assembled yeast chromatin in which histone H3 was replaced by the centromere-specific histone variant, Cse4. The positions occupied by recombinant frog and native chicken histones were found to be very similar. In contrast, nucleosomes containing the canonical yeast octamer or, in particular, the Cse4 octamer were assembled at distinct populations of locations, a property that was more apparent on particular genomic DNA fragments. The factors that may contribute to this variation in nucleosome positioning and the implications of the behavior are discussed.

Phenotypic consequences of variation across the aldosterone synthase and 11-beta hydroxylase locus in a hypertensive cohort: data from the MRC BRIGHT Study.

Background  Aldosterone is an important cardiovascular hormone; 15% of hypertensive subjects have alteration in aldosterone regulation, defined by a raised ratio of aldosterone to renin (ARR). Studies of the aldosterone synthase gene (CYP11B2) have focused on a single nucleotide polymorphism in the 5′promoter region (–344 C/T). In normotensive subjects, the T allele associates with raised levels of the 11‐deoxysteroids, deoxycorticosterone and 11‐deoxycortisol which are substrates for 11β‐hydroxylase, encoded by the adjacent and homologous gene, CYP11B1. We have speculated that this altered 11β‐hydroxylase efficiency leads to increased ACTH drive to the adrenal gland to maintain cortisol production and reported herein the association between the –344 C/T single nucleotide polymorphism (SNP) and adrenal steroid production in subjects with essential hypertension.