Richard W. Francis

Predominance of nontypeable Haemophilus influenzae in children with otitis media following introduction of a 3 + 0 pneumococcal conjugate vaccine schedule

In Australia the 7-valent pneumococcal conjugate vaccine (PCV7) is administered at 2, 4 and 6 months of age, with no booster dose. Information on bacterial carriage and the aetiology of recurrent acute otitis media (rAOM) after introduction of PCV7 using the 3+0 schedule is required to evaluate the potential impact of second generation pneumococcal vaccines. We found that 2-4 years after introduction of PCV7 in the National Immunisation Program, nontypeable Haemophilus influenzae (NTHi) was the predominant pathogen isolated from the nasopharynx and middle ear of children with a history of rAOM. Compared with healthy controls (n=81), NTHi and Streptococcus pneumoniae carriage rates were significantly higher in children with a history of rAOM (n=186) (19% vs. 56% p<0.0001 and 26% vs. 41%, p=0.02, respectively). Carriage of PCV7 pneumococcal serotypes was rare, whereas PCV7-related and non-PCV7 serotypes were isolated of 38% of cases and 24% of controls. Serotype 19A was the most common serotype isolated from the nasopharynx and middle ear and accounted for 36% (14/39) of total pneumococcal isolates with reduced susceptibility to cotrimoxazole. Of the 119 children carrying NTHi, 17% of isolates were β-lactamase positive. The scarcity of PCV7 serotypes in children with and without a history of rAOM indicates that the 3+0 PCV7 schedule is preventing carriage and rAOM from PCV7 serotypes. Introduction of new vaccines in Australia with increased pneumococcal serotype and pathogen coverage, including 19A and NTHi, should decrease the circulation of antibiotic-resistant bacteria and reduce the burden of rAOM.

Autism risk associated with parental age and with increasing difference in age between the parents.

Advancing paternal and maternal age have both been associated with risk for autism spectrum disorders (ASD). However, the shape of the association remains unclear, and results on the joint associations is lacking. This study tests if advancing paternal and maternal ages are independently associated with ASD risk and estimates the functional form of the associations. In a population-based cohort study from five countries (Denmark, Israel, Norway, Sweden and Western Australia) comprising 5 766 794 children born 1985–2004 and followed up to the end of 2004–2009, the relative risk (RR) of ASD was estimated by using logistic regression and splines. Our analyses included 30 902 cases of ASD. Advancing paternal and maternal age were each associated with increased RR of ASD after adjusting for confounding and the other parents age (mothers 40–49 years vs 20–29 years, RR=1.15 (95% confidence interval (CI): 1.06–1.24), P-value<0.001; fathers⩾50 years vs 20–29 years, RR=1.66 (95% CI: 1.49–1.85), P-value<0.001). Younger maternal age was also associated with increased risk for ASD (mothers <20 years vs 20–29 years, RR=1.18 (95% CI: 1.08–1.29), P-value<0.001). There was a joint effect of maternal and paternal age with increasing risk of ASD for couples with increasing differences in parental ages. We did not find any support for a modifying effect by the sex of the offspring. In conclusion, as shown in multiple geographic regions, increases in ASD was not only limited to advancing paternal or maternal age alone but also to differences parental age including younger or older similarly aged parents as well as disparately aged parents.

A screen for endocytic motifs.

Sorting signals for cargo selection into coated vesicles are usually in the form of short linear motifs. Three motifs for clathrin‐mediated endocytosis have been identified: YXXΦ, [D/E]XXXL[L/I] and FXNPXY. To search for new endocytic motifs, we made a library of CD8 chimeras with random sequences in their cytoplasmic tails, and used a novel fluorescence‐activated cell sorting (FACS)‐based assay to select for endocytosed constructs. Out of the five tails that were most efficiently internalized, only one was found to contain a conventional motif. Two contain dileucine‐like sequences that appear to be variations on the [D/E]XXXL[L/I] motif. Another contains a novel internalization signal, YXXXΦN, which is able to function in cells expressing a mutant µ2 that cannot bind YXXΦ, indicating that it is not a variation on the YXXΦ motif. Similar sequences are present in endogenous proteins, including a functional YXXXΦN (in addition to a classical YXXΦ) in cytotoxic T‐lymphocyte‐associated protein 4 (CTLA‐4). Thus, the repertoire of endocytic motifs is more extensive than the three well‐characterized sorting signals.

DataSHIELD: taking the analysis to the data, not the data to the analysis

