John Broxholme

Positional cloning of a quantitative trait locus on chromosome 13q14 that influences immunoglobulin E levels and asthma

Atopic or immunoglobulin E (IgE)-mediated diseases include the common disorders of asthma, atopic dermatitis and allergic rhinitis. Chromosome 13q14 shows consistent linkage to atopy and the total serum IgE concentration. We previously identified association between total serum IgE levels and a novel 13q14 microsatellite (USAT24G1; ref. 7) and have now localized the underlying quantitative-trait locus (QTL) in a comprehensive single-nucleotide polymorphism (SNP) map. We found replicated association to IgE levels that was attributed to several alleles in a single gene, PHF11. We also found association with these variants to severe clinical asthma. The gene product (PHF11) contains two PHD zinc fingers and probably regulates transcription. Distinctive splice variants were expressed in immune tissues and cells.

Positional cloning of a novel gene influencing asthma from chromosome 2q14

Asthma is a common disease in children and young adults. Four separate reports have linked asthma and related phenotypes to an ill-defined interval between 2q14 and 2q32 (refs. 1–4), and two mouse genome screens have linked bronchial hyper-responsiveness to the region homologous to 2q14 (refs. 5,6). We found and replicated association between asthma and the D2S308 microsatellite, 800 kb distal to the IL1 cluster on 2q14. We sequenced the surrounding region and constructed a comprehensive, high-density, single-nucleotide polymorphism (SNP) linkage disequilibrium (LD) map. SNP association was limited to the initial exons of a solitary gene of 3.6 kb (DPP10), which extends over 1 Mb of genomic DNA. DPP10 encodes a homolog of dipeptidyl peptidases (DPPs) that cleave terminal dipeptides from cytokines and chemokines, and it presents a potential new target for asthma therapy.

A Fine-Scale Chimpanzee Genetic Map from Population Sequencing

Going Ape Over Genetic Maps Recombination is an important process in generating diversity and producing selectively advantageous genetic combinations. Thus, changes in recombination hotspots may influence speciation. To investigate the variation in recombination processes in humans and their closest existing relatives, Auton et al. (p. 193, published online 15 March) prepared a fine-scale genetic map of the Western chimpanzee and compared it with that of humans. While rates of recombination are comparable between humans and chimpanzees, the locations and genetic motifs associated with recombination differ between the species. Chimpanzees show similar genetic recombination rates as humans but differ in the genomic regions involved. To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

Clinical whole-genome sequencing in severe early-onset epilepsy reveals new genes and improves molecular diagnosis

In severe early-onset epilepsy, precise clinical and molecular genetic diagnosis is complex, as many metabolic and electro-physiological processes have been implicated in disease causation. The clinical phenotypes share many features such as complex seizure types and developmental delay. Molecular diagnosis has historically been confined to sequential testing of candidate genes known to be associated with specific sub-phenotypes, but the diagnostic yield of this approach can be low. We conducted whole-genome sequencing (WGS) on six patients with severe early-onset epilepsy who had previously been refractory to molecular diagnosis, and their parents. Four of these patients had a clinical diagnosis of Ohtahara Syndrome (OS) and two patients had severe non-syndromic early-onset epilepsy (NSEOE). In two OS cases, we found de novo non-synonymous mutations in the genes KCNQ2 and SCN2A. In a third OS case, WGS revealed paternal isodisomy for chromosome 9, leading to identification of the causal homozygous missense variant in KCNT1, which produced a substantial increase in potassium channel current. The fourth OS patient had a recessive mutation in PIGQ that led to exon skipping and defective glycophosphatidyl inositol biosynthesis. The two patients with NSEOE had likely pathogenic de novo mutations in CBL and CSNK1G1, respectively. Mutations in these genes were not found among 500 additional individuals with epilepsy. This work reveals two novel genes for OS, KCNT1 and PIGQ. It also uncovers unexpected genetic mechanisms and emphasizes the power of WGS as a clinical tool for making molecular diagnoses, particularly for highly heterogeneous disorders.

Mutations in TCF12 , encoding a basic helix-loop-helix partner of TWIST1, are a frequent cause of coronal craniosynostosis

