Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Investigating ANKH and ENPP1 in Slovakian families with chondrocalcinosis

Familial articular chondrocalcinosis (CC) was Wrst reported in 1963. It is characterised by multiple calciWcations of hyaline and Wbrous cartilage in the joints and intervertebral discs. Mutations in ANKH have been identified in several pedigrees as a monogenic cause for this disorder. ANKH is a key protein in pyrophosphate metabolism and is involved in pyrophosphate transport across the cell membrane. The objective of this work was to screen ANKH and ENPP1, two key genes in pyrophosphate metabolism, in Slovakian kindreds with familial CC. DNA samples from 25 individuals (10 aVected, 15 unaVected) from 8 families were obtained. The promoter, coding regions and intron-exon boundaries of ANKH and ENPP1 were sequenced. Twelve DNA sequence variants, six in each gene, were identiWed. All the variants had been previously identified. None segregated with the disease. Our results suggest that neither ANKH nor ENPP1 mutations are the cause of CC in these families, indicating that possibly other major genes are involved in the aethiopathogenesis of this condition in these families.

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10-10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10-9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10-14, OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10-14, OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma. © 2011 Nature America, Inc. All rights reserved.

Robust face clustering for real-world data

This thesis has investigated how to cluster a large number of faces within a multi-media corpus in the presence of large session variation. Quality metrics are used to select the best faces to represent a sequence of faces; and session variation modelling improves clustering performance in the presence of wide variations across videos. Findings from this thesis contribute to improving the performance of both face verification systems and the fully automated clustering of faces from a large video corpus.

Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 103). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10-9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 × 10-7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10-9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression. © 2009 The Author(s).

Partial epilepsy syndrome in a Gypsy family linked to 5q31.3-q32

PURPOSE The restricted genetic diversity and homogeneous molecular basis of Mendelian disorders in isolated founder populations have rarely been explored in epilepsy research. Our long-term goal is to explore the genetic basis of epilepsies in one such population, the Gypsies. The aim of this report is the clinical and genetic characterization of a Gypsy family with a partial epilepsy syndrome. METHODS Clinical information was collected using semistructured interviews with affected subjects and informants. At least one interictal electroencephalography (EEG) recording was performed for each patient and previous data obtained from records. Neuroimaging included structural magnetic resonance imaging (MRI). Linkage and haplotype analysis was performed using the Illumina IVb Linkage Panel, supplemented with highly informative microsatellites in linked regions and Affymetrix SNP 5.0 array data. RESULTS We observed an early-onset partial epilepsy syndrome with seizure semiology strongly suggestive of temporal lobe epilepsy (TLE), with mild intellectual deficit co-occurring in a large proportion of the patients. Psychiatric morbidity was common in the extended pedigree but did not cosegregate with epilepsy. Linkage analysis definitively excluded previously reported loci, and identified a novel locus on 5q31.3-q32 with an logarithm of the odds (LOD) score of 3 corresponding to the expected maximum in this family. DISCUSSION The syndrome can be classified as familial temporal lobe epilepsy (FTLE) or possibly a new syndrome with mild intellectual deficit. The linked 5q region does not contain any ion channel-encoding genes and is thus likely to contribute new knowledge about epilepsy pathogenesis. Identification of the mutation in this family and in additional patients will define the full phenotypic spectrum.

Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.

Investigation of the role of ENPP1 and TNAP genes in chondrocalcinosis

Objectives. Extracellular inorganic pyrophosphate (ePPi) inhibits certain forms of pathological mineralization while promoting others. Three molecules involved in ePPi regulation are important candidates for the development of calcium pyrophosphate dihydrate chondrocalcinosis (CPPD CC). These include ANKH, ectonucleotide pyrophosphatase (ENPP1) and TNAP. We have previously showed that genetic variation in ANKH is a cause of autosomal dominant familial CC and also some sporadic cases of CPPD CC. We now investigate the possible role of ENPP1 and TNAP in CPPD CC. Methods. Exons, untranslated regions (UTR) and exon-intron boundaries of ENPP1 and TNAP were sequenced using ABI Big Dye chemistry on automated sequencers. Sixteen variants were identified (3 in ENPP1 and 13 in TNAP) and were subsequently genotyped in 128 sporadic Caucasian CPPD CC patients and 600 healthy controls using a combination of polymerase chain reaction/restriction fragment-length polymorphism analysis or using Taqman. Allele and genotype frequencies were compared between cases and controls using the χ 2 test. Linkage disequilibrium, haplotype and the single nucleotide polymorphism-specific analyses were also performed. This study had 80% power to detect an odds ratio of 2.2 or more at these loci. Results. No difference was observed in the allele or genotype frequencies between patients and controls at either ENPP1 or TNAP. Conclusions. Polymorphisms of ENPP1 and TNAP are not major determinants of susceptibility to CC in the population studied. Further studies of the aetiology of sporadic CPPD CC are required to determine its causes.

Genetic, functional, and histopathological evaluation of two C-terminal BRCA1 missense variants

