Robert C. Free

GWAS Central: a comprehensive resource for the comparison and interrogation of genome-wide association studies.

To facilitate broad and convenient integrative visualization of and access to GWAS data, we have created the GWAS Central resource (http://www.gwascentral.org). This database seeks to provide a comprehensive collection of summary-level genetic association data, structured both for maximal utility and for safe open access (i.e., non-directional signals to fully preclude research subject identification). The resource emphasizes on advanced tools that allow comparison and discovery of relevant data sets from the perspective of genes, genome regions, phenotypes or traits. Tested markers and relevant genomic features can be visually interrogated across up to 16 multiple association data sets in a single view, starting at a chromosome-wide view and increasing in resolution down to individual bases. In addition, users can privately upload and view their own data as temporary files. Search and display utility is further enhanced by exploiting phenotype ontology annotations to allow genetic variants associated with phenotypes and traits of interest to be precisely identified, across all studies. Data submissions are accepted from individual researchers, groups and consortia, whereas we also actively gather data sets from various public sources. As a result, the resource now provides over 67 million P-values for over 1600 studies, making it the world’s largest openly accessible online collection of summary-level GWAS association information.

HGVbaseG2P: a central genetic association database

The Human Genome Variation database of Genotype to Phenotype information (HGVbaseG2P) is a new central database for summary-level findings produced by human genetic association studies, both large and small. Such a database is needed so that researchers have an easy way to access all the available association study data relevant to their genes, genome regions or diseases of interest. Such a depository will allow true positive signals to be more readily distinguished from false positives (type I error) that fail to consistently replicate. In this paper we describe how HGVbaseG2P has been constructed, and how its data are gathered and organized. We present a range of user-friendly but powerful website tools for searching, browsing and visualizing G2P study findings. HGVbaseG2P is available at http://www.hgvbaseg2p.org.

DNA analysis of outdoor air reveals a high degree of fungal diversity, temporal variability, and genera not seen by spore morphology.

Fungi are ubiquitous with many capable of causing disease by direct infection, toxicoses, or allergy. Fungal spores are present in outdoor air throughout the year, yet airborne diversity is poorly characterised. Airborne fungal spores are routinely counted by microscopy, enabling identification to genera at best. We generated traditional microscopic counts over a year, then used environmental sequencing techniques to assess and compare 3 d selected from the main fungal spore season. The days selected corresponded to one with a high quantity of spores unidentifiable by microscopy, and two representing dry and wet summer periods. Over 86 % of genera detected by sequencing were not routinely identifiable by microscopy. A high degree of temporal variability was detected, with the percentage of clones attributed to Basidiomycota or Ascomycota, and composition of genera within each phylum varying greatly between days. Throughout the year Basidiomycota spores were found at higher levels than Ascomycota, but levels fluctuated daily with Ascomycota comprising 11-84 % of total spores and Basidiomycota 7-81 %. No significant difference was found between the proportion of clones attributed to each morphological group detected by sequencing to that counted by microscopy (P = 0.477, 0.985, and 0.561). The majority of abundant genera detected by DNA analysis are not routinely identified by microscopy (>75 %). Of those, several are known human and plant pathogens, and may represent unrecognised aeroallergens.

ciliaFA: a research tool for automated, high-throughput measurement of ciliary beat frequency using freely available software.

BackgroundAnalysis of ciliary function for assessment of patients suspected of primary ciliary dyskinesia (PCD) and for research studies of respiratory and ependymal cilia requires assessment of both ciliary beat pattern and beat frequency. While direct measurement of beat frequency from high-speed video recordings is the most accurate and reproducible technique it is extremely time consuming. The aim of this study was to develop a freely available automated method of ciliary beat frequency analysis from digital video (AVI) files that runs on open-source software (ImageJ) coupled to Microsoft Excel, and to validate this by comparison to the direct measuring high-speed video recordings of respiratory and ependymal cilia. These models allowed comparison to cilia beating between 3 and 52 Hz.MethodsDigital video files of motile ciliated ependymal (frequency range 34 to 52 Hz) and respiratory epithelial cells (frequency 3 to 18 Hz) were captured using a high-speed digital video recorder. To cover the range above between 18 and 37 Hz the frequency of ependymal cilia were slowed by the addition of the pneumococcal toxin pneumolysin. Measurements made directly by timing a given number of individual ciliary beat cycles were compared with those obtained using the automated ciliaFA system.ResultsThe overall mean difference (± SD) between the ciliaFA and direct measurement high-speed digital imaging methods was −0.05 ± 1.25 Hz, the correlation coefficient was shown to be 0.991 and the Bland-Altman limits of agreement were from −1.99 to 1.49 Hz for respiratory and from −2.55 to 3.25 Hz for ependymal cilia.ConclusionsA plugin for ImageJ was developed that extracts pixel intensities and performs fast Fourier transformation (FFT) using Microsoft Excel. The ciliaFA software allowed automated, high throughput measurement of respiratory and ependymal ciliary beat frequency (range 3 to 52 Hz) and avoids operator error due to selection bias. We have included free access to the ciliaFA plugin and installation instructions in Additional file 1 accompanying this manuscript that other researchers may use.

Single-step QuantiFERON screening of adult contacts: a prospective cohort study of tuberculosis risk

Background The efffectiveness of tuberculosis (TB) contact screening programmes using interferon γ release assays remains uncertain as prospective contact TB risk is not well characterised. Objectives To quantify 2-year TB risk and evaluate screening performance with single-step QuantiFERON TB Gold-In Tube (QFT) in adult contacts. To compare TB risk between QFT tested subgroups stratified by exposure type (smear positive pulmonary (SP) versus non-smear positive (NSP) TB) and age (younger (16–35 years) versus older (≥36 years)). Methods Screening involved QFT testing in older contacts of SP and all younger contacts, 8–12 weeks after index notification. Chemoprevention (3RH) was offered to QFT positive (+) younger adults. TB risk was determined in a prospective cohort study. Results 43 TB events occurred in 1769 adult contacts observed for median 717 days (2-year rate (95% CI)=2·5% (1.7 to 3.2)). Index-contact strain matching was demonstrable for 18 of 22 (82%) paired samples. No contacts (0/98) receiving 3RH developed TB. 215 of 817 appropriately tested adults (26.3%) were QFT+. 14 of 112 untreated QFT+ adults developed TB (2-year rate (95% CI)=13·4% (7.7 to 21.1)). The model required 35 contacts screened with QFT to identify one contact developing TB at 2 years. TB rates were comparable in QFT+ contacts of SP and NSP (rate ratio (RR)=0.98, p=0·962). For QFT+ older contacts, the disease rate was lower (8.9% (3.3 to 19.1)) and similar to the overall group rate (RR=1.4, p=0.503). Conclusions QFT based single-step contact screening is effective in young adults.

A Two-Tube Combined TaqMan/SYBR Green Assay to Identify Mycobacteria and Detect Single Global Lineage-Defining Polymorphisms in Mycobacterium tuberculosis

We have developed a novel real-time PCR assay to identify and perform preliminary genotyping of mycobacteria in a manner tailored to our local service. Within a single thermocycler run, mycobacterial 16S rDNA and the Mycobacterium tuberculosis global lineage-defining RD750 polymorphism are targeted in separate reaction tubes, each of which includes both TaqMan and SYBR Green chemistries. The results of this 16S-RD assay differentiate M. tuberculosis complex (MTBC) from nontuberculous mycobacteria (NTM) and recognize whether or not MTBC isolates belong to the East African-Indian lineage, the single most frequently isolated global MTBC lineage in our service. If required, NTM amplicons may be sequenced to provide more specific identities. We report the technical performance of this assay on 88 mycobacteria-positive cultures and discuss its use in the initial management of mycobacterial infections. The 16S-RD assay correctly identified all 70 MTBC-positive cultures and 17 NTM-positive cultures while contemporaneously recognizing 26 MTBC isolates as within and 44 outside the East African-Indian lineage. In artificial samples, the combined assay also showed limited potential to detect mixed mycobacterial infections (MTBC/NTM) and tuberculosis infections involving more than one global MTBC lineage. The approach we have established can be readily tailored to targets of particular value for any mycobacterial diagnostic service, thereby optimizing the value of the results for local clinical and public health management of mycobacterial infections.

Automated design of paralogue ratio test assays for the accurate and rapid typing of copy number variation

Motivation: Genomic copy number variation (CNV) can influence susceptibility to common diseases. High-throughput measurement of gene copy number on large numbers of samples is a challenging, yet critical, stage in confirming observations from sequencing or array Comparative Genome Hybridization (CGH). The paralogue ratio test (PRT) is a simple, cost-effective method of accurately determining copy number by quantifying the amplification ratio between a target and reference amplicon. PRT has been successfully applied to several studies analyzing common CNV. However, its use has not been widespread because of difficulties in assay design. Results: We present PRTPrimer (www.prtprimer.org) software for automated PRT assay design. In addition to stand-alone software, the web site includes a database of pre-designed assays for the human genome at an average spacing of 6 kb and a web interface for custom assay design. Other reference genomes can also be analyzed through local installation of the software. The usefulness of PRTPrimer was tested within known CNV, and showed reproducible quantification. This software and database provide assays that can rapidly genotype CNV, cost-effectively, on a large number of samples and will enable the widespread adoption of PRT. Availability: PRTPrimer is available in two forms: a Perl script (version 5.14 and higher) that can be run from the command line on Linux systems and as a service on the PRTPrimer web site (www.prtprimer.org). Contact: cjt14@le.ac.uk Supplementary Information: Supplementary data are available at Bioinformatics online.

Semantically enabling a genome-wide association study database

BackgroundThe amount of data generated from genome-wide association studies (GWAS) has grown rapidly, but considerations for GWAS phenotype data reuse and interchange have not kept pace. This impacts on the work of GWAS Central – a free and open access resource for the advanced querying and comparison of summary-level genetic association data. The benefits of employing ontologies for standardising and structuring data are widely accepted. The complex spectrum of observed human phenotypes (and traits), and the requirement for cross-species phenotype comparisons, calls for reflection on the most appropriate solution for the organisation of human phenotype data. The Semantic Web provides standards for the possibility of further integration of GWAS data and the ability to contribute to the web of Linked Data.ResultsA pragmatic consideration when applying phenotype ontologies to GWAS data is the ability to retrieve all data, at the most granular level possible, from querying a single ontology graph. We found the Medical Subject Headings (MeSH) terminology suitable for describing all traits (diseases and medical signs and symptoms) at various levels of granularity and the Human Phenotype Ontology (HPO) most suitable for describing phenotypic abnormalities (medical signs and symptoms) at the most granular level. Diseases within MeSH are mapped to HPO to infer the phenotypic abnormalities associated with diseases. Building on the rich semantic phenotype annotation layer, we are able to make cross-species phenotype comparisons and publish a core subset of GWAS data as RDF nanopublications.ConclusionsWe present a methodology for applying phenotype annotations to a comprehensive genome-wide association dataset and for ensuring compatibility with the Semantic Web. The annotations are used to assist with cross-species genotype and phenotype comparisons. However, further processing and deconstructions of terms may be required to facilitate automatic phenotype comparisons. The provision of GWAS nanopublications enables a new dimension for exploring GWAS data, by way of intrinsic links to related data resources within the Linked Data web. The value of such annotation and integration will grow as more biomedical resources adopt the standards of the Semantic Web.

The relationship between the Leicester cough questionnaire, eosinophilic airway inflammation and asthma patient related outcomes in severe adult asthma

BackgroundSevere asthma is characterised by a variety of symptoms, which include chronic cough, however the mechanisms responsible for cough reflex hypersensitivity in asthma remain poorly elucidated. Current asthma patient-related outcome instruments such as the six-point Juniper Asthma Control Score (ACQ-6) and Asthma Quality of Life Questionnaire (AQLQ) were not primarily designed to capture cough and its related morbidity in asthma. The Leicester Cough Questionnaire (LCQ) is a patient-related outcome instrument designed to capture the health-related quality of life associated with cough. To date the LCQ has not been evaluated in a severe asthma population.MethodsWe evaluated 262 extensively characterised adult patients with severe asthma attending the Leicester Severe Asthma Service. All patients had a clinician diagnosis of asthma and objective physiological evidence and met the ATS/ERS criterion for servere asthma. In all patients we evaluated a) the LCQ distribution and b) the relationships between the LCQ and ACQ-6, AQLQ, airway inflammation in sputum.ResultsThe LCQ demonstrated the following properties; mean: 15.0, standard deviation: 4.54, median: 15.48, and range: 11.6–19.2. We found a moderate correlation between LCQ and ACQ-6 (r = − 0.605, p < 0.0001) and a LCQ and AQLQ (r = 0.710, p < 0.0001). There was no relationship between LCQ and log10 sputum percentage eosinophils (%).ConclusionA proportion of patients with severe asthma have a significant degree of cough-related morbidity that appears independent of eosinophilic airway inflammation and is not captured fully by existing asthma patient-reported outcome instruments. Our preliminary findings suggest that further research is now required to validate the LCQ and its responsiveness in severe asthma populations to capture cough-related morbidity and response to specific interventions.

Cohort Profile: Extended Cohort for E-health, Environment and DNA (EXCEED)

EXCEED is a longitudinal population-based cohort which facilitates investigation of genetic, environmental and lifestyle-related determinants of a broad range of diseases and of multiple morbidity through data collected at baseline and via electronic healthcare record linkage. Recruitment has taken place in Leicester, Leicestershire and Rutland since 2013 and is ongoing, with 9,840 participants aged 30-69 to date. The population of Leicester is diverse and additional recruitment from the local South Asian community is ongoing. Participants have consented to follow-up for up to 25 years through electronic health records and additional bespoke data collection is planned. Data available includes baseline demographics, anthropometry, spirometry, lifestyle factors (smoking and alcohol use) and longitudinal health information from primary care records, with additional linkage to other EHR datasets planned. Patients have consented to be contacted for recall-by-genotype and recall-by-phenotype sub-studies, providing an important resource for precision medicine research. We welcome requests for collaboration and data access by contacting the study management team via exceed@le.ac.uk.