Craniosynostosis, the premature fusion of the cranial sutures, is a heterogeneous disorder with a prevalence of ∼1 in 2,200 (refs. 1,2). A specific genetic etiology can be identified in ∼21% of cases, including mutations of TWIST1, which encodes a class II basic helix-loop-helix (bHLH) transcription factor, and causes Saethre-Chotzen syndrome, typically associated with coronal synostosis. Using exome sequencing, we identified 38 heterozygous TCF12 mutations in 347 samples from unrelated individuals with craniosynostosis. The mutations predominantly occurred in individuals with coronal synostosis and accounted for 32% and 10% of subjects with bilateral and unilateral pathology, respectively. TCF12 encodes one of three class I E proteins that heterodimerize with class II bHLH proteins such as TWIST1. We show that TCF12 and TWIST1 act synergistically in a transactivation assay and that mice doubly heterozygous for loss-of-function mutations in Tcf12 and Twist1 have severe coronal synostosis. Hence, the dosage of TCF12-TWIST1 heterodimers is critical for normal coronal suture development.

SW-ARRAY: a dynamic programming solution for the identification of copy-number changes in genomic DNA using array comparative genome hybridization data

Comparative genome hybridization (CGH) to DNA microarrays (array CGH) is a technique capable of detecting deletions and duplications in genomes at high resolution. However, array CGH studies of the human genome noting false negative and false positive results using large insert clones as probes have raised important concerns regarding the suitability of this approach for clinical diagnostic applications. Here, we adapt the Smith–Waterman dynamic-programming algorithm to provide a sensitive and robust analytic approach (SW-ARRAY) for detecting copy-number changes in array CGH data. In a blind series of hybridizations to arrays consisting of the entire tiling path for the terminal 2 Mb of human chromosome 16p, the method identified all monosomies between 267 and 1567 kb with a high degree of statistical significance and accurately located the boundaries of deletions in the range 267–1052 kb. The approach is unique in offering both a nonparametric segmentation procedure and a nonparametric test of significance. It is scalable and well-suited to high resolution whole genome array CGH studies that use array probes derived from large insert clones as well as PCR products and oligonucleotides.

Next-generation sequencing (NGS) as a diagnostic tool for retinal degeneration reveals a much higher detection rate in early-onset disease.

Inherited retinal degeneration (IRD) is a common cause of visual impairment (prevalence ∼1/3500). There is considerable phenotype and genotype heterogeneity, making a specific diagnosis very difficult without molecular testing. We investigated targeted capture combined with next-generation sequencing using Nimblegen 12plex arrays and the Roche 454 sequencing platform to explore its potential for clinical diagnostics in two common types of IRD, retinitis pigmentosa and cone-rod dystrophy. 50 patients (36 unknowns and 14 positive controls) were screened, and pathogenic mutations were identified in 25% of patients in the unknown, with 53% in the early-onset cases. All patients with new mutations detected had an age of onset <21 years and 44% had a family history. Thirty-one percent of mutations detected were novel. A de novo mutation in rhodopsin was identified in one early-onset case without a family history. Bioinformatic pipelines were developed to identify likely pathogenic mutations and stringent criteria were used for assignment of pathogenicity. Analysis of sequencing metrics revealed significant variability in capture efficiency and depth of coverage. We conclude that targeted capture and next-generation sequencing are likely to be very useful in a diagnostic setting, but patients with earlier onset of disease are more likely to benefit from using this strategy. The mutation-detection rate suggests that many patients are likely to have mutations in novel genes.

Host Genetic Factors and Vaccine-Induced Immunity to Hepatitis B Virus Infection