Background: The vast majority of BRCA1 missense sequence variants remain uncharacterised for their possible effect on protein expression and function, and therefore are unclassified in terms of their pathogenicity. BRCA1 plays diverse cellular roles and it is unlikely that any single functional assay will accurately reflect the total cellular implications of missense mutations in this gene. Objective: To elucidate the effect of two BRCA1 variants, 5236G>C (G1706A) and 5242C>A (A1708E) on BRCA1 function, and to survey the relative usefulness of several assays to direct the characterisation of other unclassified variants in BRCA genes. Methods and Results: Data from a range of bioinformatic, genetic, and histopathological analyses, and in vitro functional assays indicated that the 1708E variant was associated with the disruption of different cellular functions of BRCA1. In transient transfection experiments in T47D and 293T cells, the 1708E product was mislocalised to the cytoplasm and induced centrosome amplification in 293T cells. The 1708E variant also failed to transactivate transcription of reporter constructs in mammalian transcriptional transactivation assays. In contrast, the 1706A variant displayed a phenotype comparable to wildtype BRCA1 in these assays. Consistent with functional data, tumours from 1708E carriers showed typical BRCA1 pathology, while tumour material from 1706A carriers displayed few histopathological features associated with BRCA1 related tumours. Conclusions: A comprehensive range of genetic, bioinformatic, and functional analyses have been combined for the characterisation of BRCA1 unclassified sequence variants. Consistent with the functional analyses, the combined odds of causality calculated for the 1706A variant after multifactorial likelihood analysis (1:142) indicates a definitive classification of this variant as "benign". In contrast, functional assays of the 1708E variant indicate that it is pathogenic, possibly through subcellular mislocalisation. However, the combined odds of 262:1 in favour of causality of this variant does not meet the minimal ratio of 1000:1 for classification as pathogenic, and A1708E remains formally designated as unclassified. Our findings highlight the importance of comprehensive genetic information, together with detailed functional analysis for the definitive categorisation of unclassified sequence variants. This combination of analyses may have direct application to the characterisation of other unclassified variants in BRCA1 and BRCA2.

In vivo positive contrast IRON sequence and quantitative T2* measurement confirms inflammatory burden in a patient with asymptomatic carotid atheroma after USPIO-enhanced MR Imaging

The authors report an in vivo human examination of carotid atheroma by using the inversion-recovery ON resonance (IRON) sequence, which is able to produce positive contrast after the infusion of an ultrasmall super paramagnetic iron oxide (USPIO) contrast medium. This technique provides a method of potentially identifying inflammatory burden within carotid atheroma. This may be particularly useful in patients who currently do not meet criteria for intervention (ie, moderate symptomatic stenosis or <70% asymptomatic stenosis) to further risk-stratify this important patient cohort. A 63-year-old man was imaged at 1.5 T before and 36 hours after USPIO infusion by using the IRON sequence. Regions of interest showing profound signal loss at T2*-weighted imaging corresponded well with regions of positive contrast at IRON imaging after the administration of USPIO. These regions also showed a profound decrease in T2* measurements after USPIO infusion, whereas surrounding tissue did not. It has been shown that such strong signal loss on T2*-weighted images after USPIO infusion is indicative of USPIO uptake.

Analysis of word embeddings and sequence features for clinical information extraction

This study investigates the use of unsupervised features derived from word embedding approaches and novel sequence representation approaches for improving clinical information extraction systems. Our results corroborate previous findings that indicate that the use of word embeddings significantly improve the effectiveness of concept extraction models; however, we further determine the influence that the corpora used to generate such features have. We also demonstrate the promise of sequence-based unsupervised features for further improving concept extraction.

Accounting for individual variability in the von Bertalanffy growth model

Estimation of von Bertalanffy growth parameters has received considerable attention in fisheries research. Since Sainsbury (1980, Can. J. Fish. Aquat. Sci. 37: 241-247) much of this research effort has centered on accounting for individual variability in the growth parameters. In this paper we demonstrate that, in analysis of tagging data, Sainsbury's method and its derivatives do not, in general, satisfactorily account for individual variability in growth, leading to inconsistent parameter estimates (the bias does not tend to zero as sample size increases to infinity). The bias arises because these methods do not use appropriate conditional expectations as a basis for estimation. This bias is found to be similar to that of the Fabens method. Such methods would be appropriate only under the assumption that the individual growth parameters that generate the growth increment were independent of the growth parameters that generated the initial length. However, such an assumption would be unrealistic. The results are derived analytically, and illustrated with a simulation study. Until techniques that take full account of the appropriate conditioning have been developed, the effect of individual variability on growth has yet to be fully understood.

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

Effect of individual variability on estimation of population parameters from length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock when there is individual variability in the von Bertalanffy growth parameter L-infinity and investigate the possible bias in the estimates when the individual variability is ignored. Three methods are examined: (i) the regression method based on the Beverton and Holt's (1956, Rapp. P.V. Reun. Cons. Int. Explor. Mer, 140: 67-83) equation; (ii) the moment method of Powell (1979, Rapp. PV. Reun. Int. Explor. Mer, 175: 167-169); and (iii) a generalization of Powell's method that estimates the individual variability to be incorporated into the estimation. It is found that the biases in the estimates from the existing methods are, in general, substantial, even when individual variability in growth is small and recruitment is uniform, and the generalized method performs better in terms of bias but is subject to a larger variation. There is a need to develop robust and flexible methods to deal with individual variability in the analysis of length-frequency data.