Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

Purpose:The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.
Methods:Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.
Results:Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.
Conclusion:Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(^-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Partitioning Heritability of Regulatory and Cell-Type-Specific Variants across 11 Common Diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs () across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of  from imputed SNPs (5.1× enrichment; p = 3.7 × 10⁻¹⁷) and 38% (SE = 4%) of  from genotyped SNPs (1.6× enrichment, p = 1.0 × 10⁻⁴). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of  despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Complete nucleotide sequence of φCHU: A Luz24likevirus infecting Pseudomonas aeruginosa and displaying a unique host range

A complete nucleotide sequence of the new Pseudomonas aeruginosa Luz24likevirus phiCHU was obtained. This virus was shown to have a unique host range whereby it grew poorly on the standard laboratory strain PAO1, but infected 26 of 46 clinical isolates screened, and strains harboring IncP2 plasmid pMG53. It was demonstrated that phiCHU has single strand interruptions in its genome. Analysis of the phiCHU genome also suggested that recombination event(s) participated in the evolution of the leftmost portion of the genome, presumably encoding early genes.

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD)

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined P ADJ = 1.15 × 10(-6), OR 1.332 (1.187-1.496)]. This association was characterized by a highly significant sex difference (P diff = 0.0032) in that it was clearly confined to females with genome-wide significance [P ADJ = 2.62 × 10(-8), OR 1.541 (1.324-1.796); males: P ADJ = 0.382, OR 1.084 (0.905-1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMSaveD pathogenesis.

New insights into sequence variation in the IGS region of 21 cyathostomin species and the implication for molecular identification

Cyathostomins comprise a group of 50 species of parasitic nematodes that infect equids. Ribosomal DNA sequences, in particular the intergenic spacer (IGS) region, have been utilized via several methodologies to identify pre-parasitic stages of the commonest species that affect horses. These methods rely on the availability of accurate sequence information for each species, as well as detailed knowledge of the levels of intra- and inter-specific variation. Here, the IGS DNA region was amplified and sequenced from 10 cyathostomin species for which sequence was not previously available. Also, additional IGS DNA sequences were generated from individual worms of 8 species already studied. Comparative analysis of these sequences revealed a greater range of intra-specific variation than previously reported (up to 23%); whilst the level of inter-specific variation (3-62%) was similar to that identified in earlier studies. The reverse line blot (RLB) method has been used to exploit the cyathostomin IGS DNA region for species identification. Here, we report validation of novel and existing DNA probes for identification of cyathostomins using this method and highlight their application in differentiating life-cycle stages such as third-stage larvae that cannot be identified to species by morphological means.

The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution

BACKGROUND: The liver fluke Fasciola hepatica is a major pathogen of livestock worldwide, causing huge economic losses to agriculture, as well as 2.4 million human infections annually.

RESULTS: Here we provide a draft genome for F. hepatica, which we find to be among the largest known pathogen genomes at 1.3 Gb. This size cannot be explained by genome duplication or expansion of a single repeat element, and remains a paradox given the burden it may impose on egg production necessary to transmit infection. Despite the potential for inbreeding by facultative self-fertilisation, substantial levels of polymorphism were found, which highlights the evolutionary potential for rapid adaptation to changes in host availability, climate change or to drug or vaccine interventions. Non-synonymous polymorphisms were elevated in genes shared with parasitic taxa, which may be particularly relevant for the ability of the parasite to adapt to a broad range of definitive mammalian and intermediate molluscan hosts. Large-scale transcriptional changes, particularly within expanded protease and tubulin families, were found as the parasite migrated from the gut, across the peritoneum and through the liver to mature in the bile ducts. We identify novel members of anti-oxidant and detoxification pathways and defined their differential expression through infection, which may explain the stage-specific efficacy of different anthelmintic drugs.

CONCLUSIONS: The genome analysis described here provides new insights into the evolution of this important pathogen, its adaptation to the host environment and external selection pressures. This analysis also provides a platform for research into novel drugs and vaccines.

Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor

BACKGROUND: Klebsiella pneumoniae strains are pathogenic to animals and humans, in which they are both a frequent cause of nosocomial infections and a re-emerging cause of severe community-acquired infections. K. pneumoniae isolates of the capsular serotype K2 are among the most virulent. In order to identify novel putative virulence factors that may account for the severity of K2 infections, the genome sequence of the K2 reference strain Kp52.145 was determined and compared to two K1 and K2 strains of low virulence and to the reference strains MGH 78578 and NTUH-K2044.

RESULTS: In addition to diverse functions related to host colonization and virulence encoded in genomic regions common to the four strains, four genomic islands specific for Kp52.145 were identified. These regions encoded genes for the synthesis of colibactin toxin, a putative cytotoxin outer membrane protein, secretion systems, nucleases and eukaryotic-like proteins. In addition, an insertion within a type VI secretion system locus included sel1 domain containing proteins and a phospholipase D family protein (PLD1). The pld1 mutant was avirulent in a pneumonia model in mouse. The pld1 mRNA was expressed in vivo and the pld1 gene was associated with K. pneumoniae isolates from severe infections. Analysis of lipid composition of a defective E. coli strain complemented with pld1 suggests an involvement of PLD1 in cardiolipin metabolism.

CONCLUSIONS: Determination of the complete genome of the K2 reference strain identified several genomic islands comprising putative elements of pathogenicity. The role of PLD1 in pathogenesis was demonstrated for the first time and suggests that lipid metabolism is a novel virulence mechanism of K. pneumoniae.

Characterization of nontypable Haemophilus influenzae isolates recovered from adult patients with underlying chronic lung disease reveals genotypic and phenotypic traits associated with persistent infection

Nontypable Haemophilus influenzae (NTHi) has emerged as an important opportunistic pathogen causing infection in adults suffering obstructive lung diseases. Existing evidence associates chronic infection by NTHi to the progression of the chronic respiratory disease, but specific features of NTHi associated with persistence have not been comprehensively addressed. To provide clues about adaptive strategies adopted by NTHi during persistent infection, we compared sequential persistent isolates with newly acquired isolates in sputa from six patients with chronic obstructive lung disease. Pulse field gel electrophoresis (PFGE) identified three patients with consecutive persistent strains and three with new strains. Phenotypic characterisation included infection of respiratory epithelial cells, bacterial self-aggregation, biofilm formation and resistance to antimicrobial peptides (AMP). Persistent isolates differed from new strains in showing low epithelial adhesion and inability to form biofilms when grown under continuous-flow culture conditions in microfermenters. Self-aggregation clustered the strains by patient, not by persistence. Increasing resistance to AMPs was observed for each series of persistent isolates; this was not associated with lipooligosaccharide decoration with phosphorylcholine or with lipid A acylation. Variation was further analyzed for the series of three persistent isolates recovered from patient 1. These isolates displayed comparable growth rate, natural transformation frequency and murine pulmonary infection. Genome sequencing of these three isolates revealed sequential acquisition of single-nucleotide variants in the AMP permease sapC, the heme acquisition systems hgpB, hgpC, hup and hxuC, the 3-deoxy-D-manno-octulosonic acid kinase kdkA, the long-chain fatty acid transporter ompP1, and the phosphoribosylamine glycine ligase purD. Collectively, we frame a range of pathogenic traits and a repertoire of genetic variants in the context of persistent infection by NTHi.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Whole genome protein microarrays for serum profiling of immunodominant antigens of Bacillus anthracis

A commercial Bacillus anthracis (Anthrax) whole genome protein microarray has been used to identify immunogenic Anthrax proteins (IAP) using sera from groups of donors with (a) confirmed B. anthracis naturally acquired cutaneous infection, (b) confirmed B. anthracis intravenous drug use-acquired infection, (c) occupational exposure in a wool-sorters factory, (d) humans and rabbits vaccinated with the UK Anthrax protein vaccine and compared to naïve unexposed controls. Anti-IAP responses were observed for both IgG and IgA in the challenged groups; however the anti-IAP IgG response was more evident in the vaccinated group and the anti-IAP IgA response more evident in the B. anthracis-infected groups. Infected individuals appeared somewhat suppressed for their general IgG response, compared with other challenged groups. Immunogenic protein antigens were identified in all groups, some of which were shared between groups whilst others were specific for individual groups. The toxin proteins were immunodominant in all vaccinated, infected or other challenged groups. However, a number of other chromosomally-located and plasmid encoded open reading frame proteins were also recognized by infected or exposed groups in comparison to controls. Some of these antigens e.g., BA4182 are not recognized by vaccinated individuals, suggesting that there are proteins more specifically expressed by live Anthrax spores in vivo that are not currently found in the UK licensed Anthrax Vaccine (AVP). These may perhaps be preferentially expressed during infection and represent expression of alternative pathways in the B. anthracis "infectome." These may make highly attractive candidates for diagnostic and vaccine biomarker development as they may be more specifically associated with the infectious phase of the pathogen. A number of B. anthracis small hypothetical protein targets have been synthesized, tested in mouse immunogenicity studies and validated in parallel using human sera from the same study.

Integrative analysis of chemo-transcriptomic profiles for drug-feature specific gene expression signatures

One of the major challenges in systems biology is to understand the complex responses of a biological system to external perturbations or internal signalling depending on its biological conditions. Genome-wide transcriptomic profiling of cellular systems under various chemical perturbations allows the manifestation of certain features of the chemicals through their transcriptomic expression profiles. The insights obtained may help to establish the connections between human diseases, associated genes and therapeutic drugs. The main objective of this study was to systematically analyse cellular gene expression data under various drug treatments to elucidate drug-feature specific transcriptomic signatures. We first extracted drug-related information (drug features) from the collected textual description of DrugBank entries using text-mining techniques. A novel statistical method employing orthogonal least square learning was proposed to obtain drug-feature-specific signatures by integrating gene expression with DrugBank data. To obtain robust signatures from noisy input datasets, a stringent ensemble approach was applied with the combination of three techniques: resampling, leave-one-out cross validation, and aggregation. The validation experiments showed that the proposed method has the capacity of extracting biologically meaningful drug-feature-specific gene expression signatures. It was also shown that most of signature genes are connected with common hub genes by regulatory network analysis. The common hub genes were further shown to be related to general drug metabolism by Gene Ontology analysis. Each set of genes has relatively few interactions with other sets, indicating the modular nature of each signature and its drug-feature-specificity. Based on Gene Ontology analysis, we also found that each set of drug feature (DF)-specific genes were indeed enriched in biological processes related to the drug feature. The results of these experiments demonstrated the pot- ntial of the method for predicting certain features of new drugs using their transcriptomic profiles, providing a useful methodological framework and a valuable resource for drug development and characterization.

The mitochondrial and autosomal mutation landscapes of prostate cancer

BACKGROUND: Prostate cancer (PCa) is the most common cancer in men. PCa is strongly age associated; low death rates in surveillance cohorts call into question the widespread use of surgery, which leads to overtreatment and a reduction in quality of life. There is a great need to increase the understanding of tumor characteristics in the context of disease progression.

OBJECTIVE: To perform the first multigenome investigation of PCa through analysis of both autosomal and mitochondrial DNA, and to integrate exome sequencing data, and RNA sequencing and copy-number alteration (CNA) data to investigate how various different tumor characteristics, commonly analyzed separately, are interconnected.

DESIGN, SETTING, AND PARTICIPANTS: Exome sequencing was applied to 64 tumor samples from 55 PCa patients with varying stage and grade. Integrated analysis was performed on a core set of 50 tumors from which exome sequencing, CNA, and RNA sequencing data were available.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Genes, mutated at a significantly higher rate relative to a genomic background, were identified. In addition, mitochondrial and autosomal mutation rates were correlated to CNAs and proliferation, assessed as a cell cycle gene expression signature.

RESULTS AND LIMITATIONS: Genes not previously reported to be significantly mutated in PCa, such as cell division cycle 27 homolog (Saccharomyces cerevisiae) (CDC27), myeloid/lymphoid or mixed-lineage leukemia 3 (MLL3), lysine (K)-specific demethylase 6A (KDM6A), and kinesin family member 5A (KIF5A) were identified. The mutation rate in the mitochondrial genome was 55 times higher than that of the autosomes. Multilevel analysis demonstrated a tight correlation between high reactive-oxygen exposure, chromosomal damage, high proliferation, and in parallel, a transition from multiclonal indolent primary PCa to monoclonal aggressive disease. As we only performed targeted sequence analysis; copy-number neutral rearrangements recently described for PCa were not accounted for.

CONCLUSIONS: The mitochondrial genome displays an elevated mutation rate compared to the autosomal chromosomes. By integrated analysis, we demonstrated that different tumor characteristics are interconnected, providing an increased understanding of PCa etiology.

Global identification of androgen response elements

Chromatin immunoprecipitation (ChIP) is an invaluable tool in the study of transcriptional regulation. ChIP methods require both a priori knowledge of the transcriptional regulators which are important for a given biological system and high-quality specific antibodies for these targets. The androgen receptor (AR) is known to play essential roles in male sexual development, in prostate cancer and in the function of many other AR-expressing cell types (e.g. neurons and myocytes). As a ligand-activated transcription factor the AR also represents an endogenous, inducible system to study transcriptional biology. Therefore, ChIP studies of the AR can make use of treatment contrast experiments to define its transcriptional targets. To date several studies have mapped AR binding sites using ChIP in combination with genome tiling microarrays (ChIP-chip) or direct sequencing (ChIP-seq), mainly in prostate cancer cell lines and with varying degrees of genomic coverage. These studies have provided new insights into the DNA sequences to which the AR can bind, identified AR cooperating transcription factors, mapped thousands of potential AR regulated genes and provided insights into the biological processes regulated by the AR. However, further ChIP studies will be required to fully characterise the dynamics of the AR-regulated transcriptional programme, to map the occupancy of different AR transcriptional complexes which result in different transcriptional output and to delineate the transcriptional networks downstream of the AR.

A stacked Late Quaternary fluvio-periglacial sequence from the Axe valley, southern England with implications for landscape evolution and Palaeolithic archaeology

The current model of mid-latitude late Quaternary terrace sequences, is that they are uplift-driven but climatically controlled terrace staircases, relating to both regional-scale crustal and tectonic factors, and palaeohydrological variations forced by quasi-cyclic climatic conditions in the 100 K world (post Mid Pleistocene Transition). This model appears to hold for the majority of the river valleys draining into the English Channel which exhibit 8–15 terrace levels over approximately 60–100 m of altitudinal elevation. However, one valley, the Axe, has only one major morphological terrace and has long-been regarded as anomalous. This paper uses both conventional and novel stratigraphical methods (digital granulometry and terrestrial laser scanning) to show that this terrace is a stacked sedimentary sequence of 20–30 m thickness with a quasi-continuous (i.e. with hiatuses) pulsed, record of fluvial and periglacial sedimentation over at least the last 300–400 K yrs as determined principally by OSL dating of the upper two thirds of the sequence. Since uplift has been regional, there is no evidence of anomalous neotectonics, and climatic history must be comparable to the adjacent catchments (both of which have staircase sequences) a catchment-specific mechanism is required. The Axe is the only valley in North West Europe incised entirely into the near-horizontally bedded chert (crypto-crystalline quartz) and sand-rich Lower Cretaceous rocks creating a buried valley. Mapping of the valley slopes has identified many large landslide scars associated with past and present springs. It is proposed that these are thaw-slump scars and represent large hill-slope failures caused by Vauclausian water pressures and hydraulic fracturing of the chert during rapid permafrost melting. A simple 1D model of this thermokarstic process is used to explore this mechanism, and it is proposed that the resultant anomalously high input of chert and sand into the valley during terminations caused pulsed aggradation until the last termination. It is also proposed that interglacial and interstadial incision may have been prevented by the over-sized and interlocking nature of the sub-angular chert clasts until the Lateglacial when confinement of the river overcame this immobility threshold. One result of this hydrogeologically mediated valley evolution was to provide a sequence of proximal Palaeolithic archaeology over two MIS cycles. This study demonstrates that uplift tectonics and climate alone do not fully determine Quaternary valley evolution and that lithological and hydrogeological conditions are a fundamental cause of variation in terrestrial Quaternary records and landform evolution.