Background Vaccination against hepatitis B virus infection (HBV) is safe and effective; however, vaccine-induced antibody level wanes over time. Peak vaccine-induced anti-HBs level is directly related to antibody decay, as well as risk of infection and persistent carriage despite vaccination. We investigated the role of host genetic factors in long-term immunity against HBV infection based on peak anti-HBs level and seroconversion to anti-HBc. Methods We analyzed 715 SNP across 133 candidate genes in 662 infant vaccinees from The Gambia, assessing peak vaccine-induced anti-HBs level and core antibody (anti-HBc) status, whilst adjusting for covariates. A replication study comprised 43 SNPs in a further 393 individuals. Results In our initial screen we found variation in IFNG, MAPK8, and IL10RA to affect peak anti-HBs level (GMTratio of <0.6 or >1.5 and P≤0.001) and lesser associations in other genes. Odds of core-conversion was associated with variation in CD163. A coding change in ITGAL (R719V) with likely functional relevance showed evidence of association with increased peak anti-HBs level in both screens (1st screen: s595_22 GMTratio 1.71, P = 0.013; 2nd screen: s595_22 GMTratio 2.15, P = 0.011). Conclusion This is to our knowledge the largest study to date assessing genetic determinants of HBV vaccine-induced immunity. We report on associations with anti-HBs level, which is directly related to durability of antibody level and predictive of vaccine efficacy long-term. A coding change in ITGAL, which plays a central role in immune cell interaction, was shown to exert beneficial effects on induction of peak antibody level in response to HBV vaccination. Variation in this gene does not appear to have been studied in relation to immune responses to viral or vaccine challenges previously. Our findings suggest that genetic variation in loci other than the HLA region affect immunity induced by HBV vaccination.

Missense mutation in sterile alpha motif of novel protein SamCystin is associated with polycystic kidney disease in (cy/+)rat

Autosomal dominant polycystic kidney disease (PKD) is the most common genetic disease that leads to kidney failure in humans. In addition to the known causative genes PKD1 and PKD2, there are mutations that result in cystic changes in the kidney, such as nephronophthisis, autosomal recessive polycystic kidney disease, or medullary cystic kidney disease. Recent efforts to improve the understanding of renal cystogenesis have been greatly enhanced by studies in rodent models of PKD. Genetic studies in the (cy/+) rat showed that PKD spontaneously develops as a consequence of a mutation in a gene different from the rat orthologs of PKD1 and PKD2 or other genes that are known to be involved in human cystic kidney diseases. This article reports the positional cloning and mutation analysis of the rat PKD gene, which revealed a C to T transition that replaces an arginine by a tryptophan at amino acid 823 in the protein sequence. It was determined that Pkdr1 is specifically expressed in renal proximal tubules and encodes a novel protein, SamCystin, that contains ankyrin repeats and a sterile alpha motif. The characterization of this protein, which does not share structural homologies with known polycystins, may give new insights into the pathophysiology of renal cyst development in patients.

Pervasive haplotypic variation in the spliceo-transcriptome of the human major histocompatibility complex

The human major histocompatibility complex (MHC) on chromosome 6p21 is a paradigm for genomics, showing remarkable polymorphism and striking association with immune and non-immune diseases. The complex genomic landscape of the MHC, notably strong linkage disequilibrium, has made resolving causal variants very challenging. A promising approach is to investigate gene expression levels considered as tractable intermediate phenotypes in mapping complex diseases. However, how transcription varies across the MHC, notably relative to specific haplotypes, remains unknown. Here, using an original hybrid tiling and splice junction microarray that includes alternate allele probes, we draw the first high-resolution strand-specific transcription map for three common MHC haplotypes (HLA-A1-B8-Cw7-DR3, HLA-A3-B7-Cw7-DR15, and HLA-A26-B18-Cw5-DR3-DQ2) strongly associated with autoimmune diseases including type 1 diabetes, systemic lupus erythematosus, and multiple sclerosis. We find that haplotype-specific differences in gene expression are common across the MHC, affecting 96 genes (46.4%), most significantly the zing finger protein gene ZFP57. Differentially expressed probes are correlated with polymorphisms between haplotypes, consistent with cis effects that we directly demonstrate for ZFP57 in a cohort of healthy volunteers (P = 1.2 × 10(-14)). We establish that alternative splicing is significantly more frequent in the MHC than genome-wide (72.5% vs. 62.1% of genes, P ≤ 1 × 10(-4)) and shows marked haplotypic differences. We also unmask novel and abundant intergenic transcription involving 31% of transcribed blocks identified. Our study reveals that the renowned MHC polymorphism also manifests as transcript diversity, and our novel haplotype-based approach marks a new step toward identification of regulatory variants involved in the control of MHC-associated phenotypes and diseases.