Background: Research in modern biomedicine and social science requires sample sizes so large that they can often only be achieved through a pooled co-analysis of data from several studies. But the pooling of information from individuals in a central database that may be queried by researchers raises important ethico-legal questions and can be controversial. In the UK this has been highlighted by recent debate and controversy relating to the UK’s proposed ‘care.data’ initiative, and these issues reflect important societal and professional concerns about privacy, confidentiality and intellectual property. DataSHIELD provides a novel technological solution that can circumvent some of the most basic challenges in facilitating the access of researchers and other healthcare professionals to individual-level data. Methods: Commands are sent from a central analysis computer (AC) to several data computers (DCs) storing the data to be co-analysed. The data sets are analysed simultaneously but in parallel. The separate parallelized analyses are linked by non-disclosive summary statistics and commands transmitted back and forth between the DCs and the AC. This paper describes the technical implementation of DataSHIELD using a modified R statistical environment linked to an Opal database deployed behind the computer firewall of each DC. Analysis is controlled through a standard R environment at the AC. Results: Based on this Opal/R implementation, DataSHIELD is currently used by the Healthy Obese Project and the Environmental Core Project (BioSHaRE-EU) for the federated analysis of 10 data sets across eight European countries, and this illustrates the opportunities and challenges presented by the DataSHIELD approach. Conclusions: DataSHIELD facilitates important research in settings where: (i) a co-analysis of individual-level data from several studies is scientifically necessary but governance restrictions prohibit the release or sharing of some of the required data, and/or render data access unacceptably slow; (ii) a research group (e.g. in a developing nation) is particularly vulnerable to loss of intellectual property—the researchers want to fully share the information held in their data with national and international collaborators, but do not wish to hand over the physical data themselves; and (iii) a data set is to be included in an individual-level co-analysis but the physical size of the data precludes direct transfer to a new site for analysis.

FusionFinder: a software tool to identify expressed gene fusion candidates from RNA-Seq data.

The hallmarks of many haematological malignancies and solid tumours are chromosomal translocations, which may lead to gene fusions. Recently, next-generation sequencing techniques at the transcriptome level (RNA-Seq) have been used to verify known and discover novel transcribed gene fusions. We present FusionFinder, a Perl-based software designed to automate the discovery of candidate gene fusion partners from single-end (SE) or paired-end (PE) RNA-Seq read data. FusionFinder was applied to data from a previously published analysis of the K562 chronic myeloid leukaemia (CML) cell line. Using FusionFinder we successfully replicated the findings of this study and detected additional previously unreported fusion genes in their dataset, which were confirmed experimentally. These included two isoforms of a fusion involving the genes BRK1 and VHL, whose co-deletion has previously been associated with the prevalence and severity of renal-cell carcinoma. FusionFinder is made freely available for non-commercial use and can be downloaded from the project website (http://bioinformatics.childhealthresearch.org.au/software/fusionfinder/).

Novel BRD4-NUT fusion isoforms increase the pathogenic complexity in NUT midline carcinoma.

Nuclear protein in testis (NUT)-midline carcinoma (NMC) is a rare, aggressive disease typically presenting with a single t(15;19) translocation that results in the generation of a bromodomain-containing protein 4 (BRD4)–NUT fusion. PER-624 is a cell line generated from an NMC patient with an unusually complex karyotype that gave no initial indication of the involvement of the NUT locus. Analysis of PER-624 next-generation transcriptome sequencing (RNA-Seq) using the algorithm FusionFinder identified a novel transcript in which Exon 15 of BRD4 was fused to Exon 2 of NUT, therefore differing from all published NMC fusion transcripts. The three additional exons contained in the PER-624 fusion encode a series of polyproline repeats, with one predicted to form a helix. In the NMC cell line PER-403, we identified the ‘standard’ NMC fusion and two novel isoforms. Knockdown by small interfering RNA in either cell line resulted in decreased proliferation, increased cell size and expression of cytokeratins consistent with epithelial differentiation. These data demonstrate that the novel BRD4–NUT fusion in PER-624 encodes a functional protein that is central to the oncogenic mechanism in these cells. Genomic PCR indicated that in both PER-624 and PER-403, the translocation fuses an intron of BRD4 to a region upstream of the NUT coding sequence. Thus, the generation of BRD4–NUT fusion transcripts through post-translocation RNA-splicing appears to be a common feature of these carcinomas that has not previously been appreciated, with the mechanism facilitating the expression of alternative isoforms of the fusion. Finally, ectopic expression of wild-type NUT, a protein normally restricted to the testis, could be demonstrated in PER-403, indicating additional pathways for aberrant cell signaling in NMC. This study contributes to our understanding of the genetic diversity of NMC, an important step towards finding therapeutic targets for a disease that is refractory to current treatments.

First Genome-Wide Association Study in an Australian Aboriginal Population Provides Insights into Genetic Risk Factors for Body Mass Index and Type 2 Diabetes

A body mass index (BMI) >22kg/m2 is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,075,436 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 P genotyped = 1.50x10-6; rs11140653 Pimputed_1000G = 2.90x10-7) for BMI lies 5’ of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 Pgenotyped = 8.06x10-6; rs10841048 Pimputed_1000G = 6.28x10-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 Pgenotyped = 4.65x10-5; rs13225016 Pimputed_1000G = 6.57x10-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 Pgenotyped = 5.59x10-6, Pimputed_1000G = 5.73x10-6) for T2D lies 5’ of BCL9 that, along with TCF7L2, promotes beta-catenin’s transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (Pgenotyped = 4.49x10-4) for BMI and IGF2BP2 Pimputed_1000G = 2.55x10-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

The International Collaboration for Autism Registry Epidemiology (iCARE): Multinational Registry-Based Investigations of Autism Risk Factors and Trends

The International Collaboration for Autism Registry Epidemiology (iCARE) is the first multinational research consortium (Australia, Denmark, Finland, Israel, Norway, Sweden, USA) to promote research in autism geographical and temporal heterogeneity, phenotype, family and life course patterns, and etiology. iCARE devised solutions to challenges in multinational collaboration concerning data access security, confidentiality and management. Data are obtained by integrating existing national or state-wide, population-based, individual-level data systems and undergo rigorous harmonization and quality control processes. Analyses are performed using database federation via a computational infrastructure with a secure, web-based, interface. iCARE provides a unique, unprecedented resource in autism research that will significantly enhance the ability to detect environmental and genetic contributions to the causes and life course of autism.

ViPAR: a software platform for the Virtual Pooling and Analysis of Research Data

Abstract Background: Research studies exploring the determinants of disease require sufficient statistical power to detect meaningful effects. Sample size is often increased through centralized pooling of disparately located datasets, though ethical, privacy and data ownership issues can often hamper this process. Methods that facilitate the sharing of research data that are sympathetic with these issues and which allow flexible and detailed statistical analyses are therefore in critical need. We have created a software platform for the Virtual Pooling and Analysis of Research data (ViPAR), which employs free and open source methods to provide researchers with a web-based platform to analyse datasets housed in disparate locations. Methods: Database federation permits controlled access to remotely located datasets from a central location. The Secure Shell protocol allows data to be securely exchanged between devices over an insecure network. ViPAR combines these free technologies into a solution that facilitates ‘virtual pooling’ where data can be temporarily pooled into computer memory and made available for analysis without the need for permanent central storage. Results: Within the ViPAR infrastructure, remote sites manage their own harmonized research dataset in a database hosted at their site, while a central server hosts the data federation component and a secure analysis portal. When an analysis is initiated, requested data are retrieved from each remote site and virtually pooled at the central site. The data are then analysed by statistical software and, on completion, results of the analysis are returned to the user and the virtually pooled data are removed from memory. Conclusions: ViPAR is a secure, flexible and powerful analysis platform built on open source technology that is currently in use by large international consortia, and is made publicly available at [ http://bioinformatics.childhealthresearch.org.au/software/vipar/ ].

Official ERS technical standards: Global Lung Function Initiative reference values for the carbon monoxide transfer factor for Caucasians

There are numerous reference equations available for the single-breath transfer factor of the lung for carbon monoxide (T LCO); however, it is not always clear which reference set should be used in clinical practice. The aim of the study was to develop the Global Lung Function Initiative (GLI) all-age reference values for T LCO. Data from 19 centres in 14 countries were collected to define T LCO reference values. Similar to the GLI spirometry project, reference values were derived using the LMS (lambda, mu, sigma) method and the GAMLSS (generalised additive models for location, scale and shape) programme in R. 12 660 T LCO measurements from asymptomatic, lifetime nonsmokers were submitted; 85% of the submitted data were from Caucasians. All data were uncorrected for haemoglobin concentration. Following adjustments for elevation above sea level, gas concentration and assumptions used for calculating the anatomic dead space volume, there was a high degree of overlap between the datasets. Reference values for Caucasians aged 5–85 years were derived for T LCO, transfer coefficient of the lung for carbon monoxide and alveolar volume. This is the largest collection of normative T LCO data, and the first global reference values available for T LCO. This is the largest collection of normative TLCO data and represents a step towards standardised interpretation http://ow.ly/4PcZ30dB1tn